Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “network reduction”

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 1,477 records · Page 82

Digital Twin Framework for PIP-II Linac: AI-Driven Multi-Scale Modeling from Ion Source to 800 MeV

The PIP-II linac will enable >1.2 MW beam power for DUNE, requiring unprecedented operational reliability across its warm front-end (RFQ, MEBT) and five distinct SRF sections operating at 162.5/325/650 MHz. We present a comprehensive digital twin framework uniquely combining a fully differentiable fast beam transport code with neural network surrogates trained on high-fidelity PIC simulations, capturing space charge and nonlinear dynamics beyond traditional envelope codes while achieving 10⁴× speedup at <1% accuracy. End-to-end differentiability enables gradient-based optimization across 500+ parameters simultaneously—previously impossible with conventional tools—while the model incorporates static/dynamic errors and serves as a virtual commissioning platform for diverse hardware integration. The framework facilitates reinforcement learning for pulsed/CW mode transitions, predictive maintenance through anomaly detection, and autonomous tuning algorithm development with real-time execution capability. Validation against physics simulations shows excellent agreement for the front-end, with initial results demonstrating potential for 30% commissioning time reduction and proactive fault mitigation, providing a scalable blueprint for operating next-generation high-intensity accelerators.

Pathak, Abhishek [Fermilab] (ORCID:000000021704208↗

Forward sensitivity approach for estimating eddy viscosity closures in nonlinear model reduction

In this paper, we propose a variational approach to estimate eddy viscosity using forward sensitivity method (FSM) for closure modeling in nonlinear reduced order models. FSM is a data assimilation technique that blends model's predictions with noisy observations to correct initial state and/or model parameters. We apply this approach on a projection based reduced order model (ROM) of the one-dimensional viscous Burgers equation with a square wave defining a moving shock, and the two-dimensional vorticity transport equation formulating a decay of Kraichnan turbulence. We investigate the capability of the approach to approximate an optimal value for eddy viscosity with different measurement configurations. Specifically, we show that our approach can sufficiently assimilate information either through full field or sparse noisy measurements to estimate eddy viscosity closure to cure standard Galerkin reduced order model (GROM) predictions. Therefore, our approach provides a modular framework to correct forecasting error from a sparse observational network on a latent space. We highlight that the proposed GROM-FSM framework is promising for emerging digital twin applications, where real-time sensor measurements can be used to update and optimize surrogate model's parameters.

42 ENGINEERING↗

Training a Neural Network on Analog TaO x ReRAM Devices Irradiated With Heavy Ions: Effects on Classification Accuracy Demonstrated With CrossSim

The image classification accuracy of a TaO x ReRAM-based neuromorphic computing accelerator is evaluated after intentionally inducing a displacement damage up to a fluence of 10 14 2.5-MeV Si ions/cm 2 on the analog devices that are used to store weights. Results are consistent with a radiation-induced oxygen vacancy production mechanism. When the device is in the high-resistance state during heavy ion radiation, the device resistance, linearity, and accuracy after training are only affected by high fluence levels. Here, the findings in this paper are in accordance with the results of previous studies on TaO x -based digital resistive random access memory. When the device is in the low-resistance state during irradiation, no resistance change was detected, but devices with a 4-kΩ inline resistor did show a reduction in accuracy after training at 10 14 2.5-MeV Si ions/cm 2 . This indicates that changes in resistance can only be somewhat correlated with changes to devices’ analog properties. This paper demonstrates that TaO x devices are radiation tolerant not only for high radiation environment digital memory applications but also when operated in an analog mode suitable for neuromorphic computation and training on new data sets.

42 ENGINEERING↗

Valence State Driven Site Preference in the Quaternary Compound Ca5MgAgGe5: An Electron-Deficient Phase with Optimized Bonding

The quaternary phase Ca5Mg0.95Ag1.05(1)Ge5 (3) was synthesized by high-temperature solid-state techniques, and its crystal structure was determined by single-crystal diffraction methods in the orthorhombic space group Pnma – Wyckoff sequence c12 with a = 23.1481(4) Å, b = 4.4736(1) Å, c = 11.0128(2) Å, V = 1140.43(4) Å3, Z = 4. The crystal structure can be described as linear intergrowths of slabs cut from the CaGe (CrB-type) and the CaMGe (TiNiSi-type; M = Mg, Ag) structures. Hence, 3 is a hettotype of the hitherto missing n = 3 member of the structure series with the general formula R2+nT2X2+n, previously described with n = 1, 2, and 4. The member with n = 3 was predicted in the space group Cmcm – Wyckoff sequence f5c2. The experimental space group Pnma (in the nonstandard setting Pmcn) corresponds to a klassengleiche symmetry reduction of index two of the predicted space group Cmcm. This transition originates from the switching of one Ge and one Ag position in the TiNiSi-related slab, a process that triggers an uncoupling of each of the five 8f sites in Cmcm into two 4c sites in Pnma. The Mg/Ag site preference was investigated using VASP calculations and revealed a remarkable example of an intermetallic compound for which the electrostatic valency principle is a critical structure-directing force. The compound is deficient by one valence electron according to the Zintl concept, but LMTO electronic structure calculations indicate electronic stabilization and overall bonding optimization in the polyanionic network. Other stability factors beyond the Zintl concept that may account for the electronic stabilization are discussed.

36 MATERIALS SCIENCE↗

Surrogate models for plasma displacement and current in 3D perturbed magnetohydrodynamic equilibria in tokamaks

Abstract A numerical database of over one thousand perturbed three-dimensional (3D) equilibria has been generated, constructed based on the MARS-F (Liu et al 2000 Phys. Plasmas 7 3681) computed plasma response to the externally applied 3D field sources in multiple tokamak devices. Perturbed 3D equilibria with the n = 1–4 ( n is the toroidal mode number) toroidal periodicity are computed. Surrogate models are created for the computed perturbed 3D equilibrium utilizing model order reduction (MOR) techniques. In particular, retaining the first few eigenstates from the singular value decomposition (SVD) of the data is found to produce reasonably accurate MOR-representations for the key perturbed quantities, such as the perturbed parallel plasma current density and the plasma radial displacement. SVD also helps to reveal the core versus edge plasma response to the applied 3D field. For the database covering the conventional aspect ratio devices, about 95% of data can be represented by the truncated SVD-series with inclusion of only the first five eigenstates, achieving a relative error (RE) below 20%. The MOR-data is further utilized to train neural networks (NNs) to enable fast reconstruction of perturbed 3D equilibria, based on the two-dimensional equilibrium input and the 3D source field. The best NN-training is achieved for the MOR-data obtained with a global SVD approach, where the full set of samples used for NN training and testing are stretched and form a large matrix which is then subject to SVD. The fully connected multi-layer perceptron, with one or two hidden layers, can be trained to predict the MOR-data with less than 10% RE. As a key insight, a better strategy is to train separate NNs for the plasma response fields with different toroidal mode numbers. It is also better to apply MOR and to subsequently train NNs separately for conventional and low aspect ratio devices, due to enhanced toroidal coupling of Fourier spectra in the plasma response in the latter case.

3D equilibrium↗

Engineering Inertial and Primary-Frequency Response for Distributed Energy Resources

We propose a framework to engineer synthetic-inertia and droop-control parameters for distributed energy resources (DERs) so that the system frequency in a network composed of DERs and synchronous generators conforms to prescribed transient and steady-state performance specifications. Our approach is grounded in a second-order lumped-parameter model that captures the dynamics of synchronous generators and frequency-responsive DERs endowed with inertial and droop control. A key feature of this reduced-order model is that its parameters can be related to those of the originating higherorder dynamical model. This allows one to systematically design the DER inertial and droop-control coefficients leveraging classical frequency-domain response characteristics of second-order systems. Time-domain simulations validate the accuracy of the model-reduction method and demonstrate how DER controllers can be designed to meet steady-state-regulation and transient-performance specifications.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Quantifying the value of redundant measurements at GCOS Reference Upper-Air Network sites

The potential for measurement redundancy to reduce uncertainty in atmospheric variables has not been investigated comprehensively for climate observations. We evaluated the usefulness of entropy and mutual correlation concepts, as defined in information theory, for quantifying random uncertainty and redundancy in time series of the integrated water vapour (IWV) and water vapour mixing ratio profiles provided by five highly instrumented GRUAN (GCOS, Global Climate Observing System, Reference Upper-Air Network) stations in 2010–2012. Results show that the random uncertainties on the IWV measured with radiosondes, global positioning system, microwave and infrared radiometers, and Raman lidar measurements differed by less than 8%. Comparisons of time series of IWV content from ground-based remote sensing instruments with in situ soundings showed that microwave radiometers have the highest redundancy with the IWV time series measured by radiosondes and therefore the highest potential to reduce the random uncertainty of the radiosondes time series. Moreover, the random uncertainty of a time series from one instrument can be reduced by ~ 60% by constraining the measurements with those from another instrument. The best reduction of random uncertainty is achieved by conditioning Raman lidar measurements with microwave radiometer measurements. In conclusion, specific instruments are recommended for atmospheric water vapour measurements at GRUAN sites. This approach can be applied to the study of redundant measurements for other climate variables.

54 ENVIRONMENTAL SCIENCES↗

Optimizing Traffic Signal Control to Enhance Transportation Efficiency and Maximize Pedestrian Benefits in the Road Network

Increasing urban mobility requirements demand efficient transportation system strategies for both vehicular and pedestrian movement. This study enhances the Decentralized Graph-based Multi-Agent Reinforcement Learning (DGMARL) approach, originally tailored for vehicular traffic signal timing, to incorporate pedestrian traffic dynamics. The improved algorithm considers crucial metrics such as Eco_PI, assesses vehicle fuel consumption by factoring in stops and delays, and addresses pedestrian waiting time, crucial for system efficiency while acknowledging driver waiting time impact. Utilizing Digital Twin simulation along the MLK Smart Corridor in Chattanooga, Tennessee, the algorithm's performance is compared for various pedestrian control scenarios. To evaluate the effectiveness of DGMARL, this study compared DGMARL-enabled signal management with automated pedestrian traffic detection and an actuated signal management system (real-word baseline) with pedestrian recall, which predetermingly enforces a pedestrian phase every cycle. Findings indicate substantial improvements with DGMARL, showing a 28.29% enhancement in vehicle Eco_PI, a 60.55 % reduction in pedestrian waiting time, and a 55.74% decrease in driver stop delay, on average, compared to the baseline actuated signal timing plan.

Kumarasamy, Vijayalakshmi K [The University of Ten↗

Structure-guided approach to modulate small molecule binding to a promiscuous ligand-activated protein

Ligand-binding promiscuity in detoxification systems protects the body from toxicological harm but is a roadblock to drug development due to the difficulty in optimizing small molecules to both retain target potency and avoid metabolic events. Immense effort is invested in evaluating metabolism of molecules to develop safer, more effective treatments, but engineering specificity into or out of promiscuous proteins and their ligands is a challenging task. To better understand the promiscuous nature of detoxification networks, we have used X-ray crystallography to characterize a structural feature of pregnane X receptor (PXR), a nuclear receptor that is activated by diverse molecules (with different structures and sizes) to up-regulate transcription of drug metabolism genes. We found that large ligands expand PXR’s ligand-binding pocket, and the ligand-induced expansion occurs through a specific unfavorable compound-protein clash that likely contributes to reduced binding affinity. Removing the clash by compound modification resulted in more favorable binding modes with significantly enhanced binding affinity. We then engineered the unfavorable ligand-protein clash into a potent, small PXR ligand, resulting in marked reduction in PXR binding and activation. Structural analysis showed that PXR is remodeled, and the modified ligands reposition in the binding pocket to avoid clashes, but the conformational changes result in less favorable binding modes. Thus, ligand-induced binding pocket expansion increases ligand-binding potential of PXR but is an unfavorable event; therefore, drug candidates can be engineered to expand PXR’s ligand-binding pocket and reduce their safety liability due to PXR binding.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

On the Strong Scaling of the Spectral Element Solver Nek5000 on Petascale Systems. In: EASC '16 Proceedings of the Exascale Applications and Software Conference 2016, Article No. 5

The present work is targeted at performing a strong scaling study of the high-order spectral element fluid dynamics solver Nek5000. Prior studies such as [5] indicated a recommendable metric for strong scalability from a theoretical viewpoint, which we test here extensively on three parallel machines with different performance characteristics and interconnect networks, namely Mira (IBM Blue Gene/Q), Beskow (Cray XC40) and Titan (Cray XK7). The test cases considered for the simulations correspond to a turbulent flow in a straight pipe at four different friction Reynolds numbers Reτ = 180, 360, 550 and 1000. Considering the linear model for parallel communication we quantify the machine characteristics in order to better assess the scaling behaviors of the code. Subsequently sampling and profiling tools are used to measure the computation and communication times over a large range of compute cores. We also study the effect of the two coarse grid solvers XXT and AMG on the computational time. Super-linear scaling due to a reduction in cache misses is observed on each computer. The strong scaling limit is attained for roughly 5000-10,000 degrees of freedom per core on Mira, 30,000 - 50,0000 on Beskow, with only a small impact of the problem size for both machines, and ranges between 10,000 and 220,000 depending on the problem size on Titan. This work aims at being a reference for Nek5000 users and also serves as a basis for potential issues to address as the community heads towards exascale supercomputers.

Offermans, Nicolas↗

Integration of LIBS with Machine Learning for Real-Time Monitoring of Feedstock in H 2 Gasification Applications

This project, funded by the U.S. Department of Energy (DOE) – Office of Fossil Energy under Award Number DE-FE0032177, aimed to assess the feasibility of an integrated Laser-Induced Breakdown Spectroscopy (LIBS) system with advanced machine learning (ML) models for real-time characterization and potential control of hydrogen gasifiers running on waste materials as feedstocks. This was a multidisciplinary effort that encompassed the acquisition and standardized analysis of individual and blended feedstocks—comprising biomass, coal waste, and plastic waste, followed by the development of a dynamic LIBS bench system for material sample analysis and development of predictive ML models. Comprehensive laboratory testing enabled the creation of a robust elemental dataset that served as the foundation for ML model training. Techniques such as Random Forest, Gradient Boosting, Support Vector Regression, and Neural Networks were employed to predict key feedstock properties, including higher heating value (HHV), moisture content, thermal conductivity, and ash composition with high accuracy. The results were validated against experimental data and demonstrated strong potential for real-time application in gasifier control systems. The project concluded with a study on the integration of the LIBS+ML approach for gasifier control and a techno-economic analysis of the implementation of the approach into hydrogen (H 2 ) gasification systems. Dissemination of results was carried out at a DOE meeting. This work establishes a scalable framework for automated, in-line feedstock quality assessment, offering significant implications for process optimization and emissions reduction in hydrogen production.

01 COAL, LIGNITE, AND PEAT↗

Bacterial Community Dynamics in Dichloromethane-Contaminated Groundwater Undergoing Natural Attenuation

The uncontrolled release of the industrial solvent methylene chloride, also known as dichloromethane (DCM), has resulted in widespread groundwater contamination in the United States. Here we investigate the role of groundwater bacterial communities in the natural attenuation of DCM at an undisclosed manufacturing site in New Jersey. Here, we investigate the bacterial community structure of groundwater samples differentially contaminated with DCM to better understand the biodegradation potential of these autochthonous bacterial communities. Bacterial community analysis was completed using high-throughput sequencing of the 16S rRNA gene of groundwater samples (n = 26) with DCM contamination ranging from 0.89 to 9,800,000 μg/L. Significant DCM concentration-driven shifts in overall bacterial community structure were identified between samples, including an increase in the abundance of Firmicutes within the most contaminated samples. And across all samples, a total of 6,134 unique operational taxonomic units (OTUs) were identified, with 16 taxa having strong correlations with increased DCM concentration. Putative DCM degraders such as Pseudomonas, Dehalobacterium and Desulfovibrio were present within groundwater across all levels of DCM contamination. Interestingly, each of these taxa dominated specific DCM contamination ranges respectively. Potential DCM degrading lineages yet to be cited specifically as a DCM degrading organisms, such as the Desulfosporosinus, thrived within the most heavily contaminated groundwater samples. Co-occurrence network analysis revealed aerobic and anaerobic bacterial taxa with DCM-degrading potential were present at the study site. Our 16S rRNA gene survey serves as the first in situ bacterial community assessment of contaminated groundwater harboring DCM concentrations ranging over seven orders of magnitude. Diversity analyses revealed known as well as potentially novel DCM degrading taxa within defined DCM concentration ranges, indicating niche-specific responses of these autochthonous populations. Altogether, our findings suggest that monitored natural attenuation is an appropriate remediation strategy for DCM contamination, and that high-throughput sequencing technologies are a robust method for assessing the potential role of biodegrading bacterial assemblages in the apparent reduction of DCM concentrations in environmental scenarios.

16S rRNA↗

Modelling of long-range transport of Southeast Asia biomass-burning aerosols to Taiwan and their radiative forcings over East Asia

Biomass burning is a major source of aerosols and air pollutants during the springtime in Southeast Asia. At Lulin mountain background station (elevation 2862 m) in Taiwan, the concentrations of carbon monoxide (CO), ozone (O3) and particulate matter particles with diameter less than 10 μm (PM10), were measured around 150-250 ppb, 40-60 ppb, and 10-30μg/m3, respectively at spring time (February-April) during 2006 and 2009, which are about 2~3 times higher than those in other seasons. Observations and simulation results indicate that the higher concentrations during the spring time are clearly related to biomass burning plumes transported from the Indochina Peninsula of Southeast Asia. The spatial distribution of high aerosols optical depth (AOD) were identified by the satellite measurement and Aerosol Robotic Network (AERONET) ground observation, and could be reasonably captured by the WRF-Chem model during the study period of 15-18 March, 2008. AOD reached as high as 0.8-1.0 in Indochina ranging from 10 to 22°N and 95 to 107°E. Organic carbon (OC) is a major contributor of AOD over Indochina according to simulation results. The contributor of AOD from black carbon (BC) is minor when compared with OC over the Indochina. However, the direct absorption radiative forcing of BC in the atmosphere could reach 35-50 W m-2, which is about 8-10 times higher than that of OC. The belt shape of radiation reduction at surface from Indochina to Taiwan could be as high 20-40 W m-2 during the study period. The implication of the radiative forcing from biomass burning aerosols and their impact on the regional climate in East Asia is our major concern.

aerosols optical depth↗

Biomass-burning Aerosols in South East-Asia: Smoke Impact Assessment (BASE-ASIA)

Biomass burning has been a regular practice for land clearing and land conversion in many countries, especially those in Africa, South America, and Southeast Asia. However, the unique climatology of Southeast Asia is very different than that of Africa and South America, such that large-scale biomass burning causes smoke to interact extensively with clouds during the peak-burning season of March to April. Significant global sources of greenhouse gases (e.g., CO2, CH4), chemically active gases (e.g., NO, CO, HC, CH3Br), and atmospheric aerosols are produced by biomass burning processes. These gases influence the Earth-atmosphere system, impacting both global climate and tropospheric chemistry. Some aerosols can serve as cloud condensation nuclei, which play an important role in determining cloud lifetime and precipitation, hence, altering the earth's radiation and water budget. Biomass burning also affects the biogeochemical cycling of nitrogen and carbon compounds from the soil to the atmosphere; the hydrological cycle (i.e., run off and evaporation); land surface reflectivity and emissivity; as well as ecosystem biodiversity and stability. Analyses from satellite measurements reveal that smoke is frequently present solar (emitted thermal) radiation from clouds due to smoke aerosols can be reduced (enhanced) by as much as 100 (20) W/sq m over the month of March 2000. In addition, the reduction in cloud spectral reflectance at 670 run is large enough to lead to significant errors in retrieving cloud properties (e.g., optical thickness and effective radius) from satellite measurements. The fresh water distribution in this region is highly dependent on monsoon rainfall; in fact, the predictability of the tropical climate system is much reduced during the boreal spring. Estimating the burning fuel (e.g., bark, branches, and wood), an important part of studying regional carbon cycle, may rely on utilizing a wide range of distinctive spectral features in the shortwave and longwave regions. Therefore, to accurately assess the impact of smoke aerosols in this region requires continuous observations from satellites, aircraft, networks of ground-based instruments and dedicated field experiments. A new initiative will be proposed and discussed.

Tsay, Si-Chee↗

Crossing the bridge from molecular catalysis to a heterogenous electrode in electrocatalytic water oxidation

Significant progress has been made in designing single-site molecular Ru(II)-polypyridyl-aqua catalysts for homogenous catalytic water oxidation. Surface binding and transfer of the catalytic reactivity onto conductive substrates provides a basis for heterogeneous applications in electrolytic cells and dye-sensitized photoelectrosynthesis cells (DSPECs). Earlier efforts have focused on phosphonic acid (-PO 3 H 2 ) or carboxylic acid (-CO 2 H) bindings on oxide surfaces. However, issues remain with limited surface stabilities, especially in aqueous solutions at higher pH under conditions that favor water oxidation by reducing the thermodynamic barrier and accelerating the catalytic rate using atom-proton transfer (APT) pathways. Here, we address the problem by combining silane surface functionalization and surface reductive electropolymerization on mesoporous, nanofilms of indium tin oxide (ITO) on fluorine-doped tin oxide (FTO) substrates (FTO| nano ITO). FTO| nano ITO electrodes were functionalized with vinyltrimethoxysilane (VTMS) to introduce vinyl groups on the electrode surfaces by silane attachment, followed by surface electropolymerization of the vinyl-derivatized complex, [Ru II (Mebimpy)(dvbpy)(OH 2 )] 2+ (1 2+ ; Mebimpy: 2,6-bis(1-methyl-1 H -benzo[ d ]imidazol-2-yl)pyridine; dvbpy: 5,5'-divinyl-2,2'-bipyridine), in a mechanism dominated by a grafting-through method. The surface coverage of catalyst 1 2+ was controlled by the number of electropolymerization cycles. The combined silane attachment/cross-linked polymer network stabilized 1 2+ on the electrode surface under a variety of conditions especially at pH > ~6. Surface-grafted poly1 2+ was stable toward redox cycling at pH ~ 7.5 over an ~4-h period. Finally, sustained heterogeneous electrocatalytic water oxidation by the electrode gave steady-state currents for at least ~6 h with a Faradaic efficiency of ~68% for O 2 production.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

PySFD: comprehensive molecular insights from significant feature differences detected among many simulated ensembles

Motivation Many modeling analyses of molecular dynamics (MD) simulations are based on a definition of states that can be (groups of) clusters of simulation frames in a feature space composed of molecular coordinates. With increasing dimension of this feature space (due to the increasing size or complexity of a simulated molecule), it becomes very difficult to cluster the underlying MD data and estimate a statistically robust model. To mitigate this “curse of dimensionality”, one can reduce the feature space, e.g., with principal component or time-lagged independent component analysis transformations, focusing the analysis on the most important modes of transitions. In practice, however, all these reduction strategies may neglect important molecular details that are susceptible to experimental verification. Results To recover such molecular details, I have developed PySFD (Significant Feature Differences analyzer for Python), a multi-processing software package that efficiently selects significantly different features of any user-defined feature type among potentially many different simulated state ensembles, such as meta-stable states of a Markov State Model (MSM). Applying PySFD on MSMs of an aggregate of 300 microseconds MD simulations recently performed on the major histocompatibility complex class II (MHCII) protein, I demonstrate how this toolkit can extract and visualize valuable mechanistic information from big MD simulation data, e.g., in form of networks of dynamic interaction changes connecting functionally relevant sites of a protein complex.

Biochemistry & Molecular Biology↗

Charging-management And Infrastructure-planning (cmip) Model

CMIP model explores various charging infrastructure network designs to serve a free-floating car-sharing fleet and determine the charging downtime experienced by the fleet for each design. Development of the CMIP model had two major steps: (1) describing modeling assumptions and (2) developing an integer program (IP) that jointly optimizes decisions about locations to install DC fast chargers and EV-to-charger assignments. The CMIP model integrates an EV charging model, EV energy consumption model, and heterogeneous, real-world vehicle use data with an integer programming optimization model to identify optimal location of new charging stations and calculate vehicle downtime for charging. The CMIP model can be applied to understand: (a) the reduction of EV fleet downtime if an additional fast-charging station is added to the current infrastructure and (b) to what extent total vehicle downtime would be sensitive to additional charging infrastructure.

Roni, MohammadS↗

Audit of the global carbon budget: estimate errors and their impact on uptake uncertainty

Over the last 5 decades monitoring systems have been developed to detect changes in the accumulation of carbon (C) in the atmosphere and ocean; however, our ability to detect changes in the behavior of the global C cycle is still hindered by measurement and estimate errors. Here we present a rigorous and flexible framework for assessing the temporal and spatial components of estimate errors and their impact on uncertainty in net C uptake by the biosphere. We present a novel approach for incorporating temporally correlated random error into the error structure of emission estimates. Based on this approach, we conclude that the 2&sigma; uncertainties of the atmospheric growth rate have decreased from 1.2 Pg C yr ₋1 in the 1960s to 0.3 Pg C yr ₋1 in the 2000s due to an expansion of the atmospheric observation network. The 2&sigma; uncertainties in fossil fuel emissions have increased from 0.3 Pg C yr ₋1 in the 1960s to almost 1.0 Pg C yr ₋1 during the 2000s due to differences in national reporting errors and differences in energy inventories. Lastly, while land use emissions have remained fairly constant, their errors still remain high and thus their global C uptake uncertainty is not trivial. Currently, the absolute errors in fossil fuel emissions rival the total emissions from land use, highlighting the extent to which fossil fuels dominate the global C budget. Because errors in the atmospheric growth rate have decreased faster than errors in total emissions have increased, a ~20% reduction in the overall uncertainty of net C global uptake has occurred. Given all the major sources of error in the global C budget that we could identify, we are 93% confident that terrestrial C uptake has increased and 97% confident that ocean C uptake has increased over the last 5 decades. Thus, it is clear that arguably one of the most vital ecosystem services currently provided by the biosphere is the continued removal of approximately half of atmospheric CO 2 emissions from the atmosphere, although there are certain environmental costs associated with this service, such as the acidification of ocean waters.

54 ENVIRONMENTAL SCIENCES↗