The geographic distribution and seasonal to inter‐annual variability of the particulate backscattering coefficient, bbp, is described and compared with that of chlorophyll concentration, Chl. Both ...variables are obtained from satellite‐derived ocean color and inverse modeling. In general, bbp exhibits larger seasonal variations than Chl, and the bbp values are more evenly distributed in space, especially at low and middle latitudes. A phase shift between the annual cycles of bbp and Chl is evidenced, and is attributed to the presence of an important pool of non‐pigmented particles. Converting bbp to particulate organic carbon, POC, the annual mean surface POC pool between 60°S and 60°N is about 19 Tg C per meter.
We recently found a significant bias between spectral diffuse attenuation coefficient (K
(λ)) retrievals by common ocean color algorithms and measurements from profiling floats Remote. Sens.14, 4500 ...(2022)10.3390/rs14184500. Here we show, using a multi-satellite match-up dataset, that the bias is markedly reduced by simple "tuning" of the algorithm's empirical coefficients. However, while the float dataset encompasses a larger proportion of the ocean's variability than previously used datasets, it does not cover the whole range of variability of observed remote sensing reflectance (R
). Thus, using algorithms tuned to this more comprehensive dataset may still result in a temporal and/or geographical bias in global application. To address this generalization issue, we evaluated a variety of analytical algorithms based on radiative transfer theory and settled on a specific one. This algorithm computes K
(λ) from inherent optical properties (IOPs) obtained from an R
inversion and information about the angular distribution of the radiance transmitted through the air/ocean interface. The resulting K
(λ) estimates at 412 and 490 nm were not appreciably biased against the float measurements. Evaluation using other in-situ datasets and radiative transfer simulations was also satisfactory. Statistical performance was good in both clear and turbid waters. Further work should be conducted to examine whether the tuned algorithms and/or the new analytical algorithm demonstrate adequate hyperspectral performance.
Two satellite-based methods to estimate daily averaged photosynthetically available radiation (PAR) at the ocean surface are evaluated at high northern latitudes. The first method employs a ...precomputed Look-Up-Table (LUT) generated from radiative transfer simulations. The LUT associates spectral irradiance reaching the surface to a given set of input parameters, namely solar zenith angle, cloud optical thickness, cloud fraction, ozone concentration, and surface albedo. The second method, as implemented by NASA's Ocean Biology Processing Group (OBPG) in the standard Ocean Color data processing chain, expresses the energy budget of the atmosphere-surface-ocean system via a simple radiative transfer model. The performance of these methods is evaluated using an extensive in situ PAR dataset collected in the Arctic Ocean between 1998 and 2014, with daily values ranging from 0.08 to 61.07Em−2d−1. A methodology is developed to compare in situ measurements and satellite products of different spatial and temporal resolution. Agreement is generally good between satellite-derived estimates and ship-based data and between methods. Specifically, both methods yield a small positive bias of 6% and 2% and a relative uncertainty larger than that observed at low latitude, with a root mean squared error (RMSE) of 33% and 20% for the LUT and OPBG methods, respectively. This is attributed to the peculiar environmental conditions encountered in the Arctic, namely low solar elevation, changing surface albedo due to sea ice, and persistent cloudiness. The RMSE difference among methods is due to the high temporal resolution (3h) of the International Satellite Cloud Climatology Project (ISCCP) LUT input not fully compensating for its low spatial resolution (280km grid cells). The LUT method has the major advantage of providing PAR estimates in all conditions, including ice-covered regions, while the OBPG method is currently limited to open waters and a solar zenith angle lower than 83°. Consequently, the OBPG method may not account for as much as 38% of PAR reaching the Arctic ocean surface annually. Both methods have the potential to provide useful PAR estimates just below the ice, by including information about ice transmittance.
•Performances of two satellite-based methods to estimate daily PAR are examined.•Over 700days of quality controlled Arctic in situ PAR dataset are gathered.•Errors are larger than those found at lower latitudes.•Each method is recommended for a particular spatial scale.
•A Monte Carlo radiative transfer code, accelerated with graphics card is described.•It simulates polarized light propagation in the ocean atmosphere coupled system.•It is fast and accurate: ...comparisons with reference codes are presented.•Impacts of atmospheric sphericity and adjacency effects are examined.•Perspectives for a fast, high spectral resolution version of the code are introduced.
SMART-G (Speed-Up Monte-Carlo Advanced Radiative Transfer code with GPU) is a radiative transfer solver for the coupled ocean-atmosphere system with a wavy interface. It is based on the Monte-Carlo technique, works in either plane-parallel or spherical-shell geometry, and accounts for polarization. The vector code is written in CUDA (Compute Unified Device Architecture) and runs on GPUs (Graphic Processing Units). For typical simulations, the GPU-based code running on the Nvidia GTX 1070 card is shown to be 100 time faster than a state of the art CPU-based code running on an AMD PhenomIIx4 965 at 3.4GHz. This makes SMART-G competitive, in terms of computational efficiency, with codes based on other techniques (e.g., discrete ordinate, doubling-adding, matrix-operator, and successive-orders-of-scattering), while allowing maximum flexibility regarding the scope of the simulations. The monochromatic version of the code without trans-spectral processes (Raman scattering, fluorescence) is described, including the treatment of photon propagation and interactive processes (elastic scattering, absorption, reflection, refraction, but not thermal emission) and variance reduction (local estimate). Benchmark values are accurately reproduced for clear and cloudy atmospheres over a wavy reflecting surface and a black ocean. Results obtained for a diffusely reflecting ocean agree with those from a discrete ordinate code. SMART-G may be used, not only as a reference code, but also to simulate the signal/imagery observed/produced by optical sensors, create accurate look-up tables, and investigate new remote sensing techniques.
The operational MEdium Resolution Imaging Spectrometer (MERIS) daily mean photosynthetically available radiation (PAR) product generated by the NASA Ocean Biology Processing Group (OBPG) was ...evaluated in clear sky conditions against in-situ measurements at various sites in the northwestern Mediterranean Sea (BOUSSOLE buoy), the northwestern Pacific (CCE-1 and -2 moorings), and the northeastern Atlantic (COVE platform). The measurements were first checked and corrected for calibration errors and uncertainties in data processing by comparing daily means for clear days (i.e., no clouds from sunrise to sunset and low aerosol abundance) with theoretical values from an accurate Monte Carlo radiative transfer code. The OBPG algorithm performed well when sky was completely cloudless during daytime, with a bias of 0.26 E/m 2 /d (0.6%) and a RMS difference of 1.7 E/m 2 /d (4.0%). Using satellite-derived aerosol optical thickness (AOT) and Angström coefficient instead of climatology slightly degraded the results, which was likely due to uncertainties in the aerosol retrievals. A sensitivity study to aerosol properties indicated that climatology may not work in some situations (e.g., episodic dust, pollution, or biomass burning events), suggesting that it is best to use actual aerosol estimates in clear sky conditions. The analysis also revealed that specifying aerosol properties, therefore atmospheric transmittance, from AOT and Angström coefficient, even retrieved from the satellite imagery, may not be sufficient in the presence of absorbing aerosols, especially when loadings are important. Performance was degraded when including situations of clear sky at the time of the MERIS observation but cloudy sky before and/or after overpass, resulting in a bias (overestimation) of 2.8 E/m 2 /d (7.3%) and a RMS difference of 6.0 E/m 2 /d (15.8%). The relatively large overestimation was due to the inability of the OBPG PAR algorithm to detect cloudiness at times other than the time of satellite overpass. The key to improving the daily mean PAR estimates in such situations does not reside so much in improving the radiative transfer treatment or specifying more accurately aerosol properties, but rather in accounting properly for the diurnal variability of cloudiness. To this end, a methodology that utilized Modern Era Retrospective Reanalysis for Research and Applications, Version 2 (MERRA-2) hourly cloud data (fractional coverage, optical thickness) was proposed and tested, reducing the bias to 1.6 E/m 2 /d (4.2%). Improvement was not sufficient in some situations, due to the coarse resolution and uncertainties of the MERRA-2 products, which could not describe properly the cloud properties at the local scale (MERIS pixel). The treatment is applicable to any cloud situation and should be considered in a future version of the of OBPG PAR algorithm. This would require, however, refreshing the standard OBPG PAR products generated as part of the ocean-color processing line according to MERRA-2 data availability.
The 2017–2027 National Academies' Decadal Survey, Thriving on Our Changing Planet, recommended Surface Biology and Geology (SBG) as a “Designated Targeted Observable” (DO). The SBG DO is based on the ...need for capabilities to acquire global, high spatial resolution, visible to shortwave infrared (VSWIR; 380–2500 nm; ~30 m pixel resolution) hyperspectral (imaging spectroscopy) and multispectral midwave and thermal infrared (MWIR: 3–5 μm; TIR: 8–12 μm; ~60 m pixel resolution) measurements with sub-monthly temporal revisits over terrestrial, freshwater, and coastal marine habitats. To address the various mission design needs, an SBG Algorithms Working Group of multidisciplinary researchers has been formed to review and evaluate the algorithms applicable to the SBG DO across a wide range of Earth science disciplines, including terrestrial and aquatic ecology, atmospheric science, geology, and hydrology. Here, we summarize current state-of-the-practice VSWIR and TIR algorithms that use airborne or orbital spectral imaging observations to address the SBG DO priorities identified by the Decadal Survey: (i) terrestrial vegetation physiology, functional traits, and health; (ii) inland and coastal aquatic ecosystems physiology, functional traits, and health; (iii) snow and ice accumulation, melting, and albedo; (iv) active surface composition (eruptions, landslides, evolving landscapes, hazard risks); (v) effects of changing land use on surface energy, water, momentum, and carbon fluxes; and (vi) managing agriculture, natural habitats, water use/quality, and urban development. We review existing algorithms in the following categories: snow/ice, aquatic environments, geology, and terrestrial vegetation, and summarize the community-state-of-practice in each category. This effort synthesizes the findings of more than 130 scientists.
•The 2017 Decadal Survey recommended Surface Biology and Geology mission•Visible to shortwave infrared hyperspectral and multi-band thermal data•Global high resolution measurements at sub-monthly temporal resolution•Applications in snow/ice, aquatic environment, geology, and terrestrial vegetation•We review existing relevant algorithms and community-state-of-practice
•Vector radiative transfer computations were performed for atmosphere-ocean models.•Four models, four wavelengths, two altitudes, and 100 geometries were considered.•Three radiative transfer codes ...were used to validate the accuracy of computations.•Tabulated testbed values are presented for Stokes parameters I, Q, and u.•The accuracy of the testbed tables is at least 10−5 and mostly better than 10−6.
We generate and tabulate reflectance values of the Stokes parameters I, Q, and U of upwelling radiance just above a rough ocean surface and at the top of the atmosphere (TOA) for 100 scattering geometries, four atmosphere-ocean systems, and four wavelengths. The atmosphere-ocean systems increase in complexity from (a) a molecular atmosphere above a rough ocean surface (AOS-I model); to (b) a pure water body below a rough ocean surface (AOS-II model); to (c) a fully-coupled simple atmosphere-ocean system (AOS-III model) containing a molecular atmosphere, rough ocean surface, and pure water; to (d) a fully-coupled complex atmosphere-ocean system (AOS-IV model) that includes scattering by molecules, rough ocean surface, pure water, and hydrosols. Our wavelengths (350, 450, 550, and 650 nm) capture the ultraviolet-visible range. Our tables provide radiative transfer (RT) testbed results for atmosphere-ocean systems with an accuracy that surpasses the measurement accuracy of state-of-the-art polarimeters. To validate the accuracy of these tables we performed computations using three independent RT codes that provide deterministic numerical solutions for the RT equation. The agreement is 10–5 for AOS-IV model, and 10–6 for the other models. The degree of linear polarization computed by these RT codes differs by ≤0.2% for 15 isolated cases of tabulated reflectance values, and by ≤0.1% for all remaining cases. We also provide comparisons with results obtained by a stochastic RT code for AOS-I model. The agreement between the deterministic and stochastic results for this model is 10–5 at TOA, and 10–6 above the ocean surface.
Abstract
In situ observations and output from a numerical model are utilized to examine three dust outbreaks that occurred in the northwestern Sonoran Desert. Via analysis of these events, it is ...shown that trapped waves generated in the lee of an upwind mountain range produced high surface wind speeds along the desert floor and the observed dust storms. Based on analysis of observational and model output, general characteristics of dust outbreaks generated by trapped waves are suggested, including dust-layer depths and concentrations that are dependent upon wave phase and height above the surface, emission and transport associated with the presence of a low-level jet, and wave-generated high wind speeds and thus emission that occurs far downwind of the wave source. Trapped lee waves are ubiquitous in Earth’s atmosphere and thus it is likely that the meteorological aspects of the dust storms examined here are also relevant to understanding dust in other regions. These dust outbreaks occurred near the Salton Sea, an endorheic inland body of water that is rapidly drying due to changes in water-use management. As such, these findings are also relevant in terms of understanding how future changes in size of the Salton Sea will impact dust storms and air quality there.
Significance Statement
Dust storms are ubiquitous in Earth’s atmosphere, yet the physical processes underlying dust emission and subsequent transport are not always understood, in part due to the wide variety of meteorological processes that can generate high winds and dust. Here we use in situ measurements and numerical modeling to demonstrate that vertically trapped atmospheric waves generated by air flowing over a mountain are one such mechanism that can produce dust storms. We suggest several features of these dust outbreaks that are specific to their production by trapped waves. As the study area is a region undergoing rapid environmental change, these results are relevant in terms of predicting future dust there.
Vicarious calibration coefficients (
k
v
) of Second-generation GLobal Imager (SGLI) for ocean color processing were derived using in-situ radiometric buoy measurements from the Marine Optical BuoY ...(MOBY) and the BOUée pour l'acquiSition d'une Série Optique à Long termE (BOUSSOLE). Two aerosol-model look up tables (LUTs) used in the GCOM-C aerosol retrieval algorithm (LUT-A) and in the previous version of ocean color atmospheric correction algorithm (LUT-B) were tested in the procedures to calculate
k
v
and retrieve remote sensing reflectance (
R
rs
) and aerosol optical thickness (AOT). Bias of the processed
R
rs
compared to AERONET-OC
R
rs
was reduced by applying the determined
k
v
(i.e., corrected SGLI radiance = original SGLI radiance/
k
v
). LUT-A yielded smaller AOT bias compared to AERONET-OC AOT; on the other hand, LUT-B gave smaller
R
rs
noise due to gentle slope of the aerosol reflectance even though it caused AOT overestimation. When
k
v
was derived by adjusting to the AOT measurements,
k
v
was about 1.1 by LUT-A and 1.2 by LUT-B in the near-infrared (NIR) channel. However, the
k
v
in the NIR channel was close to 1.0 when AOT and land surface reflectance measurements of Radiometric Calibration Network (RadCalNet) were used. The LUT-A with
k
v
from MOBY and BOUSSOLE are currently adopted for the SGLI standard ocean color processing. Improvement is needed, however, to design an optimal LUT suitable for both aerosol and ocean color purposes.