Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “transport models”

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,171 records · Page 65

Investigating the Sensitivity of Surface-Level Nitrate Seasonality in Antarctica to Primary Sources Using a Global Model

Determining the sources of total nitrate (TNIT NO-3 + HNO3) reaching Antarctica is a long-standing challenge. Here we analyze the monthly sensitivity of surface-level TNIT in Antarctica to primary sources using a global 3-D chemical transport model, GEOS-Chem, and its adjoint. Modeled seasonal variation of TNIT concentrations shows good agreement with several measurement studies, given that the lack of post-depositional processing in the model leads to an expected underestimate of maximum values in November through January. Remote NOx sources have the greatest impact May-July, during when the model background concentrations are sensitive to NOx emissions from fossil fuel combustion, soil, and lightning originating from 25S to 65S. In this season, NOx is transported to Antarctica as TNIT, which is formed above continental source regions at an altitude of 5-11 km. In other seasons, more NOx is transported as a reservoir species (e.g., peroxyacetyl nitrate, PAN) through the free troposphere, transforming into TNIT within a cone of influence that extends to 35S and above 4 km altitude. Photolysis of PAN over Antarctica is the main driver of modeled NOx seasonality. Stratospheric production and loss of tracers are relatively unimportant in monthly sensitivities in GEOS-Chem, driving only a few percent of surface level variability of TNIT. A small peak concentration in August is captured by the model, although some measured values in August fall outside the range of simulated concentrations. Modifications to the model to represent sedimentation of polar stratospheric clouds (PSCs) lead to increased surface level August TNIT concentrations. However, this simple representation does not explicitly account for PSC particle deposition or disappearance of the tropopause in the middle of winter, and thus the influence of stratospheric nitrate sources estimated in this study is likely a lower bound.

stratosphere↗

The solar cycle; Proceedings of the National Solar Observatory/Sacramento Peak 12th Summer Workshop, Sunspot, NM, Oct. 15-18, 1991

Attention is given to a flux-transport model, the effect of fractal distribution on the evolution of solar surface magnetic fields, active nests on the sun, magnetic flux transport in solar active regions, recent advances in stellar cycle research, magnetic intermittency on the sun, a search for existence of large-scale motions on the sun, and new solar cycle data from the NASA/NSO spectromagnetograph. Attention is also given to the solar cycle variation of coronal temperature during cycle 22, the distribution of the north-south asymmetry for the various activity cycles, solar luminosity variation, a two-parameter model of total solar irradiance variation over the solar cycle, the origin of the solar cycle, nonlinear feedbacks in the solar dynamo, and long-term dynamics of the solar cycle.

Harvey, Karen L.↗

Wind tunnel investigation of a large-scale upper surface blown-flap model having four engines

Investigations were conducted in the Ames 40- by 80-Foot Wind Tunnel to determine the aerodynamic characteristics of a large-scale subsonic jet transport model with an upper surface blown flap system. The model had a 25 deg swept wing of aspect ratio 7.28 and four turbofan engines. The lift of the flap system was augmented by turning the turbofan exhaust over the Coanda surface. Results were obtained for several flap deflections with several wing leading-edge configurations at jet momentum coefficients from 0 to 4.0. Three-component longitudinal data are presented with four engines operating. In addition, longitudinal and lateral data are presented with an engine out. The maximum lift and stall angle of the four engine model were lower than those obtained with a two engine model that was previously investigated. The addition of the outboard nacelles had an adverse effect on these values. Efforts to improve these values were successful. A maximum lift of 8.8 at an angle-of-attack of 27 deg was obtained with a jet thrust coefficient of 2 for the landing flap configuration.

Aoyagi, K.↗

Comparison of distributed reacceleration and leaky-box models of cosmic-ray abundances (Z = 3-28)

A large collection of elemental and isotopic cosmic-ray data has been analyzed using the leaky-box transport model with and without reacceleration in the interstellar medium. Abundances of isotopes and elements with charges Z = 3-28 and energies E = 10 MeV/nucleon-1 TeV/nucleon were explored. Our results demonstrate that reacceleration models make detailed and accurate predictions with the same number of parameters or fewer as standard leaky-box models. Ad hoc fitting parameters in the standard model are replaced by astrophysically significant reacceleration parameters. Distributed reacceleration models explain the peak in secondary-to-primary ratios around 1 GeV/nucleon. They diminish the discrepancy between rigidity-dependent leakage and energy-independent anisotropy. They also offer the possibility of understanding isotopic anomalies at low energy.

Letaw, John R.↗

Source and Variability of Formaldehyde (HCHO) at Northern High Latitude: an Integrated Satellite, Aircraft, and Model Study

Here we use satellite observations of formaldehyde (HCHO)vertical column densities (VCD) from the TROPOspheric Monitoring Instrument (TROPOMI), aircraft measurements, combined with a nested regional chemical transport model (GEOS-Chem at 0.5°×0.625° resolution), to understand the variability and sources of summertime HCHO better in Alaska. We first evaluate GEOS-Chem with in-situ airborne measurements during Atmospheric Tomography Mission 1 (ATom-1) aircraft campaign. We show reasonable agreement between observed and modeled HCHO, isoprene, monoterpenes and the sum of methyl vinyl ketone and methacrolein (MVK+MACR) in continental boundary layer. In particular, HCHO profiles show spatial homogeneity in Alaska, suggesting a minor contribution of biogenic emissions to HCHO VCD. We further examine the TROPOMI HCHO product in Alaska summer, reprocessed by GEOS-Chem model output for a priori profiles and shape factors. For the year with low wildfire activity (e.g.,2018), we find that HCHO VCDs are largely dominated by background HCHO (58-71%), with minor contributions from wildfires (20-32%) and biogenic VOC emissions (8-10%). For the year with intense wildfires (e.g.,2019), summertime HCHO VCD is dominated by wildfire emissions (50-72%), with minor contributions from background (22-41%) and biogenic VOCs (6-10%). In particular, the model indicates a major contribution of wildfires from direct emissions of HCHO, instead of secondary production of HCHO from oxidation of larger VOCs. We find that the column contributed by biogenic VOC is often small and below the TROPOMI detection limit, in part due to the slow HCHO production from isoprene oxidation under low NOx conditions. This work highlights challenges for quantifying HCHO and its precursors in remote pristine regions.

HCHO↗

Quantum bath augmented stochastic nonequilibrium atomistic simulations for molecular heat conduction

Classical molecular dynamics (MD) has been shown to be effective in simulating heat conduction in certain molecular junctions since it inherently takes into account some essential methodological components which are lacking in the quantum Landauer-type transport model, such as many-body full force-field interactions, anharmonicity effects and nonlinear responses for large temperature biases. However, the classical MD reaches its limit in the environments where the quantum effects are significant (e.g. with low-temperatures substrates, presence of extremely high frequency molecular modes). Here, we present an atomistic simulation methodology for molecular heat conduction that incorporates the quantum Bose–Einstein statistics into an “effective temperature” in the form of a modified Langevin equation. We show that the results from such a quasi-classical effective temperature MD method deviates drastically when the baths temperature approaches zero from classical MD simulations and the results converge to the classical ones when the bath approaches the high-temperature limit, which makes the method suitable for full temperature range. In addition, we show that our quasi-classical thermal transport method can be used to model the conducting substrate layout and molecular composition (e.g. anharmonicities, high-frequency modes). Anharmonic models are explicitly simulated via the Morse potential and compared to pure harmonic interactions to show the effects of anharmonicities under quantum colored bath setups. Finally, the chain length dependence of heat conduction is examined for one-dimensional polymer chains placed in between quantum augmented baths.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Distribution Development for the RDX Regional Model at Los Alamos National Laboratory - 20385

Representing uncertainty in model inputs often means finding a balance between uncertainty and physical reality. Developing wide distributions may seem conservative in principle, but this approach may lead to unrealistic model results. Characterizing the current state of knowledge of stochastic inputs presents many challenges, especially if the data available are limited or have limited relevance to the site. Relationships among these inputs may also be important to represent but are typically complex or difficult to define. Often special adjustments must be made to account for reduced credibility in particular data. If parameters are strongly related to one another, a correlation structure may be developed for input to the model. Other techniques such as regression models may be used to incorporate relationships between the information available and the desired parameters. The process of developing distributions must consider details of the model in terms of what the distribution is meant to represent. This paper uses the example of a probabilistic fate and transport model for hexahydro-1,3,5-trinitro-1,3,5-triazine (RDX) in the regional aquifer at Los Alamos National Laboratory (LANL). For many parameters, a single draw is applied to all space and time over which the model is run, for a single iteration. This simplification is often made for many reasons, and can often be beneficial, but also adds additional complexity in the distribution development process. Defining the distributional goals as they relate to the modeling process is an important step, which should take place prior to evaluation of the data. Usually, the distributions developed are meant to characterize the average value of the parameter over the spatial and temporal domain of the model. Distribution development requires consideration of many sources of information on the parameter where available, ideally from multiple references. Examples of different sources include data from different references but also from different conditions, measurement methods, or experimental types. Depending on these conditions and the reliability or relevance of particular references, different sources of data may each contribute valuable information but have varying relevance to the site. In these cases, weighting data unequally is a useful way to incorporate this information. As an example, aqueous dispersivity data are available for a variety of rock types. Only a few values are available for the desired rock type, and this is not enough to develop a distribution. Therefore, dispersivity values from other rock types are included in distribution development but are down-weighted such that the best data have the most influence on the distribution developed. In another case, K{sub d} distributions in the model are meant to represent a known composition of multiple soil types. Data from these materials are weighted accordingly to develop a distribution for the weighted average K{sub d} across soil types. Other cases include varying reliability of different sources, and weighting data according to the confidence in these sources. In some cases, input parameters are correlated with one another. An example is advective porosity, which is positively correlated with total porosity and must be less than total porosity. Paired data with both parameters must exist to discern the relationship between the two parameters and if it is necessary to build a correlation structure into the model. In general, correlated parameters may be represented in the model by a multivariate distribution, or perhaps more desirably, capturing the correlation within the developed distributions which may be treated as independent from one another. In the example of porosity, this can be done by transforming advective porosity into a proportion of total porosity which may be drawn independently from the distribution of total porosity. This paper explores how complex data and correlations can be incorporated to meet the distributional goals of the model, using the RDX regional model as a detailed example. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Measurements of H Λ 3 and H Λ 4 Lifetimes and Yields in Au + Au Collisions in the High Baryon Density Region

We report precision measurements of hypernuclei $^3_ΛH$ and $^4_ΛH$ lifetimes obtained from Au + Au collisions at $\sqrt{sNN}$ = 3.0 GeV and 7.2 GeV collected by the STAR experiment at the Relativistic Heavy Ion Collider, and the first measurement of $^3_ΛH$ and $^4_ΛH$ midrapidity yields in Au + Au collisions at $\sqrt{sNN}$ = 3.0 GeV. $^3_ΛH$ and $^4_ΛH$, being the two simplest bound states composed of hyperons and nucleons, are cornerstones in the field of hypernuclear physics. Their lifetimes are measured to be 221 ± 15 (stat) ± 19 (syst) ps for $^3_ΛH$ and 218 ± 6 (stat) ± 13 (syst) ps for $^4_ΛH$. The p T -integrated yields of $^3_ΛH$ and $^4_ΛH$ are presented in different centrality and rapidity intervals. It is observed that the shape of the rapidity distribution of $4_ΛH$ is different for 0%–10% and 10%–50% centrality collisions. Thermal model calculations, using the canonical ensemble for strangeness, describes the $^3_ΛH$ yield well, while underestimating the $^4_ΛH$ yield. Transport models, combining baryonic mean-field and coalescence (jam) or utilizing dynamical cluster formation via baryonic interactions (phqmd) for light nuclei and hypernuclei production, approximately describe the measured $^3_ΛH$ and $^4_ΛH$ yields. Our measurements provide means to precisely assess our understanding of the fundamental baryonic interactions with strange quarks, which can impact our understanding of more complicated systems involving hyperons, such as the interior of neutron stars or exotic hypernuclei.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Electrochemical Measurement of Water Transport Numbers in Anion-Exchange Membranes

Anion-exchange membranes (AEMs) are a possible replacement for perfluorosulfonic-acid membranes in energy-conversion devices, primarily due to the hydroxide mobile ion allowing the devices to operate in alkaline conditions with less expensive electrocatalysts. However, the transport properties of AEMs remain understudied, especially electro-osmosis. In this work, an electrochemical technique, where the open-circuit voltage is measured between two ends of a membrane maintained at different relative humidities, is used to determine the water transport number of various ionomers, including Versogen and Sustainion AEMs and Nafion cation-exchange membrane (CEM), as a function of water content and temperature. In addition, the CEMs and AEMs are examined in differing single-ion forms, specifically proton and sodium (CEM) and hydroxide and carbonate (AEM). Carbonate-form AEMs have the highest transport number (∼11), followed by sodium-form CEMs (∼8), hydroxide-form AEMs (∼6), and proton-form CEMs (∼3). Finally, a multicomponent transport model based on the Stefan-Maxwell-Onsager framework of binary interactions is used to develop a link between water transport number and water-transport properties, extracting a range for the unmeasured membrane water permeability of Versogen as a function of water content.

30 DIRECT ENERGY CONVERSION↗

Tracer Evolution at High Latitudes in Winter (SOLVE)

Prior to the SOLVE deployment, we modified and ran a chemical transport model (CTM) to set up a stratospheric climatology for 1 September 1999. During and immediately following the SOLVE deployments, we ran the model in near-real-time and used these results in conjunction with SOLVE observations to draw conclusions about stratospheric transport.

Plumb, R. Alan↗

Planetary and Interplanetary Environmental Models for Radiation Analysis

The essence of environmental modeling is presented as suited for radiation analysis purposes. The variables of fundamental importance for radiation environmental assessment are discussed. The characterization is performed by dividing modeling into three areas, namely the interplanetary medium, the circumplanetary environment, and the planetary or satellite surface. In the first area, the galactic cosmic rays (GCR) and their modulation by the heliospheric magnetic field as well as and solar particle events (SPE) are considered, in the second area the magnetospheres are taken into account, and in the third area the effect of the planetary environment is also considered. Planetary surfaces and atmospheres are modeled based on results from the most recent targeted spacecraft. The results are coupled with suited visualization techniques and radiation transport models in support of trade studies of health risks for future exploration missions.

DeAngelis, G.↗

Aeolian Sediment Transport through Large Patches of Roughness in the Atmospheric Inertial Sublayer

Roughness influences the flux of wind driven sand transport. In this paper, we report on sediment transport measurements for four different surface roughness configurations composed of the same size (solid) roughness elements in the atmospheric inertial sublayer (ISL). Results of these tests indicate that sediment transport rates through patches of roughness in the atmospheric inertial sublayer are to a large extent controlled and scale proportionally with the roughness density (lambda = n b h/S, where n is number of elements of breadth b and height h in area S) of the surface. However, element size apparently increases the magnitude of the reduction beyond that attributable to lambda. A sediment transport model that incorporates the effect of shear stress partitioning appears to predict reasonably well the effect of roughness on sand transport in the cases where the roughness elements are less than or equal to 0.10 m in height. However, when the dimensions of the roughness itself are equivalent to or are greater than the range of saltation lengths (vertical and horizontal), additional interactions of the elements with the saltation cloud appear to reduce the transport efficiency.

Gillies, John A.↗

Modern Era Retrospective Restrospective-Analysis for Research and Applications (MERRA) Data and Services at the GES DISC

The Modern Era Retrospective-analysis for Research and Applications (MERRA) dataset is a NASA satellite era, 30 year (1979 - present), reanalysis using the Goddard Earth Observing System Data Assimilation System, Version 5 (GEOS-5). The project, run out of NASA's Global Modeling and Assimilation Office at Goddard Space Flight Center, provides the science and application communities with a state-of-the-art global analysis with emphasis on improved estimates of the hydrological cycle over a broad range of weather and climate time scales. MERRA products are generated as a long-term synthesis that places the NASA EOS suite of observations in a climate context. The MERRA analysis is performed at a horizontal resolution of 2/3 longitude x 1/2 latitude (540x361 global gridpoints) with observational analyses every 6 hours. The MERRA output data will include 3 dimensional state fields for every 6 hourly analysis cycle on 42 pressure levels (or 72 terrain following model coordinate levels) from the surface through the stratosphere. Several data products are specifically designed to support chemistry and stratosphere transport modeling. The 2 dimensional surface and atmospheric diagnostics (numbering 259) are being stored on the native grid at 1 hourly intervals. These include radiation and vertical integrals of the atmosphere for water and energy budget studies and also surface diagnostics where the diurnal cycle is important. The one hourly surface and near surface data product will also facilitate research on the integrated analysis of Earth system observations in the land, ocean and cryosphere. The MERRA products are archived and distributed by the Goddard Earth Sciences Data and Information Services Center (GES DISC) through its Modeling DISC Web (MDISC) portal. Multiple data access methods and services are available for MERRA data through MDISC: (1) Mirador offers a quick, comprehensive search of MERRA and all GES DISC archived data holdings, allowing searches on keywords, location names or latitude/longitude box, and date/time, with responses within a few seconds. (2) Giovanni is a GES DISC developed Web application that provides data visualization and analysis online. Giovanni features popular visualizations such as latitude-longitude maps, animations, cross sections, profiles, time series, etc. and some basic statistical analysis functions such as scatter plots and correlation coefficient maps. Users are able to download results in several different formats, including Google Earth. (3) On-the-fly parameter subsetting of data within a spatial/temporal window is provided through a simple select and click Web page. (4) MERRA data are also available via OPeNDAP, GrADS Data Server (GDS) and can be converted to netCDF on the fly.

Berrick, Stephen W.↗

A Satellite-Based Estimate of Combustion Aerosol Cloud Microphysical Effects over the Arctic Ocean

Climate predictions for the rapidly changing Arctic are highly uncertain, largely due to a poor understanding of the processes driving cloud properties. In particular, cloud fraction (CF) and cloud phase (CP) have major impacts on energy budgets, but are poorly represented in most models, often because of uncertainties in aerosol–cloud interactions. Here, we use over 10 million satellite observations coupled with aerosol transport model simulations to quantify large scale microphysical effects of aerosols on CF and CP over the Arctic Ocean during polar night, when direct and semi-direct aerosol effects are minimal. Combustion aerosols over sea ice are associated with very large (~10W/sq.m) differences in longwave cloud radiative effects at the sea ice surface. However, co-varying meteorological changes on factors such as CF likely explain the majority of this signal. For example, combustion aerosols explain at most 40% of the CF differences between the full dataset and the clean-condition subset, compared to between 57% and 91% of the differences that can be predicted by co-varying meteorology. After normalizing for meteorological regime, aerosol microphysical effects have small but significant impacts on CF, CP, and precipitation frequency on an Arctic-wide scale. These effects indicate that dominant aerosol–cloud microphysical mechanisms are related to the relative fraction of liquid-containing clouds, with implications for a warming Arctic.

Zamora, Lauren M.↗

Static Aeroelastic Modeling of a Sub-Scale Wind Tunnel Model with Novel Flap Concept

This paper examines the static aeroelastic modeling of a flexible wind tunnel model equipped with a novel control effector known as the Variable Camber Continuous Trailing Edge Flap (VCCTEF) system. The wind tunnel model is an approximately 10 percent sub-scale version of a wing designed to be equipped on a commerical transport aircraft such as the full-scale NASA Generic Transport Model (GTM). The structure of the model is made highly flexible such that a 10 percent of semi-span wing tip deflection is expected at a design lift condition, and a representation of the VCCTEF concept is incorporated on the model. Static aeroelastic modeling is conducted by using a representative single beam structural finite-element model coupled to a vortex-lattice aerodynamic model. The resulting aeroelastic model of the flexible wind tunnel model is compared against experimental wind tunnel test results from the actual model tested at the University of Washington Aeronautical Laboratory (UWAL). Comparison cases are made using reference and VICON measurements of VCCTEF flap deflections from the wind tunnel test data, and percent errors between the lift curve parameters are presented. The results show up to a maximum percent of error of <10% with regards to lift curve slope values between the aeroelastic model and the UWAL test data. The agreement of lift curve slope values from the aeroelastic model and experimental test results serves as validation for the coupled vortex-lattice finite-element static aeroelastic model.

Aeroelasticity↗

Biogeochemistry simulations for the Salt Marsh Accretion Response to Temperature eXperiment (SMARTX)

Coastal ecosystems have been largely ignored in Earth system models but are important zones for carbon and nutrient processing. Interactions between water, microbes, soil, sediments, and vegetation are important for mechanistic representation of coastal processes and ecosystem function. To investigate the role of these feedbacks, we used a reactive transport model (PFLOTRAN) that has the capability to be connected to the Energy Exascale Earth System Model (E3SM). PFLOTRAN was used to incorporate redox reactions and track chemical species important for coastal ecosystems as well as define simple representations of vegetation dynamics. Our goal was to incorporate oxygen flux, salinity, pH, sulfur cycling, and methane production along with plant-mediated transport of gases and tidal flux. Using porewater profile and incubation data for model calibration and evaluation, we were able to create depth-resolved biogeochemical soil profiles for saltmarsh habitat and use this updated representation to simulate direct and indirect effects of elevated CO2 and temperature on subsurface biogeochemical cycling. We found that simply changing the partial pressure of CO2 or increasing temperature in the model did not fully reproduce observed changes in the porewater profile, but the inclusion of plant or microbial responses to CO2 and temperature manipulations was more accurate in representing porewater concentrations. This indicates the importance of characterizing tightly coupled vegetation-subsurface processes for developing predictive understanding and the need for measurement of plant-soil interactions on the same time scale to understand how hotspots or moments are generated.Included in this data package are PFLOTRAN input (PFLOTRAN input files and chemical database) files for simulating single column biogeochemistry, root, and tide interactions at the Global Change Research Wetland (Kirkpatrick Marsh; Edgewater, MD). The biogeochemical network includes soil organic matter decomposition, nitrogen, iron, and sulfur cycling, and methanogenesis. Reduced species can be oxidized and plant processes include oxygen and nutrient priming, methane release, and nutrient uptake.Inputs:TAI_database.dat - geochemical database for reactions, more information on database structure and variables can be found here https://www.pflotran.org/documentation/user_guide/cards/pages/geochemical_database.htmlswamp.in - input file for biogeochemical network in PFLOTRANswamp_eCO2.in - input file for biogeochemical network in PFLOTRAN with input gas partial pressures/concentrations adjusted for elevated CO2 treatmentsOutputs:swamp_obs_0.tec - hourly porewater concentrations from from multiple depths in the soil columnswamp_eCO2_obs_0.tec - hourly porewater concentrations from multiple depths in the soil column for elevated CO2 treatmentsPFLOTRAN code access: https://github.com/fmyuan/pflotran-elm-interface.git

54 ENVIRONMENTAL SCIENCES↗

Model Predictive Optimal Control of a Time-Delay Distributed-Parameter Systems

This paper presents an optimal control method for a class of distributed-parameter systems governed by first order, quasilinear hyperbolic partial differential equations that arise in many physical systems. Such systems are characterized by time delays since information is transported from one state to another by wave propagation. A general closed-loop hyperbolic transport model is controlled by a boundary control embedded in a periodic boundary condition. The boundary control is subject to a nonlinear differential equation constraint that models actuator dynamics of the system. The hyperbolic equation is thus coupled with the ordinary differential equation via the boundary condition. Optimality of this coupled system is investigated using variational principles to seek an adjoint formulation of the optimal control problem. The results are then applied to implement a model predictive control design for a wind tunnel to eliminate a transport delay effect that causes a poor Mach number regulation.

Nguyen, Nhan↗

Developing a Redox Network for Coastal Saltmarsh Systems in the PFLOTRAN Reaction Model

Coastal ecosystems have been largely ignored in Earth system models but are important zones for carbon and nutrient processing. Interactions between water, microbes, soil, sediments, and vegetation are important for mechanistic representation of coastal processes and ecosystem function. To investigate the role of these feedbacks, we used a reactive transport model (PFLOTRAN) that has the capability to be connected to the Energy Exascale Earth System Model (E3SM). PFLOTRAN was used to incorporate redox reactions and track chemical species important for coastal ecosystems as well as define simple representations of vegetation dynamics. Here our goal was to incorporate oxygen flux, salinity, pH, sulfur cycling, and methane production along with plant-mediated transport of gases and tidal flux. Using porewater profile and incubation data for model calibration and evaluation, we were able to create depth-resolved biogeochemical soil profiles for saltmarsh habitat and use this updated representation to simulate direct and indirect effects of elevated CO 2 and temperature on subsurface biogeochemical cycling. We found that simply changing the partial pressure of CO 2 or increasing temperature in the model did not fully reproduce observed changes in the porewater profile, but the inclusion of plant or microbial responses to CO 2 and temperature manipulations was more accurate in representing porewater concentrations. This indicates the importance of characterizing tightly coupled vegetation-subsurface processes for developing predictive understanding and the need for measurement of plant-soil interactions on the same time scale to understand how hotspots or moments are generated.

54 ENVIRONMENTAL SCIENCES↗