Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “variable flow”

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 469 records · Page 26
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Simulations of divertor designs that spatially separate power and particle exhaust using mid-leg divertor particle pumping

Predictive design modeling of a Dissipation-Focused Divertor for future operation in DIII-D reveals that increasing the poloidal distance of the pump duct entrance from the target surface along the low-field side divertor baffle increases neutral compression and modifies the spatial distribution of power dissipation. With a divertor pump located mid-leg between the target and the X-point, SOLPS-ITER boundary plasma simulations without drifts predict the formation of a dense neutral cloud near the target with > 30x higher neutral compression in detachment, a more stable detachment front located further from the target, and ~25% lower outer midplane separatrix density required for detachment onset, compared to a pump located in the scrape-off layer at the target surface. Up to 19 MW of power flowing into the divertors is modeled using the following two numerical implementations for particle pumping: a specified fraction of particles incident on variable wall sections of the plasma grid is removed from the computational domain (so-called albedo pumping), and a pump duct is modeled which includes dynamics of kinetic neutrals in the duct. The simulations show that the detachment front is located between the divertor target and the X-point and is relatively stable near the pump entrance, without a strong dependence on gas puff rate or injected power. The mid-leg pump design spatially separates the two primary functions of a divertor (power handling and particle exhaust), with the majority of power dissipation occurring near the target plate and particle exhaust taking place further upstream. The benefit of enhanced dissipation using mid-leg pumping comes at the cost of a higher outer midplane separatrix density for a given amount of particle injection.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

The influence of laser power modulation on melt pool dynamics in laser powder bed fusion

While the majority of laser powder bed fusion (LPBF) metal additive manufacturing uses a continuous wave (CW) laser heat source, some commercial applications of LPBF additive manufacturing instead involve the modulation of the laser power on tens-of-microsecond timescales as an adjustable process variable. This article reports the use of in situ, high speed x-ray and optical imaging to probe melt pool fluid flow, defect formation, and nearby powder motion during LPBF with both modulated and CW laser heat sources. We observe melt pool dynamics unique to modulated laser melting even at very high duty cycles that are related to fluctuations in vapor depression depth, complex pore formation mechanisms, and changes to denudation physics when compared to CW melting. These behaviors are present in Ti–6Al–4V, 316L stainless steel, and AL1100 alloys but vary slightly as a function of material, indicating a substantial dependence on the viscosity and surface tension of the liquid metal. While high duty cycles produce weld tracks of comparable quality to CW melting, lower duty cycles introduce substantial defect concentrations. At intermediate duty cycles, careful control of modulation parameters can repeatably and precisely yield one pore per laser pulse, suggesting a method for intentionally inserting engineered porosity at specific sites during an LPBF build.

3D printing↗
DOE OSTI2022

Isolated Building Wake Experiment At Texas Tech University’s Wind Engineering Research Field Lab

Texas Tech University’s (TTU’s) Wind Engineering Research Field Lab (WERFL) building is an idealized rectangular building on the scale of a typical suburban house. It is surrounded almost entirely by unobstructed planes and is built on a turn table. WERFL is ideal for providing measurements of an idealized building under real-world conditions including atmospheric stability and variable relative wind angles. The intent of this data set is for verification and validation of computational fluid dynamics model simulations of flow around buildings. The primary focus of this dataset features 24 3D sonic anemometers (20 downwind and 4 upwind) and 10 2D sonic anemometers (9 downwind and 1 upwind) located in the immediate vicinity of the WERFL building. It also includes a variety of instruments intended to measure the ambient meteorological parameters and the undisturbed wind profile.

17 WIND ENERGY↗
DOE OSTI2020

Metal oxide semiconductor (MOS) controlled devices and methods of making the same

Metal-Oxide-Semiconductor (MOS) controlled semiconductor devices and methods of making the devices are provided. The devices include a gate which controls current flow through channel regions positioned between source/emitter and drain regions of the device. The devices include a gate oxide layer having a variable thickness. The thickness of the gate oxide layer under the edge of the gate and over the source/emitter regions is different than the thickness over the channel regions of the device. The oxide layer thickness near the edge of the gate can be greater than the oxide layer thickness over the channel regions. The source/emitter regions can be implanted to provide enhanced oxide growth during gate oxide formation. The source/emitter region can include regions that are implanted to provide enhanced oxide growth during gate oxide formation and regions which do not provide enhanced oxide growth during gate oxide formation. The devices can be SiC devices such as SiC MOSFETs and SiC IGBTs.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Investigation of onshore wind farm wake recovery with in situ aircraft measurements during AWAKEN

The share of wind power for electricity supply is increasing worldwide. This highly variable resource requires the improved prediction of power output for network stability. The interaction between wind farm wakes and the atmospheric boundary layer (ABL) introduces uncertainties in power production that warrant detailed investigation. The flow downwind of wind farms is characterized by a reduction in wind speed and an increase in turbulence, which both vary with atmospheric conditions. During the American WAKE experimeNt (AWAKEN), the Technische Universität Braunschweig conducted measurement flights with a research aircraft upwind and downwind of onshore wind farms in the southern Great Plains in Oklahoma in the USA. This study utilizes data from 20 flights conducted at approximately hub height in September 2023 to investigate the wind field variability downwind of the wind farms and vertical profiles to observe atmospheric stratification. The flights were aligned perpendicular to the main wind direction downwind of the King Plains and Armadillo Flats wind farms. Additionally, lidar data from both upwind and downwind ground-based measurement sites and sonic anemometer data were used for comprehensive analysis. Results indicate that under stable ABL conditions, the wake persists at greater downwind distances with a higher velocity deficit in the wake relative to the undisturbed flow compared to unstable stratification. In homogeneous terrain under stable conditions, wake recovery to 95 % occurs between a distance of 4.5 and 9 km downwind of the wind farm. In the semi-complex terrain characterized by shallow hills, slopes, and valleys, the wake exhibits a higher velocity deficit compared to homogeneous terrain, while in some cases the wake was amplified by the terrain resulting in higher velocity deficit 10 km downwind of the wind farm compared to the measurements closer to the wind farm. The turbulent kinetic energy (TKE) and “TKE difference” was found to be a valuable measure in understanding wakes in a semi-complex terrain, showing a clear wake recovery and formation depending on the stratification of the ABL.

17 WIND ENERGY↗
DOE OSTI2021

Decomposition pathways and mitigation strategies for highly-stable hydroxyphenazine flow battery anolytes

Aqueous organic redox flow batteries are a promising technology for large-scale energy storage. The stability of the redox active organic molecules is increasingly being recognized as one of the major hurdles. Upon extended flow battery cycling, 7,8-dihydroxyphenazine-2 sulfonic acid (DHPS) undergoes desulfonation and reduction of a phenolic C–O bond to yield a mixture of 7/8-hydroxyphenazine-2-sulfonic acid, as well as hydrogenation of the aromatic ring system. Density functional theory (DFT) analysis of the charged DHPS, its ring-hydrogenated products, and variably substituted hydroxy phenazines has led to the development of a series of dihydroxylated phenazine isomers which provide insight into the effects of substitution pattern on solubility and stability. Seven dihydroxyphenazine (DHP) isomers were synthesized and their solubilities, electrochemical properties, flow battery cycling performance, and degradation pathways have been investigated. Based on theoretical and experimental results, hydroxyl substitution at the 1, 4, 6 and 9 positions yields highly stable derivatives, while substitution at the 2, 3, 7, and 8 positions results in unstable derivatives. Flow cells of 1,4- and 1,6-DHPs coupled with ferro/ferricyanide achieved high stabilities, with temporal capacity loss of 0.029 and 0.031% per day, respectively. Decomposition of 1,8- and 2,7-DHPs were found to arise from irreversible hydrogen rearrangement (tautomerization), yielding redox-inactive species. These results provide a detailed understanding of decomposition pathways and mitigation strategies for phenazine-based anolytes, and can provide general design guidelines for the development of stable redox-active organics.

25 ENERGY STORAGE↗
DOE OSTI2022

Respiration regimes in rivers: Partitioning source–specific respiration from metabolism time series

Respiration in streams is controlled by the timing, magnitude, and quality of organic matter (OM) inputs from internal primary production and external fluxes. Here, we estimated the contribution of different OM sources to seasonal, annual, and event-driven characteristics of whole-stream ecosystem respiration (ER) using an inverse modeling framework that accounts for possible time-lags between OM inputs and respiration. We modeled site-specific, dynamic OM stocks contributing to ER: autochthonous OM from gross primary production (GPP); allochthonous OM delivered during flow events; and seasonal pulses of leaf litter. OM stored in the sediment and dissolved organic matter (DOM) transported during baseflow were modeled as a stable stock contributing to baseline respiration. We applied this modeling framework to five streams with different catchment size, climate, and canopy cover, where multi-year time series of ER and environmental variables were available. Overall, the model explained between 53% and 74% of observed ER dynamics. Respiration of autochthonous OM tracked seasonal peaks in GPP in spring or summer. Increases in ER were often associated with high-flow events. Respiration associated with litter inputs was larger in smaller streams. Time lags between leaf inputs and respiration were longer than for other OM sources, likely due to lower biological reactivity. Model estimates of source-specific ER and OM stocks compared well with existing measures of OM stocks, inputs, and respiration or decomposition. Our modeling approach has the potential to expand the scale of comparative analyses of OM dynamics within and among freshwater ecosystems.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Sensitivity-based voltage constraints for optimal power flow in low-voltage distribution feeders

The optimal power flow (OPF) problem for distribution systems can include network details down to the low-voltage (LV) points of interconnection of individual customers. This paper addresses the implementation of voltage magnitude constraints, and sets forth a practicable approach for capturing the effects on voltage from the switching behavior of loads (e.g., heat pumps, air conditioners, water heaters, or pool pumps) and from the variability of renewable generation (e.g., rooftop solar). The proposed method adjusts the OPF voltage constraints based on forecasts of load and generation upper and lower bounds, in conjunction with sensitivity factors derived from the power flow equations. An illustrative OPF formulation is also provided, which incorporates transformer models that include core loss. We demonstrate that accurate modeling of these LV network components is critical to avoid voltage violations at customer points of interconnection. Furthermore, the ideas are validated through numerical case studies on a realistic distribution feeder.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2021

Two Dimensional Topology Optimization of Heat Exchangers with the Density and Level-Set Methods

We design heat exchangers using two topology optimization approaches: the density, i.e. volume fraction and level set methods. Our goal is to maximize the heat exchange between two fluids in separate channels while constraining the pressure drop across each channel. The heat exchanger is modeled with a coupled thermal-flow formulation. The flow is governed by an isothermal and incompressible Stokes-Brinkman equation and the heat transfer is governed by a convection-diffusion equation with high Peclet number. We solve one set of Stokes-Brinkman equations per fluid. Each Brinkman term in the flow equation serves to model the other phase as a solid, thereby preventing mixing. We first represent the solid and fluid phases using a volume fraction variable and apply a SIMP-like penalization in the Brinkman term to drive the optimization to a discrete design. The cost and constraint function derivatives are automatically calculated with the library pyadjoint and the optimization is performed by the Method of Moving Asymptotes. In a second optimization formulation, we use the level set approach to define the interface that separates the two fluids. Pyadjoint calculates the shape derivatives of the cost and constraint functions and the Hamilton-Jacobi advects the interface, allowing for topological changes. We present results in two dimensions and discuss the advantages and disadvantages of each approach.

42 ENGINEERING↗
DOE OSTI2020

Chemo‐mechanical Alterations Induced From CO 2 Injection in Carbonate‐Cemented Sandstone: An Experimental Study at 71 °C and 29 MPa

Abstract Carbon capture, utilization, and storage may lead to mechanical degradation of the subsurface reservoir from fluid‐rock interaction, which could lead to wellbore instability or reservoir compaction. To better understand potential relationship between mechanical degradation with various carbonate cement textures and compositions in sandstone reservoirs, six flow‐through experiments were conducted. Formation water (TDS = 5,390 mg/L) enriched with CO 2 flowed through two types of Pennsylvanian Morrow B Sandstone: an ankerite‐siderite‐cemented sandstone (disseminated cement texture) and a calcite‐cemented sandstone (poikilotopic cement texture). The experiments produced little change in permeability in the ankerite‐siderite‐cemented sandstone, but permeability increased up to more than 1 order of magnitude in the calcite‐cemented sandstone. Ultrasonic measurements and cylinder‐splitting tests (also known as Brazilian tests) suggested negligible mechanical degradation of the ankerite‐siderite‐cemented sandstone. Variable changes, with significant mechanical degradation in the static moduli, were observed in the calcite‐cemented sandstone. Thus, dissolution of the disseminated ankerite‐siderite cement (0.28–0.30%) had minimal impact on modifying the flow network and the mechanical integrity of the sandstone, whereas dissolution of the poikilotopic calcite cement (0.89–1.13%, quantified with fluid chemistry and visualized with X‐ray microcomputed tomography) impacted the mechanical strength of the sandstone by disconnecting framework grains. With the high water‐to‐rock mass ratios (7.3–8.2) and number of pore volumes (147–675) employed in these experiments, potential risks are most relevant to regions near injection wells. Ultimately, the chemo‐mechanical effects induced by CO 2 injection are strongly influenced by the cement texture and composition and the burial history of the reservoir rock.

Wu, Z.↗
DOE OSTI2022

Coordinated operation of pumped-storage hydropower with power and water distribution systems

Small pumped-storage hydropower (PSH) units have gained popularity as distributed energy storage options that can provide flexibility to the operation of power distribution systems. Optimal operation of small PSH units is not only dependent on the energy storage provided to power distribution system, but also on the inflow and outflow of water from and to the water distribution system. Here, in this context, this paper develops an optimization model for coordinated operation of PSH units with power and water distribution systems. The proposed model optimizes the operation of water tanks, variable-speed pumps and PSH in pumping and generating modes to minimize the operation cost of power distribution system, while respecting the power flow constraints of power distribution and hydraulic constraints of water distribution system. Appropriate electricity tariffs are implemented to avoid additional expenses in water distribution system that can be enforced by its coordinated operation in favor of power distribution system. The proposed model is implemented on a 33-bus and a 123-bus test power distribution system connected to a 16-node test water distribution system. Results demonstrate the effectiveness of proposed model in tapping PSH flexibility to reduce the operation cost of power and water distribution systems, while meeting the power and water demands.

13 HYDRO ENERGY↗
DOE OSTI2023

Fabric controls on fracture surface roughness of an architected rock material

Fluid flow through fractures is intimately linked to the fracture surfaces that define the void geometry through which fluids flow. Thus, an understanding of what controls fracture surface roughness is essential to the development of models for predicting fluid transport through fractured rock. The difficulty in predicting surface roughness arises from the complexity of rock which is inherently heterogeneous and nonuniform in composition, fabric, and structural components, even when samples are acquired from the same rock mass. Here, a benchmarked-simulation approach motivated from geo-architected 3D printed synthetic gypsum rocks is used to provide insight into the competing contributions from fabric and layering on fracture roughness formation. Simulation results from a discrete element model (Particle Flow Code, Itasca Consulting Group, Inc.) clearly indicate that the relative orientation between mineral layers and in-layer mineral fabric, and the variability in mineral bonding strengths determine whether anisotropic corrugated surfaces or isotropic surfaces are formed. Weak mineral layers oriented perpendicular to the applied load resulted in strong roughness anisotropy. Peak failure loads were found to vary up to 30% depending on the strength of the mineral fabric at the location of fracture initiation, which provides insight into the observed high variability in strength values of natural rock. The uniqueness of induced fracture roughness and peak failure load is intimately linked to layering, mineral fabric, and their distribution in the rock. These findings have important implications for any architected material fabricated through serial printing of layers with local compositional heterogeneity.

3D printed rock↗
DOE OSTI2022

Global Barotropic Tide Modeling Using Inline Self‐Attraction and Loading in MPAS‐Ocean

Abstract We examine ocean tides in the barotropic version of the Model for Prediction Across Scales (MPAS‐Ocean), the ocean component of the Department of Energy Earth system model. We focus on four factors that affect tidal accuracy: self‐attraction and loading (SAL), model resolution, details of the underlying bathymetry, and parameterized topographic wave drag. The SAL term accounts for the tidal loading of Earth's crust and the self‐gravitation of the ocean and the load‐deformed Earth. A common method for calculating SAL is to decompose mass anomalies into their spherical harmonic constituents. Here, we compare a scalar SAL approximation versus an inline SAL using a fast spherical harmonic transform package. Wave drag accounts for energy lost by breaking internal tides that are produced by barotropic tidal flow over topographic features. We compare a series of successively finer quasi‐uniform resolution meshes (62.9, 31.5, 15.7, and 7.87 km) to a variable resolution (45 to 5 km) configuration. We ran MPAS‐Ocean in a single‐layer barotropic mode forced by five tidal constituents. The 45 to 5 km variable resolution mesh obtained the best total root‐mean‐square error (5.4 cm) for the deep ocean ( 1,000 m) tide compared to TPXO8 and ran twice as fast as the quasi‐uniform 8 km mesh, which had an error of 5.8 cm. This error is comparable to those found in other forward (non‐assimilative) ocean tide models. In future work, we plan to use MPAS‐Ocean to study tidal interactions with other Earth system components, and the tidal response to climate change.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Distribution of centrality measures on undirected random networks via the cavity method

The Katz centrality of a node in a complex network is a measure of the node’s importance as far as the flow of information across the network is concerned. For ensembles of locally tree-like undirected random graphs, this observable is a random variable. Its full probability distribution is of interest but difficult to handle analytically because of its “global” character and its definition in terms of a matrix inverse. Leveraging a fast Gaussian Belief Propagation-Cavity algorithm to solve linear systems on tree-like structures, we show that i) the Katz centrality of a single instance can be computed recursively in a very fast way, and ii) the probability P ( K ) that a random node in the ensemble of undirected random graphs has centrality K satisfies a set of recursive distributional equations, which can be analytically characterized and efficiently solved using a population dynamics algorithm. We test our solution on ensembles of Erdős-Rényi and Scale Free networks in the locally tree-like regime, with excellent agreement. The analytical distribution of centrality for the configuration model conditioned on the degree of each node can be employed as a benchmark to identify nodes of empirical networks with over- and underexpressed centrality relative to a null baseline. We also provide an approximate formula based on a rank- 1 projection that works well if the network is not too sparse, and we argue that an extension of our method could be efficiently extended to tackle analytical distributions of other centrality measures such as PageRank for directed networks in a transparent and user-friendly way.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Modeling of Cube Array Roughness: RANS, Large Eddy Simulation, and Direct Numerical Simulation

Abstract Flow over arrays of cubes is an extensively studied model problem for rough wall turbulent boundary layers. While considerable research has been performed in computationally investigating these topologies using direct numerical simulation (DNS) and large eddy simulation (LES), the ability of sublayer-resolved Reynolds-averaged Navier–Stokes (RANS) to predict the bulk flow phenomena of these systems is relatively unexplored, especially at low and high packing densities. Here, RANS simulations are conducted on six different packing densities of cubes in aligned and staggered configurations. The packing densities investigated span from what would classically be defined as isolated, up to those in the d-type roughness regime, filling in the gap in the present literature. Three different sublayer-resolved turbulence closure models were tested for each case: a low Reynolds number k–ϵ model, the Menter k–ω SST model, and a full Reynolds stress model. Comparisons of the velocity fields, secondary flow features, and drag coefficients are made between the RANS results and existing LES and DNS results. There is a significant degree of variability in the performance of the various RANS models across all comparison metrics. However, the Reynolds stress model demonstrated the best accuracy in terms of the mean velocity profile as well as drag partition across the range of packing densities.

Engineering↗
DOE OSTI2023

Utah FORGE: Fault Shear Reactivation Experimental Data for Fluid Injection-Rate Controls on Seismic Moment

Included are experimental data recorded from shear experiments that specifically explore the link between fluid-injection rate and seismic moment resulting from shear reactivation of laboratory faults. Raw mechanical data from three experiments are included alongside corresponding MATLAB scripts that import and plot the data, as well as use it to calculate shear and normal stress. Experiments are performed on 2.5-3 inch long granitoid cores from the Utah FORGE EGS demonstration site, containing a single inclined fracture with small-scale roughness added to the fracture surface. The raw data included here were recorded from an aluminum triaxial pressure vessel (TEMCO) configured with three independent servo-controlled pumps, with distilled water used as the working fluid. The pumps control confining pressure, upstream pore pressure, and axial pressure, with each pump connected to a LabView interface to record applied pressures, cumulative injected water volumes, and pump flow rates. The downstream outlet from the fracture is closed to allow pressurization, which is measured by an external pressure transducer. A linear variable differential transformer (LVDT) attached to the axial piston measures axial displacement, from which we calculate shear displacement along the fracture. Additionally, P-wave transducers are used to record acoustic signatures, where acoustic emission events and maximum amplitudes are compared against seismic moment and shear slip velocity. Fluid injection rates range between 0.05 mL/min, 0.25 mL/min, and 0.75 mL/min for each experiment. Along-fault pressure distributions are progressively less uniform as injection rates increase, representing a switch from steady-state to transient conditions. Triggered shear displacement is used as a proxy for seismic moment and is indexed against cumulative injection volume and rate. Each experiment is performed under constant shear stress conditions, and the sample is fully saturated with DI water. Axial and confining stresses are applied to 3 MPa through pressure-stepping in 500 kPa increments. The pore pressure is held constant at 200 kPa prior to initiating the experiment, and initial axial displacement is recorded. The axial stress is then increased to initiate shear mobilization during the loading phase (run-in) until a peak steady state is achieved. The initial shear stress is reduced to approximately 80% of the peak shear stress by decreasing the axial stress, then held constant for the duration of each experiment.

15 GEOTHERMAL ENERGY↗
DOE OSTI2022

Technoeconomic Analysis of Infrastructure Buildout Scenarios

The success of CCUS deployment in the southeast will depend on the optimized pipeline network from CO 2 point sources to geologic sinks for CO 2 storage. An optimal pipeline framework will reduce both environmental impacts and pipeline building cost. However, finding an optimal pipeline/transport infrastructure design is a non-trivial task and requires simultaneous usage of large volumes of data, computational resources, and a state-of-the-art simulator. Such a transport infrastructure model must consider point sources for CO 2 capture and associated volumes, sinks for CO 2 storage, and transportation from source to sink via pipeline networks. These considerations must be addressed simultaneously and systemically in an optimization process, in which a defined objective function (i.e., capital, variable CO 2 capture, transport, and storage cost function) is required to be minimized with the consideration of practical constraints (e.g., CO 2 flow through pipelines is less than the maximum capacity). Although optimization has been widely used in subsurface resources production and CO 2 storage, its application in large scale CCUS infrastructure design is rarely reported in the literature. SimCCS, developed by Los Alamos National Laboratory, is an open-source CCUS pipeline infrastructure design toolset that facilitates the optimization of pipeline infrastructure networks.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2020

Localized basal area affects soil respiration temperature sensitivity in a coastal deciduous forest

Abstract. Soil respiration (Rs), the flow of CO2 from the soil surface to the atmosphere, is one of the largest carbon fluxes in the terrestrial biosphere. The spatial variability of Rs is both large and poorly understood, limiting our ability to robustly scale it in space. One factor in Rs spatial variability is the autotrophic contribution from plant roots, but it is uncertain how the presence of plants affects the magnitude and temperature sensitivity of Rs. This study used 1 year of Rs measurements to examine the effect of localized basal area on Rs in the growing and dormant seasons, as well as during moisture-limited times, in a temperate, coastal, deciduous forest in eastern Maryland, USA. In a linear mixed-effects model, tree basal area within a 5 m radius (BA5) exerted a significant positive effect on the temperature sensitivity of soil respiration. Soil moisture was the dominant control on Rs during the dry portions of the year, while soil moisture, temperature, and BA5 all exerted significant effects on Rs in wetter periods. Our results suggest that autotrophic respiration is more sensitive to temperature than heterotrophic respiration at these sites, although we did not measure these source fluxes directly, and that soil respiration is highly moisture sensitive, even in a record-rainfall year. The Rs flux magnitudes (0.46–15.0 µmol m-2 s-1) and variability (coefficient of variability 10 %–23 % across plots) observed in this study were comparable to values observed in similar forests. Six Rs observations would be required in order to estimate the mean across all study sites to within 50 %, and 518 would be required in order to estimate it to within 5 %, with 95 % confidence. A better understanding of the spatial interactions between plants and microbes, as well as the strength and speed of above- and belowground coupling, is necessary to link these processes with large-scale soil-to-atmosphere C fluxes.

54 ENVIRONMENTAL SCIENCES↗