Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Archives”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 901 records · Page 50

Powerful Radio Sources in the Southern Sky. II. A Swift X-Ray Perspective

We recently constructed the G4Jy-3CRE, a catalog of extragalactic radio sources based on the GLEAM 4-Jy (G4Jy) sample, with the aim of increasing the number of powerful radio galaxies and quasars with similar selection criteria to those of the revised release of the Third Cambridge Catalog (3CR). The G4Jy-3CRE consists of a total of 264 radio sources mainly visible from the Southern Hemisphere. Here, we present an initial X-ray analysis of 89 G4Jy-3CRE radio sources with archival X-ray observations from the Neil Gehrels Swift Observatory. We reduced a total of 624 Swift observations, for about 0.9 Ms of integrated exposure time. We found X-ray counterparts for 59 radio sources belonging to the G4Jy-3CRE, nine of them showing extended X-ray emission. The remaining 30 sources do not show any X-ray emission associated with their radio cores. Our analysis demonstrates that X-ray snapshot observations, even if lacking uniform exposure times, as those carried out with Swift, allow us to (i) verify and/or refine the host galaxy identification; (ii) discover the extended X-ray emission around radio galaxies of the intracluster medium when harbored in galaxy clusters, as the case of G4Jy 1518 and G4Jy 1664; and (iii) detect X-ray radiation arising from their radio lobes, as for G4Jy 1863.

79 ASTRONOMY AND ASTROPHYSICS↗

IceCat-1: The IceCube Event Catalog of Alert Tracks

We present a catalog of likely astrophysical neutrino track-like events from the IceCube Neutrino Observatory. IceCube began reporting likely astrophysical neutrinos in 2016, and this system was updated in 2019. The catalog presented here includes events that were reported in real time since 2019, as well as events identified in archival data samples starting from 2011. We report 275 neutrino events from two selection channels as the first entries in the catalog, the IceCube Event Catalog of Alert Tracks, which will see ongoing extensions with additional alerts. The Gold and Bronze alert channels respectively provide neutrino candidates with a 50% and 30% probability of being astrophysical, on average assuming an astrophysical neutrino power-law energy spectral index of 2.19. For each neutrino alert, we provide the reconstructed energy, direction, false-alarm rate, probability of being astrophysical in origin, and likelihood contours describing the spatial uncertainty in the alert's reconstructed location. We also investigate a directional correlation of these neutrino events with gamma-ray and X-ray catalogs, including 4FGL, 3HWC, TeVCat, and Swift-BAT.

79 ASTRONOMY AND ASTROPHYSICS↗

Discovery of a Gamma-Ray Black Widow Pulsar by GPU-accelerated Einstein@Home

We report the discovery of 1.97 ms period gamma-ray pulsations from the 75 minute orbital-period binary pulsar now named PSR J1653-0158. The associated Fermi Large Area Telescope gamma-ray source 4FGL J1653.6-0158 has long been expected to harbor a binary millisecond pulsar. Despite the pulsar-like gamma-ray spectrum and candidate optical/X-ray associations—whose periodic brightness modulations suggested an orbit—no radio pulsations had been found in many searches. The pulsar was discovered by directly searching the gamma-ray data using the GPU-accelerated Einstein@Home distributed volunteer computing system. The multidimensional parameter space was bounded by positional and orbital constraints obtained from the optical counterpart. More sensitive analyses of archival and new radio data using knowledge of the pulsar timing solution yield very stringent upper limits on radio emission. Any radio emission is thus either exceptionally weak, or eclipsed for a large fraction of the time. The pulsar has one of the three lowest inferred surface magnetic-field strengths of any known pulsar with B surf ≈ 4 × 10 7 G. The resulting mass function, combined with models of the companion star's optical light curve and spectra, suggests a pulsar mass gsim 2 M ⊙ . The companion is lightweight with mass ~0.01 M ⊙ , and the orbital period is the shortest known for any rotation-powered binary pulsar. This discovery demonstrates the Fermi Large Area Telescope's potential to discover extreme pulsars that would otherwise remain undetected.

79 ASTRONOMY AND ASTROPHYSICS↗

Tucana B: A Potentially Isolated and Quenched Ultra-faint Dwarf Galaxy at $D$ ≈ 1.4 Mpc*

We report the discovery of Tucana B, an isolated ultra-faint dwarf galaxy at a distance of $D$ = 1.4 Mpc. Tucana B was found during a search for ultra-faint satellite companions to the known dwarfs in the outskirts of the Local Group, although its sky position and distance indicate the nearest galaxy to be ~500 kpc distant. Deep ground-based imaging resolves Tucana B into stars, and it displays a sparse red giant branch consistent with an old, metal-poor stellar population analogous to that seen in the ultra-faint dwarf galaxies of the Milky Way, albeit at fainter apparent magnitudes. Tucana B has a half-light radius of 80 ± 40 pc and an absolute magnitude of ${M}_{V}=-{6.9}_{-0.6}^{+0.5}$ mag (${L}_{V}=({5}_{-2}^{+4})\times {10}^{4}$ L⊙), which is again comparable to the Milky Way's ultra-faint satellites. There is no evidence for a population of young stars, either in the optical color–magnitude diagram or in GALEX archival ultraviolet imaging, with the GALEX data indicating $\mathrm{log}({\mathrm{SFR}}_{\mathrm{NUV}}/{M}_{\odot }\,{\mathrm{yr}}^{-1})\lt -5.4$ for star formation on ≲100 Myr timescales. Given its isolation and physical properties, Tucana B may be a definitive example of an ultra-faint dwarf that has been quenched by reionization, providing strong confirmation of a key driver of galaxy formation and evolution at the lowest mass scales. It also signals a new era of ultra-faint dwarf galaxy discovery at the extreme edges of the Local Group.

79 ASTRONOMY AND ASTROPHYSICS↗

Compact Steep Spectrum Radio Sources with Enhanced Star Formation Are Smaller Than 10 kpc

Compact steep spectrum (CSS) radio sources are active galactic nuclei (AGN) that have radio jets propagating only on galactic scales, defined as having projected linear size (LS) of up to 20 kpc. CSS sources are generally hosted by massive early-type galaxies with little ongoing star formation; however, a small fraction are known to have enhanced star formation. Using archival data from the Faint Images of the Radio Sky at Twenty cm survey, the Very Large Array Sky Survey, and the Sloan Digital Sky Survey, we identify a volume-limited sample of 166 CSS sources at z < 0.2 with L 1.4 GHz > 10 24 W Hz -1 . Comparing the star formation rates and linear sizes of these CSS sources, we find that the ≈14% of CSS sources with specific star formation rates above 0.01 Gyr -1 all have LS < 10 kpc. We discuss the possible mechanisms driving this result, concluding that it is likely the excess star formation in these sources occurred in multiple bursts and ceased prior to the AGN jet being triggered.

79 ASTRONOMY AND ASTROPHYSICS↗

Enabling Early Transient Discovery in LSST via Difference Imaging with DECam

We present SLIDE, a pipeline that enables transient discovery in data from the Vera C. Rubin Observatory’s Legacy Survey of Space and Time (LSST), using archival images from the Dark Energy Camera as templates for difference imaging. We apply this pipeline to the recently released Data Preview 1 (DP1; the first public release of Rubin commissioning data) and search for transients in the resulting difference images. The image subtraction, photometry extraction, and transient detection are all performed on the Rubin Science Platform. We demonstrate that SLIDE effectively extracts clean photometry by circumventing poor or missing LSST templates. We identified 29 previously unreported transients, 12 of which would not have been detected based on the DP1 DiaObject catalog. SLIDE will be especially useful for transient analysis in the early years of LSST, when template coverage will be largely incomplete or when templates may be contaminated by transients present at the time of acquisition. We present multiband light curves for a sample of known transients, along with new transient candidates identified through our search. Finally, we discuss the prospects of applying this pipeline during the main LSST survey. Our pipeline is broadly applicable and will support studies of all transients with slowly evolving phases.

Dong, Yize 一泽董 [Harvard-Smithsonian Center for Ast↗

Atomic Iron and Nickel in the Coma of C/1996 B2 (Hyakutake): Production Rates, Emission Mechanisms, and Possible Parents

Abstract Two papers recently reported the detection of gaseous nickel and iron in the comae of over 20 comets from observations collected over two decades, including interstellar comet 2I/Borisov. To evaluate the state of the laboratory data in support of these identifications, we reanalyzed archived spectra of comet C/1996 B2 (Hyakutake), one of the nearest and brightest comets of the past century, using a combined experimental and computational approach. We developed a new, many-level fluorescence model that indicates that the fluorescence emissions of Fe I and Ni I vary greatly with heliocentric velocity. Combining this model with laboratory spectra of an Fe-Ni plasma, we identified 22 lines of Fe I and 14 lines of Ni I in the spectrum of Hyakutake. Using Haser models, we estimate the nickel and iron production rates as Q Ni = (2.6–4.1) × 10 22 s −1 and Q Fe = (0.4–2.8) × 10 23 s −1 . From derived column densities, the Ni/Fe abundance ratio log 10 [Ni/Fe] = −0.15 ± 0.07 deviates significantly from solar abundance ratios, and it is consistent with the ratios observed in solar system comets. Possible production and emission mechanisms are analyzed in the context of existing laboratory measurements. Based on the observed spatial distributions, excellent fluorescence model agreement, and Ni/Fe ratio, our findings support an origin consisting of a short-lived unknown parent followed by fluorescence emission. Our models suggest that the strong heliocentric velocity dependence of the fluorescence efficiencies can provide a meaningful test of the physical process responsible for the Fe I and Ni I emission.

Bromley, S. J. (ORCID:0000000321108152)↗

Long- and short-term temporal variability in cloud condensation nuclei spectra over a wide supersaturation range in the Southern Great Plains site

Abstract. When aerosol particles seed the formation of liquid water droplets in the atmosphere, they are called cloud condensation nuclei (CCN). Different aerosols will act as CCN under different degrees of water supersaturation (relative humidity above 100 %), depending on their size and composition. In this work, we build and analyze a best-estimate CCN spectrum product, tabulated at ∼ 45 min resolution, generated using high quality data from seven independent instruments at the U.S. Department of Energy Atmospheric Radiation Measurement (ARM) Southern Great Plains site. The data product spans a large supersaturation range, from 0.0001 % to ∼ 30 %, and time period of 5 years, from 2009–2013, and is available on the ARM data archive. We leverage this added statistical power to examine relationships that are unclear in smaller datasets. Our analysis is performed in three main areas. First, probability distributions of many aerosol and CCN metrics are found to exhibit skewed log-normal distribution shapes. Second, clustering analyses of CCN spectra reveal that the primary drivers of CCN differences are aerosol number size distributions, rather than hygroscopicity or composition, especially at supersaturations above 0.2 %, while also allowing for a simplified understanding of seasonal and diurnal variations in CCN behavior. The predictive ability of using limited hygroscopicity data with accurate number size distributions to estimate CCN spectra is investigated, and the uncertainties of this approach are estimated. Third, the dynamics of CCN spectral clusters and concentrations are examined with cross-correlation and autocorrelation analyses. We find that CCN concentrations change rapidly on the timescale of 1–3 h, with some conservation beyond that which is greatest for the lower supersaturation region of the spectrum.

54 ENVIRONMENTAL SCIENCES↗

The De-Icing Comparison Experiment (D-ICE): a study of broadband radiometric measurements under icing conditions in the Arctic

Abstract. Surface-based measurements of broadband shortwave (solar) and longwave (infrared) radiative fluxes using thermopile radiometers are made regularly around the globe for scientific and operational environmental monitoring. The occurrence of ice on sensor windows in cold environments – whether snow, rime, or frost – is a common problem that is difficult to prevent as well as difficult to correct in post-processing. The Baseline Surface Radiation Network (BSRN) community recognizes radiometer icing as a major outstanding measurement uncertainty. Towards constraining this uncertainty, the De-Icing Comparison Experiment (D-ICE) was carried out at the NOAA Atmospheric Baseline Observatory in Utqiaġvik (formerly Barrow), Alaska, from August 2017 to July 2018. The purpose of D-ICE was to evaluate existing ventilation and heating technologies developed to mitigate radiometer icing. D-ICE consisted of 20 pyranometers and 5 pyrgeometers operating in various ventilator housings alongside operational systems that are part of NOAA's Barrow BSRN station and the US Department of Energy Atmospheric Radiation Measurement (ARM) program North Slope of Alaska and Oliktok Point observatories. To detect icing, radiometers were monitored continuously using cameras, with a total of more than 1 million images of radiometer domes archived. Ventilator and ventilator–heater performance overall was skillful with the average of the systems mitigating ice formation 77 % (many >90 %) of the time during which icing conditions were present. Ventilators without heating elements were also effective and capable of providing heat through roughly equal contributions of waste energy from the ventilator fan and adiabatic heating downstream of the fan. This provided ∼0.6 ∘C of warming, enough to subsaturate the air up to a relative humidity (with respect to ice) of ∼105 %. Because the mitigation technologies performed well, a near complete record of verified ice-free radiometric fluxes was assembled for the duration of the campaign. This well-characterized data set is suitable for model evaluation, in particular for the Year of Polar Prediction (YOPP) first Special Observing Period (SOP1). We used the data set to calculate short- and long-term biases in iced sensors, finding that biases can be up to +60 W m−2 (longwave) and −211 to +188 W m−2 (shortwave). However, because of the frequency of icing, mitigation of ice by ventilators, cloud conditions, and the timing of icing relative to available sunlight, the biases in the monthly means were generally less than the aggregate uncertainty attributed to other conventional sources in both the shortwave and longwave.

54 ENVIRONMENTAL SCIENCES↗

The development of rainfall retrievals from radar at Darwin

Abstract. The U.S. Department of Energy Atmospheric Radiation Measurement program Tropical Western Pacific site hosted a C-band polarization (CPOL) radar in Darwin, Australia. It provides 2 decades of tropical rainfall characteristics useful for validating global circulation models. Rainfall retrievals from radar assume characteristics about the droplet size distribution (DSD) that vary significantly. To minimize the uncertainty associated with DSD variability, new radar rainfall techniques use dual polarization and specific attenuation estimates. This study challenges the applicability of several specific attenuation and dual-polarization-based rainfall estimators in tropical settings using a 4-year archive of Darwin disdrometer datasets in conjunction with CPOL observations. This assessment is based on three metrics: statistical uncertainty estimates, principal component analysis (PCA), and comparisons of various retrievals from CPOL data. The PCA shows that the variability in R can be consistently attributed to reflectivity, but dependence on dual-polarization quantities was wavelength dependent for 1 10mmh-1. Rainfall estimates during these conditions primarily originate from deep convective clouds with median drop diameters greater than 1.5 mm. An uncertainty analysis and intercomparison with CPOL show that a Colorado State University blended technique for tropical oceans, with modified estimators developed from video disdrometer observations, is most appropriate for use in all cases, such as when 1 10mmh-1 (deeper convective rain).

54 ENVIRONMENTAL SCIENCES↗

Moderate spectral resolution solar irradiance measurements, aerosol optical depth, and solar transmission, from 360 to 1070 nm, using the refurbished rotating shadow band spectroradiometer (RSS)

Abstract. This paper reports on a third-generation rotating shadow band spectroradiometer (RSS) used to measure global and diffuse horizontal plus direct normal irradiances and transmissions at 1002 wavelengths between 360 and 1070 nm. The prism-dispersed spectral data are from the Atmospheric Radiation Measurement (ARM) Southern Great Plains site in north-central Oklahoma (36.605∘ N, 97.486∘ W) and cover dates between August 2009 and February 2014. The refurbished RSS isolates the detector in a vacuum chamber with pressures near 10−7 torr. This prevents the deposition of outgassed vapors from the interior of the spectrometer shell on the cooled detector that affected the operation of the first commercial RSS. Methods for (1) ensuring the correct wavelength registration of the data and (2) deriving extraterrestrial responses over the entire spectrum, including throughout strong water vapor and oxygen bands, are described. The resulting data produced are archived as ARM data records and include cloud-screened aerosol optical depths, spectral irradiances and direct normal solar transmission, as well as normalized diffuse and global irradiances.

54 ENVIRONMENTAL SCIENCES↗

Calibrating radar wind profiler reflectivity factor using surface disdrometer observations

Abstract. This study uses surface disdrometer reflectivity factor estimates to calibrate the vertical and off-vertical pointing radar beams produced by an ultra high frequency (UHF) band radar wind profiler (RWP) deployed at the US Department of Energy (DOE) Atmospheric Radiation Measurement (ARM) program Southern Great Plains (SGP) Central Facility in northern Oklahoma from April 2011 through July 2019. The methodology consists of five steps. First, the recorded Doppler velocity power spectra are adjusted to account for Nyquist velocity aliasing and coherent integration filtering effects. Second, the spectrum moments are calculated. The third step increases the signal-to-noise ratio (SNR) due to inflated noise power estimates during convective rain events that cause SNR to be biased low. The fourth step determines the RWP calibration constant for one radar beam (called the “reference” beam) by comparing uncalibrated RWP reflectivity factors at 500 m above the ground to 1 min resolution surface disdrometer reflectivity factors. The last step uses the calibrated reference beam reflectivity factor to calibrate the other radar beams during precipitation. There are two key findings. The RWP sensitivity decreased by approximately 3 to 4 dB yr−1 as the hardware aged. This drift was slow enough that the reference calibration constant can be estimated over 3-month intervals using episodic rain events. The calibrated moments are available on the DOE ARM data archive, and the Python processing code is available on public repositories.

54 ENVIRONMENTAL SCIENCES↗

Regridding uncertainty for statistical downscaling of solar radiation

Abstract. Initial steps in statistical downscaling involve being able to compare observed data from regional climate models (RCMs). This prediction requires (1) regridding RCM outputs from their native grids and at differing spatial resolutions to a common grid in order to be comparable to observed data and (2) bias correcting RCM data, for example via quantile mapping, for future modeling and analysis. The uncertainty associated with (1) is not always considered for downstream operations in (2). This work examines this uncertainty, which is not often made available to the user of a regridded data product. This analysis is applied to RCM solar radiation data from the NA-CORDEX (North American Coordinated Regional Climate Downscaling Experiment) data archive and observed data from the National Solar Radiation Database housed at the National Renewable Energy Lab. A case study of the mentioned methods over California is presented.

14 SOLAR ENERGY↗

Data generated during the 2018 LAPSE-RATE campaign: an introduction and overview

Unmanned aircraft systems (UASs) offer innovative capabilities for providing new perspectives on the atmosphere, and therefore atmospheric scientists are rapidly expanding their use, particularly for studying the planetary boundary layer. In support of this expansion, from 14 to 20 July 2018 the International Society for Atmospheric Research using Remotely piloted Aircraft (ISARRA) hosted a community flight week, dubbed the Lower Atmospheric Profiling Studies at Elevation – a Remotely-piloted Aircraft Team Experiment (LAPSE-RATE; de Boer et al., 2020a). This field campaign spanned a 1-week deployment to Colorado's San Luis Valley, involving over 100 students, scientists, engineers, pilots, and outreach coordinators. These groups conducted intensive field operations using unmanned aircraft and ground-based assets to develop comprehensive datasets spanning a variety of scientific objectives, including a total of nearly 1300 research flights totaling over 250 flight hours. This article introduces this campaign and lays the groundwork for a special issue on the LAPSE-RATE project. The remainder of the special issue provides detailed overviews of the datasets collected and the platforms used to collect them. All of the datasets covered by this special issue have been uploaded to a LAPSE-RATE community set up at the Zenodo data archive (https://zenodo.org/communities/lapse-rate/, last access: 3 December 2020).

47 OTHER INSTRUMENTATION↗

Measurements from mobile surface vehicles during the Lower Atmospheric Profiling Studies at Elevation – a Remotely-piloted Aircraft Team Experiment (LAPSE-RATE)

Abstract. Between 14 and 20 July 2018, small unmanned aircraft systems (UASs) were deployed to the San Luis Valley of Colorado (USA) alongside surface-based remote sensors, in situ sensors, and radiosonde systems as part of the Lower Atmospheric Profiling Studies at Elevation – a Remotely-piloted Aircraft Team Experiment (LAPSE-RATE). The measurements collected as part of LAPSE-RATE targeted quantities related to enhancing our understanding of boundary layer structure, cloud and aerosol properties and surface–atmosphere exchange and provide detailed information to support model evaluation and improvement work. Additionally, intensive intercomparison between the different unmanned aircraft platforms was completed. The current paper describes the observations obtained using three different types of surface-based mobile observing vehicles. These included the University of Colorado Mobile UAS Research Collaboratory (MURC), the National Oceanic and Atmospheric Administration National Severe Storms Laboratory Mobile Mesonet, and two University of Nebraska Combined Mesonet and Tracker (CoMeT) vehicles. Over the 1-week campaign, a total of 143 h of data were collected using this combination of vehicles. The data from these coordinated activities provide detailed perspectives on the spatial variability of atmospheric state parameters (air temperature, humidity, pressure, and wind) throughout the northern half of the San Luis Valley. These datasets have been checked for quality and published to the Zenodo data archive under a specific “community” setup for LAPSE-RATE (https://zenodo.org/communities/lapse-rate/, last access: 21 January 2021) and are accessible at no cost by all registered users. The primary dataset DOIs are https://doi.org/10.5281/zenodo.3814765 (CU MURC measurements; de Boer et al., 2020d), https://doi.org/10.5281/zenodo.3738175 (NSSL MM measurements; Waugh, 2020), and https://doi.org/10.5281/zenodo.3838724 (UNL CoMeT measurements; Houston and Erwin, 2020).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

SoDaH: the SOils DAta Harmonization database, an open-source synthesis of soil data from research networks, version 1.0

Data collected from research networks present opportunities to test theories and develop models about factors responsible for the long-term persistence and vulnerability of soil organic matter (SOM). Synthesizing datasets collected by different research networks presents opportunities to expand the ecological gradients and scientific breadth of information available for inquiry. Synthesizing these data is challenging, especially considering the legacy of soil data that have already been collected and an expansion of new network science initiatives. To facilitate this effort, here we present the SOils DAta Harmonization database (SoDaH; https://lter.github.io/som-website, last access: 22 December 2020), a flexible database designed to harmonize diverse SOM datasets from multiple research networks. SoDaH is built on several network science efforts in the United States, but the tools built for SoDaH aim to provide an open-access resource to facilitate synthesis of soil carbon data. Moreover, SoDaH allows for individual locations to contribute results from experimental manipulations, repeated measurements from long-term studies, and local- to regional-scale gradients across ecosystems or landscapes. Finally, we also provide data visualization and analysis tools that can be used to query and analyze the aggregated database. The SoDaH v1.0 dataset is archived and available at https://doi.org/10.6073/pasta/9733f6b6d2ffd12bf126dc36a763e0b4 (Wieder et al., 2020).

54 ENVIRONMENTAL SCIENCES↗

University of Colorado and Black Swift Technologies RPAS-based measurements of the lower atmosphere during LAPSE-RATE

Abstract. Between 14 and 20 July 2018, small remotely piloted aircraft systems (RPASs) were deployed to the San Luis Valley of Colorado (USA) together with a variety of surface-based remote and in situ sensors as well as radiosonde systems as part of the Lower Atmospheric Profiling Studies at Elevation – a Remotely-piloted Aircraft Team Experiment (LAPSE-RATE). The observations from LAPSE-RATE were aimed at improving our understanding of boundary layer structure, cloud and aerosol properties, and surface–atmosphere exchange and provide detailed information to support model evaluation and improvement work. The current paper describes the observations obtained using four different types of RPASs deployed by the University of Colorado Boulder and Black Swift Technologies. These included the DataHawk2, the Talon and the TTwistor (University of Colorado), and the S1 (Black Swift Technologies). Together, these aircraft collected over 30 h of data throughout the northern half of the San Luis Valley, sampling altitudes between the surface and 914 m a.g.l. Data from these platforms are publicly available through the Zenodo archive and are co-located with other LAPSE-RATE data as part of the Zenodo LAPSE-RATE community (https://zenodo.org/communities/lapse-rate/, last access: 27 May 2021). The primary DOIs for these datasets are https://doi.org/10.5281/zenodo.3891620 (DataHawk2, de Boer et al., 2020a, e), https://doi.org/10.5281/zenodo.4096451 (Talon, de Boer et al., 2020d), https://doi.org/10.5281/zenodo.4110626 (TTwistor, de Boer et al., 2020b), and https://doi.org/10.5281/zenodo.3861831 (S1, Elston and Stachura, 2020).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A daily, 250 m and real-time gross primary productivity product (2000–present) covering the contiguous United States

Abstract. Gross primary productivity (GPP) quantifies the amount of carbon dioxide (CO2) fixed by plants through photosynthesis. Although as a key quantity of terrestrial ecosystems, there is a lack of high-spatial-and-temporal-resolution, real-time and observation-based GPP products. To address this critical gap, here we leverage a state-of-the-art vegetation index, near-infrared reflectance of vegetation (NIRV), along with accurate photosynthetically active radiation (PAR), to produce a SatelLite Only Photosynthesis Estimation (SLOPE) GPP product for the contiguous United States (CONUS). Compared to existing GPP products, the proposed SLOPE product is advanced in its spatial resolution (250 m versus >500 m), temporal resolution (daily versus 8 d), instantaneity (latency of 1 d versus >2 weeks) and quantitative uncertainty (on a per-pixel and daily basis versus no uncertainty information available). These characteristics are achieved because of several technical innovations employed in this study: (1) SLOPE couples machine learning models with MODIS atmosphere and land products to accurately estimate PAR. (2) SLOPE couples highly efficient and pragmatic gap-filling and filtering algorithms with surface reflectance acquired by both Terra and Aqua MODIS satellites to derive a soil-adjusted NIRV (SANIRV) dataset. (3) SLOPE couples a temporal pattern recognition approach with a long-term Cropland Data Layer (CDL) product to predict dynamic C4 crop fraction. Through developing a parsimonious model with only two slope parameters, the proposed SLOPE product explains 85 % of the spatial and temporal variations in GPP acquired from 49 AmeriFlux eddy-covariance sites (324 site years), with a root-mean-square error (RMSE) of 1.63 gC m−2 d−1. The median R2 over C3 and C4 crop sites reaches 0.87 and 0.94, respectively, indicating great potentials for monitoring crops, in particular bioenergy crops, at the field level. With such a satisfactory performance and its distinct characteristics in spatiotemporal resolution and instantaneity, the proposed SLOPE GPP product is promising for biological and environmental research, carbon cycle research, and a broad range of real-time applications at the regional scale. The archived dataset is available at https://doi.org/10.3334/ORNLDAAC/1786 (download page: https://daac.ornl.gov/daacdata/cms/SLOPE_GPP_CONUS/data/, last access: 20 January 2021) (Jiang and Guan, 2020), and the real-time dataset is available upon request.

54 ENVIRONMENTAL SCIENCES↗