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.     

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.     

Determination of land surface temperature using precipitable water based Split-Window and Artificial Neural Network in Turkey

Land surface temperature (LST) calculation utilizing satellite thermal images is very difficult due to the great temporal variance of atmospheric water vapor in the atmosphere which strongly affects the thermal radiance incoming to satellite sensors. In this study, Split-Window (SW) and Radial Basis Function (RBF) methods were utilized for prediction of LST using precipitable water for Turkey. Coll 94 Split-Window algorithm was modified using regional precipitable water values estimated from upper-air Radiosond observations for the years 1990–2007 and Local Split-Window (LSW) algorithms were generated for the study area. Using local algorithms and Advanced Very High Resolution Radiometer (AVHRR) data, monthly mean daily sum LST values were calculated. In RBF method latitude, longitude, altitude, surface emissivity, sun shine duration and precipitable water values were used as input variables of the structure. Correlation coefficients between estimated and measured LST values were obtained as 99.23% (for RBF) and 94.48% (for LSW) at 00:00 UTC and 92.77% (for RBF) and 89.98% (for LSW) at 12:00 UTC. These meaningful statistical results suggest that RBF and LSW methods could be used for LST calculation.     

Seasonal variability of atmospheric aerosol over the North Indian region during 2005–2009

The Indo-Gangetic basin (IGB) extends 2000 km in length along NW–SE and has 400 km width, in the north the basin is bounded by towering Himalaya. High aerosol optical depth (AOD) is observed over the IGB throughout the year. The Himalaya restricts the transport of aerosols across Tibet and China. We have used ground based Kanpur and Gandhi College Aerosol Robotic Network (AERONET) stations and Multiangle Imaging SpectroRadiometer (MISR) and Moderate Resolution Imaging Spectroradiometer (MODIS) Terra level-3 AOD products for the years 2005–2009 to study the variability of aerosol over the Indo-Gangetic (IG) plains. An increase in both satellite-derived as well as ground observed aerosol loading during 2005–2009 has been found over major cities located in the IG plains. The correlation coefficients between AERONET and MISR data are found to be 0.70, 0.36 0.82, in contrast the correlation coefficients between AERONET and MODIS 0.49, 0.68, and 0.43, respectively during summer, winter and monsoon seasons. The AOD estimation using MISR is found to be close to AERONET data during summer and monsoon seasons, in contrast MODIS estimation is better during winter season.     

Exploring the effect of COVID-19 pandemic lockdowns on urban cooling: A tale of three cities

COVID-19 pandemic has had a major impact on our society, environment and public health, in both positive and negative ways. The main aim of this study is to monitor the effect of COVID-19 pandemic lockdowns on urban cooling. To do so, satellite images of Landsat 8 for Milan and Rome in Italy, and Wuhan in China were used to look at pre-lockdown and during the lockdown. First, the surface biophysical characteristics for the pre-lockdown and within-lockdown dates of COVID-19 were calculated. Then, the land surface temperature (LST) retrieved from Landsat thermal data was normalized based on cold pixels LST and statistical parameters of normalized LST (NLST) were calculated. Thereafter, the correlation coefficient (r) between the NLST and index-based built-up index (IBI) was estimated. Finally, the surface urban heat island intensity (SUHII) of different cities on the lockdown and pre-lockdown periods was compared with each other. The mean NLST of built-up lands in Milan (from 7.71 °C to 2.32 °C), Rome (from 5.05 °C to 3.54 °C) and Wuhan (from 3.57 °C to 1.77 °C) decreased during the lockdown dates compared to pre-lockdown dates. The r (absolute value) between NLST and IBI for Milan, Rome and Wuhan decreased from 0.43, 0.41 and 0.16 in the pre-lockdown dates to 0.25, 0.24, and 0.12 during lockdown dates respectively, which shows a large decrease for all cities. Analysis of SUHI for these cities showed that SUHII during the lockdown dates compared to pre-lockdown dates decreased by 0.89 °C, 1.78 °C, and 1.07 °C respectively. The results indicated a high and substantial impact of anthropogenic activities and anthropogenic heat flux (AHF) on the SUHI due to the substantial reduction of huge anthropogenic pressure in cities. Our conclusions draw attention to the contribution of COVID-19 lockdowns (reducing the anthropogenic activities) to creating cooler cities.     

基于DTW的HJ-1B卫星地表温度产品质量检验

HJ-1B卫星热红外数据应用广泛,但其地表温度反演产品的质量检验工作尚不完善。以黑河流域为研究区,利用普适性单通道算法得到HJ 1B地表温度,基于7个地面站点（下垫面为荒漠、沙漠、植被、农作物、雪地和湿地）同步观测资料和MODIS地表温度产品（MOD11A1）,引入动态时间规整方法（DTW）对站点处HJ 1B地表温度进行验证。试验结果表明：HJ 1B地表温度反演产品与地面观测值的偏差值在沙漠和荒漠站点为1K以内,均方根误差在05K左右;对于植被和农作物站点的偏差在2K以内,均方根误差为1~2K;基于DTW的验证对时序不匹配的数据评价结果与现有指标表现一致。HJ-1B地表温度反演产品与地面站点的相对偏差均低于其与MODIS地表温度反演产品的相对偏差。     

Exploring the potential of Sentinel-3 delay Doppler altimetry for enhanced detection of coastal currents along the Northwest Atlantic shelf

In this study, we evaluate Sentinel-3A satellite synthetic aperture radar (SAR) altimeter observations along the Northwest Atlantic coast, spanning the Nova Scotian Shelf, Gulf of Maine, and Mid-Atlantic Bight. Comparisons are made of altimeter sea surface height (SSH) measurements from three different altimeter data processing approaches: fully-focused synthetic aperture radar (FFSAR), un-focused SAR (UFSAR), and conventional low-resolution mode (LRM). Results show that fully-focused SAR data always outperform LRM data and are comparable or slightly better than the nominal un-focused SAR product. SSH measurement noise in both SAR-mode datasets is significantly reduced compared to LRM. FFSAR SSH 20-Hz noise levels, derived from 80-Hz FFSAR data, are lower than 20-Hz UFSAR SSH with 25% noise reduction offshore of 5 km, and 55–70% within 5 km of the coast. The offshore noise improvement is most likely due to the higher native along-track data posting rate (80 Hz for FFSAR, and 20 Hz for UFSAR), while the large coastal improvement indicates an apparent FFSAR data processing advantage approaching the coastlines. FFSAR-derived geostrophic ocean current estimates exhibit the lowest bias and noise when compared to in situ buoy-measured currents. Assessment at short spatial scales of 5–20 km reveals that Sentinel-3A SAR data provide sharper and more realistic measurement of small-scale sea surface slopes associated with expected nearshore coastal currents and small-scale gyre features that are much less well resolved in conventional altimetric LRM data.     

Impact of rain on wave retrieval from Sentinel-1 synthetic aperture radar images in tropical cyclones

The main objective of our work was to investigate the impact of rain on wave observations from C-band (~5.3 GHz) synthetic aperture radar (SAR) in tropical cyclones. In this study, 10 Sentinel-1 SAR images were available from the Satellite Hurricane Observation Campaign, which were taken under cyclonic conditions during the 2016 hurricane season. The third-generation wave model, known as Simulating WAves Nearshore (SWAN) (version 41.31), was used to simulate the wave fields corresponding to these Sentinel-1 SAR images. In addition, rainfall data from the Tropical Rainfall Measuring Mission satellite passing over the spatial coverage of the Sentinel-1 SAR images were collected. The simulated results were validated against significant wave heights (SWHs) from the Jason-2 altimeter and European Centre for Medium-Range Weather Forecasts data, revealing a root mean square error (RMSE) of ~0.5 m with a 0.25 scatter index. Winds retrieved from the VH-polarized Sentinel-1 SAR images using the Sentinel-1 Extra Wide-swath Mode Wind Speed Retrieval Model after Noise Removal were taken as prior information for wave retrieval. It was discovered that rain did indeed affect the SAR wave retrieval, as evidenced by the 3.21-m RMSE of SWHs between the SAR images and the SWAN model, which was obtained for the ~1000 match-ups with raindrops. The raindrops dampened the wave retrieval when the rain rate was < ~5 mm/hr; however, they enhanced wave retrieval for higher rain rates. It was also found that the portion of the rain-induced ring wave with a wave number > 0.05 rad/m (~125 m wavelength) was clearly observed in the SAR-derived wave spectra.     

Exploiting the Sentinel-3 tandem phase dataset and azimuth oversampling to better characterize the sensitivity of SAR altimeter sea surface height to long ocean waves

Being the very first SAR mode altimeter tandem phase, the Sentinel-3 A/B tandem phase has provided an unprecedented opportunity to better characterize the sensitivity of SAR altimetry retrievals to high-frequency processes, such as long ocean waves. In this paper, we show that for some sea-state conditions, that are still to be precisely characterized, long ocean waves are responsible for high-frequency (spatial and temporal) coherent Sea Level Anomaly (SLA) signals. It is found that the peak wavelength corresponds to the dominant swell wavelength. Furthermore, the short time lag between S3-A and S3-B acquisitions allows performing cross-spectral analyses that reveal phase shifts consistent with waves travelling according to the wave dispersion relation. It is also demonstrated that the classical 20 Hz sampling frequency is insufficient to properly sample most swell-induced SLA signals and that aliasing can generate errors over the entire frequency spectrum, including at long wavelengths. These results advocate for the use of azimuth oversampling (40 Hz or 80 Hz). Low-pass filtering should be applied prior to any down-sampling to 20 Hz, in order to prevent long-wavelength errors induced by spectral leakage.     

Validity assessment of SAMOSA retracking for fully-focused SAR altimeter waveforms

We demonstrate in this work how we can take advantage of known unfocused SAR (UF-SAR) retracking methods (e.g. the physical SAMOSA model) for retracking of fully-focused SAR (FF-SAR) waveforms. Our insights are an important step towards consistent observations of sea surface height, significant wave height and backscatter coefficient (wind speed) with both UF-SAR and FF-SAR. This is of particular interest for SAR altimetry in the coastal zone, since coastal clutter may be filtered out more efficiently in the high-resolution FF-SAR waveform data, which has the potential to improve data quality. We implemented a multi-mission FF-SAR altimetry processor for Sentinel-3 (S3) and Sentinel-6 Michael Freilich (S6), using a back-projection algorithm, and analysed ocean waveform statistics compared to multilooked UF-SAR. We find for Sentinel-3 that the averaged power waveforms of UF-SAR and FF-SAR over ocean are virtually identical, while for Sentinel-6 the FF-SAR power waveforms better resemble the UF-SAR zero-Doppler beam. We can explain and model the similarities and differences in the data via theoretical considerations of the waveform integrals. These findings suggest to use the existing UF-SAR SAMOSA model for retracking S3 FF-SAR waveforms but the SAMOSA zero-Doppler beam model for S6 FF-SAR waveforms, instead. Testing the outlined approach over short track segments, we obtain range biases between UF-SAR and FF-SAR lower than 2 mm and significant wave height biases lower than 5 cm.     

Performance assessment of RTPPP positioning with SSR corrections and PPP-AR positioning with FCB for multi-GNSS from MADOCA products

The state-space representation (SSR) product of satellite orbit and clock is one of the most essential corrections for real-time precise point positioning (RTPPP). When it comes to PPP ambiguity resolution (PPP-AR), the fractional cycle bias (FCB) matters. The Japan Aerospace Exploration Agency (JAXA) has developed a multi-GNSS (i.e., global navigation satellite system) advanced demonstration tool for orbit and clock analysis (MADOCA), providing free and precise orbit and clock products. Because of the shortage of relevant studies on performance evaluation, this paper focuses on the performance assessment of RTPPP and PPP-AR by real-time and offline MADOCA products. To begin with, the real-time MADOCA products are evaluated by comparing orbit and clock with JAXA final products, which gives an objective impression of the correction. Second, PPP tests in static and simulated kinematic mode are conducted to further verify the quality of real-time MADOCA products. Finally, the offline MADOCA products are assessed by PPP and PPP-AR comparisons. The results are as follows: (1) Orbit comparisons produced an average error of about 0.04–0.13 m for the global positioning system (GPS), 0.14–0.16 m for the global navigation satellite system (GLONASS), and 0.07–0.08 m for the quasi-zenith satellite system (QZSS). The G15 satellite had the most accurate orbit, with a difference of 0.04 m between the JAXA orbit products and MADOCA’s counterpart, while the R07 satellite had the least accurate orbit with a difference of 0.16 m. Clock products had an accuracy of 0.4–1.3 ns for GPS, 1.4–1.6 ns for GLONASS, and 0.7–0.8 ns for QZSS in general. The G15 satellite had the most accurate clock with a difference of only 0.40 ns between the JAXA clock products and MADOCA products, and the R07 satellite had the least accurate clock with a difference of 1.55 ns. The orbit and clock products for GLONASS performed worse than those of GPS and QZSS. (2) After convergence, the positioning accuracy was 3.0–8.1 cm for static PPP and 8.1–13.7 cm for kinematic PPP when using multi-GNSS observations and precise orbit and clock products. The PFRR station performed the good performance both in static and kinematic mode with an accuracy of 2.99 cm and 8.08 cm, respectively, whereas the CPNM station produced the worst static performance with an error of 8.09 cm, and the ANMG station produced the worst kinematic performance with a counterpart of 13.69 cm. (3) The PPP-AR solution was superior to the PPP solution, given that, with respect to PPP, post-processing PPP-AR improved the positioning accuracy and convergence time by 13–32 % (3–89 %) in GPS-only mode by 2–15 % (5–60 %) in GPS/QZSS mode. Thus, we conclude that the current MADOCA products can provide SSR corrections and FCB products with positioning accuracy at the decimeter or even centimeter level, which could meet the demands of the RTPPP and PPP-AR solutions.     

A simplified data assimilation method for reconstructing time-series MODIS NDVI data

The Normalized Difference Vegetation Index (NDVI) is an important vegetation index, widely applied in research on global environmental and climatic change. However, noise induced by cloud contamination and atmospheric variability impedes the analysis and application of NDVI data. In this work, a simplified data assimilation method is proposed to reconstruct high-quality time-series MODIS NDVI data. We extracted 16-Day L3 Global 1 km SIN Grid NDVI data sets for western China from MODIS vegetation index (VI) products (MOD13A2) for the period 2003–2006. NDVI data in the first three years (2003–2005) were used to generate the background field of NDVI based on a simple three-point smoothing technique, which captures annual features of vegetation change. NDVI data for 2006 were used to test our method. For every time step, the quality assurance (QA) flags of the MODIS VI products were adopted to empirically determine the weight between the background field and NDVI observations. Ultimately, more reliable NDVI data can be produced. The results indicate that the newly developed method is robust and effective in reconstructing high-quality MODIS NDVI time-series.     

Developing a dust storm detection method combining Support Vector Machine and satellite data in typical dust regions of Asia

Enhancing the dust storm detection is a key part for the environmental protection, human healthy and economic development. The goal of this paper is to propose a new Support Vector Machine (SVM)-based method to automatically detect dust storms using remote sensing data. Existing methods dealing with this problem are usually threshold-based that are of great complexity and uncertainty. In this paper we propose a simple and reliable method combining SVM with MODIS L1 data and explore the optimal band combinations used as the feature vectors of SVM. The developed method was evaluated by MODIS and OMI data qualitatively and quantitatively on three study sites located in the Arabian Desert, Gobi Desert and Taklimakan Desert, and it was also compared to three other traditional methods based on their accuracy, complexity, reliability and sensitivity to thresholds. The detection results demonstrated that the combination of (Band7 − Band3)/(Band7 + Band3) ((B7 − B3)/(B7 + B3)), Band20 − Band31 (B20 − B31), and Band31/Band32 (B31/B32) can detect the dust storms more precisely than other individual bands or their combination. The comparison among those cases indicated that the proposed automatic method exhibited an advantage of minimizing the uncertainty and complexity, which were the limits of defining thresholds based on the threshold-based methods. The conclusions can provide references for studies that focus on statistical-based dust storm detection.     

A comparison of different regression models for downscaling Landsat and MODIS land surface temperature images over heterogeneous landscape

Remotely sensed high spatial resolution thermal images are required for various applications in natural resource management. At present, availability of high spatial resolution (<200 m) thermal images are limited. The temporal resolution of such images is also low. Whereas, coarser spatial resolution (∼1000 m) thermal images with high revisiting capability (∼1 day) are freely available. To bridge this gap, present study attempts to downscale coarser spatial resolution thermal image to finer spatial resolution using relationships between land surface temperature (LST) and vegetation indices over a heterogeneous landscape of India. Five regression based models namely (i) Disaggregation of Radiometric Temperature (DisTrad), (ii) Temperature Sharpening (TsHARP), (iii) TsHARP with local variant, (iv) Least median square regression downscaling (LMS_DS) and (v) Pace regression downscaling (PR_DS) are applied to downscale LST of Landsat Thematic Mapper (TM) and Terra MODIS (Moderate Resolution Imaging Spectroradiometer) images. All the five models are first evaluated on Landsat image aggregated to 960 m resolution and downscaled to 480 m and 240 m resolution. The downscale accuracy is achieved using LMS_DS and PR_DS models at 240 m resolution at 0.61 °C and 0.75 °C respectively. MODIS data downscaled from 1000 m to 250 m spatial resolution results root mean square error (RMSE) of 1.43 °C and 1.62 °C for LMS_DS and PR_DS models, respectively. The LMS_DS model is less sensitive to outliers in heterogeneous landscape and provides higher accuracy when compared to other models. Downscaling model is found to be suitable for agricultural and vegetated landscapes up to a spatial resolution of 250 m but not applicable to water bodies, dry river bed sand sandy open areas.     

BeiDou-3 orbit and clock quality of the IGS Multi-GNSS Pilot Project

Within the Multi-GNSS Pilot Project (MGEX) of the International GNSS Service (IGS), precise orbit and clock products for the BeiDou-3 global navigation satellite system (BDS-3) are routinely generated by a total of five analysis centers. The processing standards and specific properties of the individual products are reviewed and the BDS-3 orbit and clock product performance is assessed through direct inter-comparison, satellite laser ranging (SLR) residuals, clock stability analysis, and precise point positioning solutions. The orbit consistency evaluated by the signal-in-space range error is on the level of 4–8 cm for the medium Earth orbit satellites whereas SLR residuals have RMS values between 3 and 9 cm. The clock analysis reveals sytematic effects related to the elevation of the Sun above the orbital plane for all ACs pointing to deficiencies in solar radiation pressure modeling. Nevertheless, precise point positioning with the BDS-3 MGEX orbit and clock products results in 3D RMS values between 7 and 8 mm.     

High-performance altimeter Doppler processing for measuring sea level height under varying sea state conditions

This paper describes an innovative method for processing nadir altimeter data acquired in Synthetic Aperture Radar (SAR) mode, enhancing the system performances over open ocean. Similarly to the current SAR data processing scheme, the so-called LR-RMC (Low Resolution with Range Migration Correction) method, originally designed by Phalippou and Demeester (2011), includes Doppler beam forming, Doppler shift correction and range correction. In LR-RMC, however, an alternative and less complex averaging (stacking) operation is used so that all the Doppler beams produced in a radar cycle (4 bursts of 64 beams for the open-burst Sentinel-3-mode altimeter) are incoherently combined to form a multi-beam echo. In that manner, contrarily to the narrow-band SAR technique, the LR-RMC processing enlarges the effective footprint to average out the effects of surface waves and particularly those from small sub-mesoscale structures (<1 km) that are known to impact SAR-mode performances. On the other hand, the number of averaged beams is as high as in current SAR-mode processing, thus providing a noise reduction at least equally good. The LR-RMC method has the added benefit of reducing the incoherent integration time with respect to the SAR-mode processing (50 ms compared to 2.5 s) limiting possible surface movement effects. By processing one year of Sentinel-3A SRAL SAR-mode data using the LR-RMC method, it is shown that the swell impact on the SAR altimeter performances is totally removed and that an improvement of 10–50% is obtained in the measurement noise of the sea surface height and significant wave height with respect to SAR mode. Additionally, observational capabilities over the middle scales are enhanced potentially allowing the ocean mesoscale features to be retrieved and observations assimilated more usefully in ocean models.     

Evaluation of MODIS GPP over a complex ecosystem in East Asia: A case study at Gwangneung flux tower in Korea

Moderate Resolution Imaging Radiometer (MODIS) gross primary productivity (GPP) has been used widely to study the global carbon cycle associated with terrestrial ecosystems. The retrieval of the current MODIS productivity with a 1 × 1 km² resolution has limitations when presenting subgrid scale processes in terrestrial ecosystems, specifically when forests are located in mountainous areas, and shows heterogeneity in vegetation type due to intensive land use. Here, we evaluate MODIS GPP (MOD17) at Gwangneung deciduous forest KoFlux tower (deciduous forest; GDK) for 2006–2010 in Korea, where the forests comprise heterogeneous vegetation cover over complex terrain. The monthly MODIS GPP data overestimated the GDK measurements in a range of +15% to +34% and was fairly well correlated (R = 0.88) with the monthly variability at GDK during the growing season. In addition, the MODIS data partly represented the sharp GPP reduction during the Asian summer monsoon (June–September) when intensive precipitation considerably reduces solar radiation and disturbs the forest ecosystem. To examine the influence of subgrid scale heterogeneity on GPP estimates over the MODIS scale, the individual vegetation type and its area within a corresponding MODIS pixel were identified using a national forest type map (∼71-m spatial resolution), and the annual GPP in the same area as the MODIS pixel was estimated. This resulted in a slight reduction in the positive MODIS bias by ∼10%, with a high degree of uncertainty in the estimation. The MODIS discrepancy for GDK suggests further investigation is necessary to determine the MODIS errors associated with the site-specific aerodynamic and hydrological characteristics that are closely related to the mountainous topography. The accuracy of meteorological variables and the impact of the very cloudy conditions in East Asia also need to be assessed.     

Spatial and temporal monitoring of drought conditions using the satellite rainfall estimates and remote sensing optical and thermal measurements

Drought is an important natural disaster that causes devastating impacts on the ecosystem, livestock, environment, and society. So far, various remote-sensing methods have been developed to estimate drought conditions, each of which has advantages and restrictions. This study aims to monitor the real-time drought indices at the field scales via the integration of various earth observations. Our proposed method consists of two steps. In the first step, the relationships between long-term standardized precipitation indices (SPI) derived from PERSIANN-CDR rainfall data and two drought-dependent parameters derived from MODIS products, including normalized NDVI and soil-air temperature gradient, are obtained at the spatial resolution of PERSIANN-CDR grid (approximately 25 km). As the next step, the corresponding relationships are applied to estimate the drought index maps at the spatial resolution of MODIS products (1 km). Numerous analyses are carried out to evaluate the proposed method. The results revealed that, from various drought indices, including SPIs of different timescales (1, 3, 6, and 12-months), SPI-3 and SPI-6 are more appropriate to the proposed method in terms of correlation with temperature and vegetation parameters. The findings also demonstrate the competency of the proposed method in estimating SPI indices with average RMSE 0.67 and the average correlation coefficient of 0.74.     

2021 Turkey mega forest Fires: Biodiversity measurements of the IUCN Red List wildlife mammals in Sentinel-2 based burned areas

Mega wildfires are one of the environmental disasters worldwide. This study evaluates the pre-fire species diversity and the indirect effects, including habitat loss for one of the largest wildfires in Manavgat (Antalya-Turkey) in 2021, with a two-step methodology. Here, (1) burnt areas in the Manavgat district (2021) were detected with remote sensing data from Sentinel-2A by delta Normalized Burn Ratio calculation for a selected area in Google Earth Engine, and (2) mammals' habitat vector data of International Union for Conservation of Nature (IUCN) Red List were integrated into Habitat and Biodiversity modelling of Terrset to analyze the alpha, beta, gamma diversity and Range Restriction Index for the wildfire region. In the total 4210 km² study area, 696 km² of the area was damaged by different fire severity; also, there were 56 mammal species' habitats here. These species include bats (i.e. Nyctalus leisleri), felids (i.e. Felis chaus), rodents (i.e. Rattus norvegicus) and large mammals (i.e. Ursus arctos). 88 % of these species are in IUCN's Least Concern category. The remaining classes are Near Threatened (3.7 %) and Vulnerable (7.4 %). This study also indicated that the burnt area's species richness does not totally consist of endemic species. Therefore, pre-fire species richness analyses of this study can be a base for further studies about the species' post-fire activity and occupancy.Furthermore, the 2021 mega wildfires show us the necessity of wildfire monitoring and risk studies in the entire Mediterranean ecosystem to evaluate the risks to the Sustainable Development Goals. Therefore, post-fire spatial data, fire migration monitorization, and recording of the species' activities should be performed temporally. In this way, the ability of wildlife's recovering, and the direct and indirect effects of the fire will be examined for ecosystems in the long term.     

Retrieval and validation of aerosol optical depth using GF-1 WFV cameras data

Conventional AOD (Aerosol Optical Depth) retrieval is restricted to the global and regional scale due to the limited spatial resolution of satellites. This does not allow for aerosol monitoring at the city level. The Chinese GF-1 Wide Field of View (WFV) sensors have sufficiently fine resolution as a data source for AOD retrieval with fine spatial resolution and a 4-day revisit time. In this study, principles similar to those in the Deep Blue (DB) and Dark Target (DT) algorithms were used to retrieve AOD at 100 m spatial resolution from GF-1 WFV images supported by Moderate Resolution Imaging Spectraradiometer (MODIS) surface reflectance (SR) products (MOD09A1). The derived GF-1 WFV AOD were compared with a combination of MOD04_3K DT AOD and MOD04_L2 DB AOD (MODIS AOD) to find that they yield reasonable Spearman correlations (R_S > 0.82) over Taiwan and Beijing. The derived GF-1 WFV AOD were also validated against Aerosol Robotic Network (AERONET) AOD; the Spearman correlation values were R_S = 0.911 in Beijing and R_S = 0.858 in Taiwan.