Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “compiler”

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,765 records · Page 98

Developing an Undergraduate Research Laboratory for Experimental Physical Chemistry [Book Chapter]

Embarking on a faculty career in physical chemistry at a primarily undergraduate institution (PUI) is exciting, but also overwhelming. Becoming both an effective teacher and a productive researcher requires a good plan, hard work and support from others who understand what it takes to be successful. The phrase, “building a laboratory,” calls to mind the physical construction of a research instrument, but beginning a career at a PUI involves much more than assembling equipment. Funding for the lab must be secured, students must be recruited and mentored, and the lab must be overseen to keep things running on a daily basis. While many of these tasks would also be part of starting a career at a doctoral level institution, there are many special considerations that must be taken into account when working with undergraduate students because their laboratory experience and chemical knowledge are still within those beautiful but nascent stages. Here, this unique aspect makes building a research program challenging, but also very rewarding. Exposing young minds to the wonders of physical chemistry and building personal relationships through the mentoring process can be immensely gratifying. This chapter is a compilation of holistic advice for building an experimental physical chemistry research laboratory for undergraduates, reflecting a typical experience found in the United States.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Development of a Conceptual Process for Selective CO 2 Capture from Fuel Gas Streams Using [hmim][Tf 2 N] Ionic Liquid as a Physical Solvent

The Ionic Liquid (IL) [hmim][Tf 2 N] was used as a physical solvent in an Aspen Plus simulation, employing the Peng-Robinson Equation of State (P-R EOS) with Boston-Mathias (BM) alpha function and standard mixing rules, to develop a conceptual process for CO 2 capture from a shifted warm fuel gas stream produced from Pittsburgh # 8 coal for a 400 MWe power plant. The physical properties of the IL, including density, viscosity, surface tension, vapor pressure and heat capacity were obtained from literature and modeled as a function of temperature. Also, available experimental solubility values for CO 2 , H 2 , H 2 S, CO, and CH 4 in this IL were compiled and their binary interaction parameters (Δ ij and l ij ) were optimized and correlated as functions of temperature. The Span-Wager Equation-of-State EOS was also employed to generate CO 2 solubilities in [hmim][Tf 2 N] at high pressures (up to 10 MPa) and temperatures (up to 510 K). The conceptual process developed consisted of 4 adiabatic absorbers (2.4 m ID, 30 m high) arranged in parallel and packed with Plastic Pall Rings of 0.025 m for CO 2 capture; 3 flash drums arranged in series for solvent (IL) regeneration with the pressure-swing option; and a pressure-intercooling system for separating and pumping CO 2 up to 153 bar to the sequestration sites. The compositions of all process streams, CO 2 capture efficiency, and net power were calculated using Aspen Plus simulator. The results showed that, based on the composition of the inlet gas stream to the absorbers, 95.67 mol% of CO 2 was captured and sent to sequestration sites; 99.5 mol% of H 2 was separated and sent to turbines; the solvent exhibited a minimum loss of 0.31 mol%; and the net power balance of the entire system was 30.81 MW. These results indicated that [hmim][Tf 2 N] IL could be used as a physical solvent for CO 2 capture from warm shifted fuel gas streams with high efficiency.

20 FOSSIL-FUELED POWER PLANTS↗

Substitution Reactions in the Pyrolysis of Acetone Revealed through a Modeling, Experiment, Theory Paradigm

The development of high-fidelity mechanisms for chemically reactive systems is a challenging process that requires the compilation of rate descriptions for a large and somewhat ill-defined set of reactions. The present unified combination of modeling, experiment, and theory provides a paradigm for improving such mechanism development efforts. Here we combine broadband rotational spectroscopy with detailed chemical modeling based on rate constants obtained from automated ab initio transition state theory-based master equation calculations and high-level thermochemical parametrizations. Broadband rotational spectroscopy offers quantitative and isomer-specific detection by which branching ratios of polar reaction products may be obtained. Using this technique, we observe and characterize products arising from H atom substitution reactions in the flash pyrolysis of acetone (CH 3 C(O)CH 3 ) at a nominal temperature of 1800 K. The major product observed is ketene (CH 2 CO). Minor products identified include acetaldehyde (CH 3 CHO), propyne (CH 3 CCH), propene (CH 2 CHCH 3 ), and water (HDO). Literature mechanisms for the pyrolysis of acetone do not adequately describe the minor products. The inclusion of a variety of substitution reactions, with rate constants and thermochemistry obtained from automated ab initio kinetics predictions and Active Thermochemical Tables analyses, demonstrates an important role for such processes. Furthermore, the pathway to acetaldehyde is shown to be a direct result of substitution of acetone's methyl group by a free H atom, while propene formation arises from OH substitution in the enol form of acetone by a free H atom.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

SimELIT: A Novel GUI-Based Comprehensive Ion Trajectory Simulation Software for Mass Spectrometry

Ion trajectory simulation in mass spectrometry systems from injection to detection is technically challenging but very important for better understanding the ion dynamics in instrument development. Here, in this work, we present SimELIT (Simulator of Eulerian and Lagrangian Ion Trajectories), a novel ion trajectory simulation platform. SimELIT is built upon a suite of multiphysics solvers compiled into OpenFOAM (an open-source numerical solver library particularly used for computational mechanics), with a simple web-based graphical user interface (GUI) allowing users to define the details of OpenFOAM cases and run simulations. SimELIT is a modular program and can provide extensions of physics (e.g., gas flows, electrodynamic fields) and thus enable ion trajectory simulations from the ion source to detector. The current version (SimELIT) provides two numerical solvers for ion trajectory simulations–(1) a Lagrangian particle tracker in vacuum and (2) a Eulerian ion density solver in background gas in the presence of electric fields. Here, we describe the architecture of SimELIT, including its use of Docker and the React Framework, and demonstrate the computation of ion trajectories of multiple m/z values in a static/linear voltage drop in vacuum (across a 1 m long flight tube). Further, the drift motion of ions under 1 Torr pressure conditions in a static background (N 2 ) gas through a 20 V/cm static electric field is shown. The results produced from SimELIT were compared with SIMION and theoretical estimates. In addition, we report the computation of ion trajectories in electrodynamic fields within a planar FAIMS device operating at atmospheric pressure.

97 MATHEMATICS AND COMPUTING↗

Structural Control on Downdip Locking Extent of the Himalayan Megathrust

Here, geologic reconstructions of the Main Himalayan Thrust in Nepal show a laterally extensive midcrustal ramp, hypothesized to form the downdip boundary of interseismic locking. Using a recent compilation of interseismic GPS velocities and a simplified model of fault coupling, we estimate the width of coupling across Nepal using a series of two-dimensional transects. We find that the downdip width of fault coupling increases smoothly from 70 to 90 km in eastern Nepal to 100–110 km in central Nepal, then narrows again in western Nepal. The inferred coupling transition is closely aligned with geologic reconstructions of the base of the midcrustal ramp in central and eastern Nepal, but in western Nepal, the data suggest that the location is intermediate between two proposed ramp locations. The result for western Nepal implies either an anomalous coupling transition that occurs along a shallowly dipping portion of the fault or that both ramps may be partially coupled and that a proposed crustal-scale duplexing process may be active during the interseismic period. We also find that the models require a convergence rate of 15.5 ± 2 mm/year throughout Nepal, reducing the geodetic moment accumulation rate by up to 30% compared with earlier models, partially resolving an inferred discrepancy between geodetic and paleoseismic estimates of moment release across the Himalaya.

58 GEOSCIENCES↗

Reply to Comment by Peterie Et Al. on “Accelerated Fill‐Up of the Arbuckle Group Aquifer and Links to U.S. Midcontinent Seismicity”

Abstract Peterie et al. question one observation in our paper: associating pressure increases to injection volumes at distances of up to 25 km from an injection well. In this reply, we show that the comment misunderstands our analysis and the evidence that led to this conclusion. We also show that gauge‐depth‐corrected pressures, used by the authors to produce statewide pressure maps, are discrepant with the static fluid level data, provided in our original compilation and analysis. The discrepancies are a result of the pressure correction method employed, which naïvely substitutes formation pressure for bottomhole pressure to calculate wellbore fluid density. Their linearly interpolated pressure maps, based on sparse data, contain interpolation and extrapolation artifacts that contradict injection trends in the state, the Theis solution, and the superposition principle. We reiterate that pressure and static fluid level increases in Class I wells existed prior to 2013, most notably in central Kansas, where recent earthquakes are cited in the comment as evidence of a pressure plume emanating from the Kansas‐Oklahoma border, 90 km away. We show that the space‐time pattern of seismicity in this area is inconsistent with a northward propagating pressure plume and, instead, seismicity appears to be centered on and near a cluster of high‐rate injection wells, two of which are among the highest rate wells in the state. These observations, along with recent M4 + earthquakes during continued decreases in wastewater injection in southern Kansas and northern Oklahoma, question the usefulness of the comment for understanding and managing societally significant earthquakes.

Ansari, Esmail↗

A data-driven global soil heterotrophic respiration dataset and the drivers of its inter-annual variability

Soil heterotrophic respiration (SHR), one of the primary carbon fluxes from terrestrial ecosystems to the atmosphere, is important for carbon-climate feedbacks because of its sensitivity to available litter and soil carbon, climatic conditions, and nutrient availability. However, until recently limited SHR data were available, and most published global SHR estimates have either a short time span, coarse spatial resolution, or reply on overly-simple model formulations. To better understand and quantify the global distribution of SHR and its sensitivity to climate variability, we produced a new global SHR dataset using Random Forest algorithms, up-scaling 455 point data from the Global Soil Respiration Database (SRDB 4.0) with gridded fields of climatic, edaphic and productivity as explanatory variables. We estimated a global total SHR of 46.8 Pg C yr-1 over 1985-2013 (95% confidence interval: 38.6-56.3 Pg C yr-1), with a significant increasing trend of 0.03 Pg C yr-2 during this period. We found that the choice of soil moisture datasets contributes more to the difference among these data-driven SHR members rather than that of productivity, temperature and precipitation data sources. We also analyzed the influence of climatic variables on the inter-annual variability (IAV) of our SHR product. Water availability was the dominant driver of IAV at global scales, although the inferred sensitivity depends on the choice of the soil moisture gridded dataset. At the ecosystem scale, temperature strongly controls the IAV of SHR in tropical forests, while water availability dominates in extra-tropical forest and semi-arid regions. Our machine-learning gridded SHR dataset and outputs from process-based land surface models (TRENDYv6) show agreement for a strong association between water variability and SHR IAV at the global scale, but the two approaches lead to different temporal trend globally and different controlling variables for IAV at the ecosystem scale. Our study provides evidence for the pervasive and important role of water availability in driving SHR, indicating both a direct effect limiting decomposition rates and an indirect effect through the amount of fresh organic matter made available to SHR from productivity. In consideration of potential limitations and uncertainties remaining in our data-driven SHR datasets, we call for a more scientifically designed observation network for SHR, more observation data compilation, and increased use of deep learning methods making maximum use of observation data in hand. This will benefit process-based models, and improve our understanding of SHR response to future anomalous environmental conditions.

Yao, Yitong↗

Identifying Key Environmental Factors Explaining Temporal Patterns of DOC Export From Watersheds in the Conterminous United States

The export of dissolved organic carbon (DOC) from a watershed is a key component of the terrestrial biosphere carbon cycle. There is a need to improve our understanding of how and by how much various environmental factors are driving the temporal patterns of DOC export in order to accurately model and evaluate terrestrial carbon storage and fluxes. In this synthesis, we compiled observational data sets from 14 watersheds in the conterminous United States spanning the time period from 1981 to 2017. We used these data sets to examine the relative impacts of various climate, atmospheric deposition, and land cover factors on the temporal patterns of DOC export across watersheds of different sizes and landscape conditions, as well as the time-series autocorrelation of DOC export. Overall, our results suggest that the dominant factor on an annual scale was the amount of precipitation, which had a positive correlation with DOC export from a watershed. Overall increasing nitrogen deposition was coincident with increasing DOC export, and increasing sulfur deposition was coincident with declining DOC export. The seasonal pattern of DOC export was strongly regulated by air temperature, the long-term trend was negatively influenced by increasing sulfur deposition, and no obvious autocorrelation was detected in DOC export. In addition, higher rates of DOC export were positively correlated with greater area of wetlands within a watershed, but was not found to be strongly related to any of the other land cover types.

54 ENVIRONMENTAL SCIENCES↗

Dynamics of Fungal and Bacterial Biomass Carbon in Natural Ecosystems: Site-level Applications of the CLM-Microbe Model

Explicitly representing microbial processes has been recognized as a key improvement to Earth system models for the realistic projections of soil carbon (C) and climate dynamics. The CLM-Microbe model builds upon the CLM4.5 and explicitly represents two major soil microbial groups, fungi and bacteria. Based on the compiled time-series data of fungal (FBC) and bacterial (BBC) biomass C from nine biomes, we parameterized and validated the CLM-Microbe model, and further conducted sensitivity analysis and uncertainty analysis for simulating C cycling. The model performance was evaluated with mean absolute error (MAE), root mean square error (RMSE), and coefficient of determination (R2) for relative change in FBC and BBC. The CLM-Microbe model is able to reasonably capture the seasonal dynamics of FBC and BBC across biomes, particularly for tropical/subtropical forest, temperate broadleaf forest, and grassland, with MAE < 0.49 for FBC and <0.36 for BBC and RMSE <0.52 FBC and <0.39 for BBC, while R2 values are relatively smaller in some biomes (e.g., shrub) due to small sample sizes. We found good consistencies between simulated and observed FBC (R 2 =0.70, P<0.001) and BBC (R 2 =0.26, P<0.05) on average across biomes, but the model is not able to fully capture the large variation in observed FBC and BBC. Sensitivity analysis shows the most critical parameters are turnover rate, carbon-to-nitrogen ratio of fungi and bacteria, and microbial assimilation efficiency. This study confirms that the explicit representation of soil microbial mechanisms enhances model performance in simulating C variables such as heterotrophic respiration and soil organic C density. The further application of the CLM-Microbe model would deepen our understanding of microbial contributions to the global C cycle.

54 ENVIRONMENTAL SCIENCES↗

Subduction, Underplating, and Return Flow Recorded in the Cycladic Blueschist Unit Exposed on Syros, Greece

Exhumed high-pressure/low-temperature (HP/LT) metamorphic rocks provide insights into deep (~20–70 km) subduction interface dynamics. On Syros Island (Cyclades, Greece), the Cycladic Blueschist Unit preserves blueschist-to-eclogite facies oceanic- and continental-affinity rocks that record the structural and thermal evolution linked to Eocene subduction. Despite decades of research, the metamorphic and deformation history (P-T-D) and timing of subduction and exhumation are matters of ongoing discussion. We suggest that Syros comprises three coherent tectonic slices and that each slice underwent subduction, underplating, and syn-subduction return flow along similar P-T trajectories, but at progressively younger times. Subduction and exhumation are distinguished by lineations and ductile fold axis orientations, and are kinematically consistent with previous studies that document top-to-the-S-SW shear (prograde-to-peak subduction), top-to-the-NE shear (blueschist facies exhumation), and then E-W coaxial stretching (greenschist facies exhumation). Amphibole zonations record cooling during decompression, indicating return flow above a cold slab. Multi-mineral Rb-Sr isochrons and compiled metamorphic geochronology show that the three slices record distinct stages of peak subduction (53–52, ~50, and 45 Ma) that young with structural depth. Retrograde blueschist and greenschist facies fabrics span ~50–40 and ~43–20 Ma, respectively, and also young with structural depth. Synthesized data sets support a revised tectonic framework for Syros, involving subduction of structurally distinct coherent slices and simultaneous return flow of previously accreted tectonic slices in the subduction channel shear zone. Distributed, ductile, dominantly coaxial return flow in an Eocene-Oligocene subduction channel proceeded at rates of ~1.5–5 mm/yr and accommodated ~80% of the total exhumation of this HP/LT complex.

58 GEOSCIENCES↗

Precambrian Crystalline Basement Properties From Pressure History Matching and Implications for Induced Seismicity in the US Midcontinent

Injection-induced seismicity across the US midcontinent has almost exclusively occurred in the crystalline basement that underlies the Arbuckle Group aquifer and its equivalents, the primary wastewater disposal zone in this region. However, the properties of the basement are not well known. Newly compiled data, from Class I wells in Kansas, provide a unique record of pressures in the Arbuckle and an opportunity to constrain the reservoir-scale properties of the basement such as permeability, diffusivity, and specific storage. Constraints on these parameters are critical for modeling fluid flow and pressures across the entire Arbuckle-basement system, and are necessary for accurate evaluation and prediction of injection-induced earthquakes. Here, we present a detailed, three-dimensional geological and pressure history-matched numerical model for the Arbuckle and basement, based on data from >400 wells covering a large region in south-central Kansas, where injection-induced seismicity has been concentrated since 2014. Simulations of dynamic data from 319 wells indicate that Arbuckle pressures have increased by 1.1 MPa in high injection rate areas and an overpressure of <0.1 MPa may be the cause of seismicity in the basement. Pressure-history matching also yields the likely range in porosity (0.3%–7%), permeability (0.1–0.7 mD), and diffusivity (0.004–0.07 m2 /s) for the basement. The resulting estimates suggest reservoir-scale properties of the basement are enhanced by faults and fractures. Importantly, the diffusivities determined in this study are lower than estimates derived from Kansas earthquake triggering fronts, and suggest that such seismicity-based techniques may have limitations, particularly where spacetime patterns between injection and seismicity are complex.

58 GEOSCIENCES↗

The Lithophile Element Budget of Earth's Core

Abstract The relative composition of Earth's core and mantle were set during core formation. By determining how elements partition between metal and silicate at high pressures and temperatures, measurements of the mantle composition and geophysical observations of the core can be used to understand the mechanisms by which Earth formed. Here we present the results of metal‐silicate partitioning experiments for a range of nominally lithophile elements (Al, Ca, K, Mg, O, Si, Th, and U) and S to 85 GPa and up to 5400 K. With our results and a compilation of literature data, we developed a parameterization for partitioning that accounts for compositional dependencies in both the metal and silicate phases. Using this parameterization in a range of planetary growth models, we find that, in general, lithophile element partitioning into the metallic phase is enhanced at high temperatures. The relative abundances of FeO, SiO 2 , and MgO in the mantle vary significantly between planetary growth models, and the mantle abundances of these elements can be used to provide important constraints on Earth's accretion. To match Earth's core mass and mantle composition, Earth's building blocks must have been enriched in Fe and depleted in Si compared with CI chondrites. Finally, too little Mg, Si, and O are partitioned into the core for precipitation of oxides to be a major source of energy for the geodynamo. In contrast, several ppb of U can be partitioned into the core at high temperatures, and this energy source must be accounted for in thermal evolution models.

Chidester, B. A.↗

The NPF Effect on CCN Number Concentrations: A Review and Re-Evaluation of Observations From 35 Sites Worldwide

The new particle formation (NPF) effect on cloud condensation nuclei (CCN) differs from region to region and remains highly uncertain. In this work, we re-evaluate the NPF-initiated enhancements in CCN number concentrations (E_N CCN ) by compiling results from 35 field sites worldwide and five global models. We show the E_N CCN at sites close to the anthropogenic sources are more than 2-fold as large as that at remote sites with supersaturations of 0.1%–1.0%. The enhancements anti-correlate with the time scales for the new particles to grow to CCN size, which are 5–10 and 20–30 hr at polluted and remote sites, respectively. We further reveal that under polluted environments the condensational growth of new particles may be accelerated, leading to larger E_N CCN . The models fail to simulate the E_NCCN in polluted regions. Our study highlights the importance of parameterizing the time scale of NPF to affect CCN differently for remote and polluted regions in models.

54 ENVIRONMENTAL SCIENCES↗

Late Paleozoic Evolution of the Anadarko Basin: Implications for Laurentian Tectonics and the Assembly of Pangea

Abstract Interplay between marginal and intraplate late Paleozoic tectonics in southeastern Laurentia is a subject of debate. In southern Laurentia, intraplate deformation has been attributed to combinations of compressive forces from orogenic belts in the west, southwest, and southeast. While left‐lateral transpression is generally an accepted mechanism for formation of uplifts in the Southern Oklahoma Aulacogen (SOA), to date no studies have explained the driving mechanisms of these kinematics. Here, we analyze the Southern Oklahoma Transpressional System (SOTS), a ∼500 km long ∼50 km wide fault zone that extends from the southeastern Laurentian margin into the plate interior along Cambrian lithospheric weaknesses inherited from the SOA. We use a compilation of published studies to present detailed, time‐integrated sedimentary thickness maps of SOTS‐related basins. We combine these maps with a new interpretation of three 2D seismic reflection lines to constrain timing and kinematics of all major late Paleozoic SOTS fault zones. We show that the SOTS kinematics evolved from contraction to left‐lateral strike‐slip during Mississippian‐Pennsylvanian time. Tectonic activity and fault kinematics in the SOTS are explained by a model in which a remnant ocean basin closed diachronously from northeast to southwest along the Southern Appalachian and Ouachita‐Marathon orogens. Left‐lateral transpression during diachronous closure of the Rheic oceanic basin along the Ouachita‐Marathon margin is related to rotation of the stress field due to changes in slab‐pull and eastward drift of Laurentia toward Gondwana. We speculate that the SOTS acted in part as a STEP (subduction‐transfer edge propagator) fault that followed the pre‐existing SOA.

Hobbs, Noah. F.↗

Evaluation of Nominal Energy Storage at Existing Hydropower Reservoirs in the US

Long-term planning and operation of hydropower reservoirs require an understanding of both water and energy storage. As energy storage needs of the evolving grid increase, we must account for the water and energy storage potential of these reservoirs. Given the limitations of current data on existing hydropower, we compile statistics related to storage volume and hydraulic head from publicly available data sets and examine differences in descriptions of US hydropower storage. Assembled characteristics are used to calculate nominal energy storage capacity, a simple measure of potential to generate power from a given volume of water, not factoring in detailed constraints. Inventory-based estimates of energy storage are calculated at 2,075 dams, which helps put the potential for US hydropower to support energy storage in context with similar evaluations in other regions and with other energy storage technologies. The national energy storage capacity ranges between 34.5 and 45.1 TWh depending on the information used, with 52% of energy storage located at the 10 largest reservoirs in the US. Energy storage capacities are also calculated at 236 dams with historical volume and elevation data. Finally, reservoir inflows provide context for the storage volumes and sensitivities to hydrologic variability. Larger reservoirs with greater storage volume to inflow ratios are concentrated in the Western US, but the majority of hydropower reservoirs store less than the annual inflow. We address several infrastructure and water resource informatics challenges and highlight remaining issues, including representing seasonal or shorter variability in water volumes and representing connected hydropower facilities.

13 HYDRO ENERGY↗

Characterizing the Atmospheric Mn Cycle and Its Impact on Terrestrial Biogeochemistry

Abstract The role of manganese (Mn) in ecosystem carbon (C) biogeochemical cycling is gaining increasing attention. While soil Mn is mainly derived from bedrock, atmospheric deposition could be a major source of Mn to surface soils, with implications for soil C cycling. However, quantification of the atmospheric Mn cycle, which comprises emissions from natural (desert dust, sea salts, volcanoes, primary biogenic particles, and wildfires) and anthropogenic sources (e.g., industrialization and land‐use change due to agriculture), transport, and deposition, remains uncertain. Here, we use compiled emission data sets for each identified source to model and quantify the atmospheric Mn cycle by combining an atmospheric model and in situ atmospheric concentration measurements. We estimated global emissions of atmospheric Mn in aerosols (<10 μm in aerodynamic diameter) to be 1,400 Gg Mn year −1 . Approximately 31% of the emissions come from anthropogenic sources. Deposition of the anthropogenic Mn shortened Mn “pseudo” turnover times in 1‐m‐thick surface soils (ranging from 1,000 to over 10,000,000 years) by 1–2 orders of magnitude in industrialized regions. Such anthropogenic Mn inputs boosted the Mn‐to‐N ratio of the atmospheric deposition in non‐desert dominated regions (between 5 × 10 −5 and 0.02) across industrialized areas, but that was still lower than soil Mn‐to‐N ratio by 1–3 orders of magnitude. Correlation analysis revealed a negative relationship between Mn deposition and topsoil C density across temperate and (sub)tropical forests, consisting with atmospheric Mn deposition enhancing carbon respiration as seen in in situ biogeochemical studies.

Environmental Sciences & Ecology↗

Permafrost Region Greenhouse Gas Budgets Suggest a Weak CO 2 Sink and CH 4 and N 2 O Sources, But Magnitudes Differ Between Top-Down and Bottom-Up Methods

Large stocks of soil carbon (C) and nitrogen (N) in northern permafrost soils are vulnerable to remobilization under climate change. However, there are large uncertainties in present-day greenhouse gas (GHG) budgets. We compare bottom-up (data-driven upscaling and process-based models) and top-down (atmospheric inversion models) budgets of carbon dioxide (CO 2 ), methane (CH 4 ) and nitrous oxide (N 2 O) as well as lateral fluxes of C and N across the region over 2000–2020. Bottom-up approaches estimate higher land-to-atmosphere fluxes for all GHGs. Both bottom-up and top-down approaches show a sink of CO 2 in natural ecosystems (bottom-up: -29 (-709, 455), top-down: -587 (-862, -312) Tg CO 2 -C yr -1 ) and sources of CH 4 (bottom-up: 38 (22, 53), top-down: 15 (11, 18) Tg CH 4 -C y -1 ) and N 2 O (bottom-up: 0.7 (0.1, 1.3), top-down: 0.09 (-0.19, 0.37) Tg N 2 O-N yr -1 ). The combined global warming potential of all three gases (GWP-100) cannot be distinguished from neutral. Over shorter timescales (GWP-20), the region is a net GHG source because CH 4 dominates the total forcing. The net CO 2 sink in Boreal forests and wetlands is largely offset by fires and inland water CO 2 emissions as well as CH 4 emissions from wetlands and inland waters, with a smaller contribution from N 2 O emissions. Priorities for future research include the representation of inland waters in process-based models and the compilation of process-model ensembles for CH 4 and N 2 O. Discrepancies between bottom-up and top-down methods call for analyses of how prior flux ensembles impact inversion budgets, more and well-distributed in situ GHG measurements and improved resolution in upscaling techniques.

54 ENVIRONMENTAL SCIENCES↗