We estimate global terrestrial gross primary production (GPP) based on models that use satellite data within a simplified light-use efficiency framework that does not rely upon other meteorological inputs. Satellite-based geometry-adjusted reflectances are from the MODerate-resolution Imaging Spectroradiometer (MODIS) and provide information about vegetation structure and chlorophyll content at both high temporal (daily to monthly) and spatial (∼1 km) resolution. We use satellite-derived solar-induced fluorescence (SIF) to identify regions of high productivity crops and also evaluate the use of downscaled SIF to estimate GPP. We calibrate a set of our satellite-based models with GPP estimates from a subset of distributed eddy covariance flux towers (FLUXNET 2015). The results of the trained models are evaluated using an independent subset of FLUXNET 2015 GPP data. We show that variations in light-use efficiency (LUE) with incident PAR are important and can be easily incorporated into the models. Unlike many LUE-based models, our satellite-based GPP estimates do not use an explicit parameterization of LUE that reduces its value from the potential maximum under limiting conditions such as temperature and water stress. Even without the parameterized downward regulation, our simplified models are shown to perform as well as or better than state-of-the-art satellite data-driven products that incorporate such parameterizations. A significant fraction of both spatial and temporal variability in GPP across plant functional types can be accounted for using our satellite-based models. Our results provide an annual GPP value of ∼140 Pg C year ${}^{-1}$ for 2007 that is within the range of a compilation of observation-based, model, and hybrid results, but is higher than some previous satellite observation-based estimates.

While solar-induced fluorescence (SIF) shows promise as a remotely-sensed measurement directly related to photosynthesis, interpretation and validation of satellite-based SIF retrievals remains a challenge. SIF is influenced by the fraction of absorbed photosynthetically-active radiation at the canopy level that depends upon illumination geometry as well as the escape of SIF through the canopy that depends upon the viewing geometry. Several approaches to estimate the effects of sun-sensor geometry on satellite-based SIF have been proposed, and some have been implemented, most relying upon satellite reflectance measurements and/or other ancillary data sets. These approaches, designed to ultimately estimate intrinsic or physiological components of SIF related to photosynthesis, have not generally been applied globally to satellite measurements. Here, we examine in detail how SIF and related reflectance-based indices from wide swath polar orbiting satellites in low Earth orbit vary systematically due to the host satellite orbital characteristics. We compare SIF and reflectance-based parameters from the Global Ozone Mapping Experiment 2 (GOME-2) on the MetOp-B platform and from the TROPOspheric Monitoring Instrument (TROPOMI) on the Sentinel 5 Precursor satellite with a focus on high northern latitudes in summer where observations at similar geometries and local times occur. We show that GOME-2 and TROPOMI SIF observations agree nearly to within estimated uncertainties when they are compared at similar observing geometries. We show that the cross-track dependence of SIF normalized by PAR and related reflectance-based indices are highly correlated for dense canopies, but diverge substantially as the vegetation within a field-of-view becomes more sparse. This has implications for approaches that utilize reflectance measurements to help account for SIF geometrical dependences in satellite measurements. To further help interpret the GOME-2 and TROPOMI SIF observations, we simulated cross-track dependences of PAR normalized SIF and reflectance-based indices with the one dimensional Soil-Canopy Observation Photosynthesis and Energy fluxes (SCOPE) canopy radiative transfer model at sun–satellite geometries that occur across the wide swaths of these instruments and examine the geometrical dependencies of the various components (e.g., fraction of absorbed PAR, SIF yield, and escape of SIF from the canopy) of the observed SIF signal. The simulations show that most of the cross-track variations in SIF result from the escape of SIF through the scattering canopy and not the illumination.

The increasing frequency and intensity of climate extremes and complex ecosystem responses motivate the need for integrated observational studies at low latency to determine biosphere responses and carbon‐climate feedbacks. Here, we develop a satellite‐based rapid attribution workflow and demonstrate its use at a 1–2‐month latency to attribute drivers of the carbon cycle feedbacks during the 2020–2021 Western US drought and heatwave. In the first half of 2021, concurrent negative photosynthesis anomalies and large positive column CO₂ anomalies were detected with satellites. Using a simple atmospheric mass balance approach, we estimate a surface carbon efflux anomaly of 132 TgC in June 2021, a magnitude corroborated independently with a dynamic global vegetation model. Integrated satellite observations of hydrologic processes, representing the soil–plant–atmosphere continuum (SPAC), show that these surface carbon flux anomalies are largely due to substantial reductions in photosynthesis because of a spatially widespread moisture‐deficit propagation through the SPAC between 2020 and 2021. A causal model indicates deep soil moisture stores partially drove photosynthesis, maintaining its values in 2020 and driving its declines throughout 2021. The causal model also suggests legacy effects may have amplified photosynthesis deficits in 2021 beyond the direct effects of environmental forcing. The integrated, observation framework presented here provides a valuable first assessment of a biosphere extreme response and an independent testbed for improving drought propagation and mechanisms in models. The rapid identification of extreme carbon anomalies and hotspots can also aid mitigation and adaptation decisions.

Space-based quantitative passive optical remote sensing of the Earth's surface typically involves the detection and elimination of cloud-contaminated pixels as an initial processing step. We explore a fundamentally different approach; we use machine learning with cloud contaminated satellite multi-spectral data to estimate underlying terrestrial surface reflectances at red, green, and blue (RGB) wavelengths. The NN reproduces land RGB reflectances with high fidelity even in scenes with moderate to high cloud optical thicknesses. This implies that spectral features of the Earth's surface can be detected and distinguished in the presence of clouds, even when they are partially obscured by clouds; the NN is able to separate the spectral fingerprint of the Earth's surface from that of the clouds, aerosols, gaseous absorption, and Rayleigh scattering, provided that there are adequately different spectral features and that the clouds are not completely opaque. Once trained, the NN enables rapid estimates of RGB reflectances with little computational cost. Aside from the training data, there is no requirement of prior information regarding the land surface spectral reflectance, nor is there need for radiative transfer calculations. We test different wavelength windows for reconstruction of surface reflectance. This work provides an initial example of a general approach that has many potential applications in land and ocean remote sensing as well as other practical uses such as in search and rescue, precision agriculture, and change detection.

We present a new method to infer nitrogen oxides (NO_x) emissions and lifetimes based on tropospheric nitrogen dioxide (NO₂) observations together with reanalysis wind fields for cities located in polluted backgrounds. Since the accuracy of the method is difficult to assess due to lack of “true values” that can be used as a benchmark, we apply the method to synthetic NO₂ observations derived from the NASA-Unified Weather Research and Forecasting (NU-WRF) model at a high horizontal spatial resolution of 4 km × 4 km for cities over the continental US. We compare the inferred emissions and lifetimes with the values given by the NU-WRF model to evaluate the method. The method is applicable to 26 US cities. The derived results are generally in good agreement with the values given by the model, with the relative differences of 2 % ± 17 % (mean ± standard deviation) and 15 % ± 25 % for lifetimes and emissions, respectively. Our investigation suggests that the use of wind data prior to satellite overpass time improves the performance of the method. The correlation coefficients between inferred and NU-WRF lifetimes increase from 0.56 to 0.79 and for emissions increase from 0.88 to 0.96 when comparing results based on wind fields sampled simultaneously with satellite observations and averaged over 9 hours data prior to satellite observations, respectively. We estimate that uncertainties in NO_x lifetime and emissions arising from the method are approximately 15 % and 20 %, respectively, for typical (US) cities. We expect this new method to be applicable to NO₂ observations from the TROPOspheric Monitoring Instrument (TROPOMI) and geostationary satellites, such as Geostationary Environment Monitoring Spectrometer (GEMS) or the Tropospheric Emissions: Monitoring Pollution (TEMPO) instrument, to estimate urban NO_x emissions and lifetimes globally.