Parameterization of air temperature in high temporal and spatial resolution from a combination of the SEVIRI and MODIS instruments

While data like HJ-1 CCD images have advantageous spatial characteristics for describing crop properties, the temporal resolution of the data is rather low, which can be easily made worse by cloud contamination. In contrast, although Moderate Resolution Imaging Spectroradiometer (MODIS) can only achieve a spatial resolution of 250 m in its normalised difference vegetation index (NDVI) product, it has a high temporal resolution, covering the Earth up to multiple times per day. To combine the high spatial resolution and high temporal resolution of different data sources, a new method (Spatial and Temporal Adaptive Vegetation index Fusion Model [STAVFM]) for blending NDVI of different spatial and temporal resolutions to produce high spatial–temporal resolution NDVI datasets was developed based on Spatial and Temporal Adaptive Reflectance Fusion Model (STARFM). STAVFM defines a time window according to the temporal variation of crops, takes crop phenophase into consideration and improves the temporal weighting algorithm. The result showed that the new method can combine the temporal information of MODIS NDVI and spatial difference information of HJ-1 CCD NDVI to generate an NDVI dataset with both high spatial and high temporal resolution. An application of the generated NDVI dataset in crop biomass estimation was provided. An average absolute error of 17.2% was achieved. The estimated winter wheat biomass correlated well with observed biomass (R 2 of 0.876). We conclude that the new dataset will improve the application of crop biomass estimation by describing the crop biomass accumulation in detail. There is potential to apply the approach in many other studies, including crop production estimation, crop growth monitoring and agricultural ecosystem carbon cycle research, which will contribute to the implementation of Digital Earth by describing land surface processes in detail.

Satellite data holds considerable potential as a source of information on rice crop growth which can be used to inform agronomy. However, given the typical field sizes in many rice-growing countries such as China, data from coarse spatial resolution satellite systems such as the Moderate Resolution Imaging Spectroradiometer (MODIS) are inadequate for resolving crop growth variability at the field scale. Nevertheless, systems such as MODIS do provide images with sufficient frequency to be able to capture the detail of rice crop growth trajectories throughout a growing season. In order to generate high spatial and temporal resolution data suitable for mapping rice crop phenology, this study fused MODIS data with lower frequency, higher spatial resolution Landsat data. An overall workflow was developed which began with image preprocessing, calculation of multi-temporal normalized difference vegetation index (NDVI) images, and spatiotemporal fusion of data from the two sensors. The Spatial and Temporal Adaptive Reflectance Fusion Model was used to effectively downscale the MODIS data to deliver a time-series of 30 m spatial resolution NDVI data at 8-day intervals throughout the rice-growing season. Zonal statistical analysis was used to extract NDVI time-series for individual fields and signal filtering was applied to the time-series to generate rice phenology curves. The downscaled MODIS NDVI products were able to characterize the development of paddy rice at fine spatial and temporal resolutions, across wide spatial extents over multiple growing seasons. These data permitted the extraction of key crop seasonality parameters that quantified inter-annual growth variability for a whole agricultural region and enabled mapping of the variability in crop performance between and within fields. Hence, this approach can provide rice crop growth data that is suitable for informing agronomic policy and practice across a wide range of scales.

The occurrence of cloud directly affects the spatial and temporal continuity of surface temperature inverted by satellite remote sensing. In this study, a framework for reconstructing surface temperature with high spatial and temporal resolution based on data assimilation is constructed on the basis of multiple subsurface validation. (1) The accuracy of MYD11A1 LST data varies with terrain and land cover characteristics. High-altitude alpine terrains (Kalasai and Arou) and undulating desert terrains (Tazhong A-E Sites) show high precision and less error, while agricultural fields (Daman) and desert transitional zones (Bajitan) exhibit more variability and larger errors. This suggests that the uniformity and stability of certain terrains, coupled with minimal atmospheric interference, enhance the accuracy of remote sensing observations. (2) A systematic bias, indicating a consistent underestimation of LST by the MYD11A1 product compared to ground-based observations, is observed across all sites. This bias is particularly pronounced in the presence of a sanding phenomenon, which results in a mixture of sand and air near the surface, leading to a lower station observation and a significant bias. (3) The Land Surface Temperature (LST) simulated by noah-MP exhibits a high degree of consistency with the LST observed through remote sensing. The significant correlation between the simulated LST and MODIS observations at the Kalsai and Arou stations indicates that noah-MP is highly applicable to mountain grassland surfaces. (4) A framework has been developed for reconstructing surface temperature with high temporal and spatial resolution, based on data assimilation. This method can generate all-weather, hourly surface temperature data.

Land Surface Temperature (LST) and Surface Net Radiation (SNR) are vital inputs for many land surface and hydrological models. However, current remote sensing datasets of these variables come mostly at coarse resolutions. Although high-resolution LST and SNR retrievals are available, they have large gaps due to cloud-cover that hinder their use as input in models. Here, we present a downscaled and continuous daily LST and SNR product across Europe for 2018–2019. The LST product is based on all-sky LST retrievals from the Spinning Enhanced Visible and InfraRed Imager (SEVIRI) onboard the geostationary Meteosat Second Generation (MSG) satellite, and clear-sky LST retrievals from the Sea and Land Surface Temperature Radiometer (SLSTR) onboard the polar-orbiting Sentinel 3 satellites. The product combines the medium spatial (approx. 5–7 km) but high temporal (30 minute) resolution, gap-free data from MSG, with the low temporal (2–3 days) but high spatial (1 km) resolution of the Sentinel 3 LST retrievals. The resulting 1 km and daily LST dataset is based on an hourly merging of both datasets through bias-correction and Kalman Filter assimilation. Longwave outgoing radiation is computed from the merged LST product in combination with MSG-based emissivity data. Shortwave outgoing radiation is computed from the incoming shortwave radiation from MSG and downscaled albedos using 1 km PROBA-V data. MSG incoming shortwave and longwave radiation and the outgoing radiation components at 1 km spatial resolution are used together to compute the final daily SNR dataset in a consistent manner. Validation results indicate an improvement of the root mean squared error by ca. 8 % with a substantial increase in spatial detail compared to the original MSG product. The resulting pan-European LST and SNR dataset can be used for hydrological modelling and as input to models dedicated to estimating evaporation and surface turbulent heat fluxes and will be regularly updated in the future.

<strong class="journal-contentHeaderColor">Abstract.</strong> Land Surface Temperature (LST) and Surface Net Radiation (SNR) are vital inputs for many land surface and hydrological models. However, current remote sensing datasets of these variables come mostly at coarse resolutions. Although high-resolution LST and SNR retrievals are available, they have large gaps due to cloud-cover that hinder their use as input in models. Here, we present a downscaled and continuous daily LST and SNR product across Europe for 2018&ndash;2019. The LST product is based on all-sky LST retrievals from the Spinning Enhanced Visible and InfraRed Imager (SEVIRI) onboard the geostationary Meteosat Second Generation (MSG) satellite, and clear-sky LST retrievals from the Sea and Land Surface Temperature Radiometer (SLSTR) onboard the polar-orbiting Sentinel 3 satellites. The product combines the medium spatial (approx. 5&ndash;7 km) but high temporal (30 minute) resolution, gap-free data from MSG, with the low temporal (2&ndash;3 days) but high spatial (1 km) resolution of the Sentinel 3 LST retrievals. The resulting 1 km and daily LST dataset is based on an hourly merging of both datasets through bias-correction and Kalman Filter assimilation. Longwave outgoing radiation is computed from the merged LST product in combination with MSG-based emissivity data. Shortwave outgoing radiation is computed from the incoming shortwave radiation from MSG and downscaled albedos using 1 km PROBA-V data. MSG incoming shortwave and longwave radiation and the outgoing radiation components at 1 km spatial resolution are used together to compute the final daily SNR dataset in a consistent manner. Validation results indicate an improvement of the root mean squared error by ca. 8 % with a substantial increase in spatial detail compared to the original MSG product. The resulting pan-European LST and SNR dataset can be used for hydrological modelling and as input to models dedicated to estimating evaporation and surface turbulent heat fluxes and will be regularly updated in the future.

The impacts of global change, including extreme heat and water scarcity, are threatening an ever-growing urban world population. Evapotranspiration (ET) mitigates the urban heat island, reducing the effect of heat waves. It can also be used as a proxy for vegetation water use, making it a crucial tool to plan resilient green cities. To optimize the trade-off between urban greening and water security, reliable and up-to-date maps of ET for cities are urgently needed. Despite its importance, few studies have mapped urban ET accurately for an entire city in high spatial and temporal resolution. We mapped the ET of Berlin, Germany in high spatial (10-m) and temporal (hourly) resolution for the year of 2019. A novel machine learning (ML) approach combining Sentinel-2 time series, open geodata, and flux footprint modeling was applied. Two eddy flux towers with contrasting surrounding land cover provided the training and testing data. Flux footprint modeling allowed us to incorporate comprehensive land cover types in training the ML models. Open remote sensing and geodata used as model inputs included Normalized Difference Vegetation Index (NDVI) from Sentinel-2, building height, impervious surface fraction, vegetation fraction, and vegetation height. NDVI was used to indicate vegetation phenology and health, as plant transpiration contributes to the majority of terrestrial ET. Hourly reference ET (RET) was calculated and used as input to capture the temporal dynamics of the meteorological conditions. Predictions were carried out using Random Forest (RF) regression. Weighted averages extracted from hourly ET maps using flux footprints were compared to measured ET from the two flux towers. Validation showed that the approach is reliable for mapping urban ET, with a mean R2 of 0.76 and 0.56 and a mean RMSE of 0.0289 mm and 0.0171 mm at the more vegetated site and the city-center site, respectively. Lastly, the variation of ET between Local Climate Zones (LCZs) was analyzed to support urban planning. This study demonstrated the capacity to map urban ET at an unprecedented high spatial and temporal resolution with a novel methodology, which can be used to support the sustainable management of green infrastructure and water resources in an urbanizing world facing climate change.

The tradeoff between spatial and temporal resolution in satellite observations substantially restrains the potential applications of Land Surface Temperature (LST) products. So far, many spatiotemporal fusion models have been developed to address the issue and a unified comparison in LST data fusion is still required. In this paper, four popular spatiotemporal fusion algorithms including Spatial and Temporal Adaptive Reflectance Fusion Model (STARFM), Unmixing-based data fusion method, Flexible Spatiotemporal Data Fusion (FSDAF) and Spatio-Temporal Integrated Temperature Fusion Model (STITFM) were adopted to generate high spatial resolution LST using polar-orbiting and geostationary satellite data. The predicted LST was evaluated by the actual LST product and the result indicates that the overall accuracy of FSDAF is satisfied (about 2.87K) and the FSDAF algorithm is recommended to generate LSTs at high spatial and temporal resolution in heterogeneous area.

Due to the limitations of hardware technology, satellite sensors cannot obtain images with high temporal, spatial, and spectral resolutions at the same time. Current spatiotemporal fusion methods try to solve the contradiction between temporal resolution and spatial resolution, which cannot achieve good reconstruction accuracy partly because the data sources are from heterogeneous platforms with long chains difficult to be modeled. Different from the crossing-platform fusion, this work proposes to improve the spatial and temporal resolutions on a single platform. For the 2-m panchromatic images, 8-m multispectral images, and 16-m wide-field-view images captured by the Gaofen-1 satellite, our goal is to produce 2-m multispectral images with high temporal resolutions. Two convolutional neural networks are built to solve this spatiotemporal-spectral fusion issue with pansharpening and spatiotemporal fusion in serial. In the validation stage, the 2-m multispectral images are built and evaluated with the panchromatic images and 8-m multispectral images. The digital and visual evaluations show that our method can produce visually acceptable fusion quality, which may enhance the feasibility of the Gaofen-1 data.

Progress in modern neuroscience critically depends on our ability to observe the activity of large neuronal populations with cellular spatial and high temporal resolution. However, two bottlenecks constrain efforts towards fast imaging of large populations. First, the resulting large video data is challenging to analyze. Second, there is an explicit tradeoff between imaging speed, signal-to-noise, and field of view: with current recording technology we cannot image very large neuronal populations with simultaneously high spatial and temporal resolution. Here we describe multi-scale approaches for alleviating both of these bottlenecks. First, we show that spatial and temporal decimation techniques based on simple local averaging provide order-of-magnitude speedups in spatiotemporally demixing calcium video data into estimates of single-cell neural activity. Second, once the shapes of individual neurons have been identified at fine scale (e.g., after an initial phase of conventional imaging with standard temporal and spatial resolution), we find that the spatial/temporal resolution tradeoff shifts dramatically: after demixing we can accurately recover denoised fluorescence traces and deconvolved neural activity of each individual neuron from coarse scale data that has been spatially decimated by an order of magnitude. This offers a cheap method for compressing this large video data, and also implies that it is possible to either speed up imaging significantly, or to “zoom out” by a corresponding factor to image order-of-magnitude larger neuronal populations with minimal loss in accuracy or temporal resolution.

Land Surface Temperature (LST) is an important resource for a variety of tasks. The data are mostly free of charge and combine high spatial and temporal resolution with reliable data collection over a historical timeframe. When remote sensing is used to provide LST data, such as the MODA11 product using information from the MODIS sensors attached to NASA satellites, data acquisition can be hindered by clouds or cloud shadows, occluding the sensors’ view on different areas of the world. This makes it difficult to take full advantage of the high resolution of the data. A common solution to interpolating LST data is statistical interpolation methods, such as fitting polynomials or thin plate spine interpolation. These methods have difficulties in incorporating additional knowledge about the research area and learning local dependencies that can help with the interpolation process. We propose a novel approach to interpolating remote sensing LST data in a fixed research area considering local ground-site air temperature measurements. The two-step approach consists of learning the LST from air temperature measurements, where the ground-site weather stations are located, and interpolating the remaining missing values with partial convolutions within a U-Net deep learning architecture. Our approach improves the interpolation of LST for our research area by 44% in terms of RMSE, when compared to state-of-the-art statistical methods. Due to the use of air temperature, we can provide coverage of 100%, even when no valid LST measurements were available. The resulting gapless coverage of high resolution LST data will help unlock the full potential of remote sensing LST data.

Downscaling AVHRR land surface temperatures for improved surface urban heat island intensity estimation

Assessment of multi-resolution and multi-sensor data for urban surface temperature retrieval

Land surface temperature (LST) images retrieved from the thermal infrared (TIR) band data of Moderate Resolution Imaging Spectroradiometer (MODIS) have much lower spatial resolution than the MODIS visible and near-infrared (VNIR) band data. The coarse pixel scale of MODIS LST images (1000 m under nadir) have limited their capability in applying to many studies required high spatial resolution in comparison of the MODIS VNIR band data with pixel scale of 250–500 m. In this paper we intend to develop an efficient approach for pixel decomposition to increase the spatial resolution of MODIS LST image using the VNIR band data as assistance. The unique feature of this approach is to maintain the thermal radiance of parent pixels in the MODIS LST image unchanged after they are decomposed into the sub-pixels in the resulted image. There are two important steps in the decomposition: initial temperature estimation and final temperature determination. Therefore the approach can be termed double-step pixel decomposition (DSPD). Both steps involve a series of procedures to achieve the final result of decomposed LST image, including classification of the surface patterns, establishment of LST change with normalized difference of vegetation index (NDVI) and building index (NDBI), reversion of LST into thermal radiance through Planck equation, and computation of weights for the sub-pixels of the resulted image. Since the Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) with much higher spatial resolution than MODIS data was on-board the same platform (Terra) as MODIS for Earth observation, an experiment had been done in the study to validate the accuracy and efficiency of our approach for pixel decomposition. The ASTER LST image was used as the reference to compare with the decomposed LST image. The result showed that the spatial distribution of the decomposed LST image was very similar to that of the ASTER LST image with a root mean square error (RMSE) of 2.7 K for entire image. Comparison with the evaluation DisTrad (E-DisTrad) and re-sampling methods for pixel decomposition also indicate that our DSPD has the lowest RMSE in all cases, including urban region, water bodies, and natural terrain. The obvious increase in spatial resolution remarkably uplifts the capability of the coarse MODIS LST images in highlighting the details of LST variation. Therefore it can be concluded that, in spite of complicated procedures, the proposed DSPD approach provides an alternative to improve the spatial resolution of MODIS LST image hence expand its applicability to the real world.

Real-time imaging of optic nerve head collagen microstructure and biomechanics using instant polarized light microscopy

&amp;lt;p&amp;gt;As the Arctic&amp;amp;#8217;s sea ice extent continues to decline, remote sensing observations are becoming even more vital for the monitoring and understanding of this process.&amp;amp;#160; Recently, the sea ice community has entered a new era of synthetic aperture radar (SAR) satellites operating at C-band with the launch of Sentinel-1A in 2014, Sentinel-1B in 2016 and the RADARSAT Constellation Mission (RCM) in 2019. These missions represent a collection of 5 spaceborne SAR sensors that together can routinely cover Arctic sea ice with a high spatial resolution (20-90 m) but also with a high temporal resolution (1-7 days) typically associated with passive microwave sensors. Here, we used ~28,000 SAR image pairs from Sentinel-1AB together with ~15,000 SAR images pairs from RCM to generate high spatiotemporal large-scale sea ice motion products across the pan-Arctic domain for 2020. The combined Sentinel-1AB and RCM sea ice motion product provides almost complete 7-day coverage over the entire pan-Arctic domain that also includes the pole-hole. Compared to the National Snow and Ice Data Center (NSIDC) Polar Pathfinder and Ocean and Sea Ice-Satellite Application Facility (OSI-SAF) sea ice motion products, ice speed was found to be faster with the Senintel-1AB and RCM product which is attributed to the higher spatial resolution of SAR imagery. More sea ice motion vectors were detected from the Sentinel-1AB and RCM product in during the summer months and within the narrow channels and inlets compared to the NSIDC Polar Pathfinder and OSI-SAF sea ice motion products. Overall, our results demonstrate that sea ice geophysical variables across the pan-Arctic domain can now be retrieved from multi-sensor SAR images at both high spatial and temporal resolution.&amp;lt;/p&amp;gt;