Clear-sky Scenes Research Articles

Improved convective cloud differential (CCD) tropospheric ozone from S5P-TROPOMI satellite data using local cloud fields

Abstract. We present the CHORA (Cloud Height Ozone Reference Algorithm) for retrieving tropospheric-ozone columns from S5P-TROPOMI (Sentinel-5 Precursor–TROPOspheric Monitoring Instrument). The method uses a local-cloud reference sector (CLC – CHORA Local Cloud) to determine the stratospheric (above-cloud) column, which is subtracted from the total column in clear-sky scenes in the same zonal band to retrieve the tropospheric column. The standard CCD (convective cloud differential) approach uses cloud data from the Pacific region (CPC – CHORA Pacific Cloud) instead. An important assumption for the standard method is the zonal invariance of stratospheric ozone. The local-cloud approach is the first step to diminish this constraint in order to extend the CCD method to mid-latitudes, where stratospheric-ozone variability is larger. An iterative approach has been developed for the automatic selection of an optimal local-cloud reference sector around each retrieval grid box varying latitudinally by ± 1° and longitudinally between ± 5 and ± 50°. The optimised CLCT (CHORA Local Cloud Theil–Sen) algorithm, a follow-up from the CLC, employs a homogeneity criterion for total ozone from the cloud reference sector in order to overcome the inhomogeneities in stratospheric ozone. It directly estimates the above-cloud column ozone for a common reference altitude of 270 hPa using the Theil–Sen regression. The latter allows for the combination of the CCD method with the cloud-slicing algorithm that retrieves upper-tropospheric ozone volume mixing ratios. Monthly averaged tropospheric-column ozone (TCO) using the Pacific cloud reference sector (CPC) and the local-cloud reference sector (CLC, CLCT) has been determined over the tropics and subtropics (26° S–22° N) using TROPOMI for the time period from 2018 to 2022. The accuracy of the various methods was investigated by means of comparisons with spatially collocated NASA/GSFC SHADOZ (Southern Hemisphere Additional Ozonesondes) measurements and the ESA TROPOMI level-2 tropospheric-ozone product. At eight out of nine tropical stations, tropospheric-ozone columns using the CLCT yield better agreement with ozonesondes than the CPC. In the tropical region (20° S–20° N), the CLCT shows a significantly lower overall mean bias and dispersion of 1 ± 7 %, outperforming both the CPC (12 ± 10 %) and CCD-ESA (22 ± 10 %). The CLCT surpasses the ESA operational product, providing more accurate tropospheric-ozone retrievals at eight out of nine stations in the tropics. For the Hilo station, with a larger stratospheric-ozone variability due to its proximity to the subtropics, the bias of +30 % (CPC) is effectively reduced to −5 % (CLCT). Similarly, in the subtropics (Reunion, Irene, Hanoi, and King's Park), the CLCT algorithm provides an overall bias and scatter of −11 ± 9 % with respect to sondes. The CLCT effectively reduces the impact of stratospheric-ozone inhomogeneity, typically at higher latitudes. These results demonstrate the advantage of the local-cloud reference sector in the subtropics. The algorithm is therefore an important basis for subsequent systematic applications in current and future missions of geostationary satellites, like GEMS (Geostationary Environment Monitoring Spectrometer, Korea), ESA Sentinel-4, and NASA TEMPO (Tropospheric Emissions: Monitoring of POllution), predominantly covering the middle latitudes.

Tropospheric NO2 retrieval algorithm for geostationary satellite instruments: applications to GEMS

Abstract. In this study, we develop an advanced retrieval algorithm for tropospheric nitrogen dioxide (NO2) from the geostationary satellite instruments and apply it to Geostationary Environment Monitoring Spectrometer (GEMS) observations. Overall, the algorithm follows previous heritage for the polar-orbiting satellites Global Ozone Monitoring Experiment-2 (GOME-2) and Tropospheric Monitoring Instrument (TROPOMI), but several improvements are implemented to account for specific features of geostationary satellites. The DLR GEMS NO2 retrieval employs an extended fitting window compared to the current fitting window used in GEMS operational v2.0 NO2 retrieval, which results in improved spectral fit quality and lower uncertainties. For the stratosphere–troposphere separation in GEMS measurements, two methods are developed and evaluated: (1) STRatospheric Estimation Algorithm from Mainz (STREAM) as used in the DLR TROPOMI NO2 retrieval and adapted to GEMS and (2) estimation of stratospheric NO2 columns from the Copernicus Atmosphere Monitoring Service (CAMS) Integrated Forecast System (IFS) cycle 48R1 model data, which introduce full stratospheric chemistry as it will be used in the operational Sentinel-4 NO2 retrieval. While STREAM provides hourly estimates of stratospheric NO2, it has limitations in describing small-scale variations and exhibits systematic biases near the boundary of the field of view. In this respect, the use of estimated stratospheric NO2 columns from the CAMS forecast model profile demonstrates better applicability by describing not only diurnal variation but also small-scale variations. For the improved air mass factor (AMF) calculation, sensitivity tests are performed using different input data. In our algorithm, cloud fractions retrieved from the Optical Cloud Recognition Algorithm (OCRA) adapted to GEMS level 1 data are applied instead of the GEMS v2.0 cloud fraction. OCRA is used operationally in TROPOMI and Sentinel-4. Compared to the GEMS level 2 cloud fraction which is typically set to around 0.1 for clear-sky scenes, OCRA sets cloud fractions close to or at 0. The OCRA-based cloud corrections result in increased tropospheric AMFs and decreased tropospheric NO2 vertical columns, leading to better agreement with results from existing TROPOMI observations. The effects of surface albedo on GEMS tropospheric NO2 retrievals are assessed by comparing the GEMS v2.0 background surface reflectance (BSR) and TROPOMI Lambertian-equivalent reflectivity (LER) climatology v2.0 product. The differences between the two surface albedo products and their impact on tropospheric AMF are particularly pronounced over snow/ice scenes during winter. A priori NO2 profiles from the CAMS forecast model, applied in the DLR GEMS algorithm, effectively capture variations in NO2 concentrations throughout the day with high spatial resolution and the advanced chemical mechanism, which demonstrates its suitability for geostationary satellite measurements. The retrieved DLR GEMS tropospheric NO2 columns show good capability for capturing hotspot signals at the scale of city clusters and describe spatial gradients from city centres to surrounding areas. Diurnal variations of tropospheric NO2 columns over Asia are well described through hourly sampling of GEMS. Evaluation of DLR GEMS tropospheric NO2 columns against TROPOMI v2.4 and GEMS v2.0 operational products shows overall good agreement. The uncertainty of DLR GEMS tropospheric NO2 vertical columns varies based on observation scenarios. In regions with low pollution levels such as open-ocean and remote rural areas, retrieval uncertainties typically range from 10 % to 50 %, primarily due to uncertainties in slant columns. For heavily polluted regions, uncertainties in tropospheric NO2 columns are mainly driven by errors in tropospheric AMF calculations. Notably, the total uncertainty in GEMS tropospheric NO2 columns is most significant in winter, particularly over heavily polluted regions with low-level clouds below or near the NO2 peak.

Generating daily 100 m resolution land surface temperature estimates continentally using an unbiased spatiotemporal fusion approach

Fine spatial resolution (i.e., ≤ 100 m) land surface temperature (LST) data are crucial to study heterogeneous landscapes (e.g., agricultural and urban). Some well-known spatiotemporal fusion methods like the Spatial and Temporal Adaptive Reflectance Fusion Model (STARFM) and the Enhanced STARFM (ESTARFM), which were originally developed to fuse surface reflectance data, may not be suitable for direct application in LST studies due to the high sub-diurnal dynamics of LST. Furthermore, the effectiveness of spatiotemporal fusion methods for LST data has not been thoroughly evaluated in previous studies that only focused on relatively small spatiotemporal extents. To address these limitations, we proposed a variant of ESTARFM, referred to as the unbiased ESTARFM (ubESTARFM), specifically designed to accommodate the high temporal dynamics of LST to generate fine-resolution LST estimates. We evaluated ubESTARFM and ESTARFM against in-situ LST and the ECOsystem Spaceborne Thermal Radiometer Experiment on Space Station (ECOSTRESS) LST across 12 regions throughout Australia, encompassing various land covers and environments. Independent validation showed that ubESTARFM had a bias of 2.55 K, unbiased root mean squared error (ubRMSE) of 2.57 K, and Pearson correlation coefficient (R) of 0.95 against the in-situ LST over 11,290 observations at the 12 sites, all of which were considerably better than those calculated for ESTARFM, being a bias of 4.73 K, ubRMSE of 3.80 K and R of 0.92. When compared to ECOSTRESS data, ubESTARFM LST had a bias of −1.69 K, ubRMSE of 2.00 K, and R of 0.70 over 43 near clear-sky scenes, while ESTARFM LST had a bias of 1.79 K, ubRMSE of 2.68 K, and R of 0.59. Overall, our results demonstrated that ubESTARFM can avoid systematic bias accumulation, substantially reduce uncertainty deviation, and maintain a good level of correlation with validation datasets when compared to ESTARFM. A further assessment underscored the potential of ubESTARFM for application using LST data acquired from geostationary platforms (e.g., Himawari-8), with a mean ubRMSE (R) of 2.22 K (0.97) against in-situ LST over 1327 observations at 3 sites from southeast Australia at the overpass time of MODIS/Terra. This promising method leverages reliable numeric values from coarse-resolution LST while borrowing spatial heterogeneity from fine-resolution LST and has the potential to be coupled with energy balance and/or radiative transfer models thus enabling better farm and/or regional-scale water management strategies to be implemented. Furthermore, both the input and generated LST data, encompassing a comprehensive spatial extent over diverse land covers and climatic conditions, are publicly available for benchmarking future algorithmic refinements.

The Far-Infrared Radiation Mobile Observation System (FIRMOS) for spectral characterization of the atmospheric emission

Abstract. The Far-Infrared Radiation Mobile Observation System (FIRMOS) is a Fourier transform spectroradiometer developed to support the Far-infrared Outgoing Radiation Understanding and Monitoring (FORUM) satellite mission by validating measurement methods and instrument design concepts, both in the laboratory and in field campaigns. FIRMOS is capable of measuring the downwelling spectral radiance emitted by the atmosphere in the spectral band from 100 to 1000 cm−1 (10–100 µm in wavelength), with a maximum spectral resolution of 0.25 cm−1. We describe the instrument design and its characterization and discuss the geophysical products obtained by inverting the atmospheric spectral radiance measured during a campaign from the high-altitude location of Mount Zugspitze in Germany, beside the Extended-range Atmospheric Emitted Radiance Interferometer (E-AERI), which is permanently installed at the site. Following the selection of clear-sky scenes, using a specific algorithm, the water vapour and temperature profiles were retrieved from the FIRMOS spectra by applying the Kyoto protocol and Informed Management of the Adaptation (KLIMA) code. The profiles were found in very good agreement with those provided by radiosondes and by the Raman lidar operating from the Zugspitze Schneefernerhaus station. In addition, the retrieval products were validated by comparing the retrieved integrated water vapour values with those obtained from the E-AERI spectra.

How adequately are elevated moist layers represented in reanalysis and satellite observations?

Abstract. We assess the representation of elevated moist layers (EMLs) in ERA5 reanalysis, the Infrared Atmospheric Sounding Interferometer (IASI) L2 retrieval Climate Data Record (CDR) and the Atmospheric Infrared Sounder (AIRS)-based Community Long-term Infrared Microwave Combined Atmospheric Product System (CLIMCAPS)-Aqua L2 retrieval. EMLs are free-tropospheric moisture anomalies that typically occur in the vicinity of deep convection in the tropics. EMLs significantly affect the spatial structure of radiative heating, which is considered a key driver for meso-scale dynamics, in particular convective aggregation. To our knowledge, the representation of EMLs in the mentioned data products has not been explicitly studied – a gap we start to address in this work. We assess the different datasets' capability of capturing EMLs by collocating them with 2146 radiosondes launched from Manus Island within the western Pacific warm pool, a region where EMLs occur particularly often. We identify and characterise moisture anomalies in the collocated datasets in terms of moisture anomaly strength, vertical thickness and altitude. By comparing the distributions of these characteristics, we deduce what specific EML characteristics the datasets are capturing well and what they are missing. Distributions of ERA5 moisture anomaly characteristics match those of the radiosonde dataset quite well, and remaining biases can be removed by applying a 1 km moving average to the radiosonde profiles. We conclude that ERA5 is a suitable reference dataset for investigating EMLs. We find that the IASI L2 CDR is subject to stronger smoothing than ERA5, with moisture anomalies being on average 13 % weaker and 28 % thicker than collocated ERA5 anomalies. The CLIMCAPS L2 product shows significant biases in its mean vertical humidity structure compared to the other investigated datasets. These biases manifest as an underestimation of mean moist layer height of about 1.3 km compared to the three other datasets that yields a general mid-tropospheric moist bias and an upper-tropospheric dry bias. Aside from these biases, the CLIMCAPS L2 product shows a similar, if not better, capability of capturing EMLs compared to the IASI L2 CDR. More nuanced evaluations of CLIMCAPS' capabilities may be possible once the underlying cause for the identified biases has been found and fixed. Biases found in the all-sky scenes do not change significantly when limiting the analysis to clear-sky scenes. We calculate radiatively driven vertical velocities ωrad derived from longwave heating rates to estimate the dynamical effect of the moist layers. Moist-layer-associated ωrad values derived from Global Climate Observing System Reference Upper-Air Network (GRUAN) soundings range between 2 and 3 hPa h−1, while mean meso-scale pressure velocities from the EUREC4A (Elucidating the Role of Clouds-Circulation Coupling in Climate) field campaign range between 1 and 2 hPa h−1, highlighting the dynamical significance of EMLs. Subtle differences in the representation of moisture and temperature structures in ERA5 and the satellite datasets create large relative errors in ωrad on the order of 40 % to 80 % with reference to GRUAN, indicating limited usefulness of these datasets to assess the dynamical impact of EMLs.

Data assimilation for the Model for Prediction Across Scales – Atmosphere with the Joint Effort for Data assimilation Integration (JEDI-MPAS 1.0.0): EnVar implementation and evaluation

Abstract. On 24 September 2021, JEDI-MPAS 1.0.0, a new data assimilation (DA) system for the Model Prediction Across Scales – Atmosphere (MPAS-A) built on the software framework of the Joint Effort for Data assimilation Integration (JEDI) was publicly released for community use. Operating directly on the native MPAS unstructured mesh, JEDI-MPAS capabilities include three-dimensional variational (3DVar) and ensemble–variational (EnVar) schemes as well as the ensemble of DA (EDA) technique. On the observation side, one advanced feature in JEDI-MPAS is the full all-sky approach for satellite radiance DA with the introduction of hydrometeor analysis variables. This paper describes the formulation and implementation of EnVar for JEDI-MPAS. JEDI-MPAS 1.0.0 is evaluated with month-long cycling 3DEnVar experiments with a global 30–60 km dual-resolution configuration. The robustness and credible performance of JEDI-MPAS are demonstrated by establishing a benchmark non-radiance DA experiment, then incrementally adding microwave radiances from three sources: Advanced Microwave Sounding Unit-A (AMSU-A) temperature sounding channels in clear-sky scenes, AMSU-A window channels in all-sky scenes, and Microwave Humidity Sounder (MHS) water vapor channels in clear-sky scenes. JEDI-MPAS 3DEnVar behaves well with a substantial and significant positive impact obtained for almost all aspects of forecast verification when progressively adding more microwave radiance data. In particular, the day 5 forecast of the best-performing JEDI-MPAS experiment yields an anomaly correlation coefficient (ACC) of 0.8 for 500 hPa geopotential height, a gap of roughly a half day when compared to cold-start forecasts initialized from operational analyses of the National Centers for Environmental Prediction, whose ACC does not drop to 0.8 until a lead time of 5.5 d. This indicates JEDI-MPAS's great potential for both research and operations.

Assimilation of AMSU-A in All-Sky Conditions

Abstract Radiances from microwave temperature sounders have been assimilated operationally at ECMWF for two decades, but observations significantly affected by clouds and precipitation have been screened out. Extending successful assimilation beyond clear-sky scenes is a challenge that has taken several years of development to achieve. In this paper we describe the all-sky treatment of AMSU-A, which enables greater numbers of temperature sounding radiances to be used in meteorologically active parts of the troposphere. Successful all-sky assimilation required combining lessons learned from the clear-sky assimilation of AMSU-A with the approach initially developed for humidity-sensitive microwave radiances. This concerned particularly observation thinning, error modeling, and variational quality control. As a result of the move to all-sky assimilation, the forecast impact of AMSU-A now replicates and exceeds that of the previous clear-sky usage. This is shown via trials in comparison to the current ECMWF assimilation system, judged with respect to forecast scores and background fits to independent observations. Persistently cloudy regions and phenomena such as tropical cyclones are better sampled when assimilating AMSU-A in all-sky conditions, causing an increase of about 13% in used channel-5 radiances globally. These impacts are explored, with an emphasis on tropical cyclones in the 2019 season. Independent observations provide consistent evidence that representation of humidity is improved, for example, while extratropical Z500 forecasts are improved by about 0.5% out to at least day 2. On the strength of these results, assimilation of AMSU-A moved to all-sky conditions with the upgrade to IFS cycle 47R3 in October 2021.

Retrieving H<sub>2</sub>O/HDO columns over cloudy and clear-sky scenes from the Tropospheric Monitoring Instrument (TROPOMI)

Abstract. This paper presents an extended scientific HDO/H2O total column data product from short-wave infrared (SWIR) measurements by the Tropospheric Monitoring Instrument (TROPOMI) including clear-sky and cloudy scenes. The retrieval employs a forward model which accounts for scattering, and the algorithm infers the trace gas column information, surface properties, and effective cloud parameters from the observations. Compared to the previous clear-sky-only data product, coverage is greatly enhanced by including scenes over low clouds, particularly enabling data over oceans as the albedo of water in the SWIR spectral range is too low to retrieve under cloud-free conditions. The new dataset is validated against co-located ground-based Fourier transform infrared (FTIR) observations by the Total Carbon Column Observing Network (TCCON). The median bias for clear-sky scenes is 1.4×1021 molec cm−2 (2.9 %) in H2O columns and 1.1×1017 molec cm−2 (−0.3 %) in HDO columns, which corresponds to −17 ‰ (9.9 %) in a posteriori δD. The bias for cloudy scenes is 4.9×1021 molec cm−2 (11 %) in H2O, 1.1×1018 molec cm−2 (7.9 %) in HDO, and −20 ‰ (9.7 %) in a posteriori δD. At low-altitude stations, the bias is small at low and middle latitudes and has a larger value at high latitudes. At high-altitude stations, an altitude correction is required to compensate for different partial columns seen by the station and the satellite. The bias in a posteriori δD after altitude correction depends on sensitivity due to shielding by clouds and on realistic a priori profile shapes for both isotopologues. Cloudy scenes generally involve low sensitivity below the clouds, and since the information is filled up by the prior, a realistic shape of the prior is important for realistic total column estimation in these cases. Over oceans, aircraft measurements with the Water Isotope System for Precipitation and Entrainment Research (WISPER) instrument from a field campaign in 2018 are used for validation, yielding biases of −3.9 % in H2O and −3 ‰ in δD over clouds. To demonstrate the added value of the new dataset, a short case study of a cold air outbreak over the Atlantic Ocean in January 2020 is presented, showing the daily evolution of the event with single-overpass results.

Total water vapour columns derived from Sentinel 5P using the AMC-DOAS method

Abstract. Water vapour is the most abundant natural greenhouse gas in the Earth's atmosphere, and global data sets are required for meteorological applications and climate research. The Tropospheric Monitoring Instrument (TROPOMI) on board Sentinel-5 Precursor (S5P) launched on 13 October 2017 has a high spatial resolution of around 5 km and a daily global coverage. Currently, there is no operational total water vapour product for S5P measurements. Here, we present first results of a new scientific total column water vapour (TCWV) product for S5P using the so-called air-mass-corrected differential optical absorption spectroscopy (AMC-DOAS) scheme. This method analyses spectral data between 688 and 700 nm and has already been successfully applied to measurements from the Global Monitoring Experiment (GOME) on ERS-2, the SCanning Imaging Absorption spectroMeter for Atmospheric CHartographY (SCIAMACHY) on Envisat and GOME-2 on MetOp. The adaptation of the AMC-DOAS method to S5P data requires an additional post-processing procedure to correct the influences of surface albedo, cloud height and cloud fraction. The quality of the new AMC-DOAS S5P water vapour product is assessed by comparisons with data from GOME-2 on MetOp-B retrieved also with the AMC-DOAS algorithm and with four independent data sets, namely reanalysis data from the European Centre for Medium range Weather Forecast (ECMWF ERA5), data obtained by the Special Sensor Microwave Imager and Sounder (SSMIS) flown on the Defense Meteorological Satellite Program (DMSP) platform 16 and two scientific S5P TCWV products derived from TROPOMI measurements. Both are recently published TCWV products for S5P provided by the Max Planck Institute for Chemistry (MPIC) in Mainz and the Netherlands Institute for Space Research (SRON), Utrecht. The SRON TCWV is limited to clear-sky scenes over land. These comparisons reveal a good agreement between the various data sets but also some systematic differences between all of them. On average, the daily derived offset between AMC-DOAS S5P TCWV and AMC-DOAS GOME-2B TCWV is negative (around −1.5 kg m−2) over land and positive over ocean surfaces (more than 1.5 kg m−2). In contrast, SSMIS TCWV is on average lower than AMC-DOAS S5P TCWV by about 3 kg m−2. Monthly averaged ERA5 TCWV and AMC-DOAS S5P TCWV comparison shows spatial features over both land and water surface. Over land, there are systematic spatial structures. There are larger differences between AMC-DOAS S5P TCWV and ERA5 TCWV in tropical regions. Over sea, AMC-DOAS S5P TCWV is slightly lower than ERA5 TCWV by around 2 kg m−2. The AMC-DOAS S5P TCWV and S5P TCWV from MPIC agree on average within 1 kg m−2 over both land and ocean. TCWV from SRON shows daily global averaged differences to AMC-DOAS S5P TCWV of around 1.2 kg m−2. All of these differences are in line with the accuracy of these products and with the typical range of differences of 5 kg m−2 obtained when comparing different TCWV data sets. The AMC-DOAS TCWV product for S5P provides therefore a valuable new and independent data set for atmospheric applications which also has a higher spatial coverage than the other S5P TCWV products.

Impacts of MOPITT cloud detection revisions on observation frequency and mapping of highly polluted scenes

Measurements made by the MOPITT (“Measurements of Pollution in the Troposphere”) instrument on the NASA Terra polar-orbiting platform enable the retrieval of tropospheric carbon monoxide (CO) concentrations. As determined by the Terra orbit and MOPITT swath width, the frequency of MOPITT observations at a specific location, or measurement sampling frequency, is typically about once every three to four days. However, because the MOPITT retrieval algorithm is only applicable to clear-sky scenes, MOPITT retrieval sampling frequency strongly depends on regional cloudiness and can be much smaller than the measurement sampling frequency. Moreover, highly polluted scenes, characterized by high aerosol optical depths, can be confused with cloudy scenes and thus be discarded unnecessarily by the MOPITT cloud detection algorithm. Herein are described revisions to this algorithm which substantially increase retrieval sampling over land in varying pollution conditions. The performance of the revised cloud detection algorithm is evaluated through validation, case studies, and continental-scale maps of retrieval sampling frequency. Presented case studies illustrate (1) why the current operational MOPITT cloud detection algorithm excludes extended areas of potentially valuable cloud-free MOPITT observations and (2) how, for the same scenes, improved retrieval coverage benefits analyses of regional CO variability. Maps of retrieval sampling frequency for South America and Asia exhibit well-defined improvements, especially in regions with poor sampling frequency in the current product.

Wetland Hydroperiod Change Along the Upper Columbia River Floodplain, Canada, 1984 to 2019

Increasing air temperatures and changing hydrological conditions in the mountainous Kootenay Region of British Columbia, Canada are expected to affect floodplain wetland extent and function along the Columbia River. The objective of this study was to determine the seasonally inundated hydroperiod for a floodplain section (28.66 km2) of the Upper Columbia River wetlands complex using time series satellite image observations and binary open water mask extraction. A mid pixel resolution (30 m) optical satellite image time series of 61 clear sky scenes from the Landsat Thematic Mapper (TM) and Operational Land Imager (OLI) sensors were used to map temporal variations in floodplain open water wetland extent during the April to October hydrologically active season from 1984 to 2019 (35 years). The hydroperiod from the first 31 scenes (T1: 18 years) was compared to the second 30 (T2: 16 years) to identify changes in the permanent and seasonal open water bodies. The seasonal variation in open water extent and duration was similar across the two time periods but the permanent water body extent diminished by ~16% (or ~3.5% of the floodplain). A simple linear model (r2 = 0.87) was established to predict floodplain open water extent as a function of river discharge downstream of the case study area. Four years of Landsat Multi-Spectral Scanner (MSS) data from 1992 to 1995 (12 scenes) were examined to evaluate the feasibility of extending the hydroperiod record back to 1972 using lower resolution (60 m) archive data. While the MSS hydroperiod produced a similar pattern of open water area to duration to the TM/OLI hydroperiod, small open water features were omitted or expanded due to the lower resolution. While MSS could potentially extend the TM/OLI hydroperiod record, this was not performed as the loss of features like the river channel diminished its value for change detection purposes. Radarsat 2 scenes from 2015 to 2019 were examined to evaluate the feasibility of continued mountain valley hydroperiod monitoring using higher spatial and temporal resolution sensors like the Radarsat Constellation Mission (RCM). From the available horizontal transmit/receive (HH) single polarization sample set (8 scenes), the hydroperiod pattern of open water extent to duration was similar to the longer Landsat time series and possessed greater feature detail, but it was significantly reduced in seasonal inundation area due to the systematic omission of open water areas containing emergent vegetation. However, accepting that differences exist in sensor-based hydroperiod attributes, the higher temporal resolution of RCM will be suited to mountain floodplain inundation monitoring and open water hydroperiod analysis.

Artificial neural networks model based on remote sensing to retrieve evapotranspiration over the Brazilian Pampa

Evapotranspiration (ET) quantification improves the comprehension of the water, heat, and carbon interactions and the feedback to the climate, which is essential for global change research. We aimed to model ET using artificial neural networks (ANNs) based on Landsat-8 and reanalysis data from the National Centers for Environmental Prediction over the grasslands of the Pampa biome. The output variable was the ET trained by eddy covariance (EC) measurements acquired from a flux tower located in Santa Maria, Brazil. ANN was performed using the backpropagation algorithm with four remote sensing input variables (albedo, normalized difference vegetation index, land surface temperature, and surface net radiation). In addition, four meteorological variables from the Environmental Prediction Climate Forecast System Version 2 hourly product were included in the model (air temperature, atmospheric pressure, relative humid, and wind speed). We analyzed 67 clear-sky scenes between 2014 and 2019. Results produced very robust daily ET estimates. ANN exhibited a correlation of 0.88 relative to in situ EC data, demonstrating a good linear relationship between ET estimated and measured and producing a root-mean-square error (mean absolute error) of 0.75 (0.58) mm/day. The ANN model was also compared with the widely known simplified surface energy balance index (S-SEBI) model. S-SEBI exhibited lower correlation with the ET in situ compared to the ANN model. Furthermore, the ANN model had a superior performance in summer and winter seasons in which S-SEBI was found to outperform the ET in situ. The model developed in our research is an alternative to approaches that need a great number of input variables or in situ data since it is only dependent on freely available data. Therefore, it should support future integrated strategies of water resources allocation over the natural grasslands of the Brazilian Pampa.

Improved SIFTER v2 algorithm for long-term GOME-2A satellite retrievals of fluorescence with a correction for instrument degradation

Abstract. Solar-induced fluorescence (SIF) data from satellites are increasingly used as a proxy for photosynthetic activity by vegetation and as a constraint on gross primary production. Here we report on improvements in the algorithm to retrieve mid-morning (09:30 LT) SIF estimates on the global scale from the GOME-2 sensor on the MetOp-A satellite (GOME-2A) for the period 2007–2019. Our new SIFTER (Sun-Induced Fluorescence of Terrestrial Ecosystems Retrieval) v2 algorithm improves over a previous version by using a narrower spectral window that avoids strong oxygen absorption and being less sensitive to water vapour absorption, by constructing stable reference spectra from a 6-year period (2007–2012) of atmospheric spectra over the Sahara and by applying a latitude-dependent zero-level adjustment that accounts for biases in the data product. We generated stable, good-quality SIF retrievals between January 2007 and June 2013, when GOME-2A degradation in the near infrared was still limited. After the narrowing of the GOME-2A swath in July 2013, we characterised the throughput degradation of the level-1 data in order to derive reflectance corrections and apply these for the SIF retrievals between July 2013 and December 2018. SIFTER v2 data compare well with the independent NASA v2.8 data product. Especially in the evergreen tropics, SIFTER v2 no longer shows the underestimates against other satellite products that were seen in SIFTER v1. The new data product includes uncertainty estimates for individual observations and is best used for mostly clear-sky scenes and when spectral residuals remain below a certain spectral autocorrelation threshold. Our results support the use of SIFTER v2 data being used as an independent constraint on photosynthetic activity on regional to global scales.

Determining the drivers of suspended sediment dynamics in tidal marsh-influenced estuaries using high-resolution ocean color remote sensing

Sediment budgets are a critical metric to assess coastal marsh vulnerability to sea level rise and declining riverine sediment inputs. However, calculating accurate sediment budgets is challenging in tidal marsh-influenced estuaries where suspended sediment concentrations (SSC) typically vary on scales of hours and hundreds of meters, and where SSC dynamics are driven by a complex and often site-specific interplay of hydrodynamic and meteorological conditions. The mapping of SSC using ocean-color remote sensing is well established and can help capture the spatio-temporal variability of SSC and determine the dominant drivers regulating sediment budgets. However, the coarse spatial resolution of traditional ocean-color sensors (1-km) generally precludes their use in coastal-marsh estuaries. Here, using the Plum Island Estuary (Massachusetts, USA) as an example, we demonstrate that high-spatial-resolution maps of SSC derived from Landsat-8 Operational Land Imager (OLI) and Sentinel-2A/B Multispectral Instruments (MSI) can be used to determine the main drivers of SSC dynamics in tidal marsh-influenced estuaries, despite the long revisit time of these sensors. Local empirical algorithms between SSC and remote-sensing reflectance were derived and applied to a total of 46 clear-sky scenes collected by the OLI and the MSI between 2013 and 2018. The analysis revealed that this 5-year record was sufficient to capture a representative range of meteorological and tidal conditions required to determine the main drivers of SSC dynamics in this mid-latitude system. The interplay between river and tidal flows dominated SSC dynamics in this estuary, whereas wind-driven resuspension had a more moderate effect. The SSC was higher during spring because of increased river discharge due to snowmelt. Tidal asymmetry also enhanced sediment resuspension during flood tides, possibly favoring deposition on marsh platforms. Together, water level, water-level rate of change, river discharge and wind speed were able to explain >60% of the variability in the main channel SSC, thereby facilitating future prediction of SSC from these readily available variables. This study demonstrates that the existing multi-year records of high-resolution remote sensing can provide a representative depiction of SSC dynamics in hydrodynamically-complex and small-scale estuaries that moderate-resolution ocean color remote sensing and in situ measurements are unable to capture.

An Efficient Framework for Producing Landsat-Based Land Surface Temperature Data Using Google Earth Engine

A long time-series land surface temperature (LST) product is useful for ecological and environmental studies. However, current LST products cannot provide a global coverage at a fine spatial resolution (∼100 m) over a long period (>30 years). Landsat series satellites that have been launched since 1972 provide a unique opportunity to fill the gap. Here, we proposed a single-channel framework for producing global long time-series Landsat LST retrievals on a Google earth engine (GEE) cloud computing platform. This framework unifies the LST, land surface emissivity (LSE) and atmospheric water vapor (AWV) estimation algorithms, as well as the emissivity and atmospheric input data for the Landsat LST retrievals from the entire Landsat thermal infrared image archive. In situ LST measurements and the MODIS LST products were employed to evaluate Landsat LST retrievals using the proposed framework over land and water surfaces, respectively. In total, 1317 clear-sky LST samples were collected from the Landsat 5–8 series after spatiotemporal registration with seven sites, and the average bias and root-mean-square error (RMSE) were 0.33 and 2.01 K, respectively. Intercomparison between Landsat and MODIS LST retrievals based on 100 clear-sky scenes over 12 inland lakes showed an average bias of 0.17 K and RMSE of 1.11 K. We conclude that the proposed single-channel framework can produce Landsat LST with high accuracy following a simple yet robust way. Implementation of the single-channel method on GEE shows promise in providing the community with freely accessible and global long time-series (>30 years) LST data.

An experimental 2D-Var retrieval using AMSR2

Abstract. A two-dimensional variational retrieval (2D-Var) is presented for a passive microwave imager. The overlapping antenna patterns of all frequencies from the Advanced Microwave Scanning Radiometer 2 (AMSR2) are explicitly simulated to attempt retrieval of near-surface wind speed and surface skin temperature at finer spatial scales than individual antenna beams. This is achieved, with the effective spatial resolution of retrieved parameters judged by analysis of 2D-Var averaging kernels. Sea surface temperature retrievals achieve about 30 km resolution, with wind speed retrievals at about 10 km resolution. It is argued that multi-dimensional optimal estimation permits greater use of total information content from microwave sensors than other methods, with no compromises on target resolution needed; instead, various targets are retrieved at the highest possible spatial resolution, driven by the channels' sensitivities. All AMSR2 channels can be simulated within near their published noise characteristics for observed clear-sky scenes, though calibration and emissivity model errors are key challenges. This experimental retrieval shows the feasibility of 2D-Var for cloud-free retrievals and opens the possibility of stand-alone 3D-Var retrievals of water vapour and hydrometeor fields from microwave imagers in the future. The results have implications for future satellite missions and sensor design, as spatial oversampling can somewhat mitigate the need for larger antennas in the push for higher spatial resolution.

Investigating the liquid water path over the tropical Atlantic with synergistic airborne measurements

Abstract. Liquid water path (LWP) is an important quantity to characterize clouds. Passive microwave satellite sensors provide the most direct estimate on a global scale but suffer from high uncertainties due to large footprints and the superposition of cloud and precipitation signals. Here, we use high spatial resolution airborne microwave radiometer (MWR) measurements together with cloud radar and lidar observations to better understand the LWP of warm clouds over the tropical North Atlantic. The nadir measurements were taken by the German High Altitude and LOng range research aircraft (HALO) in December 2013 (dry season) and August 2016 (wet season) during two Next-generation Advanced Remote sensing for VALidation (NARVAL) campaigns. Microwave retrievals of integrated water vapor (IWV), LWP, and rainwater path (RWP) are developed using artificial neural network techniques. A retrieval database is created using unique cloud-resolving simulations with 1.25 km grid spacing. The IWV and LWP retrievals share the same eight MWR frequency channels in the range from 22 to 31 GHz and at 90 GHz as their sole input. The RWP retrieval combines active and passive microwave observations and is able to detect drizzle and light precipitation. The comparison of retrieved IWV with coincident dropsondes and water vapor lidar measurements shows root-mean-square deviations below 1.4 kg m−2 over the range from 20 to 60 kg m−2. This comparison raises the confidence in LWP retrievals which can only be assessed theoretically. The theoretical analysis shows that the LWP error is constant with 20 g m−2 for LWP below 100 g m−2. While the absolute LWP error increases with increasing LWP, the relative one decreases from 20 % at 100 g m−2 to 10 % at 500 g m−2. The identification of clear-sky scenes by ancillary measurements, here backscatter lidar, is crucial for thin clouds (LWP < 12 g m−2) as the microwave retrieved LWP uncertainty is higher than 100 %. The analysis of both campaigns reveals that clouds were more frequent (47 % vs. 30 % of the time) in the dry than in the wet season. Their average LWP (63 vs. 40 g m−2) and RWP (6.7 vs. 2.7 g m−2) were higher as well. Microwave scattering of ice, however, was observed less frequently in the dry season (0.5 % vs. 1.6 % of the time). We hypothesize that a higher degree of cloud organization on larger scales in the wet season reduces the overall cloud cover and observed LWP. As to be expected, the observed IWV clearly shows that the dry season is on average less humid than the wet season (28 vs. 41 kg m−2). The results reveal that the observed frequency distributions of IWV are substantially affected by the choice of the flight pattern. This should be kept in mind when using the airborne observations to carefully mediate between long-term ground-based and spaceborne measurements to draw statistically sound conclusions.

A full-mission data set of H<sub>2</sub>O and HDO columns from SCIAMACHY 2.3 µm reflectance measurements

Abstract. A new data set of vertical column densities of the water vapour isotopologues H2O and HDO from the SCanning Imaging Absorption spectroMeter for Atmospheric CHartographY (SCIAMACHY) instrument for the whole of the mission period from January 2003 to April 2012 is presented. The data are retrieved from reflectance measurements in the spectral range 2339 to 2383 nm with the Shortwave Infrared CO Retrieval (SICOR) algorithm, ignoring atmospheric light scattering in the measurement simulation. The retrievals are validated with ground-based Fourier transform infrared measurements obtained within the Multi-platform remote Sensing of Isotopologues for investigating the Cycle of Atmospheric water (MUSICA) project. A good agreement for low-altitude stations is found with an average bias of −3.6×1021 for H2O and −1.0×1018 molec cm−2 for HDO. The a posteriori computed δD shows an average bias of −8 ‰, even though polar stations have a larger negative bias. The latter is due to the large amount of sensor noise in SCIAMACHY in combination with low albedo and high solar zenith angles. To demonstrate the benefit of accounting for light scattering in the retrieval, the quality of the data product fitting effective cloud parameters simultaneously with trace gas columns is evaluated in a dedicated case study for measurements round high-altitude stations. Due to a large altitude difference between the satellite ground pixel and the mountain station, clear-sky scenes yield a large bias, resulting in a δD bias of 125 ‰. When selecting scenes with optically thick clouds within 1000 m above or below the station altitude, the bias in a posteriori δD is reduced from 125 to 44 ‰. The insights from the present study will also benefit the analysis of the data from the new Sentinel-5 Precursor mission.

Deriving clear-sky longwave spectral flux from spaceborne hyperspectral radiance measurements: a case study with AIRS observations

Abstract. Previous studies have shown that longwave (LW) spectral fluxes have unique merit in climate studies. Using Atmospheric Infrared Sounder (AIRS) radiances as a case study, this study presents an algorithm to derive the entire LW clear-sky spectral fluxes from spaceborne hyperspectral observations. No other auxiliary observations are needed in the algorithm. A clear-sky scene is identified using a three-step detection method. The identified clear-sky scenes are then categorized into different sub-scene types using information about precipitable water, lapse rate and surface temperature inferred from the AIRS radiances at six selected channels. A previously established algorithm is then used to invert AIRS radiances to spectral fluxes over the entire LW spectrum at 10 cm−1 spectral interval. Accuracy of the algorithms is evaluated against collocated Clouds and the Earth's Radiant Energy System (CERES) observations. For nadir-view observations, the mean difference between outgoing longwave radiation (OLR) derived by this algorithm and the collocated CERES OLR is 1.52 Wm−2 with a standard deviation of 2.46 Wm−2. When the algorithm is extended for viewing zenith angle up to 45°, the performance is comparable to that for nadir-view results.

Evaluation of two Vaisala RS92 radiosonde solar radiative dry bias correction algorithms

Abstract. Solar heating of the relative humidity (RH) probe on Vaisala RS92 radiosondes results in a large dry bias in the upper troposphere. Two different algorithms (Miloshevich et al., 2009, MILO hereafter; and Wang et al., 2013, WANG hereafter) have been designed to account for this solar radiative dry bias (SRDB). These corrections are markedly different with MILO adding up to 40 % more moisture to the original radiosonde profile than WANG; however, the impact of the two algorithms varies with height. The accuracy of these two algorithms is evaluated using three different approaches: a comparison of precipitable water vapor (PWV), downwelling radiative closure with a surface-based microwave radiometer at a high-altitude site (5.3 km m.s.l.), and upwelling radiative closure with the space-based Atmospheric Infrared Sounder (AIRS). The PWV computed from the uncorrected and corrected RH data is compared against PWV retrieved from ground-based microwave radiometers at tropical, midlatitude, and arctic sites. Although MILO generally adds more moisture to the original radiosonde profile in the upper troposphere compared to WANG, both corrections yield similar changes to the PWV, and the corrected data agree well with the ground-based retrievals. The two closure activities – done for clear-sky scenes – use the radiative transfer models MonoRTM and LBLRTM to compute radiance from the radiosonde profiles to compare against spectral observations. Both WANG- and MILO-corrected RHs are statistically better than original RH in all cases except for the driest 30 % of cases in the downwelling experiment, where both algorithms add too much water vapor to the original profile. In the upwelling experiment, the RH correction applied by the WANG vs. MILO algorithm is statistically different above 10 km for the driest 30 % of cases and above 8 km for the moistest 30 % of cases, suggesting that the MILO correction performs better than the WANG in clear-sky scenes. The cause of this statistical significance is likely explained by the fact the WANG correction also accounts for cloud cover – a condition not accounted for in the radiance closure experiments.