Current Atmospheric Correction Research Articles

Atmospheric correction under cloud edge effects for Geostationary Ocean Color Imager through deep learning

The Geostationary Ocean Color Imager (GOCI) is widely employed in tracking diurnal dynamics of oceanic conditions. However, the current atmospheric correction (AC) algorithms for GOCI often mask pixels around cloud edges to exclude pixels contaminated by cloud edge effects (CEEs; including stray light, cloud shadows, and cloud adjacent effects (AEs)), which results in massive data loss. In this paper, we propose a novel AC algorithm to correct these CEE-affected pixels based on deep learning (namely, DLACC) to achieve a comparable quality level to those pixels far away from cloud edges. We used the standard near-infrared iterative AC algorithm in SeaDAS (namely, NIR) to obtain Rayleigh-corrected reflectance (Rrc) and remote sensing reflectance (Rrs). Then, high-quality Rrs values were extracted using a quality assurance (QA) algorithm. Next, we obtained 2,821,668 matchups by matching CEE-free noontime Rrs values with CEE-affected Rrc values in a time window of 1 h. These matchups were used to develop DLACC. Validations using in situ data from Aerosol Robotic Network-Ocean Color (AERONET-OC) stations showed that more matchups were obtained with the DLACC model than with the NIR algorithm and Korea Ocean Satellite Center standard AC algorithm in GDPS 2.0 (KOSC), and the accuracies were similar. More importantly, the DLACC algorithm is more tolerant of AEs and reduces AEs in a 2-pixel distance from cloud edges by 10% in the blue bands and by 50% in the green band. As a result, more valid observations are obtained with the DLACC algorithm, with the daily percentage of valid observations (DPVOs) increasing by 71% and 62% over those with the NIR algorithm and the KOSC algorithm, respectively. These increases in valid observations could further result in more consistent patterns in terms of space and time. Our DLACC algorithm can not only be used to process GOCI images but also provide an open framework to develop corresponding deep-learning models for other geostationary satellites.

The Inversion of HY-1C-COCTS Ocean Color Remote Sensing Products from High-Latitude Seas

China’s first operational ocean color satellite sensor, named the Chinese Ocean Color and Temperature Scanner (HY-1C-COCTS), was launched in September 2018 and began to provide operational data in June 2019. However, as a polar orbiting ocean color satellite sensor, HY-1C-COCTS would inevitably encounter regions impacted by large solar zenith angles when observing the high-latitude seas, especially during the winter. The current atmospheric correction algorithm used by ocean color satellite data processing software cannot effectively process observation data with solar zenith angles greater than 70°. This results in a serious lack of effective ocean color product data from high-latitude seas in winter. To solve this problem, this study developed an atmospheric correction algorithm based on a neural network model for use with HY-1C-COCTS data. The new algorithm used HY-1C-COCTS satellite data collected from latitudes greater than 50°N and between April 2020 and April 2021 to establish a direct relationship between the total radiance received by the satellite and the remote sensing reflectance products. The evaluation using the test dataset shows that the inversion accuracy of the new algorithm is relatively high under different solar zenith angles and different HY-1C-COCTS bands (the relative deviation is 3.37%, 7.05%, 5.10%, 5.29%, and 10.06% at 412 nm, 443 nm, 490 nm, 520 nm, and 565 nm, respectively; the relative deviation is 1.07% when the solar zenith angle is large (70~90°)). Cross comparison with MODIS Aqua satellite products shows that the inversion results are consistent. After verifying the accuracy and stability of the algorithm, we reconstructed the remote sensing reflectance dataset from the Arctic Ocean and surrounding high-latitude seas (latitude greater than 50°N) and successfully retrieved chlorophyll-a data and information on other marine ecological parameters.

Using Machine Learning for Timely Estimates of Ocean Color Information From Hyperspectral Satellite Measurements in the Presence of Clouds, Aerosols, and Sunglint

Retrievals of ocean color from space are important for better understanding of the ocean ecosystem but can be limited under conditions such as clouds, aerosols, and sunglint. Many ocean color algorithms use a few selected spectral bands to perform an atmospheric correction and then derive the upwelling radiance from the ocean. The limitations in the atmospheric correction under certain conditions lead to many gaps in daily spatial coverage of ocean color retrievals. To address these limitations, we introduce a new approach that uses machine learning to estimate ocean color from top of atmosphere radiances or reflectance measurements. In this approach, a principal component analysis is used to decompose the hyperspectral measurements into spectral features that describe the scattering and absorption of the atmosphere and the underlying surface. The coefficients of the principal components are then used to train a neural network to predict ocean color properties derived from the MODIS atmospheric correction algorithm. This machine learning approach is independent of a priori information and does not rely on any radiative transfer modeling. We apply the approach to two hyperspectral UV/VIS instruments, the ozone monitoring instrument (OMI) and the TROPOspheric Monitoring Instrument (TROPOMI), using measurements from 320–500 nm to show that it can be used to reproduce ocean color properties in less-than-ideal conditions. This machine learning approach complements the current atmospheric correction ocean color retrievals by filling in the gaps resulting from cloud, aerosol, and sunglint contamination. This method can be applied to the future hyperspectral Ocean Color Instrument (OCI), which will be onboard NASA’s Plankton, Aerosol Cloud, ocean Ecosystem (PACE) ocean color satellite set to launch in 2024.

GNSS-corrected InSAR displacement time-series spanning the 2019 Ridgecrest, CA earthquakes

SUMMARYInSAR displacement time-series are emerging as a valuable product to study a number of Earth processes. One challenge to current time-series processing methods, however, is that when large earthquakes occur, they can leave sharp coseismic steps in the time-series. These discontinuities can cause current atmospheric correction and noise smoothing algorithms to break down, as these algorithms commonly assume that deformation is steady through time. Here, we aim to remedy this by exploring two methods for correcting earthquake offsets in InSAR time-series: a simple difference offset estimate (SDOE) process and a multiparameter offset estimate (MPOE) parametric time-series inversion technique. We apply these methods to a 2-yr time-series of Sentinel-1 interferograms spanning the 2019 Ridgecrest, CA earthquake sequence. Descending track results indicate that the SDOE method precisely corrects for only 20 per cent of the coseismic offsets at 62 study locations included in our scene and only partially corrects or sometimes overcorrects for the rest of our study sites. On the other hand, the MPOE estimate method successfully corrects the coseismic offset for the majority of sites in our analysis. This MPOE method allows us to produce InSAR time-series and data-derived estimates of deformation during each phase of the earthquake cycle. In order to better isolate and estimate the signal of post-seismic lithospheric deformation in the InSAR time-series, we apply a GNSS-based correction to our interferograms. This correction ties the interferograms to median-filtered weekly GNSS displacements and removes additional atmospheric artefacts. We present InSAR-based estimates of post-seismic deformation for the area around the Ridgecrest rupture, as well as a 2-yr coseismic-corrected, GNSS-corrected InSAR time-series data set. This GNSS-corrected InSAR time-series will enable future modelling of post-seismic processes such as afterslip in the near field of the rupture, poroelastic deformation at intermediate distances and viscoelastic deformation at longer timescales in the far field.

Sensitivity of Satellite Ocean Color Data to System Vicarious Calibration of the Long Near Infrared Band

Satellite ocean color missions require accurate system vicarious calibrations (SVC) to retrieve the relatively small remote-sensing reflectance ( $R_{\mathrm {rs}}$ , sr−1) from the at-sensor radiance. However, the current atmospheric correction and SVC procedures do not include calibration of the “long” near infrared band (NIRL—869 nm for MODIS), partially because earlier studies, based primarily on simulations, indicate that accuracy in the retrieved $R_{\mathrm {rs}}$ is insensitive to moderate changes in the NIRL vicarious gain ( $g$ ). However, the sensitivity of ocean color data products to $g$ (NIRL) has not been thoroughly examined. Here, we first derive 10 SVC “gain configurations” (vicarious gains for all visible and NIR bands) for MODIS/Aqua using current operational NASA protocols, each time assuming a different $g$ (869). From these, we derive a suite of ~1.4E6 unique gain configurations with $g$ (869) ranging from 0.85 to 1.2. All MODIS/A data for 25 locations within each of five ocean gyres were then processed using each of these gain configurations. Resultant time series show substantial variability in dominant $R_{\mathrm {rs}}$ (547) patterns in response to changes in $g$ (869) (and associated gain configurations). Overall, mean $R_{\mathrm {rs}}$ (547) values generally decrease with increasing $g$ (869), while the standard deviations around those means show gyre-specific minima for $0.97 (869) $g$ (869) using such time series, finding $g$ (869) = 1.025 most closely comports with expectations. This approach is broadly applicable to other ocean color sensors, and highlights the importance of rigorous cross-sensor calibration of the NIRL bands, with implications on consistency of merged-sensor data sets.

Multiband Atmospheric Correction Algorithm for Ocean Color Retrievals

NASA’s current Atmospheric Correction (AC) algorithm for ocean color utilizes two bands and their ratio in the Near Infrared (NIR) to estimate aerosol reflectance and aerosol type. The algorithm then extrapolates the spectral dependence of aerosol reflectance to the visible wavelengths based on modeled spectral dependence of the identified aerosol type. Future advanced ocean color sensors, such as the Ocean Color Instrument (OCI) that will be carried on the Plankton, Aerosol, Cloud, and ocean Ecosystem (PACE) satellite, will be capable of measuring the hyperspectral radiance from 340 to 890 nm at 5-nm spectral resolution and at 7 discrete Short-wave Infrared (SWIR) channels: 940, 1038, 1250, 1378, 1615, 2130, and 2260 nm. To optimally employ this unprecedented instrument capability, we propose an improved AC algorithm that utilizes all atmospheric-window channels in the NIR to SWIR spectral range to reduce the uncertainty in the AC process. A theoretical uncertainty analysis of this, namely Multi-Band AC (MBAC), indicates that the algorithm can reduce the uncertainty in remote sensing reflectance (Rrs) retrievals of the ocean caused by sensor random noise. Furthermore, in optically complex waters, where the NIR signal is affected by contributions from highly-reflective turbid waters, the MBAC algorithm can be adaptively weighted to the strongly-absorbing SWIR channels to enable improved ocean color retrievals in coastal waters. We provide here a description of the algorithm and demonstrate the improved performance in ocean color retrievals, relative to the current NASA standard AC algorithm, through comparison with field measurements and assessment of propagated uncertainties in applying the MBAC algorithm to MODIS and simulated PACE OCI data.

Proposal for a General Atmospheric Correction Method of Breakdown and Withstand Voltages of Air-Gap Insulated Configurations Based on a Streamer–Leader Differentiated Model of the Breakdown Process

Current atmospheric correction methods for breakdown and withstand voltages in long air-gap configurations show deviations in comparison to measurement data when applied to correct for unusual or atypical atmospheric conditions (i.e., altitudes above 2000 m a.s.l.). The corresponding standards IEC 60060 and IEC 60071 show inconsistencies, although both are being based upon the same measurements and empirical models. In its first part, the concept paper at hand presents the current knowledge and problems of atmospheric correction. In the second part, the authors outline the possible steps for a necessary new investigation as well as the algorithm of a hypothetical correction method. The proposed research is based upon the recent knowledge that any breakdown process in air-gap insulation consists of the sub-processes streamer and leader. While their characteristics seem to be independent of configuration and voltage form, their ratio and proportion varies, which, in turn, defines the value of breakdown voltage. Therefore, the authors propose a sub-process differentiated, measurement-based evaluation of the atmospheric influences air pressure, air humidity, air conductivity and temperature. The main goal is to develop a generic physical model of the breakdown process.

Vicarious Calibration of sUAS Microbolometer Temperature Imagery for Estimation of Radiometric Land Surface Temperature.

In recent years, the availability of lightweight microbolometer thermal cameras compatible with small unmanned aerial systems (sUAS) has allowed their use in diverse scientific and management activities that require sub-meter pixel resolution. Nevertheless, as with sensors already used in temperature remote sensing (e.g., Landsat satellites), a radiance atmospheric correction is necessary to estimate land surface temperature. This is because atmospheric conditions at any sUAS flight elevation will have an adverse impact on the image accuracy, derived calculations, and study replicability using the microbolometer technology. This study presents a vicarious calibration methodology (sUAS-specific, time-specific, flight-specific, and sensor-specific) for sUAS temperature imagery traceable back to NIST-standards and current atmospheric correction methods. For this methodology, a three-year data collection campaign with a sUAS called “AggieAir”, developed at Utah State University, was performed for vineyards near Lodi, California, for flights conducted at different times (early morning, Landsat overpass, and mid-afternoon”) and seasonal conditions. From the results of this study, it was found that, despite the spectral response of microbolometer cameras (7.0 to 14.0 μm), it was possible to account for the effects of atmospheric and sUAS operational conditions, regardless of time and weather, to acquire accurate surface temperature data. In addition, it was found that the main atmospheric correction parameters (transmissivity and atmospheric radiance) significantly varied over the course of a day. These parameters fluctuated the most in early morning and partially stabilized in Landsat overpass and in mid-afternoon times. In terms of accuracy, estimated atmospheric correction parameters presented adequate statistics (confidence bounds under ±0.1 for transmissivity and ±1.2 W/m2/sr/um for atmospheric radiance, with a range of RMSE below 1.0 W/m2/sr/um) for all sUAS flights. Differences in estimated temperatures between original thermal image and the vicarious calibration procedure reported here were estimated from −5 °C to 10 °C for early morning, and from 0 to 20 °C for Landsat overpass and mid-afternoon times.

Performance evaluation of operational atmospheric correction algorithms over the East China Seas

To acquire high-quality operational data products for Chinese in-orbit and scheduled ocean color sensors, the performances of two operational atmospheric correction (AC) algorithms (ESA MEGS 7.4.1 and NASA SeaDAS 6.1) were evaluated over the East China Seas (ECS) using MERIS data. The spectral remote sensing reflectance R rs(λ), aerosol optical thickness (AOT), and Angstrom exponent (α) retrieved using the two algorithms were validated using in situ measurements obtained between May 2002 and October 2009. Match-ups of R rs, AOT, and α between the in situ and MERIS data were obtained through strict exclusion criteria. Statistical analysis of R rs(λ) showed a mean percentage difference (MPD) of 9%–13% in the 490–560 nm spectral range, and significant overestimation was observed at 413 nm (MPD>72%). The AOTs were overestimated (MPD>32%), and although the ESA algorithm outperformed the NASA algorithm in the blue-green bands, the situation was reversed in the red-near-infrared bands. The value of α was obviously underestimated by the ESA algorithm (MPD=41%) but not by the NASA algorithm (MPD=35%). To clarify why the NASA algorithm performed better in the retrieval of α, scatter plots of the α single scattering albedo (SSA) density were prepared. These α-SSA density scatter plots showed that the applicability of the aerosol models used by the NASA algorithm over the ECS is better than that used by the ESA algorithm, although neither aerosol model is suitable for the ECS region. The results of this study provide a reference to both data users and data agencies regarding the use of operational data products and the investigation into the improvement of current AC schemes over the ECS.

Spectral relationships for atmospheric correction I Validation of red and near infra-red marine reflectance relationships

The present study provides an extensive overview of red and near infra-red (NIR) spectral relationships found in the literature and used to constrain red or NIR-modeling schemes in current atmospheric correction (AC) algorithms with the aim to improve water-leaving reflectance retrievals, ρw(λ), in turbid waters. However, most of these spectral relationships have been developed with restricted datasets and, subsequently, may not be globally valid, explaining the need of an accurate validation exercise. Spectral relationships are validated here with turbid in situ data for ρw(λ). Functions estimating ρw(λ) in the red were only valid for moderately turbid waters (ρw(λNIR) < 3.10(-3)). In contrast, bounding equations used to limit ρw(667) retrievals according to the water signal at 555 nm, appeared to be valid for all turbidity ranges presented in the in situ dataset. In the NIR region of the spectrum, the constant NIR reflectance ratio suggested by Ruddick et al. (2006) (Limnol. Oceanogr. 51, 1167-1179), was valid for moderately to very turbid waters (ρw(λNIR) < 10(-2)) while the polynomial function, initially developed by Wang et al. (2012) (Opt. Express 20, 741-753) with remote sensing reflectances over the Western Pacific, was also valid for extremely turbid waters (ρw(λNIR) > 10(-2)). The results of this study suggest to use the red bounding equations and the polynomial NIR function to constrain red or NIR-modeling schemes in AC processes with the aim to improve ρw(λ) retrievals where current AC algorithms fail.

Correction of artifacts in the SeaWiFS atmospheric correction: removing discontinuity in the derived products

It has been found that, at the certain solar and sensor-viewing geometry and for certain atmosphere conditions, discontinuity lines appear in the Sea-viewing Wide Field-of-view Sensor (SeaWiFS) retrieved atmospheric and ocean color products. Such discontinuity lines, which do not happen very often, are apparently artifacts from the atmospheric correction scheme that uses the lookup tables to process SeaWiFS data. In this paper, a brief description of the SeaWiFS atmospheric correction algorithm, in particular, the technique that is used in retrieving the aerosol models and aerosol radiance contributions in the visible wavelengths, is provided. Results from some specific simulations that explain the causes of the discontinuity lines in the derived products are presented. We show that these discontinuities appearing in the derived SeaWiFS products are the result of the atmospheric correction due to aerosol model effects. To remove the discontinuities, a simple modification to the current atmospheric correction algorithm is proposed and tested with both simulated and real SeaWiFS data. The modification has been implemented in the SeaWiFS fourth data reprocessing in August 2002.

Remote sensing of ocean color and aerosol properties: resolving the issue of aerosol absorption

Current atmospheric correction and aerosol retrieval algorithms for ocean color sensors use measurements of the top-of-the-atmosphere reflectance in the near infrared, where the contribution from the ocean is known for case 1 waters, to assess the aerosol optical properties. Such measurements are incapable of distinguishing between weakly and strongly absorbing aerosols, and the atmospheric correction and aerosol retrieval algorithms fail if the incorrect absorption properties of the aerosol are assumed. We present an algorithm that appears promising for the retrieval of in-water biophysical properties and aerosol optical properties in atmospheres containing both weakly and strongly absorbing aerosols. By using the entire spectrum available to most ocean color instruments (412-865 nm), we simultaneously recover the ocean's bio-optical properties and a set of aerosol models that best describes the aerosol optical properties. The algorithm is applied to simulated situations that are likely to occur off the U.S. East Coast in summer when the aerosols could be of the locally generated weakly absorbing Maritime type or of the pollution-generated strongly absorbing urban-type transported over the ocean by the winds. The simulations show that the algorithm behaves well in an atmosphere with either weakly or strongly absorbing aerosol. The algorithm successfully identifies absorbing aerosols and provides close values for the aerosol optical thickness. It also provides excellent retrievals of the ocean bio-optical properties. The algorithm uses a bio-optical model of case 1 waters and a set of aerosol models for its operation. The relevant parameters of both the ocean and atmosphere are systematically varied to find the best (in a rms sense) fit to the measured top-of-the-atmosphere spectral reflectance. Examples are provided that show the algorithm's performance in the presence of errors, e.g., error in the contribution from whitecaps and error in radiometric calibration.