Impact of COVID-19 lockdown on the fine particulate matter concentration levels: Results from Bengaluru megacity,India

Leveraging the COVID-19 India-wide lockdown situation, the present study attempts to quantify the reduction in the ambient fine particulate matter concentrations during the lockdown (compared with that of the pre-lockdown period), owing to the highly reduced specific anthropogenic activities and thereby pollutant emissions. The study was conducted over Bengaluru (India), using PM_2.5 (mass concentration of particulate matter having size less than or equal to 2.5 µm) and Black Carbon mass concentration (BC) data. Open-access datasets from pollution control board (PCB) were also utilised to understand the spatial variability and region-specific reduction in PM_2.5 across the city. The highest percentage reduction was observed in BC_ff (black carbon attributable to fossil fuel combustion), followed by total BC and PM_2.5. No decrease in BC_bb (black carbon attributable to wood/biomass burning) was observed, suggesting unaltered wood-based cooking activities and biomass-burning (local/regional) throughout the study period. Results support the general understanding of multi-source (natural and anthropogenic) nature of PM_2.5 in contrast to limited-source (combustion based) nature of BC. The diurnal amplitudes in BC and BC_ff were reduced, while they remained almost the same for PM_2.5 and BC_bb. Analysis of PCB data reveal the highest reduction in PM_2.5 in an industrial cluster area. The current lockdown situation acted as a natural model to understand the role of a few major anthropogenic activities (viz., traffic, construction, industries related to non-essential goods, etc.) in enhancing the background fine particulate matter levels. Contemporary studies reporting reduction in surface fine particulate matter and satellite retrieved columnar Aerosol Optical Depth (AOD) during COVID-19 lockdown period are discussed.     

The 16-day waves in the mesosphere and lower thermosphere over Wuhan (30.6°N, 114.5°E) and Adelaide (35°S, 138°E)

Winds from a meteor radar at Wuhan (30.6°N, 114.5°E) and a MF radar at Adelaide (35°S, 138°E) are used to study the 16-day waves in the mesosphere and lower thermosphere (MLT). The height range is 78–98 km at Wuhan and 70–98 km at Adelaide. By comparison, it is found that the zonal components at both sites are generally larger than the meridional ones, and eastward motion of the zonal background winds is favorable for the 16-day waves penetration to the MLT region. The zonal maximum amplitude appears in the autumn (September–October) around 86–98 km at Wuhan and in the winter months and early spring (July–October) around 72–82 km at Adelaide. Differences are found in wave amplitudes and time of appearance between the two years of 2002 and 2003. In 2003, the intensity of the wave amplitudes is relatively smaller than that for 2002 at both sites. The summer 16-day waves are comparatively weaker at Adelaide in both years, but stronger in 2002 at Wuhan near the mesopause and the lower thermosphere (86–98 km). The strong summer waves at Wuhan may come from the winter southern hemisphere.     

Influence of land surface parameters on the spatio-seasonal land surface temperature regime in rural West Bengal,India

The magnitude and causes of changes in the land surface temperature of rural areas have not been extensively studied. The thermal band of Landsat imagery is taken to extract winter, summer, and monsoon season land surface temperature (LST) and relate it to surface parameters over a 30-year period. From the extracted parameters constructed a prospective surface temperature (PST) model using Multivariate Adaptive Regression Splines. The Chandrabhaga river basin in West Bengal of the lateritic Rarh Tract at the Chota Nagpur Plateau fringe was chosen as the study area because it is far from urban influences, to avoid the well-known heat island effect. Over the study period, summer and winter average LST increased linearly by 0.085?°C/y and 0.016?°C/y respectively. These results were validated with air temperature (RMSE?=?x and y, respectively). Over time more of the area is in the higher temperature zones, e.g., in April 2011, 4% area exceeded >32°, whereas in 2015 this proportion reached 52%. PST models of all the seasons were moderate to highly correlate (0.57–0.87) with actual LST, showing the value of this model. It also revealed the relative importance of the regional factors. Based on this information factor management is a scientific step to restrict or minimize the temperature rise effect.     

Modeling the effects of green alternative on heat island mitigation of a meso level town,West Bengal,India

Urban heat island (UHI) is emerged as a ubiquitous phenomenon that affects the outdoor thermal comfort. Hence urban microclimatic studies using numeric simulation software to find out adaption strategies are attracting the great attention of researchers and policymaker. Number of scientists carried out their research on this topic using different tools and technique whereas ENVI-met model is the most widely used simulation tools. The present study intended to examine the implication of green infrastructure of the town on enhanced microclimatic condition and to evaluate the best suitable mitigation strategy modelling cool city with ENVI-met software (V4). Five greening condition such as existing UHI condition (C1), 100% greening of roof (C2), 100% greening of roof and walls (C3), 50% greening of roof and walls (C4) and: Plantation at suitable area with 50% greening of roof and walls (C5) have been taken into consideration for the microclimate simulation for all three selected sites. Among these five conditions, C3 for the open mid-rise and compact low-rise, C5 for the open low-rise are identified as one of the most suitable strategies which can reduce the air temperature of peak hours by 2.6 °C, 1.33 °C and 1.87 °C respectively. These models are validated by simple linear regression between simulated and existing air temperature in case of existing UHI condition (C1) and in all the cases coefficient of determination value is high such as for open mid-rise, compact low-rise and open low-rise, it is 0.92, 0.92 and 0.75 respectively. Therefore it can be concluded that the application of those strategies can improve the urban thermal environment as well as the outdoor thermal comfort of English bazaar Town and its surroundings.     

Validation of non-linear split window algorithm for land surface temperature estimation using Sentinel-3 satellite imagery: Case study; Tehran Province,Iran

In recent years, land surface temperature (LST) has become critical in environmental studies and earth science. Remote sensing technology enables spatiotemporal monitoring of this parameter on large scales. This parameter can be estimated by satellite images with at least one thermal band. Sentinel-3 SLSTR data provide LST products with a spatial resolution of 1 km. In this research, direct and indirect validation procedures were employed to evaluate the Sentinel-3 SLSTR LST products over the study area in different seasons from 2018 to 2019. The validation method was based on the absolute (direct) evaluation of this product with field data and comparison (indirect) evaluation with the MODIS LST product and the estimated LST using the non-linear split-window (NSW) algorithm. Also, two emissivity estimation methods, (1) NDVI thresholding method (NDVI-THM) and (2) classification-based emissivity method (CBEM), were used to estimate the LST using the NSW method according to the two thermal bands of Sentinel-3 images. Then, the accuracy of these methods in estimating LST was evaluated using field data and temporal changes of vegetation, which the NDVI-THM method generated better results. For indirect evaluation between the Sentinel-3 LST product, MODIS LST product, and LST estimated using NSW, four filters based on spatial and temporal separates between pairs of pixels and pixel quality were used to ensure the accuracy and consistency of the compared pairs of a pixel. In general, the accuracy results of the LST products of MODIS and Sentinel-3, and LST estimated using NSW showed a similar trend for LST changes during the seasons. With respect to the two absolute and comparative validations for the Sentinel-3 LST products, summer with the highest values of bias (?1.24 K), standard deviation (StDv = 2.66 K), and RMSE (2.43 K), and winter with the lowest ones (bias of 0.14 K, StDv of 1.13 K, and RMSE of 1.12 K) provided the worst and best results for the seasons in the period of 2018–2019, respectively. According to both absolute and comparative evaluation results, the Sentinel-3 SLSTR LST products provided reliable results for all seasons on a large temporal and spatial scale over our studied area.     

Application of satellite observations to study the changes of hypoxic conditions in near-bottom water in the western part of Peter the Great Bay (the Sea of Japan)

In this paper we explore the possibilities of applying satellite ocean colour (OC) observations and SST to study the changes in the conditions of hypoxia in the near-bottom water in the western part of Peter the Great Bay. Near-bottom water hypoxia occurs in water bodies with increased organic matter influx when the dissolved oxygen (DO) consumed at its oxidation is not restored. Consumption of most DO is usually attributed to the oxidation of organic matter formed as a result of increased algae growth during water eutrophication. Satellite data on indicators of phytoplankton (chlorophyll-a concentration (Chl) and fluorescence (FLH)) allow to analyze the spatial-temporal changes of this substation. Coloured dissolved organic matter (CDOM), non-algal particles (NAP) influence on satellite Chl estimates and also on near-bottom water hypoxia formation. This study analyzes daily, seasonal, and inter-annual changes in the distributions of indicators (Chl, FLH, the coefficients of light absorption by coloured detrital matter (a_CDM) and light backscattering by suspended particles (b_bp)), based on the instant satellite OC data from MODIS-Aqua. Data on the Chl, the sea surface temperature (SST) from the MODIS-Aqua, the precipitation from the TRMM satellite and the hydrometeorological stations (HMSs), the wind speed and direction from HMS “Vladivostok” are used to study the influence of hydrometeorological conditions on the Chl values. These distributions were compared with the literary information based on field observations of the hypoxia cases in the same area and with the changes in the vertical DO, Chl, temperature, salinity distribution obtained by coastal expeditions in October-November 2010 and February-March 2011. Significant interrelations within 95% confidence level between the satellite Chl, FLH values calculated at the MUMM atmospheric correction and in situ Chl values obtained in the autumn of 2010 were reached separately for the cases with winds of northern and southern directions with the correlation coefficients of 0.71, 0.48 and 0.49, 0.71, respectively. Significant dependences of Chl on SST and Chl on wind speed explained by the influence of continental runoff and water ventilation were obtained. Therefore, the changes of Chl reflect the changes of hypoxic conditions in the near-bottom water. In Amursky Bay the onset of hypoxia was at the Chl and SST values equal to 4 mg m^?3 and 13 °C (↑ – at increasing SST); near Furugelm Island it was at 1.6 mg m^?3 and 25 °C (↑), 1 mg m^?3 and 21 °C (↓). The difference in the Chl values was reflected in the hypoxia onset timings that were the beginning of June (2011), August (2013), and September (2014), respectively. The water flow from the eastern coast of Amursky Bay in early August of 2013 recorded from the OC and SST satellite imagers appeared in an additional hypoxic zone. Decreased OC characteristics in the runoff of the Razdolnaya River in August-September of 2014 were a sign of hypoxia at its mouth. Near Furugelm Island the hypoxia destruction (increase in the DO level from 1 to 4.5 ml L^?1) was observed at the Chl of 0.9 mg m^?3 and SST = 18 °C (↓). At the autumn maximum of Chl equal to 1.7 mg m^?3 and SST = 4 °C (↓) in mid-November the DO level here increased to 8 ml L^?1. In Amursky Bay, short-term destructions/weakening of hypoxia manifested themselves in sharp increases of Chl. At that, the ratio between the Chl value and the approximation level was equal to 2 and higher for SST equal to 22–25 °C (↑), to 0.9 and higher for SST equal to 5–13 °C (↓). With the water stratification destruction in temperature and the noticeable weakening of the stratification in salinity (mid-November), the hypoxia destructed (the DO level increased from 2 to 6 ml L^?1). In this case, Chl and SST were about 3 mg m^?3 and 5 °C (↓).     

Simulation of the land surface temperature from moon-based Earth observations

The land surface temperature (LST) is a key parameter for the Earth’s energy balance. As a natural satellite of the Earth, the orbital of the moon differs from that of current Earth observation satellites. It is a new way to measure the land surface temperature from the moon and has many advantages compared with artificial satellites. In this paper, we present a new method for simulating the LST measured by moon-based Earth observations. Firstly, a modified land-surface diurnal temperature cycle (DTC) method is applied to obtain the global LST at the same coordinated universal time (UTC) using the Moderate Resolution Imaging Spectroradiometer (MODIS) LST products. The lunar elevation angles calculated using the ephemeris data (DE405) from the Jet Propulsion Laboratory (JPL) were then applied to simulate the Earth coverage observed from the moon. At the same time, the modified DTC model was validated using in situ data, MODIS LST products, and the FengYun-2F (FY-2F) LST, respectively. The results show that the fitting accuracy (root-mean-square error, RMSE) of the modified DTC model is not greater than 0.72?°C for eight in situ stations with different land cover types, and the maximum fitting RMSE of the modified model is smaller than that of current DTC models. By the comparison of the simulated LST with MODIS and FY-2F LST products, the errors of the results were feasible and accredited, and the simulated global LST has a reasonable spatiotemporal distribution and change trend. The simulated LST data can therefore be used as base datasets to simulate the thermal infrared imagery from moon-based Earth observations in future research.     

Surface temperature estimation in Singhbhum Shear Zone of India using Landsat-7 ETM+ thermal infrared data

Land surface temperature (LST) is an important factor in global change studies, heat balance and as control for climate change. A comparative study of LST over parts of the Singhbhum Shear Zone in India was undertaken using various emissivity and temperature retrieval algorithms applied on visible and near infrared (VNIR), and thermal infrared (TIR) bands of high resolution Landsat-7 ETM+ imagery. LST results obtained from satellite data of October 26, 2001 and November 2, 2001 through various algorithms were validated with ground measurements collected during satellite overpass. In addition, LST products of MODIS and ASTER were compared with Landsat-7 ETM+ and ground truth data to explore the possibility of using multi-sensor approach in LST monitoring. An image-based dark object subtraction (DOS3) algorithm, which is yet to be tested for LST retrieval, was applied on VNIR bands to obtain atmospheric corrected surface reflectance images. Normalized difference vegetation index (NDVI) was estimated from VNIR reflectance image. Various surface emissivity retrieval algorithms based on NDVI and vegetation proportion were applied to ascertain emissivities of the various land cover categories in the study area in the spectral range of 10.4–12.5 μm. A minimum emissivity value of about 0.95 was observed over the reflective rock body with a maximum of about 0.99 over dense forest. A strong correlation was established between Landsat ETM+ reflectance band 3 and emissivity. Single channel based algorithms were adopted for surface radiance and brightness temperature. Finally, emissivity correction was applied on ‘brightness temperature’ to obtain LST. Estimated LST values obtained from various algorithms were compared with field ground measurements for different land cover categories. LST values obtained after using Valor’s emissivity and single channel equations were best correlated with ground truth temperature. Minimum LST is observed over dense forest as about 26 °C and maximum LST is observed over rock body of about 38 °C. The estimated LST showed that rock bodies, bare soils and built-up areas exhibit higher surface temperatures, while water bodies, agricultural croplands and dense vegetations have lower surface temperatures during the daytime. The accuracy of the estimated LST was within ±2 °C. LST comparison of ASTER and MODIS with Landsat has a maximum difference of 2 °C. Strong correlation was found between LST and spectral radiance of band 6 of Landsat-7 ETM+. Result corroborates the fact that surface temperatures over land use/land cover types are greatly influenced by the amount of vegetation present.     

Assessing the performance of BDS-3 for multi-GNSS static and kinematic PPP-AR

As of 2021, a total of four different GNSS constellations – namely, GPS, GLONASS, Galileo, and BDS-3 – can be used with Full Operational Capability (FOC). In this work, the contribution of BDS-3 FOC to GPS + GLONASS + Galileo (GRE) PPP-AR is investigated, considering the three different cut-off angles (7°, 30°, and 45°) and different lengths of static observation sessions (24-, 12-, 6-, 3-, 1-, 0.5-, 0.25-hour). The data of 31 IGS-MGEX stations is processed with GRE PPP-AR and GREC3 (GPS + GLONASS (using float mode) + Galileo + BDS-3) PPP-AR modes. The results showed that BDS-3 degraded the horizontal (except for 24-h sessions) and vertical accuracy of static GRE PPP-AR solutions regardless of the elevation cutoff angle and observation time. The kinematic results showed that BDS-3 significantly contributed to the accuracy of GRE kinematic PPP-AR for 30° and 45° cutoff angles. The convergence time analysis showed that BDS-3 only contributes to GRE kinematic PPP-AR for the vertical component.     

Nightside 135.6 nm emission enhancements of mid/low latitudes during geomagnetic storms as observed by the ionosphere PhotoMeter (IPM) on the Chinese meteorological satellite FY-3D

IPM has detected nightside 135.6 nm emission enhancements over a wide latitude range, from the sub-auroral latitudes to the equatorial regions during geomagnetic storms. Our work, presented in this paper, uses the data of IPM to understand these 135.6 nm emission enhancements during of geomagnetic storms and studies the variations of total electron content (TEC) and the F₂ layer peak electron density (NmF2) in the region of enhanced emissions. Middle and low latitude emission enhancements are presented during several medium storms in 2018. The variations of both the integrated electron content (IEC) derived from the nighttime OI 135.6 nm emission by IPM and TEC from the International GNSS Service (IGS) relative to the daily mean of magnetically quiet days of per each latitude bin (30°≦geographic latitude < 40°, 15°≦geographic latitude < 30°, 0°≦geographic latitude < 15°, ?15°≦geographic latitude < 0°, ?30°≦geographic latitude < -15°, ?40°≦geographic latitude < -30°) are investigated and show that on magnetically storm day, IEC by IPM always increases, while TEC from IGC may increase or decrease. Even if both increase, the increase of IEC is greater than that of TEC. From the comparison of IEC and TEC during magnetic storms, it can be seen that the enhancement of the nighttime 135.6 nm emissions is not entirely due to the ionospheric change. The time of IEC enhancements at each latitude bin is in good agreement, which mainly corresponds to the main phase time of the geomagnetic storm event and lasts until the recovery phase. The available ground-based ionosonde stations provide the values of NmF2 which match the 135.6 nm emissions measured by IPM in space and time. The variations of NmF2 squared can characterize the variations of the OI 135.6 nm emissions caused by O⁺ ions and electrons radiative recombination. The study results show that the OI 135.6 nm emission enhancements caused by O⁺ ions and electrons radiative recombination (where NmF2 squared increases) are obviously a contribution to the measured 135.6 nm emission enhancements by IPM. The contribution accounts for at least one of all contributions to the measured 135.6 nm emission enhancements by IPM. However, where the NmF2 squared provided by ionosonde decrease or change little (where the OI 135.6 nm emissions cause by O⁺ ions and electrons radiative recombination also decrease or change little), the emission enhancements measured by IPM at storm-time appear to come from the contributions of other mechanisms, such as energetic neutral atoms precipitation, or the mutual neutralization emission (O+ + O-→2O + h? (135.6 nm)) which also occupies a certain proportion in 135.6 nm airglow emission at night.     

Multi-wavelength coordinated observations of ionospheric irregularity structures from an anomaly crest location during unusual solar minimum of the 24th cycle

The present paper reports coordinated ionospheric irregularity measurements at optical as well as GPS wavelengths. Optical measurements were obtained from Tiny Ionospheric Photometer (TIP) sensors installed onboard the Constellation Observing System for Meteorology, Ionosphere, and Climate (COSMIC) satellites. GPS radio signals were obtained from a dual frequency GPS receiver operational at Calcutta (22.58°N, 88.38°E geographic; geomagnetic dip: 32.96°; 13.00°N, 161.63°E geomagnetic) under the SCIntillation Network Decision Aid (SCINDA) program. Calcutta is located near the northern crest of Equatorial Ionization Anomaly (EIA) in the Indian longitude sector. The observations were conducted during the unusually low and prolonged solar minima period of 2008–2010. During this period, four cases of post-sunset GPS scintillation were observed from Calcutta. Among those cases, simultaneous fluctuations in GPS Carrier-to-Noise ratios (C/N_o) and measured radiances from TIP over a common ionospheric volume were observed only on February 2, 2008 and September 25, 2008. Fluctuations observed in measured radiances (maximum 0.95 Rayleigh) from TIP due to ionospheric irregularities were found to correspond well with C/N₀ fluctuations on the GPS links observed from Calcutta, such effects being noted even during late evening hours of 21:00–22:00 LT from locations around 40° magnetic dip. These measurements indicate the existence of electron density irregularities of scale sizes varying over several decades from 135.6 nm to 300–400 m well beyond the northern crest of the EIA in the Indian longitude sector during late evening hours even in the unusually low solar activity conditions.     

Some features of the day-to-day MLT wind variability in winter 2017–2018 as seen with a European/Siberian meteor radar network

We present results of wind measurements near the mesopause carried out with meteor radars (MRs) at Collm (51°N, 13°E), Obninsk (55°N, 37°E), Kazan (56°N, 49°E), Angarsk (52°N, 104°E) and Anadyr (65°N, 178°E) from October 1, 2017 till March 31, 2018. The Collm and Kazan MRs are SKiYMET radars with vertical transmission and radio echo height finding, while the other radars operate with horizontal transmission and without height finding. We paid particular attention to the meridional wind variability with periods of 4–6 days and 9–11 days. The waves with these periods are seen as spots of the wave activity in the wavelet spectra and include oscillations with different periods and different discrete zonal wavenumbers. These wave packets successively propagate as a group of waves from one site to another one in such a way that they are observed at one site and almost disappear at the previous one. The 4–6 wave group includes planetary-scale oscillations (individual spectral components) which have eastward phase velocities and mostly zonal wavenumbers 2 and 3, and the vertical wavelength exceeds 70 km at middle latitudes. The source of the oscillations is the polar jet instability. The wave group itself propagates westward, and the amplitudes of wind oscillations are approximately 5–6 m/s as obtained from the wind data averaged over the meteor zone. The 9–11 day wave set propagates westward as a group and mainly consists of spectral components which have westward phase velocity and zonal wavenumber 1. Amplitudes of these wind perturbations strongly vary from station to station and can reach, approximately, 8 m/s. The vertical wavenumber is 0.014 km⁻¹ as taken from the Kazan and 0.05 km⁻¹ according to the Collm data. We obtained a global view on the waves by using the AURA MLS geopotential data. We found a good correspondence between wave features obtained from the MR wind measurements and the MLS data. To our knowledge, such a wave propagation of planetary wave in the mesosphere/lower thermosphere (MLT) region has so far not obtained much attention.     

Substitution of satellite-based land surface temperature defective data using GSP method

Land surface temperature (LST) as an important environmental variable provides valuable information for earth environmental system modelling. Currently, LST is obtained through satellite thermal sensors at various spatial and temporal resolutions. Although spatially continuous satellite-based LST measurements are intended to overcome the shortcomings of sparse ground-based LST measurements, LST images often contain anomalous values due to the existence of clouds or sensor malfunctioning. The problem becomes more serious where the users deal with high spatial resolution characterized by low temporal resolution. This study examines the capability of a newly developed graph signal processing (GSP) method using two-dimensional single-date thermal data. For this purpose, four Landsat/TIRS datasets are analyzed. The data of five elliptical regions on thermal images are eliminated and then reconstructed through the GSP method and using the LST values of the enclosing rectangles containing the ellipsoids. The results indicate that the temperature variation determined by the GSP method generally conforms to the original image LST values. According to a correlation test conducted on the original image LST and those obtained through the GSP method, the values vary from 58% to 95%, which is an above-the-average rate (RMSE from 0.69 to 2.27). The statistical analysis of the original image LST in both the elliptical regions and the enclosing rectangles containing the ellipsoids indicates that an increase in the variance of LST data causes an increased error in the calculation of temperature by the GSP method, and vice versa. The results of the analysis of variance (ANOVA) and Duncan test indicated that an increase in the number of the non-zero spectral bins would result in increased RMSE values for all the dates and the regions. Moreover, the model errors were significant at the 0.05 level across all the image date and five elliptical study regions. Based on the results, the use of this method is recommended for the reconstruction of LST missing values, where dissimilarity of atmospheric conditions limits the use of other methods that depend on the time series data of various dates and a great deal of data calculation.     

Observation of aphelion cloud belt over Martian tropics,its evolution,and associated dust distribution from MCS data

The present study uses five Martian years of observations from Mars Climate Sounder onboard Mars Reconnaissance Orbiter for investigating the Aphelion Cloud Belt (ACB) over the tropics. Analysis of zonal mean water ice column opacity suggests that the spatial extension of the ACB is mainly confined over the tropics and mid-latitudes (-20 – 40°N) during L_S ~ 45 – 135° (L_S = 0° signifies northern spring equinox). The ACB is seen primarily in the nighttime only due to the truncation of the daytime profile observations at significantly higher altitudes (at ~30 km). Zonal mean ice extinction profiles show ACB’s altitudinal range within ~10 – 40 km, and the existence of a thin cloud band in the absence of a thick ACB during aphelion season. Three phases of the ACB could be identified as the formation phase during L_S = 45 – 75° (phase 1), the peak phase during L_S = 76 – 105° (phase 2), and the decaying phase during L_S = 106 – 135° (phase 3). Observation of the cloud latitude belt shows a northward movement starting from phase 2, prominent over regions nearby Lunae Planum and Xanthe Terra. During this phase, the top level of thick clouds within the ACB decreases to ~20 km in the southern hemisphere, while it increases a little over the northern hemisphere (NH). The decreasing tendency continues in phase 3 over the entire region ?10 – 10°N, and the thick cloud base moves higher over the NH, though the vertical depth of it becomes narrower than phase 2. Temperature profiles do not show any noticeable influence on the northward evolution of the ACB. However, the study at a regional level indicates a possible association of upper tropospheric dustiness with the ACB’s evolution. The mechanism is evident in the correlation analysis mostly at an altitude range of ~18 – 35 km. The migrating semidiurnal tide (SMD) as a proxy of dust or water ice forcing, and the calculated upper tropospheric dust radiative heating, shows an apparent northward movement of their peak amplitude within the three phases of the ACB. This match between the spatiotemporal variations of the SMD and the water ice was not observed previously. However, the correlating behavior seems to be prominent in the areas nearby Lunae Planum and Xanthe Terra and the upper-tropospheric region of the atmosphere.     

激光雷达观测的武汉上空对流边界层高度与夹卷层厚度

基于地基偏振激光雷达2011年1月至2017年12月对武汉中心城区的观测资料，研究了超大城市对流边界层高度与夹卷层厚度的变化特征.采用Fernald方法进行数据反演获得了1min时间分辨率、30m空间分辨率的后散比剖面，通过方差法确定了对流边界层高度和夹卷层厚度.结果表明，武汉上空对流边界层高度和夹卷层厚度具有明显的季节变化特征.对流边界层顶的均值最大处在春季为1.14km，夏季为1.25km，秋季为1.06km，冬季为0.74km；夹卷层厚度均值最大处在春季为0.40km，夏季为0.51km，秋季为0.34km，冬季为0.26km.这些特征与武汉地区地表温度的周年变化特征具有很强的相关性.      

Modeling of effect of polarization on UV sky radiance during twilight

The spatial distribution of the vector of the Stokes parameters characterizing the radiance intensity and the radiance polarization describes the radiation field in the atmosphere. A simplified treatment of light as the scalar has only restricted application. A few studies compared previously results of the vector and scalar radiative transfer models and showed that scalar models are in error by up to 10% for many cases. Though several observational conditions were exploited, an effect of polarization on modeling of UV radiance has not been investigated yet for twilight. The paper presents a preliminary study of modeled UV radiance during twilight taking into account polarization. The intensity and the degree of linear polarization of the scattered UV radiance for two cases of the ground-based observations are discussed. In the first case, radiation incoming from the zenith for the solar zenith angles (SZA) from 90° to 98° is under investigation. Radiation in the solar principal plane for the beginning of twilight (SZA = 90.1°) was calculated in the second case. The study showed that the UV radiation field in the twilight atmosphere can be handled correctly only using the vector theory. The errors of scalar radiative transfer strongly depend on wavelength, line of an observation and solar position. The revealed distortion of the zenith radiance caused by using of the scalar approximation reaches maximum of 15% at 340 nm for the solar zenith angle (SZA) equal to 98°. The shorter wavelengths have the smaller errors, about 5% at 305 nm for SZA = 98°, due to the larger part of the single scattered radiance. The error of the scalar modeling may be as large as −17% for radiance incoming from the horizon for SZA = 90.1°. Scalar radiative transfer models underestimate the integral intensity in the principal plane up to 3–4% ± 0.5% at SZA = 90.1° for wavelengths from 320 to 340 nm. This should be taken into account in problems of radiative budget estimation and remote sensing of the atmosphere exploiting the twilight period.     

Using land use change trajectories to quantify the effects of urbanization on urban heat island

This paper proposed a quantitative method of land use change trajectory, which means the succession among different land use types across time, to examine the effects of urbanization on an urban heat island (UHI). To accomplish this, multi-temporal images from Landsat 5 Thematic Mapper (TM) and Landsat 7 Enhanced Thematic Mapper Plus (ETM+) of Xiamen City in China from 1987 to 2007 were selected. First, the land use change trajectories were extracted based on the classified images from different years. Then the land surface temperatures (LST) were retrieved and the magnitudes of the UHI were evaluated using the UHI intensity (UHII) indicator. Finally, the indices of the contribution to UHI intensity (CUHII) were constructed and calculated to quantify the effects of each land use change trajectory on the UHI during urbanization. The results demonstrated that the land use change trajectories and CUHII are effective and useful in quantifying the effects of urbanization on UHI. In Xiamen City, a total of 2218 land use change trajectories were identified and 530 of them were the existing urban or urbanization trajectories. The UHII presents a trend of continuous increase from 0.83 °C in 1987 to 2.14 °C in 2007. With respect to the effects of urban growth on UHI, the contribution of existing urban area to UHI decreased during urbanization. Prior to 2007, the existing urban area of trajectory NO. 44444 had the most significant effect on UHI with the greatest CUHII, while the value has decreased from 55.00% in 1987 to 13.03% in 2007 because of the addition of new urbanized area. In 2007, the greatest CUHII was replaced by a trajectory from farmland to built-up area (NO. 22224) with the CUHII of 21.98%, followed by the existing urban area of trajectory NO. 44444 with the CUHII of 13.03%. These results provide not only a new methodology to assess the environmental effects of urbanization, but also decision-supports for the planning and management of cities.     

Application of Wuhan Ionospheric Oblique Backscattering Sounding System (WIOBSS) for the investigation of midlatitude ionospheric irregularities

An upgrade of Wuhan Ionospheric Backscattering Sounding System (WIOBSS) was developed in 2015. Based on the Universal Serial Bus (USB), and a high performance FPGA, the newly designed WIOBSS has a completely digital structure, which makes it portable and flexible. Two identical WIOBSSs, which were situated at Mile (24.31°N, 103.39°E) and Puer (22.74°N, 101.05°E) respectively, were used to investigate the ionospheric irregularities. The comparisons of group distance, Doppler shift and width between Mile-Puer and Puer-Mile VHF ionospheric propagation paths indicate that the reciprocity of the irregularities is satisfied at midlatitude region. The WIOBSS is robust in the detection of ionospheric irregularities.     

Satellite evidence of harmful algal blooms and related oceanographic features in the Bohai Sea during autumn 1998

Harmful algal blooms (HABs) are truly global marine phenomena of increasing significance. Some HAB occurrences are different to observe because of their high spatial and temporal variability and their advection, once formed, by surface currents. A serious HAB occurred in the Bohai Sea during autumn 1998, causing the largest fisheries economic loss. The present study analyzes the formation, distribution, and advection of HAB using satellite SeaWiFS ocean color data and other oceanographic data. The results show that the bloom originated in the western coastal waters of the Bohai Sea in early September, and developed southeastward when sea surface temperature (SST) increased to 25–26 °C. The bloom with a high Chl-a concentration (6.5 mg m⁻³) in center portion covered an area of 60 × 65 km². At the end of September, the bloom decayed when SST decreased to 22–23 °C. The HAB may have been initiated by a combination of the river discharge nutrients in the west coastal waters and the increase of SST; afterwards it may have been transported eastward by the local circulation that was enhanced by northwesterly winds in late September and early October.     

Semi-annual variations of chlorophyll concentration on the Eastern Tropical Pacific coast of Mexico

An analysis of the main remotly sensed oceanographic variables was conducted to characterize the dominant drivers that modulate the spatial-temporal variability of coastal phytoplankton biomass in the northern limit of the Eastern Tropical Pacific. The phytoplankton biomass was analyzed using monthly average satellite chlorophyll (Chlo) concentration data from MODIS sensors for the period 2003–2017. Oceanographic conditions were analyzed by considering (i) Sea Surface Temperature (SST) high-resolution data from the GHRSST project, (ii) wind stress calculated with data from the zonal and meridional components of the CCMP product, (iii) climatological precipitation, and (iv) climatological river flow. Chlo time series and spatial variability were analyzed using Hovmöller diagrams and Empirical Orthogonal Functions (EOF), respectively. A strong semi-annual signal in Chlo concentration along the coast was observed: the first peak occurs in winter-spring (5.0 mg·m⁻³) and a second one in summer-autumn (6.5 mg·m⁻³). A high year-round average of Chlo concentration (3.0–15.0 mg·m⁻³) was maintained by a 10 km wide along-shore coast, with maximum values spatially associated with river mouths, of which Santiago river registered the highest Chlo values (20.0 mg·m⁻³). Surface oceanographic conditions showed a marked annual cycle with warmer conditions (30–31 °C) from July to October and colder ones (23–24 °C) from December to April; during the cold period, seasonal wind stress stimulated coastal upwelling. The EOF showed that 70% of Chlo concentration variability was controlled by the semi-annual pattern, which responded to coastal upwelling conditions during cold period, while rivers outflows influenced high Chlo concentration during the warm period. These results highlighted the importance for land-ocean interface to sustain coastal ecosystems' biological production and the major role of watersheds as sources of nutrients to maintain high biological production during warm periods in transitional tropical-subtropical zones.