Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Geochemical Modelling”

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 163 records · Page 9

Time-dependent lunar density models

Simple models of the geochemical and geophysical evolution of the moon were constructed with the particular aim of investigating the mass density distribution within the moon as a function of time. The strongly inverted (densest material on top) density distribution resulting from fractional crystallization of the lunar magma ocean was found to rearrange itself into a highly stable stratification. Density stratification due to compositional variation dominated gravitational stability; thus, thermally driven convection was found to be suppressed in regions with relaxed (rearranged so as to minimize the gravitational potential energy) compositional density variation. The suppression of convection as a heat transfer process results in considerable modification of thermal evolution. The models produced sufficient density segregation from magma oceans of a few hundred km thickness, such that relaxation to the most stable rearrangement yielded a polar moment of inertia approximating that measured for the present-day moon. Thus, lunar models matching moment-of-inertia constraints but lacking a metallic core are motivated.

Herbert, F.↗

Illuminating the pathways to carbon liberation: a systems approach to characterizing the consequential unknowns of carbon transformation and loss from thawing permafrost peatlands (Final Report)

The IsoGenie3 Project delivered new systems-level insights into carbon cycling in thawing permafrost landscapes, with an emphasis on methane and carbon dioxide emissions. From >200 samples from the site collected over a decade, co-analyzed for geochemistry and microbiology, the team recovered ~1,500 assembled microbial genomes and ~1,900 viral population genomes, revealing appreciable genetic novelty - from a new highly abundant bacterial phylum, to novel methane consumers and their activities, to rampant viral novelty. IsoGenie3 linked these organisms to carbon compound transformations (which define the cycling of organic matter in soils, and the loss of the greenhouse gases carbon dioxide and methane), and saw that the microbes at each stage of permafrost thaw had different genetic potential to degrade categories of compounds, expressed that genetic potential differently, and actually transformed carbon compounds into greenhouse gases in different ways. IsoGenie 3 identified that some of the thaw-stage differences were due to plant-microbiome relationships; the plant species across the thaw gradient contributed different carbon compounds into the soil, and hosted distinct microbiota (differing among parts of plants as well as species). Lastly, microbes in the saturated post-thaw conditions appeared likely to contribute to the mobilization and toxification of mercury released during thaw. In parallel with ongoing field sampling and analysis, hypotheses arising from field observations were tested via lab incubation experiments. When communities are taken out of their native habitats, they behave differently, and the team first rigorously quantified the magnitude of this effect on microbiome composition and functional capacity, organic matter composition, and gas production; overall the main system processes were maintained in the lab incubations under the conditions tested. Further, the microbial data could inform geochemical reaction network models of those processes. Then, the team ran experiments with additions of compounds, varying temperature, and “live” vs. “dead” peat (the latter having been gamma irradiated, with a few additional variants to control for methodological artifacts). From these, we (a) determined the importance of plant-derived soluble phenolic compounds in bogs’ extraordinary recalcitrance of organic matter, and carbon gas emissions skewed to carbon dioxide; (b) proposed an abiotic ‘tanning’ mechanism, which could contribute to Sphagnum’s inhibitory effect on anaerobic decomposition through alteration of N availability. IsoGenie3 illuminated longer-term and landscape-scale interactions of permafrost thaw and carbon cycling, advancing knowledge of the drivers of methane dynamics not only across in the permafrost-associated peatland (where hydrology and plant communities dictate microbiomes) but also their interconnected lakes (where sediment carbon quality and resident microbiota are determined by position within lake, and lake features). By leveraging observations of site methane dynamics extending well before this project, the team was able to construct a 44-year portrait of the interplay of permafrost thaw, hydrology, vegetation dynamics, and carbon gas emissions, and the doubling of the fully-thawed fens over this time. From the detailed study of this focal site, IsoGenie3 also aimed to improve model representation of these kinds of sites and processes. To improve predictions of methane transformations, we incorporated acetate and isotope dynamics into the ‘DNDC’ biogeochemistry model. In addition, recovered genomes were grouped into ‘functional groups’, i.e. the genomes that perform a specific function of interest, then used to parameterize maximum growth rate and optimum growth temperature (via signatures in their sequence composition) for the BioCrunch model. The BioCrunch model was then in turn used to test the impact of increasing functional resolution of the microbes, on the carbon gas emissions. Lastly for modeling, the ecosys model was parameterized from the microbial and other data, and used to evaluate drivers of e.g. change in methane emissions. Finally, this project also led to the development of a range of new methods and tools, a new metric of organic matter decomposability, as well as a graph-database solution to multidisciplinary data storage and querying. This project’s ongoing analyses at our focal site also contributed to broader advancements in understanding elements of genetic plasticity and methane metabolism, climate change microbiology and community assembly, global peatland geochemistry and Arctic lakes’ roles in climate feedbacks.

54 ENVIRONMENTAL SCIENCES↗

UT-GOM2-2 Preliminary Report: Terrebonne Basin Northern Gulf of Mexico, 30 July-28 September 2023

In the summer and fall of 2023, the Gulf of Mexico Deepwater Hydrate Coring Expedition (UT-GOM2-2) drilled, cored, made downhole measurements, and analyzed samples from the seafloor to the base of the gas hydrate stability zone in one location (Site H, WR313) in the Terrebonne basin, deepwater Gulf of Mexico. Analyses of data and samples from the expedition will inform biological, geochemical, and geomechanical models to constrain the role of gas hydrates in the carbon cycle and the potential for gas hydrates as an energy resource. Pressure and conventional cores were collected continuously to a depth of 155.1 meters below the seafloor (mbsf). At deeper depths, cores were taken periodically from hydrate-bearing sands and their bounding muds to a total depth of 861.3 mbsf. 162.6 m of conventional core and 54.8 m of pressure core were obtained. Twelve temperature measurements were made between 27.1 and 144.5 mbsf to determine the geothermal gradient. At the seafloor, more than 4 m of sandy silt of unknown origin was encountered. Beneath this sand, to a depth of ~200 mbsf, the section was composed of interbedded mud and biogenic carbonate ooze. The biogenic ooze correlated to low density and high porosity intervals observed in the previously acquired logging while drilling (LWD) data and as measured. Calcareous nannofossil biostratigraphy constrains the entire record to the Pleistocene (< 0.91 million years) with a pronounced increase in sedimentation rate with depth. Beneath 200 mbsf, the section was predominantly composed of mud with two thicker, hydrate-bearing coarse-grained intervals, which are commonly known as the Blue and Orange sands. The dissolved gas concentration was quantified from pressure cores. In the shallow section, dissolved methane concentration increased below the sulfate-methane transition zone (SMTZ) and reaches saturation (the limit of solubility for methane) at 147 mbsf. Gas expansion was very common in conventional and depressurized pressure (conventionalized) cores below the SMTZ. At deeper depths, the methane concentration within muds bounding the Blue and Orange reservoirs was generally found to be less than saturation. The dissolved and hydrate gas composition is consistent with a microbial source, containing greater than 99.99% methane and only trace concentrations of ethane, propane, and butane. The methane to ethane ratio (C 1 /C 2 ) and the methane to ethane plus propane (C 1 /(C 2 +C 3 )) decrease with depth down to at least 678 mbsf, mainly driven by the increase in ethane with depth. It is unclear if this trend continues through the Orange sand interval. The δ 13 C isotopic signature of methane ranges between -69.9 and -78.5 ‰ Vienna Pee Dee Belemnite (VPDB). Pressure core recovery of all sandy intervals was poor. However, pressure core logs of the Orange sand show intervals of low density and high velocity, which are indicative of high hydrate saturation. One core from within the Orange sand was composed of interbedded graded sandy silt and mud. The sandy silts from this core are composed of mainly quartz and feldspar with some lithics. Most of the recovered pressure core samples are maintained at near in-situ pressure and temperature (within the hydrate stability field) at the University of Texas Pressure Core Center awaiting analysis. In the shallow section, samples will be used to determine the flux of organic carbon through the basin system, find the rate at which that carbon was consumed, and understand the microbial population responsible for these processes. In the deeper section, samples from in and around the hydrate reservoirs will be used to determine the petrophysical properties of the reservoir and bounding seals in these systems.

03 NATURAL GAS↗

Expedition UT-GOM2-2 Summary

In the summer and fall of 2023, the University of Texas (UT) Deepwater Hydrate Coring Expedition (UT-GOM2-2) drilled, cored, made downhole measurements, and analyzed samples from the seafloor to the base of the gas hydrate stability zone at Site H, in the Walker Ridge Protracted Area Block 313 (Site H, WR313), in the Terrebonne Basin, deepwater Gulf of America (Gulf of Mexico). Analyses of data and samples from the expedition will inform biological, geochemical, and geomechanical models to constrain the role of gas hydrates in the carbon cycle and the potential for gas hydrates as an energy resource.

03 NATURAL GAS↗

The atmospheres of the earth and the other planets: Origin, evolution and composition

The current understanding of the composition, chemistry, and structure of the atmospheres of the other planets and the origin, early history, and evolution of the earth's atmosphere is reviewed. The information on the atmospheres of the other planets is based on the successful Mariner, Viking, Pioneer, and Voyager missions to these planets. The information on the origin, early history, and evolution of the atmosphere, which is somewhat speculative, is largely based on numerical studies with geochemical and photochemical models.

Levine, Joel S.↗

Geochemical comparison of K-T boundaries from the Northern and Southern Hemispheres

Closely spaced (cm-scale) traverses through the K-T boundary at Stevns Klint (Denmark), Woodside Creek (New Zealand) and a new Southern Hemisphere site at Richards Bay (South Africa) were subjected to trace element and isotopic (C, O, Sr) investigation. Intercomparison between these data-sets, and correlation with the broad K-T database available in the literature, indicate that the chemistry of the boundary clays is not globally constant. Variations are more common than similarities, both of absolute concentrations, and interelement ratios. For example, the chondrite normalized platinum-group elements (PGE) patterns of Stevns Klint are not like those of Woodside Creek, with the Pt/Os ratios showing the biggest variation. These differences in PGE patterns are difficult to explain by secondary alteration of a layer that was originally chemically homogeneous, especially for elements of such dubious crustal mobility as Os and Ir. The data also show that enhanced PGE concentrations, with similar trends to those of the boundary layers, occur in the Cretaceous sediments below the actual boundary at Stevns Klint and all three the New Zealand localities. This confirms the observations of others that the geochemistry of the boundary layers apparently does not record a unique component. It is suggested that terrestrial processes, eg. an extended period of Late Cretaceous volcanism can offer a satisfactory explanation for the features of the K-T geochemical anomaly. Such models would probably be more consistent with the observed stepwise, or gradual, palaeontological changes across this boundary, than the instant catastrophe predicated by the impact theory.

Tredous, M.↗

Post-impact alteration of the Manson impact structure

Core materials from the Manson impact site (Manson, Iowa) are examined in order to evaluate post-impact alteration processes. Diagenetic interpretation of post-impact events is based on petrologic, mineralogic, and geochemical investigation of core materials including the following: target strata, disturbed and disrupted strata, ejecta, breccias, microbreccias, and impact melt. The diagenetic study utilizes research cores obtained by the continental scientific drilling project (CSDP) at the Manson structure, as well as core and cuttings of related materials. Samples include impactites (breccias, microbreccias, and melt material), crater fill material (sedimentary clast breccias), disturbed and disrupted target rocks, and reference target material (Amoco Eisheid No. 1 materials). The study of multiple cores will permit development of a regional picture of post-impact thermal history. The specific objectives are as follows: (1) provide a detailed description of authigenic and alteration mineralogy from diverse lithologies encountered in research drill cores at the Manson impact structure, and (2) identify and relate significant post-impact mineral alteration to post-impact thermal regime (extent and duration). Results will provide mineralogical and geochemical constraints on models for post-impact processes including the following: infilling of the crater depression; cooling and hydrothermal alteration of melt rocks; and subsequent long-term, low-temperature alteration of target rocks, breccias, and melt rocks. Preliminary petrologic and x-ray diffraction examination of fracture linings and void fillings from research core M1 indicate the presence of quartz, chlorite, mixed-layer clays, gypsum/anhydrite, calcite, and minor pyrite.

Crossey, L. J.↗

The Effect of Large Melt Fraction on the Deformation Behavior of Peridotite: Implications for the Rheology of Io' Mantle

One key constraint needed for refinement of the interior geochemical and geodynamic models of Io is the viscosity of the convecting partially- molten silicate mantle. To date, laboratory studies of partially molten mantle rocks have reached melt fractions up to approx.0.12, a value much smaller than thought to be appropriate for the asthenosphere of Io where the degree of partial melting may be 0.15 0.40 or higher. Therefore, we have performed a series of high temperature, triaxial compressive creep experiments on dry synthetic peridotites in a gas medium apparatus at a confining pressure of 300 MPa and temperatures from 1473 to 1573 K in order to understand the influence of large amounts of melt (0.15 < phi < 0.40) on the rheological behavior of partially molten rocks.

Scott, T.↗

Mineral induced phosphorylation of glycolate ion--a metaphor in chemical evolution

Bilateral surface-active minerals with excess positive charge concentrate glycolate and trimetaphosphate ion from l0(-3) m aqueous solution to half-saturation of the internal surface sites, and induce phosphorylation of glycolate ion in the mineral with trimetaphosphate, sorbed from l0(-2) m solution. By utilizing reactants from dilute solution at near-neutral pH, and eliminating the need for participating organic nitrogen compounds, the reaction comprises several elements considered necessary for geochemical realism in models for molecular evolution.

NASA Discipline Exobiology↗

Vadose Zone Transport of Tritium and Nitrate under Ponded Water Conditions

Vadose zone transport of tritium and nitrate can be important considerations at radioactive waste sites, landfills, or areas with industrial impacts. These contaminants are of particular concern because they typically have a relatively higher mobility in the subsurface compared to other compounds. Here, we describe a semiarid site with tritium and nitrate contamination involving a manmade ponded water source above a thick unsaturated zone at Los Alamos National Laboratory in New Mexico, USA. This study demonstrates the value of vadose zone flow and transport modeling for the development of field investigation plans (i.e., identifying optimal borehole locations and depths for contaminant characterization), and how a combination of modeling with isotope and geochemical measurements can provide insight into how tritium and nitrate transport in the vadose zone in semiarid environments. Modeling results suggest that at this location, tritium transport is well predicted by classical multiphase theory. Our work expands the demonstrated usefulness of a standard tritium conceptual model to sites with ponded surface conditions and agrees with previous results where a standard model was able to explain the evolution of a tritium plume at an arid waste disposal site. In addition, depth-based analyses of δ 18 O and δ 2 H of pore waters helped confirm the extent of pond infiltration into the vadose zone, and the δ 15 N of nitrate showed that the contaminant release history of the site was different than originally assumed.

58 GEOSCIENCES↗

Integrating Tide‐Driven Wetland Soil Redox and Biogeochemical Interactions Into a Land Surface Model

Abstract Redox processes, aqueous and solid‐phase chemistry, and pH dynamics are key drivers of subsurface biogeochemical cycling and methanogenesis in terrestrial and wetland ecosystems but are typically not included in terrestrial carbon cycle models. These omissions may introduce errors when simulating systems where redox interactions and pH fluctuations are important, such as wetlands where saturation of soils can produce anoxic conditions and coastal systems where sulfate inputs from seawater can influence biogeochemistry. Integrating cycling of redox‐sensitive elements could therefore allow models to better represent key elements of carbon cycling and greenhouse gas production. We describe a model framework that couples the Energy Exascale Earth System Model (E3SM) Land Model (ELM) with PFLOTRAN biogeochemistry, allowing geochemical processes and redox interactions to be integrated with land surface model simulations. We implemented a reaction network including aerobic decomposition, fermentation, sulfate reduction, sulfide oxidation, methanogenesis, and methanotrophy as well as pH dynamics along with iron oxide and iron sulfide mineral precipitation and dissolution. We simulated biogeochemical cycling in tidal wetlands subject to either saltwater or freshwater inputs driven by tidal hydrological dynamics. In simulations with saltwater tidal inputs, sulfate reduction led to accumulation of sulfide, higher dissolved inorganic carbon concentrations, lower dissolved organic carbon concentrations, and lower methane emissions than simulations with freshwater tidal inputs. Model simulations compared well with measured porewater concentrations and surface gas emissions from coastal wetlands in the Northeastern United States. These results demonstrate how simulating geochemical reaction networks can improve land surface model simulations of subsurface biogeochemistry and carbon cycling.

54 ENVIRONMENTAL SCIENCES↗

Metabolic skinflint or spendthrift? Insights into ground sloth integument and thermophysiology revealed by biophysical modeling and clumped isotope paleothermometry

Abstract Remains of megatheres have been known since the 18th -century and were among the first megafaunal vertebrates to be studied. While several examples of preserved integument show a thick coverage of fur for smaller ground sloths living in cold climates such as Mylodon and Nothrotheriops , comparatively very little is known about megathere skin. Assuming a typical placental mammal metabolism, it was previously hypothesized that megatheres would have had little-to-no fur as they achieved giant body sizes. Here the “hairless model of integument” is tested using geochemical analyses to estimate body temperature to generate novel models of ground sloth metabolism, fur coverage, and paleoclimate with Niche Mapper software. The simulations assuming metabolic activity akin to those of modern xenarthrans suggest that sparse fur coverage would have resulted in cold stress across most latitudinal ranges inhabited by extinct ground sloths. Specifically, Eremotherium predominantly required dense 10 mm fur with implications for seasonal changes of coat depth in northernmost latitudes and sparse fur in the tropics; Megatherium required dense 30 mm fur year-round in its exclusive range of cooler, drier climates; Mylodon and Nothrotheriops required dense 10–50 mm fur to avoid thermal stress, matching the integument remains of both genera, and further implying the use of behavioral thermoregulation. Moreover, clumped isotope paleothermometry data from the preserved teeth of four genera of ground sloth yielded reconstructed body temperatures lower than those previously reported for large terrestrial mammals (29 ± 2°–32 ± 3° C). This combination of low metabolisms and thick fur allowed ground sloths to inhabit various environments.

Deak, Michael D.↗

Cardiopulmonary Inflammatory Response to Meteorite Dust Exposures - Implications for Human Health on Earth and Beyond

This year marks the 50th anniversary of Apollo 11, the first time humans set foot on the Moon. The Apollo missions not only help answer questions related to our solar system, they also highlight many hazards associated with human space travel. One major concern is the effect of extraterrestrial dust on astronaut health. In an effort to expand upon previous work indicating lunar dust is respirable and reactive, the authors initiated an extensive study evaluating the role of a particulate’s innate geochemical features (e.g., bulk chemistry, internal composition, morphology, size, and reactivity) in generating adverse toxicological responses in vitro and in vivo. To allow for a broader planetary and geochemical assessment, seven samples were evaluated: six meteorites from either the Moon, Mars, or Asteroid 4 Vesta and a terrestrial basalt analogue. Even with the relatively small geochemical differences (all samples basaltic in nature), significant difference in cardiopulmonary inflammatory markers developed in both single exposure and multiple exposure studies. More specifically: 1) the single exposure studies reveal relationships between toxicity and a meteorite sample’s origin, its pre-ejected state (weathered versus un-weathered), and geochemical features (e.g. bulk iron content) and 2) multiple exposure studies reveal a correlation with particle derived reactive oxygen species (ROS) formation and neutrophil infiltration. Extended human exploration will further increase the probability of inadvertent and repeated exposures to extraterrestrial dusts. This comprehensive dataset allows for not only the toxicological evaluation of extraterrestrial materials but also clarifies important correlations between geochemistry and health. The utilization of an array of extraterrestrial samples from Moon, Mars, and asteroid 4Vesta will enable the development of a geochemical based toxicological hazard model that can be used for: 1) mission planning, 2) rapid risk assessment in cases of unexpected exposures, and 3) evaluation of the efficacy of various in situ techniques in gauging surface dust toxicity. Furthermore, by better understanding the importance of geochemical features on exposure related health outcomes in space, it is possible to better understand of the deleterious nature of dust exposure on Earth.

Harrington, A. D.↗

Offshore application of landslide susceptibility mapping using gradient-boosted decision trees: a Gulf of Mexico case study

Abstract Among natural hazards occurring offshore, submarine landslides pose a significant risk to offshore infrastructure installations attached to the seafloor. With the offshore being important for current and future energy production, there is a need to anticipate where future landslide events are likely to occur to support planning and development projects. Using the northern Gulf of Mexico (GoM) as a case study, this paper performs Landslide Susceptibility Mapping (LSM) using a gradient-boosted decision tree (GBDT) model to characterize the spatial patterns of submarine landslide probability over the United States Exclusive Economic Zone (EEZ) where water depths are greater than 120 m. With known spatial extents of historic submarine landslides and a Geographic Information System (GIS) database of known topographical, geomorphological, geological, and geochemical factors, the resulting model was capable of accurately forecasting potential locations of sediment instability. Results of a permutation modelling approach indicated that LSM accuracy is sensitive to the number of unique training locations with model accuracy becoming more stable as the number of training regions was increased. The influence that each input feature had on predicting landslide susceptibility was evaluated using the SHapely Additive exPlanations (SHAP) feature attribution method. Areas of high and very high susceptibility were associated with steep terrain including salt basins and escarpments. This case study serves as an initial assessment of the machine learning (ML) capabilities for producing accurate submarine landslide susceptibility maps given the current state of available natural hazard-related datasets and conveys both successes and limitations.

Dyer, Alec S. (ORCID:0000000219813904)↗

Controls on Barite (BaSO 4 ) Precipitation in Unconventional Reservoirs

Barite (BaSO 4 ) precipitation is one of the most ubiquitous examples of secondary sulfate mineral scaling in shale oil and gas reservoirs. Often, a suite of chemical additives is used during fracturing operations to inhibit the accumulation of mineral scales, though their efficacy is widely varied and poorly understood. This study combines experimental data and multi-component numerical reactive transport modeling to offer a more comprehensive understanding of the geochemical behavior of barite accumulation in shale matrices under conditions typical of fracturing operations. A variety of additives and conditions are individually tested in batch reactor experiments to identify the factors controlling barite precipitation. Our experimental results demonstrate a pH dependence in the rate of barite precipitation, which we use to develop a predictive model including a pH-dependent term that satisfactorily reproduces our observations. Further, this model is then extended to consider the behavior of three major shale samples of highly variable mineralogy (Eagle Ford, Marcellus, and Barnett). This data-validated model offers a reliable tool to predict and ultimately mitigate against secondary mineral accumulation in unconventional shale reservoirs.

54 ENVIRONMENTAL SCIENCES↗

Flow Directions and Ages of Subsurface Water in a Salt Marsh System Constrained by Isotope Tracing

Salt marshes are dynamic hydrologic systems where terrestrial groundwater, terrestrial surface water, and seawater mix due to bi-directional flows and pressure gradients. Due to the counteracting terrestrial and marine forcings that control these environments, we do not comprehensively understand water fluxes in these complex coastal systems. To understand the water sources, flow directions, and velocities in salt marsh porewater, we employed a combination of geochemical tracers and analytical models across a hillslope-to-salt marsh continuum in a salt marsh experiencing daily inundation of estuarine surface water (SW) from tides and mixing of fresh seasonal groundwater. We used tritium ( 3 H) as a hydrologic tracer to assess porewater ages and stable water isotope ($\delta$ 2 H and $\delta$ 18 O) analyses to separate isotopically distinct estuarine and terrestrial groundwater across different depths and landscape positions in the study transect. We employed electrical conductivity to constrain the role of source mixing and evapotranspiration in salt marsh hydrology. Salinity and stable isotopes revealed that transpiration, rather than evaporation, increased subsurface water salinity to concentrations above estuarine SW during summer. Elevated salinity at depth indicated that salt marsh subsurface water is recharged during the dry growing season. Seasonal recharge patterns drive long-term deep subsurface water dynamics across the salt marsh, with 3 H ages of 3–7 years, and daily tidal cycles drive short-term shallow porewater dynamics with 3 H ages of 0 ± 3.6 years. Furthermore, our conceptual understanding of the spatiotemporal changes in SW-subsurface water interactions at the terrestrial-marine interface quantifies the hydrological constraints we are missing to improve our understanding of biogeochemical cycles within the salt marsh.

54 ENVIRONMENTAL SCIENCES↗

Feasibility of enhancing carbon sequestration and stock capacity in temperate and boreal European forests via changes to management regimes

Forest management practices might act as nature-based methods to remove CO 2 from the atmosphere and slow anthropogenic climate change and thus support an EU forest-based climate change mitigation strategy. However, the extent to which diversified management actions could lead to quantitatively important changes in carbon sequestration and stocking capacity at the tree level remains to be thoroughly assessed. To that end, we used a state-of-the-science bio-geochemically based forest growth model to simulate effects of multiple forest management scenarios on net primary productivity (NPP) and potential carbon woody stocks (pCWS) under twenty scenarios of climate change in a suite of observed and virtual forest stands in temperate and boreal European forests. Previous modelling experiments indicated that the capacity of forests to assimilate and store atmospheric CO 2 in woody biomass is already being attained under business-as-usual forest management practices across a range of climate change scenarios. Nevertheless, we find that on the long-term, with increasing atmospheric CO 2 concentration and warming, managed forests show both higher productivity capacity and a larger potential pool size of stored carbon than unmanaged forests as long as thinning and tree harvesting are of moderate intensity.

54 ENVIRONMENTAL SCIENCES↗