Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “bias correction”

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 433 records · Page 24
All sourcesNASA NTRSDOE OSTI
NASA NTRS2018

WRF Hindcasts of Cold Front Passages over the ARM Eastern North Atlantic Site: A Sensitivity Study

The present study explores the ability of the Weather Research and Forecasting (WRF) Model to accurately reproduce the passage of extratropical cold fronts at the DOE-ARM Eastern North Atlantic (ENA) observation site on the Azores Islands. An analysis of three case studies is performed in which the impact of the WRF domain size, position of the model boundary relative to the ENA site, grid spacing, and spectral nudging conditions are explored. The results from these case studies indicate that model biases in the timing and duration of cold front passages change with the distance between the model domain boundary and the ENA site. For these three cases, if the western model boundary is further than 1500 km from the site, the front becomes too meridional and fails to reach the site, making 1000 km or 1500 km the optimal distances. In contrast, integrations with small distances (e.g., 500 km) between the site and domain boundaries have inadequate spatial spin-up, i.e. the domain is too small for the model to properly stabilize. For all 3 cases, regardless of domain size, the model has biases in its upper-level circulation that impact the position and timing of the front. However, this issue is most serious for 4000 km2 domains and larger. For these domains, prolonged spectral nudging can correct cold front biases. As such, this analysis provides a framework to optimize the WRF model configuration necessary for a realistic hindcast of a cold front passage at a fixed location centered in a domain as large as computationally possible.

Faycal Lamraoui↗
NASA NTRS1993

What do you gain from deconvolution? - Observing faint galaxies with the Hubble Space Telescope Wide Field Camera

We describe experiments with deconvolutions of simulations of deep HST Wide Field Camera images containing faint, compact galaxies to determine under what circumstances there is a quantitative advantage to image deconvolution, and explore whether it is (1) helpful for distinguishing between stars and compact galaxies, or between spiral and elliptical galaxies, and whether it (2) improves the accuracy with which characteristic radii and integrated magnitudes may be determined. The Maximum Entropy and Richardson-Lucy deconvolution algorithms give the same results. For medium and low S/N images, deconvolution does not significantly improve our ability to distinguish between faint stars and compact galaxies, nor between spiral and elliptical galaxies. Measurements from both raw and deconvolved images are biased and must be corrected; it is easier to quantify and remove the biases for cases that have not been deconvolved. We find no benefit from deconvolution for measuring luminosity profiles, but these results are limited to low S/N images of very compact (often undersampled) galaxies.

Schade, David J.↗
NASA NTRS2001

Application of Pressure-Based Wall Correction Methods to Two NASA Langley Wind Tunnels

This paper is a description and status report on the implementation and application of the WICS wall interference method to the National Transonic Facility (NTF) and the 14 x 22-ft subsonic wind tunnel at the NASA Langley Research Center. The method calculates free-air corrections to the measured parameters and aerodynamic coefficients for full span and semispan models when the tunnels are in the solid-wall configuration. From a data quality point of view, these corrections remove predictable bias errors in the measurement due to the presence of the tunnel walls. At the NTF, the method is operational in the off-line and on-line modes, with three tests already computed for wall corrections. At the 14 x 22-ft tunnel, initial implementation has been done based on a test on a full span wing. This facility is currently scheduled for an upgrade to its wall pressure measurement system. With the addition of new wall orifices and other instrumentation upgrades, a significant improvement in the wall correction accuracy is expected.

Iyer, V.↗
NASA NTRS2019

Performance of Landsat-8 and Sentinel-2 Surface Reflectance Products for River Remote Sensing Retrievals of Chlorophyll-A and Turbidity

Rivers and other freshwater systems play a crucial role in ecosystems, industry, transportation and agriculture. Despite the more than 40 years of inland water observations made possible by optical remote sensing, a standardized reflectance product for inland waters is yet forthcoming. The aim of this work is to compare the standard USGS land surface reflectance product to two Landsat-8 and Sentinel-2 aquatic remote sensing reflectance products over the Amazon, Columbia and Mississippi rivers. Landsat-8 reflectance products from all three routines are then evaluated for their comparative performance in retrieving chlorophyll-a and turbidity in reference to shipborne, underway in situ validation measurements. The land surface product shows the best agreement (4 percent Mean Absolute Percent Difference) with field measurements of radiometry collected on the Amazon River and generates 36 percent higher reflectance values in the visible bands compared to aquatic methods (ACOLITE (Atmospheric Correction for OLI (Operational Land Imager) 'lite') and SeaDAS (Sea-viewing Wide Field-of-View Sensor (SeaWiFS) Data Analysis System)) with larger differences between land and aquatic products observed in Sentinel-2 (0.01 per steraradian) compared to Landsat-8 (0.001 per steraradian). Choice of atmospheric correction routine can bias Landsat-8 retrievals of chlorophyll-a and turbidity by as much as 59 percent and 35 percent respectively. Using a more restrictive time window for matching in situ and satellite imagery can reduce differences by 5–31 percent depending on correction technique. This work highlights the challenges of satellite retrievals over rivers and underscores the need for future optical and biogeochemical research aimed at improving our understanding of the absorbing and scattering properties of river water and their relationships to remote sensing reflectance.

Landsat-8↗
NASA NTRS1994

Differential sea-state bias: A case study using TOPEX/POSEIDON data

We used selected data from the NASA altimeter TOPEX/POSEIDON to calculate differences in range measured by the C and Ku-band altimeters when the satellite overflew 5 to 15 m waves late at night. The range difference is due to free electrons in the ionosphere and to errors in sea-state bias. For the selected data the ionospheric influence on Ku range is less than 2 cm. Any difference in range over short horizontal distances is due only to a small along-track variability of the ionosphere and to errors in calculating the differential sea-state bias. We find that there is a barely detectable error in the bias in the geophysical data records. The wave-induced error in the ionospheric correction is less than 0.2% of significant wave height. The equivalent error in differential range is less than 1% of wave height. Errors in the differential sea-state bias calculations appear to be small even for extreme wave heights that greatly exceed the conditions on which the bias is based. The results also improved our confidence in the sea-state bias correction used for calculating the geophysical data records. Any error in the correction must influence Ku and C-band ranges almost equally.

Stewart, Robert H.↗
NASA NTRS2010

Evaluating the Impact of Vegetation Cover and Atmospheric Characteristics on the Estimation of Snow Water Equivalent from Spaceborne Microwave Radiometry

A radiative transfer model for estimating snow water equivalent (SWE, mm) from satellite-observed brightness temperature (K) at 19 and 37 GHz (respectively, T(sub B(sub, sat,19)) and T(sub B(sub, sat,37)) over partially forested area is presented, as an extension of a previously published model, by considering scattering of radiation within the canopy. For the specific case of dense vegetation covering fractional area f, the model can be written as, SWE = alpha{ A. delta (T(sub B(sub, sat)) + B - C. f}/(l f), where delta T(sub B(sub, sat)), is the difference of T(sub B(sub, sat,19)) and T(sub B(sub, sat,37)), alpha(mm/K) is the slope of SWE vs. brightness temperature difference at 19 and 37 GHz that would be obtained by ignoring the presence of atmosphere, delta(T(sub B)sub g)), for a homogeneous snow cover (which varies with grain size). The parameters A, B, and C, are determined primarily by atmospheric characteristics, and for a likely range of atmospheric conditions appear to be in the range of, respectively, 1.15-1.63, 0.69-2.84 K and 0.59-2.39 K. Ignoring atmospheric correction would introduce bias towards underestimation of SWE (and also, snow cover area and snow depth). Increasing cloud liquid water path (L) has the effect of increasing A, and ignoring this variation of A with L would have the impact of biasing the estimate of SWE (and snow extent). Such biasing is further exacerbated with increasing f, because of the appearance of term (l-f) in the denominator. The impact of ignoring the intercept parameters (B and C) would be noticeable at low values of SWE (appearing as a bias towards underestimation of SWE), which has been determined to be about 6 mm for average environmental conditions. The uncertainty in estimating SWE due to variations in the atmospheric characteristics is likely to be less than 15%, but could be up to 25% for non-vegetated snow-covered areas. Better estimates of SWE (and snow extent) would be obtained by adjusting the parameters of the above model to regional differences in the atmospheric characteristics. The biases in determining SWE arising due to variations in atmospheric conditions and due to changes in fractional forest cover are not independent, since they interact as {A/(l-f)}. The present calculations also show that improvement in determining snow cover area from the microwave data is likely to occur when these data are corrected for atmospheric effects, as demonstrated by a specific case study.

Choudhury, Bhaskar J.↗
NASA NTRS1977

Compatibility check of measured aircraft responses using kinematic equations and extended Kalman filter

An extended Kalman filter smoother and a fixed point smoother were used for estimation of the state variables in the six degree of freedom kinematic equations relating measured aircraft responses and for estimation of unknown constant bias and scale factor errors in measured data. The computing algorithm includes an analysis of residuals which can improve the filter performance and provide estimates of measurement noise characteristics for some aircraft output variables. The technique developed was demonstrated using simulated and real flight test data. Improved accuracy of measured data was obtained when the data were corrected for estimated bias errors.

Klein, V.↗
NASA NTRS2002

Implementation of the NCAR Community Land Model (CLM) in the NASA/NCAR finite-volume Global Climate Model (fvGCM)

In this study, the NCAR CLM version 2.0 land-surface model was integrated into the NASA/NCAR fvGCM. The CLM was developed collaboratively by an open interagency/university group of scientists and based on well-proven physical parameterizations and numerical schemes that combine the best features of BATS, NCAR-LSM, and IAP94. The CLM design is a one-dimensional point model with 1 vegetation layer, along with sub-grid scale tiles. The features of the CLM include 10-uneven soil layers with water, ice, and temperature states in each soil layer, and five snow layers, with water flow, refreezing, compaction, and aging allowed. In addition, the CLM utilizes two-stream canopy radiative transfer, the Bonan lake model and topographic enhanced streamflow based on TOPMODEL. The DAO fvGCM uses a genuinely conservative Flux-Form Semi-Lagrangian transport algorithm along with terrain- following Lagrangian control-volume vertical coordinates. The physical parameterizations are based on the NCAR Community Atmosphere Model (CAM-2). For our purposes, the fvGCM was run at 2 deg x 2.5 deg horizontal resolution with 55 vertical levels. The 10-year climate from the fvGCM with CLM2 was intercompared with the climate from fvGCM with LSM, ECMWF and NCEP. We concluded that the incorporation of CLM2 did not significantly impact the fvGCM climate from that of LSM. The most striking difference was the warm bias in the CLM2 surface skin temperature over desert regions. We determined that the warm bias can be partially attributed to the value of the drag coefficient for the soil under the canopy, which was too small resulting in a decoupling between the ground surface and the canopy. We also discovered that the canopy interception was high compared to observations in the Amazon region. A number of experiments were then performed focused on implementing model improvements. In order to correct the warm bias, the drag coefficient for the soil under the canopy was considered a function of LAI (Leaf Area Index). Analysis of the results revealed that there was a substantial impact, and the warm and dry bias in the CLM2 was significantly reduced. For the interception scheme, the canopy throughfall was increased to allow for more infiltration of precipitation into the soil, resulting in increased low-level moisture and a decrease in the interception loss ratio (canopy evaporation to precipitation).

Radakovich, Jon D.↗
DOE OSTI2023

The two-point correlation function covariance with fewer mocks

We present FITCOV an approach for accurate estimation of the covariance of two-point correlation functions that requires fewer mocks than the standard mock-based covariance. This can be achieved by dividing a set of mocks into jackknife regions and fitting the correction term first introduced in Mohammad & Percival (2022), such that the mean of the jackknife covariances corresponds to the one from the mocks. This extends the model beyond the shot-noise limited regime, allowing it to be used for denser samples of galaxies. We test the performance of our fitted jackknife approach, both in terms of accuracy and precision, using lognormal mocks with varying densities and approximate EZmocks mimicking the Dark Energy Spectroscopic Instrument LRG and ELG samples in the redshift range of z = [0.8, 1.1]. We find that the Mohammad–Percival correction produces a bias in the two-point correlation function covariance matrix that grows with number density and that our fitted jackknife approach does not. We also study the effect of the covariance on the uncertainty of cosmological parameters by performing a full-shape analysis. We demonstrate that our fitted jackknife approach based on 25 mocks can recover unbiased and as precise cosmological parameters as the ones obtained from a covariance matrix based on 1000 or 1500 mocks, while the Mohammad–Percival correction produces uncertainties that are twice as large.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2014

Landsat-8 Thermal Infrared Sensor (TIRS) Vicarious Radiometric Calibration

Launched in February 2013, the Landsat-8 carries on-board the Thermal Infrared Sensor (TIRS), a two-band thermal pushbroom imager, to maintain the thermal imaging capability of the Landsat program. The TIRS bands are centered at roughly 10.9 and 12 micrometers (Bands 10 and 11 respectively). They have 100 m spatial resolution and image coincidently with the Operational Land Imager (OLI), also on-board Landsat-8. The TIRS instrument has an internal calibration system consisting of a variable temperature blackbody and a special viewport with which it can see deep space; a two point calibration can be performed twice an orbit. Immediately after launch, a rigorous vicarious calibration program was started to validate the absolute calibration of the system. The two vicarious calibration teams, NASA/Jet Propulsion Laboratory (JPL) and the Rochester Institute of Technology (RIT), both make use of buoys deployed on large water bodies as the primary monitoring technique. RIT took advantage of cross-calibration opportunity soon after launch when Landsat-8 and Landsat-7 were imaging the same targets within a few minutes of each other to perform a validation of the absolute calibration. Terra MODIS is also being used for regular monitoring of the TIRS absolute calibration. The buoy initial results showed a large error in both bands, 0.29 and 0.51 W/sq m·sr·micrometers or -2.1 K and -4.4 K at 300 K in Band 10 and 11 respectively, where TIRS data was too hot. A calibration update was recommended for both bands to correct for a bias error and was implemented on 3 February 2014 in the USGS/EROS processing system, but the residual variability is still larger than desired for both bands (0.12 and 0.2 W/sq m·sr·micrometers or 0.87 and 1.67 K at 300 K). Additional work has uncovered the source of the calibration error: out-of-field stray light. While analysis continues to characterize the stray light contribution, the vicarious calibration work proceeds. The additional data have not changed the statistical assessment but indicate that the correction (particularly in band 11) is probably only valid for a subset of data. While the stray light effect is small enough in Band 10 to make the data useful across a wide array of applications, the effect in Band 11 is larger and the vicarious results suggest that Band 11 data should not be used where absolute calibration is required.

TIIRS↗
DOE OSTI2024

Radio measurements of the depth of air-shower maximum at the Pierre Auger Observatory

The (AERA), part of the Pierre Auger Observatory, is currently the largest array of radio antenna stations deployed for the detection of cosmic rays, spanning an area of 17 km 2 with 153 radio stations. It detects the radio emission of extensive air showers produced by cosmic rays in the 30–80 MHz band. Here, we report the AERA measurements of the ( X max ), a probe for mass composition, at cosmic-ray energies between 10 17.5 and 10 18.8 eV , which show agreement with earlier measurements with the fluorescence technique at the Pierre Auger Observatory. We show advancements in the method for radio X max reconstruction by comparison to dedicated sets of / air-shower simulations, including steps of reconstruction-bias identification and correction, which is of particular importance for irregular or sparse radio arrays. Using the largest set of radio air-shower measurements to date, we show the radio X max resolution as a function of energy, reaching a resolution better than 15 g cm − 2 at the highest energies, demonstrating that radio X max measurements are competitive with the established high-precision fluorescence technique. In addition, we developed a procedure for performing an extensive data-driven study of systematic uncertainties, including the effects of acceptance bias, reconstruction bias, and the investigation of possible residual biases. These results have been cross-checked with air showers measured independently with both the radio and fluorescence techniques, a setup unique to the Pierre Auger Observatory. Published by the American Physical Society 2024

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2024

Large-Volume Injection and Assessment of Reference Standards for n -Alkane δD and δ 13 C Analysis via Gas Chromatography Isotope Ratio Mass Spectrometry

Compound-specific stable isotope analysis of hydrogen (δD) and carbon (δ 13 C) in organic compounds is a valuable tool in biogeochemical research. A key limitation of this method is the relatively large amount of sample required to achieve desirable precision. We developed a large-volume (20 μL) injection method that allows for high throughput analysis of less concentrated samples and tested it for δ 13 C and δD measurements of n-alkanes. We also conducted a comparison of reference standards and assessed several methods to normalize and correct n-alkane δD and δ13C measurements. The mean precision of the δD method based on 233 environmental n-alkane samples (two to three replications per sample) is 4.0‰ (1σ, estimated from the weighted mean of the pooled unbiased standard deviations) and 0.46‰ (1σ) for δ 13 C from 37 environmental samples (two to three replications per sample). The evaluation of reference standards shows that the use of n-alkane standards with large offsets in δD values in adjacent n-alkane chains can lead to biases in measurement correction. The large-volume injection method shows good reproducibility of δ 13 C and δD measurements of n-alkanes and reduces the required sample concentration by about 80%. We propose that for δD measurements, a reference standard set should be used in which each reference standard has a limited range of δD values and no adjacent n-alkane chains, to minimize memory effects.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2019

Joint Assimilation of the Aura Microwave Limb Sounder and Ozone Mapping and Profiler Suite Limb Profiler Data: Towards a Reanalysis of Stratospheric Ozone for Trend Studies

The future trajectory of the stratospheric ozone recovery will be sensitive to greenhouse gas concentrations through thermal control of chemical loss and via stratospheric circulation changes. The latter in particular is subject to considerable uncertainty meriting continuing monitoring of the evolution of ozone throughout the depth of the stratosphere. Atmospheric reanalyses utilize the data assimilation methodology to obtain comprehensive representations of the state of the atmosphere, including its composition, on multidecadal scales by combining diverse measurements from satellite-borne and conventional data sources. Systematic biases among these various data types pose a challenge for assimilation by introducing spurious discontinuities that affect the utility of reanalyses for studies of long-term variability and trends.In this presentation we will outline an approach, developed at NASA's Global Modeling and Assimilation Office (GMAO), that allows joint assimilation of stratospheric ozone profiles from the Microwave Limb Sounder (MLS) on EOS Aura and the Ozone Mapping and Profiler Suite Limb Profiler (OMPS-LP) currently flying on the Suomi-NPP satellite with future missions projected into the 2030s. We will demonstrate that a simple offline correction significantly reduces biases between MLS and OMPS-LP ozone data providing a strategy for generating a long-term vertically resolved homogenized representation of stratospheric ozone in future reanalyses. One novel element of our approach compared to previous GMAO reanalysis is the use of a version of the Goddard Earth Observing System model with full stratospheric chemistry. We will show selected comparisons of MLS and OMPS-LP assimilation experiments with independent ozonesonde and satellite data as well as two examples of process-based evaluation focused on the 2016 QBO disruption and Arctic winter ozone loss focusing on the relative performance of the MLS and OMPS-LP analyses.

Wargan, K.↗
NASA NTRS2011

Video Meteor Fluxes

The flux of meteoroids, or number of meteoroids per unit area per unit time, is critical for calibrating models of meteoroid stream formation and for estimating the hazard to spacecraft from shower and sporadic meteors. Although observations of meteors in the millimetre to centimetre size range are common, flux measurements (particularly for sporadic meteors, which make up the majority of meteoroid flux) are less so. It is necessary to know the collecting area and collection time for a given set of observations, and to correct for observing biases and the sensitivity of the system. Previous measurements of sporadic fluxes are summarized in Figure 1; the values are given as a total number of meteoroids striking the earth in one year to a given limiting mass. The Gr n et al. (1985) flux model is included in the figure for reference. Fluxes for sporadic meteoroids impacting the Earth have been calculated for objects in the centimeter size range using Super-Schmidt observations (Hawkins & Upton, 1958); this study used about 300 meteors, and used only the physical area of overlap of the cameras at 90 km to calculate the flux, corrected for angular speed of meteors, since a large angular speed reduces the maximum brightness of the meteor on the film, and radiant elevation, which takes into account the geometric reduction in flux when the meteors are not perpendicular to the horizontal. They bring up corrections for both partial trails (which tends to increase the collecting area) and incomplete overlap at heights other than 90 km (which tends to decrease it) as effects that will affect the flux, but estimated that the two effects cancelled one another. Halliday et al. (1984) calculated the flux of meteorite-dropping fireballs with fragment masses greater than 50 g, over the physical area of sky accessible to the MORP fireball cameras, counting only observations in clear weather. In the micron size range, LDEF measurements of small craters on spacecraft have been used to estimate the flux (Love & Brownlee, 1993); here the physical area of the detector is well known, but the masses depend strongly on the unknown velocity distribution. In the same size range, Thomas & Netherway (1989) used the narrow-beam radar at Jindalee to calculate the flux of sporadics. In between these very large and very small sizes, a number of video and photographic observations were reduced by Ceplecha (2001). These fluxes were calculated (details are given in Ceplecha, 1988) taking the Halliday et al. (1984) MORP fireball fluxes, slightly corrected in mass, as a calibration, and adjusting the flux of small cameras to overlap with the number/mass relation from that work.

Campbell-Brown, M. D.↗
DOE OSTI2021

Bayesian sparse learning with preconditioned stochastic gradient MCMC and its applications

Deep neural networks have been successfully employed in an extensive variety of research areas, including solving partial differential equations. Despite its significant success, there are some challenges in effectively training DNN, such as avoiding overfitting in over-parameterized DNNs and accelerating the optimization in DNNs with pathological curvature. Here, we propose a Bayesian type sparse deep learning algorithm. The algorithm utilizes a set of spike-and-slab priors for the parameters in the deep neural network. The hierarchical Bayesian mixture will be trained using an adaptive empirical method. That is, one will alternatively sample from the posterior using preconditioned stochastic gradient Langevin Dynamics (PSGLD), and optimize the latent variables via stochastic approximation. The sparsity of the network is achieved while optimizing the hyperparameters with adaptive searching and penalizing. A popular SG-MCMC approach is Stochastic gradient Langevin dynamics (SGLD). However, considering the complex geometry in the model parameter space in nonconvex learning, updating parameters using a universal step size in each component as in SGLD may cause slow mixing. To address this issue, we apply a computationally manageable preconditioner in the updating rule, which provides a step-size parameter to adapt to local geometric properties. Moreover, by smoothly optimizing the hyperparameter in the preconditioning matrix, our proposed algorithm ensures a decreasing bias, which is introduced by ignoring the correction term in the preconditioned SGLD. According to the existing theoretical framework, we show that the proposed algorithm can asymptotically converge to the correct distribution with a controllable bias under mild conditions. Numerical tests are performed on both synthetic regression problems and learning solutions of elliptic PDE, which demonstrate the accuracy and efficiency of the present work.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS2014

An AeroCom Assessment of Black Carbon in Arctic Snow and Sea Ice

Though many global aerosols models prognose surface deposition, only a few models have been used to directly simulate the radiative effect from black carbon (BC) deposition to snow and sea ice. Here, we apply aerosol deposition fields from 25 models contributing to two phases of the Aerosol Comparisons between Observations and Models (AeroCom) project to simulate and evaluate within-snow BC concentrations and radiative effect in the Arctic. We accomplish this by driving the offline land and sea ice components of the Community Earth System Model with different deposition fields and meteorological conditions from 2004 to 2009, during which an extensive field campaign of BC measurements in Arctic snow occurred. We find that models generally underestimate BC concentrations in snow in northern Russia and Norway, while overestimating BC amounts elsewhere in the Arctic. Although simulated BC distributions in snow are poorly correlated with measurements, mean values are reasonable. The multi-model mean (range) bias in BC concentrations, sampled over the same grid cells, snow depths, and months of measurements, are −4.4 (−13.2 to +10.7) ng/g for an earlier phase of AeroCom models (phase I), and +4.1 (−13.0 to +21.4) ng/g for a more recent phase of AeroCom models (phase II), compared to the observational mean of 19.2 ng/g. Factors determining model BC concentrations in Arctic snow include Arctic BC emissions, transport of extra-Arctic aerosols, precipitation, deposition efficiency of aerosols within the Arctic, and meltwater removal of particles in snow. Sensitivity studies show that the model-measurement evaluation is only weakly affected by meltwater scavenging efficiency because most measurements were conducted in non-melting snow. The Arctic (60-90degN) atmospheric residence time for BC in phase II models ranges from 3.7 to 23.2 days, implying large inter-model variation in local BC deposition efficiency. Combined with the fact that most Arctic BC deposition originates from extra-Arctic emissions, these results suggest that aerosol removal processes are a leading source of variation in model performance. The multi-model mean (full range) of Arctic radiative effect from BC in snow is 0.15 (0.07-0.25) W/sq m and 0.18 (0.06-0.28) W/sq m in phase I and phase II models, respectively. After correcting for model biases relative to observed BC concentrations in different regions of the Arctic, we obtain a multi-model mean Arctic radiative effect of 0.17 W/sq m for the combined AeroCom ensembles. Finally, there is a high correlation between modeled BC concentrations sampled over the observational sites and the Arctic as a whole, indicating that the field campaign provided a reasonable sample of the Arctic.

carbon↗
DOE OSTI2022

NEONiso R Package, versions greater than 0.5.2

Carbon fluxes in terrestrial ecosystems and their response to environmental change are a major source of uncertainty in the modern carbon cycle. The National Ecological Observatory Network (NEON) presents the opportunity to merge eddy covariance (EC)-derived fluxes with CO2 isotope ratio measurements to gain insights into carbon cycle processes. Collected continuously and consistently across >40 sites, NEON EC and isotope data facilitate novel integrative analyses. However, currently provisioned atmospheric isotope data are uncalibrated, greatly limiting ability to perform cross-site analyses. Here, we present two approaches to calibrating NEON CO2 isotope ratios, along with an R package to calibrate NEON data. We find that calibrating CO2 isotopologues independently yields a lower d13C bias (<0.05‰) and higher precision (<0.40‰) than directly correcting d13C with linear regression (bias: <0.11‰, precision: 0.42‰), but with slightly higher error and lower precision in calibrated CO2 mole fraction. The magnitude of the corrections to d13C and CO2 mole fractions vary substantially by site, underscoring the need for users to apply a consistent calibration framework to data in the NEON archive. Post-calibration data sets show that site mean annual d13C correlates negatively with precipitation, temperature, and aridity, but positively with elevation. Forested and agricultural ecosystems exhibit larger gradients in CO2 and d13C than other sites, particularly during the summer and at night. The overview and analysis tools developed here will facilitate cross-site analysis using NEON data, provide a model for other continental-scale observational networks, and enable new advances leveraging the isotope ratios of specific carbon fluxes.

Fiorella, Richard↗