Anthropogenic heat (AH) is an important input for the urban thermal environment. While reduction in AH during the Coronavirus disease 2019 (COVID-19) pandemic may have weakened urban heat islands (UHI), quantitative assessments on this are lacking. Here, a new AH estimation method based on a remote sensing surface energy balance (RS-SEB) without hysteresis from heat storage was proposed to clarify the effects of COVID-19 control measures on AH. To weaken the impact of shadows, a simple and novel calibration method was developed to estimate the SEB in multiple regions and periods. To overcome the hysteresis of AH caused by heat storage, RS-SEB was combined with an inventory-based model and thermal stability analysis framework. The resulting AH was consistent with the latest global AH dataset and had a much higher spatial resolution, providing objective and refined features of human activities during the pandemic. Our study of four Chinese megacities (Wuhan, Shanghai, Beijing, and Guangzhou) indicated that COVID-19 control measures severely restricted human activities and notably reduced AH. The reduction was up to 50% in Wuhan during the lockdown in February 2020 and gradually decreased after the lockdown was eased in April 2020, similar to that in Shanghai during the Level 1 pandemic response. In contrast, AH was less reduced in Guangzhou during the same period and increased in Beijing owing to extended central heating use in winter. AH decreased more in urban centers and the change in AH varied in terms of urban land use between cities and periods. Although UHI changes during the COVID-19 pandemic cannot be entirely attributed to AH changes, the considerable reduction in AH is an important feature accompanying the weakening of the UHI.

Atmospheric pollutant data retrieved through satellite sensors are continually used to assess changes in air quality in the lower atmosphere. During the COVID-19 pandemic, several studies started to use satellite measurements to evaluate changes in air quality in many different regions worldwide. However, although satellite data is continuously validated, it is known that its accuracy may vary between monitored areas, requiring regionalized quality assessments. Thus, this study aimed to evaluate whether satellites could measure changes in the air quality of the state of São Paulo, Brazil, during the COVID-19 outbreak; and to verify the relationship between satellite-based data [Tropospheric NO column density and Aerosol Optical Depth (AOD)] and ground-based concentrations [NO and particulate material (PM; coarse: PM and fine: PM)]. For this purpose, tropospheric NO obtained from the TROPOMI sensor and AOD retrieved from MODIS sensor data by using the Multi-Angle Implementation of Atmospheric Correction (MAIAC) algorithm were compared with concentrations obtained from 50 automatic ground monitoring stations. The results showed low correlations between PM and AOD. For PM, most stations showed correlations lower than 0.2, which were not significant. The results for PM were similar, but some stations showed good correlations for specific periods (before or during the COVID-19 outbreak). Satellite-based Tropospheric NO proved to be a good predictor for NO concentrations at ground level. Considering all stations with NO measurements, correlations >0.6 were observed, reaching 0.8 for specific stations and periods. In general, it was observed that regions with a more industrialized profile had the best correlations, in contrast with rural areas. In addition, it was observed about 57% reductions in tropospheric NO throughout the state of São Paulo during the COVID-19 outbreak. Variations in air pollutants were linked to the region economic vocation, since there were reductions in industrialized areas (at least 50% of the industrialized areas showed >20% decrease in NO) and increases in areas with farming and livestock characteristics (about 70% of those areas showed increase in NO). Our results demonstrate that Tropospheric NO column densities can serve as good predictors of NO concentrations at ground level. For MAIAC-AOD, a weak relationship was observed, requiring the evaluation of other possible predictors to describe the relationship with PM. Thus, it is concluded that regionalized assessment of satellite data accuracy is essential for assertive estimates on a regional/local level. Good quality information retrieved at specific polluted areas does not assure a worldwide use of remote sensor data.

Frequent observations of surface water at fine spatial scales will provide critical data to support the management of aquatic habitat, flood risk and water quality. Sentinel-1 and Sentinel-2 satellites can provide such observations, but algorithms are still needed that perform well across diverse climate and vegetation conditions. We developed surface inundation algorithms for Sentinel-1 and Sentinel-2, respectively, at 12 sites across the conterminous United States (CONUS), covering a total of >536,000 km and representing diverse hydrologic and vegetation landscapes. Each scene in the 5-year (2017-2021) time series was classified into open water, vegetated water, and non-water at 20 m resolution using variables from Sentinel-1 and Sentinel-2, as well as variables derived from topographic and weather datasets. The Sentinel-1 algorithm was developed distinct from the Sentinel-2 model to explore if and where the two time series could potentially be integrated into a single high-frequency time series. Within each model, open water and vegetated water (vegetated palustrine, lacustrine, and riverine wetlands) classes were mapped. The models were validated using imagery from WorldView and PlanetScope. Classification accuracy for open water was high across the 5-year period, with an omission and commission error of only 3.1% and 0.9% for the Sentinel-1 algorithm and 3.1% and 0.5% for the Sentinel-2 algorithm, respectively. Vegetated water accuracy was lower, as expected given that the class represents mixed pixels. The Sentinel-2 algorithm showed higher accuracy (10.7% omission and 7.9% commission error) relative to the Sentinel-1 algorithm (28.4% omission and 16.0% commission error). Patterns over time in the proportion of area mapped as open or vegetated water by the Sentinel-1 and Sentinel-2 algorithms were charted and correlated for a subset of all 12 sites. Our results showed that the Sentinel-1 and Sentinel-2 algorithm open water time series can be integrated at all 12 sites to improve the temporal resolution, but sensor-specific differences, such as sensitivity to vegetation structure versus pixel color, complicate the data integration for mixed-pixel, vegetated water. The methods developed here provide inundation at 5-day (Sentinel-2 algorithm) and 12-day (Sentinel-1 algorithm) time steps to improve our understanding of the short- and long-term response of surface water to climate and land use drivers in different ecoregions.

Remote detection and monitoring of the vegetation responses to stress became relevant for sustainable agriculture. Ongoing developments in optical remote sensing technologies have provided tools to increase our understanding of stress-related physiological processes. Therefore, this study aimed to provide an overview of the main spectral technologies and retrieval approaches for detecting crop stress in agriculture. Firstly, we present integrated views on: i) biotic and abiotic stress factors, the phases of stress, and respective plant responses, and ii) the affected traits, appropriate spectral domains and corresponding methods for measuring traits remotely. Secondly, representative results of a systematic literature analysis are highlighted, identifying the current status and possible future trends in stress detection and monitoring. Distinct plant responses occurring under shortterm, medium-term or severe chronic stress exposure can be captured with remote sensing due to specific light interaction processes, such as absorption and scattering manifested in the reflected radiance, i.e. visible (VIS), near infrared (NIR), shortwave infrared, and emitted radiance, i.e. solar-induced fluorescence and thermal infrared (TIR). From the analysis of 96 research papers, the following trends can be observed: increasing usage of satellite and unmanned aerial vehicle data in parallel with a shift in methods from simpler parametric approaches towards more advanced physically-based and hybrid models. Most study designs were largely driven by sensor availability and practical economic reasons, leading to the common usage of VIS-NIR-TIR sensor combinations. The majority of reviewed studies compared stress proxies calculated from single-source sensor domains rather than using data in a synergistic way. We identified new ways forward as guidance for improved synergistic usage of spectral domains for stress detection: (1) combined acquisition of data from multiple sensors for analysing multiple stress responses simultaneously (holistic view); (2) simultaneous retrieval of plant traits combining multi-domain radiative transfer models and machine learning methods; (3) assimilation of estimated plant traits from distinct spectral domains into integrated crop growth models. As a future outlook, we recommend combining multiple remote sensing data streams into crop model assimilation schemes to build up Digital Twins of agroecosystems, which may provide the most efficient way to detect the diversity of environmental and biotic stresses and thus enable respective management decisions.

Cloud detection is a crucial step in the optical satellite image processing pipeline for Earth observation. Clouds in optical remote sensing images seriously affect the visibility of the background and greatly reduce the usability of images for land applications. Traditional methods based on thresholding, multi-temporal or multi-spectral information are often specific to a particular satellite sensor. Convolutional Neural Networks for cloud detection often require labeled cloud masks for training that are very time-consuming and expensive to obtain. To overcome these challenges, this paper presents a hybrid cloud detection method based on the synergistic combination of generative adversarial networks (GAN) and a physics-based cloud distortion model (CDM). The proposed weakly-supervised GAN-CDM method (available online https://github.com/Neooolee/GANCDM) only requires patch-level labels for training, and can produce cloud masks at pixel-level in both training and testing stages. GAN-CDM is trained on a new globally distributed Landsat 8 dataset (WHUL8-CDb, available online doi:https://doi.org/10.5281/zenodo.6420027) including image blocks and corresponding block-level labels. Experimental results show that the proposed GAN-CDM method trained on Landsat 8 image blocks achieves much higher cloud detection accuracy than baseline deep learning-based methods, not only in Landsat 8 images (L8 Biome dataset, 90.20% versus 72.09%) but also in Sentinel-2 images ("S2 Cloud Mask Catalogue" dataset, 92.54% versus 77.00%). This suggests that the proposed method provides accurate cloud detection in Landsat images, has good transferability to Sentinel-2 images, and can quickly be adapted for different optical satellite sensors.

The unprecedented availability of optical satellite data in cloud-based computing platforms, such as Google Earth Engine (GEE), opens new possibilities to develop crop trait retrieval models from the local to the planetary scale. Hybrid retrieval models are of interest to run in these platforms as they combine the advantages of physically- based radiative transfer models (RTM) with the flexibility of machine learning regression algorithms. Previous research with GEE primarily relied on processing bottom-of-atmosphere (BOA) reflectance data, which requires atmospheric correction. In the present study, we implemented hybrid models directly into GEE for processing Sentinel-2 (S2) Level-1C (L1C) top-of-atmosphere (TOA) reflectance data into crop traits. To achieve this, a training dataset was generated using the leaf-canopy RTM PROSAIL in combination with the atmospheric model 6SV. Gaussian process regression (GPR) retrieval models were then established for eight essential crop traits namely leaf chlorophyll content, leaf water content, leaf dry matter content, fractional vegetation cover, leaf area index (LAI), and upscaled leaf variables (i.e., canopy chlorophyll content, canopy water content and canopy dry matter content). An important pre-requisite for implementation into GEE is that the models are sufficiently light in order to facilitate efficient and fast processing. Successful reduction of the training dataset by 78% was achieved using the active learning technique Euclidean distance-based diversity (EBD). With the EBD-GPR models, highly accurate validation results of LAI and upscaled leaf variables were obtained against in situ field data from the validation study site Munich-North-Isar (MNI), with normalized root mean square errors (NRMSE) from 6% to 13%. Using an independent validation dataset of similar crop types (Italian Grosseto test site), the retrieval models showed moderate to good performances for canopy-level variables, with NRMSE ranging from 14% to 50%, but failed for the leaf-level estimates. Obtained maps over the MNI site were further compared against Sentinel-2 Level 2 Prototype Processor (SL2P) vegetation estimates generated from the ESA Sentinels' Application Platform (SNAP) Biophysical Processor, proving high consistency of both retrievals ( from 0.80 to 0.94). Finally, thanks to the seamless GEE processing capability, the TOA-based mapping was applied over the entirety of Germany at 20 m spatial resolution including information about prediction uncertainty. The obtained maps provided confidence of the developed EBD-GPR retrieval models for integration in the GEE framework and national scale mapping from S2-L1C imagery. In summary, the proposed retrieval workflow demonstrates the possibility of routine processing of S2 TOA data into crop traits maps at any place on Earth as required for operational agricultural applications.

Wildland fire smoke contains large amounts of PM that can traverse tens to hundreds of kilometers, resulting in significant deterioration of air quality and excess mortality and morbidity in downwind regions. Estimating PM levels while considering the impact of wildfire smoke has been challenging due to the lack of ground monitoring coverage near the smoke plumes. We aim to estimate total PM concentration during the Camp Fire episode, the deadliest wildland fire in California history. Our random forest (RF) model combines calibrated low-cost sensor data (PurpleAir) with regulatory monitor measurements (Air Quality System, AQS) to bolster ground observations, Geostationary Operational Environmental Satellite-16 (GOES-16)'s high temporal resolution to achieve hourly predictions, and oversampling techniques (Synthetic Minority Oversampling Technique, SMOTE) to reduce model underestimation at high PM levels. In addition, meteorological fields at 3 km resolution from the High-Resolution Rapid Refresh model and land use variables were also included in the model. Our AQS-only model achieved an out of bag (OOB) R (RMSE) of 0.84 (12.00 μg/m) and spatial and temporal cross-validation (CV) R (RMSE) of 0.74 (16.28 μg/m) and 0.73 (16.58 μg/m), respectively. Our AQS + Weighted PurpleAir Model achieved OOB R (RMSE) of 0.86 (9.52 μg/m) and spatial and temporal CV R (RMSE) of 0.75 (14.93 μg/m) and 0.79 (11.89 μg/m), respectively. Our AQS + Weighted PurpleAir + SMOTE Model achieved OOB R (RMSE) of 0.92 (10.44 μg/m) and spatial and temporal CV R (RMSE) of 0.84 (12.36 μg/m) and 0.85 (14.88 μg/m), respectively. Hourly predictions from our model may aid in epidemiological investigations of intense and acute exposure to PM during the Camp Fire episode.

Urbanization is the second largest mega-trend right after climate change. Accurate measurements of urban morphological and demographic figures are at the core of many international endeavors to address issues of urbanization, such as the United Nations' call for "Sustainable Cities and Communities". In many countries - particularly developing countries -, however, this database does not yet exist. Here, we demonstrate a novel deep learning and big data analytics approach to fuse freely available global radar and multi-spectral satellite data, acquired by the Sentinel-1 and Sentinel-2 satellites. Via this approach, we created the first-ever global and quality controlled urban local climate zones classification covering all cities across the globe with a population greater than 300,000 and made it available to the community (https://doi.org/10.14459/2021mp1633461). Statistical analysis of the data quantifies a global inequality problem: approximately 40% of the area defined as compact or light/large low-rise accommodates about 60% of the total population, whereas approximately 30% of the area defined as sparsely built accommodates only about 10% of the total population. Beyond, patterns of urban morphology were discovered from the global classification map, confirming a morphologic relationship to the geographical region and related cultural heritage. We expect the open access of our dataset to encourage research on the global change process of urbanization, as a multidisciplinary crowd of researchers will use this baseline for spatial perspective in their work. In addition, it can serve as a unique dataset for stakeholders such as the United Nations to improve their spatial assessments of urbanization.

Over the past decades, solar panels have been widely used to harvest solar energy owing to the decreased cost of silicon-based photovoltaic (PV) modules, and therefore it is essential to remotely map and monitor the presence of solar PV modules. Many studies have explored on PV module detection based on color aerial photography and manual photo interpretation. Imaging spectroscopy data are capable of providing detailed spectral information to identify the spectral features of PV, and thus potentially become a promising resource for automated and operational PV detection. However, PV detection with imaging spectroscopy data must cope with the vast spectral diversity of surface materials, which is commonly divided into spectral intra-class variability and inter-class similarity. We have developed an approach to detect PV modules based on their physical absorption and reflection characteristics using airborne imaging spectroscopy data. A large database was implemented for training and validating the approach, including spectra-goniometric measurements of PV modules and other materials, a HyMap image spectral library containing 31 materials with 5627 spectra, and HySpex imaging spectroscopy data sets covering Oldenburg, Germany. By normalizing the widely used Hydrocarbon Index (HI), we solved the intra-class variability caused by different detection angles, and validated it against the spectra-goniometric measurements. Knowing that PV modules are composed of materials with different transparencies, we used a group of spectral indices and investigated their interdependencies for PV detection with implementing the image spectral library. Finally, six well-trained spectral indices were applied to HySpex data acquired in Oldenburg, Germany, yielding an overall PV map. Four subsets were selected for validation and achieved overall accuracies, producer's accuracies and user's accuracies, respectively. This physics-based approach was validated against a large database collected from multiple platforms (laboratory measurements, airborne imaging spectroscopy data), thus providing a robust, transferable and applicable way to detect PV modules using imaging spectroscopy data. We aim to create greater awareness of the potential importance and applicability of airborne and spaceborne imaging spectroscopy data for PV modules identification.

Lakes and other surface fresh waterbodies provide drinking water, recreational and economic opportunities, food, and other critical support for humans, aquatic life, and ecosystem health. Lakes are also productive ecosystems that provide habitats and influence global cycles. Chlorophyll concentration provides a common metric of water quality, and is frequently used as a proxy for lake trophic state. Here, we document the generation and distribution of the complete MEdium Resolution Imaging Spectrometer (MERIS; Appendix A provides a complete list of abbreviations) radiometric time series for over 2300 satellite resolvable inland bodies of water across the contiguous United States (CONUS) and more than 5,000 in Alaska. This contribution greatly increases the ease of use of satellite remote sensing data for inland water quality monitoring, as well as highlights new horizons in inland water remote sensing algorithm development. We evaluate the performance of satellite remote sensing Cyanobacteria Index (CI)-based chlorophyll lgorithms, the retrievals for which provide surrogate estimates of phytoplankton concentrations in cyanobacteria dominated lakes. Our analysis quantifies the algorithms' abilities to assess lake trophic state across the CONUS. As a case study, we apply a bootstrapping approach to derive a new CI-to-chlorophyll relationship, Chl, which performs relatively well with a multiplicative bias of 1.11 (11%) and mean absolute error of 1.60 (60%). While the primary contribution of this work is the distribution of the MERIS radiometric timeseries, we provide this case study as a roadmap for future stakeholders' algorithm development activities, as well as a tool to assess the strengths and weaknesses of applying a single algorithm across CONUS.

Exposure to fine particulate matter (PM) has been linked to a substantial disease burden globally, yet little has been done to estimate the population health risks of PM in South Africa due to the lack of high-resolution PM exposure estimates. We developed a random forest model to estimate daily PM concentrations at 1 km resolution in and around industrialized Gauteng Province, South Africa, by combining satellite aerosol optical depth (AOD), meteorology, land use, and socioeconomic data. We then compared PM concentrations in the study domain before and after the implementation of the new national air quality standards. We aimed to test whether machine learning models are suitable for regions with sparse ground observations such as South Africa and which predictors played important roles in PM modeling. The cross-validation R and Root Mean Square Error of our model was 0.80 and 9.40 μg/m, respectively. Satellite AOD, seasonal indicator, total precipitation, and population were among the most important predictors. Model-estimated PM levels successfully captured the temporal pattern recorded by ground observations. Spatially, the highest annual PM concentration appeared in central and northern Gauteng, including northern Johannesburg and the city of Tshwane. Since the 2016 changes in national PM standards, PM concentrations have decreased in most of our study region, although levels in Johannesburg and its surrounding areas have remained relatively constant. This is anadvanced PM model for South Africa with high prediction accuracy at the daily level and at a relatively high spatial resolution. Our study provided a reference for predictor selection, and our results can be used for a variety of purposes, including epidemiological research, burden of disease assessments, and policy evaluation.

High resolution imaging spectrometers are prerequisite to address significant data gaps in inland optical water quality monitoring. In this work, we provide a data-driven alignment of chlorophyll- and turbidity derived from the Sentinel-2 MultiSpectral Imager (MSI) with corresponding Sentinel-3 Ocean and Land Colour Instrument (OLCI) products. For chlorophyll- retrieval, empirical 'ocean colour' blue-green band ratios and a near infra-red (NIR) band ratio algorithm, as well as a semi-analytical three-band NIR-red ratio algorithm, were included in the analysis. Six million co-registrations with MSI and OLCI spanning 24 lakes across five continents were analysed. Following atmospheric correction with POLYMER, the reflectance distributions of the red and NIR bands showed close similarity between the two sensors, whereas the distribution for blue and green bands was positively skewed in the MSI results compared to OLCI. Whilst it is not possible from this analysis to determine the accuracy of reflectance retrieved with either MSI or OLCI results, optimizing water quality algorithms for MSI against those previously derived for the Envisat Medium Resolution Imaging Spectrometer (MERIS) and its follow-on OLCI, supports the wider use of MSI for aquatic applications. Chlorophyll- algorithms were thus tuned for MSI against concurrent OLCI observations, resulting in significant improvements against the original algorithm coefficients. The mean absolute difference (MAD) for the blue-green band ratio algorithm decreased from 1.95 mg m to 1.11 mg m, whilst the correlation coefficient increased from 0.61 to 0.80. For the NIR-red band ratio algorithms improvements were modest, with the MAD decreasing from 4.68 to 4.64 mg m for the empirical red band ratio algorithm, and 3.73 to 3.67 for the semi-analytical 3-band algorithm. Three implementations of the turbidity algorithm showed improvement after tuning with the resulting distributions having reduced bias. The MAD reduced from 0.85 to 0.72, 1.22 to 1.10 and 1.93 to 1.55 FNU for the 665, 708 and 778 nm implementations respectively. However, several sources of uncertainty remain: adjacent land showed high divergence between the sensors, suggesting that high product uncertainty near land continues to be an issue for small water bodies, while it cannot be stated at this point whether MSI or OLCI results are differentially affected. The effect of spectrally wider bands of the MSI on algorithm sensitivity to chlorophyll- and turbidity cannot be fully established without further availability of in situ optical measurements.

Remote sensing-based measurements of solar-induced chlorophyll fluorescence (SIF) are useful for assessing plant functioning at different spatial and temporal scales. SIF is the most direct measure of photosynthesis and is therefore considered important to advance capacity for the monitoring of gross primary production (GPP) while it has also been suggested that its yield facilitates the early detection of vegetation stress. However, due to the influence of different confounding effects, the apparent SIF signal measured at canopy level differs from the fluorescence emitted at leaf level, which makes its physiological interpretation challenging. One of these effects is the scattering of SIF emitted from leaves on its way through the canopy. The escape fraction ( ) describes the scattering of SIF within the canopy and corresponds to the ratio of apparent SIF at canopy level to SIF at leaf level. In the present study, the fluorescence correction vegetation index (FCVI) was used to determine of far-red SIF for three structurally different crops (sugar beet, winter wheat, and fruit trees) from a diurnal data set recorded by the airborne imaging spectrometer HyPlant. This unique data set, for the first time, allowed a joint analysis of spatial and temporal dynamics of structural effects and thus the downscaling of far-red SIF from canopy ( ) to leaf level ( ). For a homogeneous crop such as winter wheat, it seems to be sufficient to determine once a day to reliably scale SIF from canopy to leaf level. In contrast, for more complex canopies such as fruit trees, calculating for each observation time throughout the day is strongly recommended. The compensation for structural effects, in combination with normalizing SIF to remove the effect of incoming radiation, further allowed the estimation of SIF emission efficiency ( ) at leaf level, a parameter directly related to the diurnal variations of plant photosynthetic efficiency.

A framework is proposed for assessing the physical consistency between two terrestrial Essential Climate Variables (ECVs) products retrieved from Earth Observation at global scale. The methodology assessed the level of agreement between the temporal variations of Leaf Area Index (LAI) and Fraction of Absorbed Photosynthetically Active Radiation (FAPAR). The simultaneous changes were classified according to their sign, magnitude and level of confidence, whereby the respective products uncertainties were taken into consideration. A set of proposed agreement metrics were used to identify temporal and spatial biases of non-coherency, non-significance, sensitivity and the overall level of agreement of the temporal changes between two ECVs. We applied the methodology using the Joint Research Center (JRC) Two-stream Inversion Package (TIP) products at 1 km, those provided by the Copernicus Global Land Service (CGLS) based on the SPOT/VGT and Proba-V at 1 km, and the MODIS MCD15A3 at 500 m. In addition, the same analysis was applied with aggregated products at a larger scale over Southern Africa. We found that the CGLS LAI and FAPAR products lacked consistency in their spatial and temporal changes and were severely affected by trends. The MCD15A3 products were characterized by the highest number of non-coherent changes between the two ECVs but temporal inconsistencies were mainly located over the eastern hemisphere. The JRC-TIP products were highly consistent. The results showed the advantages of physically-based retrieval algorithms, in both JRC-TIP and MODIS products, and indicated also that, except for MODIS over forests, aggregated products using an uncertainty-based weighted average led to higher agreement between the ECVs changes.

The early detection of () infections is critical to the management of this dangerous plan pathogen across the world. Recent studies with remote sensing (RS) sensors at different scales have shown that -infected olive trees have distinct spectral features in the visible and infrared regions (VNIR). However, further work is needed to integrate remote sensing in the management of plant disease epidemics. Here, we research how the spectral changes picked up by different sets of RS plant traits (i.e., pigments, structural or leaf protein content), can help capture the spatial dynamics of spread. We coupled a spatial spread model with the probability of infection predicted by a RS-driven support vector machine (RS-SVM) model. Furthermore, we analyzed which RS plant traits contribute most to the output of the prediction models. For that, in almond orchards affected by ( = 1426 trees), we conducted a field campaign simultaneously with an airborne campaign to collect high-resolution thermal images and hyperspectral images in the visible-near-infrared (VNIR, 400-850 nm) and short-wave infrared regions (SWIR, 950-1700 nm). The best performing RS-SVM model (OA = 75%; kappa = 0.50) included as predictors leaf protein content, nitrogen indices (NIs), fluorescence and a thermal indicator (T), alongside pigments and structural parameters. Leaf protein content together with NIs contributed 28% to the explanatory power of the model, followed by chlorophyll (22%), structural parameters (LAI and LIDF), and chlorophyll indicators of photosynthetic efficiency. Coupling the RS model with an epidemic spread model increased the accuracy (OA = 80%; kappa = 0.48). In the almond trees where the presence of was assayed by qPCR ( = 318 trees), the combined RS-spread model yielded an OA of 71% and kappa = 0.33, which is higher than the RS-only model and visual inspections (both OA = 64-65% and kappa = 0.26-31). Our work demonstrates how combining spatial epidemiological models and remote sensing can lead to highly accurate predictions of plant disease spatial distribution.

Data collected at large scale and low cost (e.g. satellite and street level imagery) have the potential to substantially improve resolution, spatial coverage, and temporal frequency of measurement of urban inequalities. Multiple types of data from different sources are often available for a given geographic area. Yet, most studies utilize a single type of input data when making measurements due to methodological difficulties in their joint use. We propose two deep learning-based methods for jointly utilizing satellite and street level imagery for measuring urban inequalities. We use London as a case study for three selected outputs, each measured in decile classes: income, overcrowding, and environmental deprivation. We compare the performances of our proposed multimodal models to corresponding unimodal ones using mean absolute error (MAE). First, satellite tiles are appended to street level imagery to enhance predictions at locations where street images are available leading to improvements in accuracy by 20, 10, and 9% in units of decile classes for income, overcrowding, and living environment. The second approach, novel to the best of our knowledge, uses a U-Net architecture to make predictions for all grid cells in a city at high spatial resolution (e.g. for 3 m × 3 m pixels in London in our experiments). It can utilize city wide availability of satellite images as well as more sparse information from street-level images where they are available leading to improvements in accuracy by 6, 10, and 11%. We also show examples of prediction maps from both approaches to visually highlight performance differences.

The frequent acquisitions of fine spatial resolution imagery (10 m) offered by recent multispectral satellite missions, including Sentinel-2, can resolve single agricultural fields and thus provide crop-specific phenology metrics, a crucial information for crop monitoring. However, effective phenology retrieval may still be hampered by significant cloud cover. Synthetic aperture radar (SAR) observations are not restricted by weather conditions, and Sentinel-1 thus ensures more frequent observations of the land surface. However, these data have not been systematically exploited for phenology retrieval so far. In this study, we extracted crop-specific land surface phenology (LSP) from Sentinel-1 and Sentinel-2 of major European crops (common and durum wheat, barley, maize, oats, rape and turnip rape, sugar beet, sunflower, and dry pulses) using ground-truth information from the "Copernicus module" of the Land Use/Cover Area frame statistical Survey (LUCAS) of 2018. We consistently used a single model-fit approach to retrieve LSP metrics on temporal profiles of CR (Cross Ratio, the ratio of the backscattering coefficient VH/VV from Sentinel-1) and NDVI (Normalized Difference Vegetation Index from Sentinel-2). Our analysis revealed that LSP retrievals from Sentinel-1 are comparable to those of Sentinel-2, particularly for winter crops. The start of season (SOS) timings, as derived from Sentinel-1 and -2, are significantly correlated (average r of 0.78 for winter and 0.46 for summer crops). The correlation is lower for end of season retrievals (EOS, r of 0.62 and 0.34). Agreement between LSP derived from Sentinel-1 and -2 varies among crop types, ranging from  = 0.89 and mean absolute error MAE = 10 days (SOS of dry pulses) to  = 0.15 and MAE = 53 days (EOS of sugar beet). Observed deviations revealed that Sentinel-1 and -2 LSP retrievals can be complementary; for example for winter crops we found that SAR detected the start of the spring growth while multispectral data is sensitive to the vegetative growth before and during winter. To test if our results correspond reasonably to in-situ data, we compared average crop-specific LSP for Germany to average phenology from ground phenological observations of 2018 gathered from the German Meteorological Service (DWD). Our study demonstrated that both Sentinel-1 and -2 can provide relevant and at times complementary LSP information at field- and crop-level.

Predicting long-term spatiotemporal characteristics of fine particulate matter (PM) is important in China to understand historical levels of PM, to support health effects research of both long-term and short-term exposures to PM, and to evaluate the efficacy of air pollution control policies. Satellite-retrieved aerosol optical depth (AOD) provides a unique opportunity to characterize the long-term trends of ground-level PM at high spatial resolution. However, the missing rate of AOD in Northeastern China (NEC) is very high, especially in winter, and challenges the accuracy of long-term predictions of PM if left unresolved. Using random forest algorithms, this study developed a gap-filling approach combing satellite AOD, meteorological data, land use parameters, population and visibility in the NEC during 2005-2016. The model, including all predictors, combined with a model without AOD was able to fill the gap of PM predictions caused by missing AOD at 1-km resolution. The R (RMSE) of the full-coverage predictions was 0.81 (18.5 μg/m) at the daily level. Gap-filled PM predictions on days with missing AOD reduced the relative prediction error from 28% to 2.5% in winter. The leave-one-year-out-cross-validation R (RMSE) of the full-coverage predictions was 0.65 (16.3 μg/m) at the monthly level, indicating relatively high accuracy of predicted historical PM concentrations. Our results suggested that AOD helped increase the reliability of historical PM prediction when ground PM measurements were unavailable, even though predictions from the AOD model only accounted for approximate 37% of the whole dataset. Predicted PM level in NEC have increased since 2005, reached its peak during 2013-2015, then saw a major decline in 2016. Our high-resolution predictions also showed a south to north gradient and many pollution hot spots in the city clusters surrounding provincial capitals, as well as within large cities. Overall, by combining predictions from the AOD model with higher accuracy and predictions from the non-AOD model to achieve full coverage, our modeling approach could produce long-term, full-coverage historical PM levels in high-latitude areas in China, despite the widespread and persistent AOD missingness.

Urban areas and their vertical characteristics have a manifold and far-reaching impact on our environment. However, openly accessible information at high spatial resolution is still missing at large for complete countries or regions. In this study, we combined Sentinel-1A/B and Sentinel-2A/B time series to map building heights for entire Germany on a 10 m grid resolving built-up structures in rural and urban contexts. We utilized information from the spectral/polarization, temporal and spatial dimensions by combining band-wise temporal aggregation statistics with morphological metrics. We trained machine learning regression models with highly accurate building height information from several 3D building models. The novelty of this method lies in the very fine resolution yet large spatial extent to which it can be applied, as well as in the use of building shadows in optical imagery. Results indicate that both radar-only and optical-only models can be used to predict building height, but the synergistic combination of both data sources leads to superior results. When testing the model against independent datasets, very consistent performance was achieved (frequency-weighted RMSE of 2.9 m to 3.5 m), which suggests that the prediction of the most frequently occurring buildings was robust. The average building height varies considerably across Germany with lower buildings in Eastern and South-Eastern Germany and taller ones along the highly urbanized areas in Western Germany. We emphasize the straightforward applicability of this approach on the national scale. It mostly relies on freely available satellite imagery and open source software, which potentially permit frequent update cycles and cost-effective mapping that may be relevant for a plethora of different applications, e.g. physical analysis of structural features or mapping society's resource usage.

Land surface temperature (LST) is a key diagnostic indicator of agricultural water use and crop stress. LST data retrieved from thermal infrared (TIR) band imagery, however, tend to have a coarser spatial resolution (e.g., 100 m for Landsat 8) than surface reflectance (SR) data collected from shortwave bands on the same instrument (e.g., 30 m for Landsat). Spatial sharpening of LST data using the higher resolution multi-band SR data provides an important path for improved agricultural monitoring at sub-field scales. A previously developed Data Mining Sharpener (DMS) approach has shown great potential in the sharpening of Landsat LST using Landsat SR data co-collected over various landscapes. This work evaluates DMS performance for sharpening ECOsystem Spaceborne Thermal Radiometer Experiment on Space Station (ECOSTRESS) LST (~70 m native resolution) and Visible Infrared Imaging Radiometer Suite (VIIRS) LST (375 m) data using Harmonized Landsat and Sentinel-2 (HLS) SR data, providing the basis for generating 30-m LST data at a higher temporal frequency than afforded by Landsat alone. To account for the misalignment between ECOSTRESS/VIIRS and Landsat/HLS caused by errors in registration and orthorectification, we propose a modified version of the DMS approach that employs a relaxed box size for energy conservation (EC). Sharpening experiments were conducted over three study sites in California, and results were evaluated visually and quantitatively against LST data from unmanned aerial vehicles (UAV) flights and from Landsat 8. Over the three sites, the modified DMS technique showed improved sharpening accuracy over the standard DMS for both ECOSTRESS and VIIRS, suggesting the effectiveness of relaxing EC box in relieving misalignment-induced errors. To achieve reasonable accuracy while minimizing loss of spatial detail due to the EC box size increase, an optimal EC box size of 180-270 m was identified for ECOSTRESS and about 780 m for VIIRS data based on experiments from the three sites. Results from this work will facilitate the development of a prototype system that generates high spatiotemporal resolution LST products for improved agricultural water use monitoring by synthesizing multi-source remote sensing data.

Satellite remote sensing offers an effective remedy to challenges in ground-based and aerial mapping that have previously impeded quantitative assessments of global seagrass extent. Commercial satellite platforms offer fine spatial resolution, an important consideration in patchy seagrass ecosystems. Currently, no consistent protocol exists for image processing of commercial data, limiting reproducibility and comparison across space and time. Additionally, the radiometric performance of commercial satellite sensors has not been assessed against the dark and variable targets characteristic of coastal waters. This study compared data products derived from two commercial satellites: DigitalGlobe's WorldView-2 and Planet's RapidEye. A single scene from each platform was obtained at St. Joseph Bay in Florida, USA, corresponding to a November 2010 field campaign. A reproducible processing regime was developed to transform imagery from basic products, as delivered from each company, into analysis-ready data usable for various scientific applications. Satellite-derived surface reflectances were compared against field measurements. WorldView-2 imagery exhibited high disagreement in the coastal blue and blue spectral bands, chronically overpredicting. RapidEye exhibited better agreement than WorldView-2, but overpredicted slightly across all spectral bands. A deep convolutional neural network was used to classify imagery into deep water, land, submerged sand, seagrass, and intertidal classes. Classification results were compared to seagrass maps derived from photointerpreted aerial imagery. This study offers the first radiometric assessment of WorldView-2 and RapidEye over a coastal system, revealing inherent calibration issues in shorter wavelengths of WorldView-2. Both platforms demonstrated as much as 97% agreement with aerial estimates, despite differing resolutions. Thus, calibration issues in WorldView-2 did not appear to interfere with classification accuracy, but could be problematic if estimating biomass. The image processing routine developed here offers a reproducible workflow for WorldView-2 and RapidEye imagery, which was tested in two additional coastal systems. This approach may become platform independent as more sensors become available.