Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “model parameterization”

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 541 records · Page 30
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Plant traits, productivity, biomass and soil properties from forest sites in the Pacific Northwest, 1999–2014

Plant trait measurements are needed for evaluating ecological responses to environmental conditions and for ecosystem process model development, parameterization, and testing. Here, we present a standardized dataset integrating measurements from projects conducted by the Terrestrial Ecosystem Research and Regional Analysis- Pacific Northwest (TERRA-PNW) research group between 1999 and 2014 across Oregon and Northern California, where measurements were collected for scaling and modeling regional terrestrial carbon processes with models such as Biome-BGC and the Community Land Model. The dataset contains measurements of specific leaf area, leaf longevity, leaf carbon and nitrogen for 35 tree and shrub species derived from more than 1,200 branch samples collected from over 200 forest plots, including several AmeriFlux sites. The dataset also contains plot-level measurements of forest composition, structure (e.g., tree biomass), and productivity, as well as measurements of soil structure (e.g., bulk density) and chemistry (e.g., carbon). Publically-archiving regional datasets of standardized, co-located, and geo-referenced plant trait measurements will advance the ability of earth system models to capture species-level climate sensitivity at regional to global scales.

09 BIOMASS FUELS↗
DOE OSTI2022

Thermodynamic modeling of MDEA(aq)-NH 3 (aq)-K 2 CO 3 (aq)-CO 2 (aq) using the Extended UNIQUAC model

Here, the development of advanced processes for carbon dioxide capture requires reliable thermodynamic models to ensure accurate process design. This work presents an optimized parameter set of the Extended UNIQUAC model for advanced mixed solvent solutions containing CO 2 -MDEA-NH 3 -K 2 CO 3 in water. The model was parameterized using vapor–liquid equilibrium, pure component saturation pressure, excess enthalpy, heat of absorption, molar heat capacity, apparent molar heat capacity, and solid–liquid equilibrium. These thermal and equilibrium properties cover temperatures between 253.6 K to 600 K and pressures up to 75.6 bar. The current model represents an extension of previous versions that were suitable for modeling aqueous NH 3 and K 2 CO 3 solutions. This parametrization accounts for molecular interactions between MDEA and these two components, which allows for the design of advanced capture processes. This paper demonstrates that the new model parameters are valid within the range of operational conditions typically found in CO 2 capture processes.

42 ENGINEERING↗
DOE OSTI2020

A long short-term memory embedding for hybrid uplifted reduced order models

In this paper, we introduce an uplifted reduced order modeling (UROM) approach through the integration of standard projection based methods with long short-term memory (LSTM) embedding. Our approach has three modeling layers or components. In the first layer, we utilize an intrusive projection approach to model dynamics represented by the largest modes. The second layer consists of an LSTM model to account for residuals beyond this truncation. This closure layer refers to the process of including the residual effect of the discarded modes into the dynamics of the largest scales. However, the feasibility of generating a low rank approximation tails off for higher Kolmogorov n -width systems due to the underlying nonlinear processes. The third uplifting layer, called super-resolution, addresses this limited representation issue by expanding the span into a larger number of modes utilizing the versatility of LSTM. Therefore, our model integrates a physics-based projection model with a memory embedded LSTM closure and an LSTM based super-resolution model. In several applications, we exploit the use of Grassmann manifold to construct UROM for unseen conditions. We performed numerical experiments by using the Burgers and Navier-Stokes equations with quadratic nonlinearity. Finally, our results show robustness of the proposed approach in building reduced order models for parameterized systems and confirm the improved trade-off between accuracy and efficiency.

42 ENGINEERING↗
DOE OSTI2015

Representing northern peatland microtopography and hydrology within the Community Land Model

Predictive understanding of northern peatland hydrology is a necessary precursor to understanding the fate of massive carbon stores in these systems under the influence of present and future climate change. Current models have begun to address microtopographic controls on peatland hydrology, but none have included a prognostic calculation of peatland water table depth for a vegetated wetland, independent of prescribed regional water tables. We introduce here a new configuration of the Community Land Model (CLM) which includes a fully prognostic water table calculation for a vegetated peatland. Our structural and process changes to CLM focus on modifications needed to represent the hydrologic cycle of bogs environment with perched water tables, as well as distinct hydrologic dynamics and vegetation communities of the raised hummock and sunken hollow microtopography characteristic of peatland bogs. The modified model was parameterized and independently evaluated against observations from an ombrotrophic raised-dome bog in northern Minnesota (S1-Bog), the site for the Spruce and Peatland Responses Under Climatic and Environmental Change experiment (SPRUCE). Simulated water table levels compared well with site-level observations. The new model predicts hydrologic changes in response to planned warming at the SPRUCE site. At present, standing water is commonly observed in bog hollows after large rainfall events during the growing season, but simulations suggest a sharp decrease in water table levels due to increased evapotranspiration under the most extreme warming level, nearly eliminating the occurrence of standing water in the growing season. Simulated soil energy balance was strongly influenced by reduced winter snowpack under warming simulations, with the warming influence on soil temperature partly offset by the loss of insulating snowpack in early and late winter. Furthermore, the new model provides improved predictive capacity for seasonal hydrological dynamics in northern peatlands, and provides a useful foundation for investigation of northern peatland carbon exchange.

58 GEOSCIENCES↗
DOE OSTI2015

Representing northern peatland microtopography and hydrology within the Community Land Model

Predictive understanding of northern peatland hydrology is a necessary precursor to understanding the fate of massive carbon stores in these systems under the influence of present and future climate change. Current models have begun to address microtopographic controls on peatland hydrology, but none have included a prognostic calculation of peatland water table depth for a vegetated wetland, independent of prescribed regional water tables. We introduce here a new configuration of the Community Land Model (CLM) which includes a fully prognostic water table calculation for a vegetated peatland. Our structural and process changes to CLM focus on modifications needed to represent the hydrologic cycle of bogs environment with perched water tables, as well as distinct hydrologic dynamics and vegetation communities of the raised hummock and sunken hollow microtopography characteristic of peatland bogs. The modified model was parameterized and independently evaluated against observations from an ombrotrophic raised-dome bog in northern Minnesota (S1-Bog), the site for the Spruce and Peatland Responses Under Climatic and Environmental Change experiment (SPRUCE). Simulated water table levels compared well with site-level observations. The new model predicts significant hydrologic changes in response to planned warming at the SPRUCE site. At present, standing water is commonly observed in bog hollows after large rainfall events during the growing season, but simulations suggest a sharp decrease in water table levels due to increased evapotranspiration under the most extreme warming level, nearly eliminating the occurrence of standing water in the growing season. Simulated soil energy balance was strongly influenced by reduced winter snowpack under warming simulations, with the warming influence on soil temperature partly offset by the loss of insulating snowpack in early and late winter. The new model provides improved predictive capacity for seasonal hydrological dynamics in northern peatlands, and provides a useful foundation for investigation of northern peatland carbon exchange.

58 GEOSCIENCES↗
DOE OSTI2014

Neutron star masses and radii from quiescent low-mass x-ray binaries

We perform a systematic analysis of neutron star radius constraints from five quiescent low-mass X-ray binaries and examine how they depend on measurements of their distances and amounts of intervening absorbing material, as well as their assumed atmospheric compositions. We construct and calibrate to published results a semi-analytic model of the neutron star atmosphere which approximates these effects for the predicted masses and radii. Starting from mass and radius probability distributions established from hydrogen-atmosphere spectral fits of quiescent sources, we apply this model to compute alternate sets of probability distributions. We perform Bayesian analyses to estimate neutron star mass-radius curves and equation of state (EOS) parameters that best-fit each set of distributions, assuming the existence of a known low-density neutron star crustal EOS, a simple model for the high-density EOS, causality, and the observation that the neutron star maximum mass exceeds 2 M {sub ☉}. We compute the posterior probabilities for each set of distance measurements and assumptions about absorption and composition. We find that, within the context of our assumptions and our parameterized EOS models, some absorption models are disfavored. We find that neutron stars composed of hadrons are favored relative to those with exotic matter with strong phase transitions. In addition, models in which all five stars have hydrogen atmospheres are found to be weakly disfavored. Our most likely models predict neutron star radii that are consistent with current experimental results concerning the nature of the nucleon-nucleon interaction near the nuclear saturation density.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

ARM data-oriented metrics and diagnostics package for climate model evaluation (ARM-DIAGS-V3) version 3

A Python-based metrics and diagnostics package is currently being developed by the ARM Infrastructure Team at Lawrence Livermore National Laboratory to facilitate the use of long-term high frequency measurements from the ARM program in evaluating the regional climate simulation of clouds, radiation, precipitation, and aerosols. This metrics and diagnostics package computes climatological means of targeted climate model simulation and generates tables and plots for comparing the model simulation with ARM observational data. The CMIP model data sets are also included in the package to enable model inter-comparison as demonstrated in Zhang et al. (2018) and Zhang et al. (2020). The mean of the CMIP model can be served as a reference for individual models. Basic performance metrics are computed to measure the accuracy of mean state and variability of climate models. The evaluated physical quantities include cloud fraction, temperature, relative humidity, cloud liquid water path, total column water vapor, precipitation, sensible and latent heat fluxes, aerosol optical depth, and radiative fluxes, with plan to extend to more fields, such as the evaluation of model simulated aerosol physicochemical properties and cloud microphysics properties. Process-oriented diagnostics focusing on aerosol, cloud, and precipitation-related phenomena are also being developed for the evaluation and development of specific model physical parameterizations. In addition to the Southern Great Plains (SGP), North Slope of Alaska (NSA) and Tropical Western Pacific (TWP) atmospheric observatories in the ARMDIAGS version 2.0, the version 3.0 package have extended to the data collected at the ARM Eastern North Atlantic (ENA) atmospheric observatory and the Observation and Modeling of the Green Ocean Amazon (GOAMAZON) field campaign. The metrics and diagnostics package are currently built upon standard Python libraries and additional Python packages developed by DOE (CDAT). The ARM metrics and diagnostic package is available publicly with the hope that it can serve as an easy entry point for climate modelers to compare their models with ARM data. In this report, we first provide an overview of major metrics in section 2. The input data, which constitutes the core content of the metrics and diagnostics package, is summarized in section 3. A user's guide documenting the workflow/structure of the version 3.0 codes and including step-by-step instruction for running the package is described in section 4.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2019

Ensemble-based flood vulnerability assessment for probable maximum flood in a changing environment

The magnitude and frequency of hydro-meteorological extremes are expected to increase in a changing environment in ways that threaten the security of US energy-water assets. These include probable maximum precipitation (PMP) and probable maximum flood (PMF), which are used as hydraulic design standards for highly sensitive infrastructures such as nuclear power plants and main dams. To assess the flood vulnerability due to PMP/PMF, an integrated high-resolution process-based hydro-meteorologic modeling framework was used to develop ensemble-based probabilistic flood maps based on best-available historic observations and future climate projections. A graphics processing unit–accelerated 2-dimensional hydrodynamic model was used to simulate the surface inundation areas corresponding to a total of 120 PMF hydrographs. These ensemble-based PMF maps were compared with flood maps obtained from the conventional deterministic PMP/PMF approach, revealing added information about conditional probability of flooding. Further, a relative sensitivity test was conducted to explore the effects of various factors in the framework, such as meteorological forcings, antecedent hydrologic conditions, reservoir storage, and flood model input resolution and parameters. In conclusion, the proposed framework better illustrates the uncertainties associated with model inputs, parameterization, and hydro-meteorological factors, allowing more informed decision-making for future emergency preparation.

58 GEOSCIENCES↗
DOE OSTI2018

Role of Troposphere-Convection-Land Coupling in the Southwestern Amazon Precipitation Bias of the Community Earth System Model version 1 (CESM1)

The Amazon forest plays an important role in the tropical circulation and global water and biogeochemical cycles. However, the common Amazon precipitation biases in global climate models undermine our ability to understand future changes in the Amazon and their global consequences. The Community Earth System Model exhibits a persistent wet season dry bias over the southwestern Amazon. The bias is also characterized by inaccurate diurnal timing and intensity of rainfall events. Analysis of a diverse set of simulations suggests that the bias originates from local processes. A known land model bias in the partitioning of surface energy fluxes is found not to contribute significantly to the rainfall bias. Key model processes are identified to be the low-level divergent circulation controlling horizontal moisture flow and the sensitivity of the modeled deep convection to the lower tropospheric moisture. The lower tropospheric moisture divergence dictates the spatiotemporal distributions of convective precipitation through parameterized entrainment mixing, which generally improves the rainfall diurnal cycle. However, it also limits convective precipitation over the areas with frequent and large-scale divergence in the lower troposphere. The southwestern Amazon is such an area, particularly in the late afternoon, although uncertainty exists in the low-level flow in the simulation and reanalyses. One way to sustain rainfall in the presence of low-level divergence is a convective trigger based on both the moisture and temperature changes by large-scale processes, which is shown to help in reducing the dry bias.

Sakaguchi, Koichi↗
DOE OSTI2021

Composited structure of non-precipitating shallow cumulus clouds

The normalized distributions of thermodynamic and dynamical variables both within and outside shallow clouds are investigated through a composite algorithm using large-eddy simulations of oceanic and continental cases. The normalized magnitude is maximum near the cloud centre and decreases outwards. While relative humidity (RH) and cloud liquid water ( q l ) decrease smoothly to match the environment, the vertical velocity, virtual potential temperature (θ v ), and potential temperature (θ) perturbations have more complicated behaviour towards the cloud boundary. Below the inversion layer, θ' v becomes negative before the vertical velocity has turned from an updraft to a subsiding shell outside the cloud, indicating the presence of a transition zone where the updraft is negatively buoyant. Due to the downdraft outside the cloud and enhanced horizontal turbulent mixing across the edge, the normalized turbulent kinetic energy (TKE) and horizontal turbulent kinetic energy (HTKE) decrease more slowly from the cloud centre outwards than the thermodynamic variables. The distributions all present asymmetric structures in response to the vertical wind shear, with more negatively buoyant air, stronger downdrafts, and larger TKE on the downshear side. We discuss several implications of the distributions for theoretical models and parameterizations. Positive buoyancy near the cloud base is mostly due to the virtual effect of water vapour, emphasizing the role of moisture in triggering. The mean vertical velocity is found to be approximately half the maximum vertical velocity within each cloud, providing a constraint to achieve possible power-law distributions for some models. Finally, the normalized distributions for different variables are used to estimate the vertical heat and moisture fluxes within clouds. The results suggest that distributions near the cloud edge and variability of maximum perturbations need careful treatment. The fluxes are underestimated in the inversion layer because cloud-top downdrafts cannot be captured well.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

How Low Can You Go? Using Synthetic 3D Imagery to Drastically Reduce Real-World Training Data for Object Detection

Deep convolutional neural networks (DCNNs) currently provide state-of-the-art performance on image classification and object detection tasks, and there are many global security mission areas where such models could be extremely useful. Crucially, the success of these models is driven in large part by the widespread availability of high-quality open source data sets such as Image Net, Common Objects in Context (COCO), and KITTI, which contain millions of images with thousands of unique labels. However, global security relevant objects-of-interest can be difficult to obtain: relevant events are low frequency and high consequence; the content of relevant images is sensitive; and adversaries and proliferators seek to obscure their activities. For these cases where exemplar data is hard to come-by, even fine-tuning an existing model with available data can be effectively impossible. Recent work demonstrated that models can be trained using a combination of real-world and synthetic images generated from 3D representations; that such models can exceed the performance of models trained using real-world data alone; and that the generated images need not be perfectly realistic (Tremblay, et al., 2018). However, this approach still required hundreds to thousands of real-world images for training and fine tuning, which for sparse, global security-relevant datasets can be an unrealistic hurdle. In this research, we validate the performance and behavior of DCNN models as we drive the number of real-world images used for training object detection tasks down to a minimal set. We perform multiple experiments to identify the best approach to train DCNNs from an extremely small set of real-world images. In doing so, we: Develop state-of-the-art, parameterized 3D models based on real-world images and sample from their parameters to increase the variance in synthetic image training data; Use machine learning explainability techniques to highlight and correct through targeted training the biases that result from training using completely synthetic images; and Validate our results by comparing the performance of the models trained on synthetic data to one another, and to a control model created by fine-tuning an existing ImageNet-trained model with a limited number (hundreds) of real-world images.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Quantifying global N 2 O emissions from natural ecosystem soils using trait-based biogeochemistry models

A group of soil microbes plays an important role in nitrogen cycling and N 2 O emissions from natural ecosystem soils. Wedeveloped a trait-based biogeochemical model based on an extant process-based biogeochemistry model, the Terrestrial Ecosystem Model (TEM), by incorporating the detailed microbial physiological processes of nitrification. The effect of ammonia-oxidizing Archaea (AOA),ammonia-oxidizing bacteria (AOB), and nitrite-oxidizing bacteria (NOB) was considered in modeling nitrification. Microbial traits, including microbial biomass and density, were explicitly considered. In addition,nitrogen cycling was coupled with carbon dynamics based on stoichiometry theory between carbon and nitrogen. The model was parameterized using observational data and then applied to quantifying global N 2 O emissions from global terrestrial ecosystem soils from 1990 to 2000. Our estimates of 8.7±1.6 Tg N yr –1 generally agreed with previous estimates during the study period. Tropical forests are a major emitter,accounting for 42 % of the global emissions. The model was more sensitive to temperature and precipitation and less sensitive to soil organic carbon and nitrogen contents. Compared to the model without considering the detailed microbial activities, the new model shows more variations in response to seasonal changes in climate. Our study suggests that further information on microbial diversity and ecophysiology features is needed. The more specific guilds and their traits shall be considered in future soil N 2 O emission quantifications.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2014

Relative Biological Effectiveness Variation Along Monoenergetic and Modulated Bragg Peaks of a 62-MeV Therapeutic Proton Beam: A Preclinical Assessment

Purpose: The biological optimization of proton therapy can be achieved only through a detailed evaluation of relative biological effectiveness (RBE) variations along the full range of the Bragg curve. The clinically used RBE value of 1.1 represents a broad average, which disregards the steep rise of linear energy transfer (LET) at the distal end of the spread-out Bragg peak (SOBP). With particular attention to the key endpoint of cell survival, our work presents a comparative investigation of cell killing RBE variations along monoenergetic (pristine) and modulated (SOBP) beams using human normal and radioresistant cells with the aim to investigate the RBE dependence on LET and intrinsic radiosensitvity. Methods and Materials: Human fibroblasts (AG01522) and glioma (U87) cells were irradiated at 6 depth positions along pristine and modulated 62-MeV proton beams at the INFN-LNS (Catania, Italy). Cell killing RBE variations were measured using standard clonogenic assays and were further validated using Monte Carlo simulations and the local effect model (LEM). Results: We observed significant cell killing RBE variations along the proton beam path, particularly in the distal region showing strong dose dependence. Experimental RBE values were in excellent agreement with the LEM predicted values, indicating dose-averaged LET as a suitable predictor of proton biological effectiveness. Data were also used to validate a parameterized RBE model. Conclusions: The predicted biological dose delivered to a tumor region, based on the variable RBE inferred from the data, varies significantly with respect to the clinically used constant RBE of 1.1. The significant RBE increase at the distal end suggests also a potential to enhance optimization of treatment modalities such as LET painting of hypoxic tumors. The study highlights the limitation of adoption of a constant RBE for proton therapy and suggests approaches for fast implementation of RBE models in treatment planning.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2022

Multilayer interface tracking model of pure tungsten oxidation

Here, we present a numerical model to predict oxide scale growth on tungsten surfaces under exposure to oxygen at high temperatures. The model captures the formation of four thermodynamically-compatible oxide sublayers, WO 2 , WO 2.72 , WO 2.9 , and WO 3 , on top of the metal substrate. Oxide layer growth is simulated by tracking the oxide/oxide and oxide/metal interfaces using a sharp-interface Stefan model coupled to diffusion kinetics. The model is parameterized using selected experimental measurements and electronic structure calculations of the diffusivities of all the oxide subphases involved. We simulate oxide growth at temperatures of 600°C and above, extracting the power law growth exponents in each case, which we find to deviate from classical parabolic growth in several cases. We conduct a comparison of the model predictions with an extensive experimental data set, with reasonable agreement at most temperatures. While many gaps in our understanding still exist, this work is a first attempt at embedding the thermodynamic and kinetic complexity of tungsten oxide growth into a comprehensive mesoscale kinetic model that attempts to capture the essential features of tungsten oxidation to fill existing knowledge gaps and guide and enhance future tungsten oxidation models.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Simultaneous numerical representation of soil microsite production and consumption of carbon dioxide, methane, and nitrous oxide using probability distribution functions

Production and consumption of nitrous oxide (N 2 O), methane (CH 4 ), and carbon dioxide (CO 2 ) are affected by complex interactions of temperature, moisture, and substrate supply, which are further complicated by spatial heterogeneity of the soil matrix. This microsite heterogeneity is often invoked to explain non-normal distributions of greenhouse gas (GHG) fluxes, also known as hot spots and hot moments. To advance numerical simulation of these belowground processes, we expanded the Dual Arrhenius and Michaelis-Menten (DAMM) model, to apply it consistently for all three GHGs with respect to the biophysical processes of production, consumption, and diffusion within the soil, including the contrasting effects of oxygen (O 2 ) as substrate or inhibitor for each process. High-frequency chamber-based measurements of all three GHGs at the Howland Forest (ME, USA) were used to parameterize the model using a multiple constraint approach. The area under a soil chamber is partitioned according to a bivariate lognormal probability distribution function (PDF) of carbon (C) and water content across a range of microsites, which leads to a PDF of heterotrophic respiration and O 2 consumption among microsites. Linking microsite consumption of O 2 with a diffusion model generates a broad range of microsite concentrations of O 2 , which then determines the PDF of microsites that produce or consume CH 4 and N 2 O, such that a range of microsites occurs with both positive and negative signs for net CH 4 and N 2 O flux. Results demonstrate that it is numerically feasible for microsites of N 2 O reduction and CH 4 oxidation to co-occur under a single chamber, thus explaining occasional measurement of simultaneous uptake of both gases. Simultaneous simulation of all three GHGs in a parsimonious modeling framework is challenging, however it increases confidence that agreement between simulations and measurements is based on skillful numerical representation of processes across a heterogeneous environment.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Local-scale Arctic tundra heterogeneity affects regional-scale carbon dynamics

In northern Alaska nearly 65% of the terrestrial surface is composed of polygonal ground, where geomorphic tundra landforms disproportionately influence carbon and nutrient cycling over fine spatial scales. Process-based biogeochemical models used for local to Pan-Arctic projections of ecological responses to climate change typically operate at coarse-scales (1km 2 –0.5°) at which fine-scale (<1km 2 ) tundra heterogeneity is often aggregated to the dominant land cover unit. Here, we evaluate the importance of tundra heterogeneity for representing soil carbon dynamics at fine to coarse spatial scales. We leveraged the legacy of data collected near Utqiagvik, Alaska between 1973 and 2016 for model initiation, parameterization, and validation. Simulation uncertainty increased with a reduced representation of tundra heterogeneity and coarsening of spatial scale. Hierarchical cluster analysis of an ensemble of 21 st -century simulations reveals that a minimum of two tundra landforms (dry and wet) and a maximum of 4km 2 spatial scale is necessary for minimizing uncertainties (<10%) in regional to Pan-Arctic modeling applications.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2016

Intra-binary shock heating of black widow companions

The low-mass companions of evaporating binary pulsars (black widows and similar) are strongly heated on the side facing the pulsar. However, in high-quality photometric and spectroscopic data, the heating pattern does not match that expected for direct pulsar illumination. Here we explore a model where the pulsar power is intercepted by an intra-binary shock (IBS) before heating the low-mass companion. We develop a simple analytic model and implement it in the popular “ICARUS” light curve code. The model is parameterized by the wind momentum ratio β and the companion wind speed f{sub v}v{sub orb}, and assumes that the reprocessed pulsar wind emits prompt particles or radiation to heat the companion surface. We illustrate an interesting range of light curve asymmetries controlled by these parameters. The code also computes the IBS synchrotron emission pattern, and thus can model black widow X-ray light curves. As a test, we apply the results to the high-quality asymmetric optical light curves of PSR J2215+5135; the resulting fit gives a substantial improvement upon direct heating models and produces an X-ray light curve consistent with that seen. The IBS model parameters imply that at the present loss rate, the companion evaporation has a characteristic timescale of τ{sub evap}≈150 Myr. Still, the model is not fully satisfactory, indicating that there are additional unmodeled physical effects.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Extending a land-surface model with <i>Sphagnum</i> moss to simulate responses of a northern temperate bog to whole ecosystem warming and elevated CO<sub>2</sub>

Abstract. Mosses need to be incorporated into Earth system models to better simulate peatland functional dynamics under the changing environment. Sphagnum mosses are strong determinants of nutrient, carbon, and water cycling in peatland ecosystems. However, most land-surface models do not include Sphagnum or other mosses as represented plant functional types (PFTs), thereby limiting predictive assessment of peatland responses to environmental change. In this study, we introduce a moss PFT into the land model component (ELM) of the Energy Exascale Earth System Model (E3SM) by developing water content dynamics and nonvascular photosynthetic processes for moss. The model was parameterized and independently evaluated against observations from an ombrotrophic forested bog as part of the Spruce and Peatland Responses Under Changing Environments (SPRUCE) project. The inclusion of a Sphagnum PFT with some Sphagnum-specific processes in ELM allows it to capture the observed seasonal dynamics of Sphagnum gross primary production (GPP) albeit with an underestimate of peak GPP. The model simulated a reasonable annual net primary production (NPP) for moss but with less interannual variation than observed, and it reproduced aboveground biomass for tree PFTs and stem biomass for shrubs. Different species showed highly variable warming responses under both ambient and elevated atmospheric CO2 concentrations, and elevated CO2 altered the warming response direction for the peatland ecosystem. Microtopography is critical: Sphagnum mosses on hummocks and hollows were simulated to show opposite warming responses (NPP decreasing with warming on hummocks but increasing in hollows), and hummock Sphagnum was modeled to have a strong dependence on water table height. The inclusion of this new moss PFT in global ELM simulations may provide a useful foundation for the investigation of northern peatland carbon exchange, enhancing the predictive capacity of carbon dynamics across the regional and global scales.

54 ENVIRONMENTAL SCIENCES↗