Stoy Lab
Permanent URI for this collectionhttps://scholarworks.montana.edu/handle/1/14931
In the Stoy Lab, we study the role of vegetation in the climate system. To do so we measure and model the exchange of water, heat, and trace gases like carbon dioxide and methane between the terrestrial surface and the atmosphere. Recent efforts seek to understand feedbacks between land management and precipitation processes.
Browse
17 results
Search Results
Item The Effects of Elevated Atmospheric CO2 and Nitrogen Amendments on Subsurface CO2 Production and Concentration Dynamics in a Maturing Pine Forest(2009-05) Daly, Edoardo; Palmroth, Sari; Stoy, Paul C.; Siqueira, Mario B. S.; Oishi, A. Christopher; Juang, Jehn-Yih; Oren, Ram; Porporato, Amilcare; Katul, Gabriel G.Profiles of subsurface soil CO2 concentration, soil temperature, and soil moisture, and throughfall were measured continuously during the years 2005 and 2006 in 16 locations at the free air CO2 enrichment facility situated within a temperate loblolly pine (Pinus taeda L.) stand. Sampling at these locations followed a 4 by 4 replicated experimental design comprised of two atmospheric CO2 concentration levels (ambient [CO2]a, ambient + 200 ppmv, [CO2]e) and two soil nitrogen (N) deposition levels (ambient, ambient + fertilization at 11.2 gN m−2 year−1). The combination of these measurements permitted indirect estimation of below ground CO2 production and flux profiles in the mineral soil. Adjacent to the soil CO2 profiles, direct (chamber-based) measurements of CO2 fluxes from the soil–litter complex were simultaneously conducted using the automated carbon efflux system. Based on the measured soil CO2 profiles, neither [CO2]e nor N fertilization had a statistically significant effect on seasonal soil CO2, CO2 production, and effluxes from the mineral soil over the study period. Soil moisture and temperature had different effects on CO2 concentration depending on the depth. Variations in CO2 were mostly explained by soil temperature at deeper soil layers, while water content was an important driver at the surface (within the first 10 cm), where CO2 pulses were induced by rainfall events. The soil effluxes were equal to the CO2 production for most of the time, suggesting that the site reached near steady-state conditions. The fluxes estimated from the CO2 profiles were highly correlated to the direct measurements when the soil was neither very dry nor very wet. This suggests that a better parameterization of the soil CO2 diffusivity is required for these soil moisture extremes.Item The relationship between reference canopy conductance and simplified hydraulic architecture(2009-06) Novick, Kimberly A.; Oren, Ram; Stoy, Paul C.; Juang, Jehn-Yih; Siqueira, Mario B. S.; Katul, Gabriel G.Terrestrial ecosystems are dominated by vascular plants that form a mosaic of hydraulic conduits to water movement from the soil to the atmosphere. Together with canopy leaf area, canopy stomatal conductance regulates plant water use and thereby photosynthesis and growth. Although stomatal conductance is coordinated with plant hydraulic conductance, governing relationships across species has not yet been formulated at a practical level that can be employed in large-scale models. Here, combinations of published conductance measurements obtained with several methodologies across boreal to tropical climates were used to explore relationships between canopy conductance rates and hydraulic constraints. A parsimonious hydraulic model requiring sapwood-to-leaf area ratio and canopy height generated acceptable agreement with measurements across a range of biomes (r2 = 0.75) . The results suggest that, at long time scales, the functional convergence among ecosystems in the relationship between water-use and hydraulic architecture eclipses inter-specific variation in physiology and anatomy of the transport system. Prognostic applicability of this model requires independent knowledge of sapwood-to-leaf area. In this study, we did not find a strong relationship between sapwood-to-leaf area and physical or climatic variables that are readily determinable at coarse scales, though the results suggest that climate may have a mediating influence on the relationship between sapwood-to-leaf area and height. Within temperate forests, canopy height alone explained a large amount of the variance in reference canopy conductance (r2 = 0.68) and this relationship may be more immediately applicable in the terrestrial ecosystem models.Item Investigating a Hierarchy of Eulerian Closure Models for Scalar Transfer Inside Forested Canopies(2008-04) Juang, Jehn-Yih; Katul, Gabriel G.; Siqueira, Mario B. S.; Stoy, Paul C.; McCarthy, Heather R.Modelling the transfer of heat, water vapour, and CO2 between the biosphere and the atmosphere is made difficult by the complex two-way interaction between leaves and their immediate microclimate. When simulating scalar sources and sinks inside canopies on seasonal, inter-annual, or forest development time scales, the so-called well-mixed assumption (WMA) of mean concentration (i.e. vertically constant inside the canopy but dynamically evolving in time) is often employed. The WMA eliminates the need to model how vegetation alters its immediate microclimate, which necessitates formulations that utilize turbulent transport theories. Here, two inter-related questions pertinent to the WMA for modelling scalar sources, sinks, and fluxes at seasonal to inter-annual time scales are explored: (1) if the WMA is to be replaced so as to resolve this two-way interaction, how detailed must the turbulent transport model be? And (2) what are the added predictive skills gained by resolving the two-way interaction vis-à-vis other uncertainties such as seasonal variations in physiological parameters. These two questions are addressed by simulating multi-year mean scalar concentration and eddy-covariance scalar flux measurements collected in a Loblolly pine (P. taeda L.) plantation near Durham, North Carolina, U.S.A. using turbulent transport models ranging from K-theory (or first-order closure) to third-order closure schemes. The multi-layer model calculations with these closure schemes were contrasted with model calculations employing the WMA. These comparisons suggested that (i) among the three scalars, sensible heat flux predictions are most biased with respect to eddy-covariance measurements when using the WMA, (ii) first-order closure schemes are sufficient to reproduce the seasonal to inter-annual variations in scalar fluxes provided the canonical length scale of turbulence is properly specified, (iii) second-order closure models best agree with measured mean scalar concentration (and temperature) profiles inside the canopy as well as scalar fluxes above the canopy, (iv) there are no clear gains in predictive skills when using third-order closure schemes over their second-order closure counterparts. At inter-annual time scales, biases in modelled scalar fluxes incurred by using the WMA exceed those incurred when correcting for the seasonal amplitude in the maximum carboxylation capacity (V cmax, 25) provided its mean value is unbiased. The role of local thermal stratification inside the canopy and possible computational simplifications in decoupling scalar transfer from the generation of the flow statistics are also discussed. “The tree, tilting its leaves to capture bullets of light; inhaling, exhaling; its many thousand stomata breathing, creating the air”. Ruth Stone, 2002, In the Next Galaxy"Item Role of vegetation in determining carbon sequestration along ecological succession in the southeastern United States(2008-06) Stoy, Paul C.; Katul, Gabriel G.; Siqueira, Mario B. S.; Juang, Jehn-Yih; Novick, Kimberly A.; McCarthy, Heather R.; Oishi, A. Christopher; Oren, RamVegetation plays a central role in controlling terrestrial carbon (C) exchange, but quantifying its impacts on C cycling on time scales of ecological succession is hindered by a lack of long‐term observations. The net ecosystem exchange of carbon (NEE) was measured for several years in adjacent ecosystems that represent distinct phases of ecological succession in the southeastern USA. The experiment was designed to isolate the role of vegetation – apart from climate and soils – in controlling biosphere–atmosphere fluxes of CO2 and water vapor. NEE was near zero over 5 years at an early successional old‐field ecosystem (OF). However, mean annual NEE was nearly equal, approximately −450 g C m−2 yr−1, at an early successional planted pine forest (PP) and a late successional hardwood forest (HW) due to the sensitivity of the former to drought and ice storm damage. We hypothesize that these observations can be explained by the relationships between gross ecosystem productivity (GEP), ecosystem respiration (RE) and canopy conductance, and long‐term shifts in ecosystem physiology in response to climate to maintain near‐constant ecosystem‐level water‐use efficiency (EWUE). Data support our hypotheses, but future research should examine if GEP and RE are causally related or merely controlled by similar drivers. At successional time scales, GEP and RE observations generally followed predictions from E. P. Odum's ‘Strategy of Ecosystem Development’, with the surprising exception that the relationship between GEP and RE resulted in large NEE at the late successional HW. A practical consequence of this research suggests that plantation forestry may confer no net benefit over the conservation of mature forests for C sequestration.Item Eco-hydrological controls on summertime convective rainfall triggers(2007-01) Juang, Jehn-Yih; Katul, Gabriel G.; Porporato, Amilcare; Stoy, Paul C.; Siqueira, Mario B. S.; Detto, Matteo; Kim, Hyun-Seok; Oren, RamTriggers of summertime convective rainfall depend on numerous interactions and feedbacks, often compounded by spatial variability in soil moisture and its impacts on vegetation function, vegetation composition, terrain, and all the complex turbulent entrainment processes near the capping inversion. To progress even within the most restricted and idealized framework, many of the governing processes must be simplified and parameterized. In this work, a zeroth‐order representation of the dynamical processes that control convective rainfall triggers – namely land surface fluxes of heat and moisture – is proposed and used to develop a semianalytical model to explore how differential sensitivities of various ecosystems to soil moisture states modify convective rainfall triggers. The model is then applied to 4 years (2001–2004) of half‐hourly precipitation, soil moisture, environmental, and eddy‐covariance surface heat flux data collected at a mixed hardwood forest (HW), a maturing planted loblolly pine forest (PP), and an abandoned old field (OF) experiencing the same climatic and edaphic conditions. We found that the sensitivity of PP to soil moisture deficit enhances the trigger of convective rainfall relative to HW and OF, with enhancements of about 25% and 30% for dry moisture states, and 5% and 15% for moist soil moisture states, respectively. We discuss the broader implications of these findings on potential modulations of convective rainfall triggers induced by projected large‐scale changes in timberland composition within the Southeastern United States.Item An evaluation of methods for partitioning eddy covariance-measured net ecosystem exchange into photosynthesis and respiration(2006-12) Stoy, Paul C.; Katul, Gabriel G.; Siqueira, Mario B. S.; Juang, Jehn-Yih; Novick, Kimberly A.; Uebelherr, Joshua M.; Oren, RamWe measured net ecosystem CO2 exchange (NEE) using the eddy covariance (EC) technique for 4 years at adjoining old field (OF), planted pine (PP) and hardwood forest (HW) ecosystems in the Duke Forest, NC. To compute annual sums of NEE and its components – gross ecosystem productivity (GEP) and ecosystem respiration (RE) – different ‘flux partitioning’ models (FPMs) were tested and the resulting C flux estimates were compared against published estimates from C budgeting approaches, inverse models, physiology-based forward models, chamber respiration measurements, and constraints on assimilation based on sapflux and evapotranspiration measurements. Our analyses demonstrate that the more complex FPMs, particularly the ‘non-rectangular hyperbolic method’, consistently produced the most reasonable C flux estimates. Of the FPMs that use nighttime data to estimate RE, one that parameterized an exponential model over short time periods generated predictions that were closer to expected flux values. To explore how much ‘new information’ was injected into the data by the FPMs, we used formal information theory methods and computed the Shannon entropy for: (1) the probability density, to assess alterations to the flux measurement distributions, and (2) the wavelet energy spectra, to assess alterations to the internal autocorrelation within the NEE time series. Based on this joint analysis, gap-filling had little impact on the IC of daytime data, but gap-filling significantly altered nighttime data in both the probability and wavelet spectral domains.Item Are ecosystem carbon inputs and outputs coupled at short time scales? A case study from adjacent pine and hardwood forests using impulse-response analysis(2007-06) Stoy, Paul C.; Palmroth, Sari; Oishi, A. Christopher; Siqueira, Mario B. S.; Juang, Jehn-Yih; Novick, Kimberly A.; Ward, Eric J.; Katul, Gabriel G.; Oren, RamA number of recent studies have attributed a large proportion of soil respiration (Rsoil) to recently photoassimilated carbon (C). Time lags (τPR) associated with these pulses of photosynthesis and responses of Rsoil have been found on time scales of hours to weeks for different ecosystems, but most studies find evidence for τPR on the order of 1–5 d. We showed that such time scales are commensurate with CO2 diffusion time scales from the roots to the soil surface, and may thus be independent from photosynthetic pulses. To further quantify the role of physical (i.e. edaphic) and biological (i.e. vegetative) controls on such lags, we investigated τPR at adjacent planted pine (PP) and hardwood (HW) forest ecosystems over six and four measurement years, respectively, using both autocorrelation analysis on automated soil surface flux measurements and their lagged cross‐correlations with drivers for and surrogates of photosynthesis. Evidence for τPR on the order of 1–3 d was identified in both ecosystems and using both analyses, but this lag could not be attributed to recently photoassimilated C because the same analysis yielded comparable lags at HW during leaf‐off periods. Future efforts to model ecosystem C inputs and outputs in a pulse–response framework must combine measurements of transport in the physical and biological components of terrestrial ecosystems.Item Carbon dioxide and water vapor exchange in a warm temperate grassland(2004-01) Novick, Kimberly A.; Stoy, Paul C.; Katul, Gabriel G.; Ellsworth, D. S.; Siqueira, Mario B. S.; Juang, Jehn-Yih; Oren, RamGrasslands cover about 40% of the ice-free global terrestrial surface, but their contribution to local and regional water and carbon fluxes and sensitivity to climatic perturbations such as drought remains uncertain. Here, we assess the direction and magnitude of net ecosystem carbon exchange (NEE) and its components, ecosystem carbon assimilation (A c) and ecosystem respiration (R E), in a southeastern United States grassland ecosystem subject to periodic drought and harvest using a combination of eddy-covariance measurements and model calculations. We modeled A c and evapotranspiration (ET) using a big-leaf canopy scheme in conjunction with ecophysiological and radiative transfer principles, and applied the model to assess the sensitivity of NEE and ET to soil moisture dynamics and rapid excursions in leaf area index (LAI) following grass harvesting. Model results closely match eddy-covariance flux estimations on daily, and longer, time steps. Both model calculations and eddy-covariance estimates suggest that the grassland became a net source of carbon to the atmosphere immediately following the harvest, but a rapid recovery in LAI maintained a marginal carbon sink during summer. However, when integrated over the year, this grassland ecosystem was a net C source (97 g C m−2 a−1) due to a minor imbalance between large A c (−1,202 g C m−2 a−1) and R E (1,299 g C m−2 a−1) fluxes. Mild drought conditions during the measurement period resulted in many instances of low soil moisture (θ<0.2 m3m−3), which influenced A c and thereby NEE by decreasing stomatal conductance. For this experiment, low θ had minor impact on R E. Thus, stomatal limitations to A c were the primary reason that this grassland was a net C source. In the absence of soil moisture limitations, model calculations suggest a net C sink of −65 g C m−2 a−1 assuming the LAI dynamics and physiological properties are unaltered. These results, and the results of other studies, suggest that perturbations to the hydrologic cycle are key determinants of C cycling in grassland ecosystems.Item Variability in net ecosystem exchange from hourly to inter-annual time scales at adjacent pine and hardwood forests: a wavelet analysis(2005-07) Stoy, Paul C.; Katul, Gabriel G.; Siqueira, Mario B. S.; Juang, Jehn-Yih; McCarthy, Heather R.; Kim, Hyun-Seok; Oishi, A. Christopher; Oren, RamOrthonormal wavelet transformation (OWT) is a computationally efficient technique for quantifying underlying frequencies in nonstationary and gap-infested time series, such as eddy-covariance-measured net ecosystem exchange of CO2 (NEE). We employed OWT to analyze the frequency characteristics of synchronously measured and modeled NEE at adjacent pine (PP) and hardwood (HW) ecosystems. Wavelet cospectral analysis showed that NEE at PP was more correlated to light and vapor pressure deficit at the daily time scale, and NEE at HW was more correlated to leaf area index (LAI) and temperature, especially soil temperature, at seasonal time scales. Models were required to disentangle the impacts of environmental drivers on the components of NEE, ecosystem carbon assimilation (Ac) and ecosystem respiration (RE). Sensitivity analyses revealed that using air temperature rather than soil temperature in RE models improved the modeled wavelet spectral frequency response on time scales longer than 1 day at both ecosystems. Including LAI improved RE model fit on seasonal time scales at HW, and incorporating parameter variability improved the RE model response at annual time scales at both ecosystems. Resolving variability in canopy conductance, rather than leaf-internal CO2, was more important for modeling Ac at both ecosystems. The PP ecosystem was more sensitive to hydrologic variables that regulate canopy conductance: vapor pressure deficit on weekly time scales and soil moisture on seasonal to interannual time scales. The HW ecosystem was sensitive to water limitation on weekly time scales. A combination of intrinsic drought sensitivity and non-conservative water use at PP was the basis for this response. At both ecosystems, incorporating variability in LAI was required for an accurate spectral representation of modeled NEE. However, nonlinearities imposed by canopy light attenuation were of little importance to spectral fit. The OWT revealed similarities and differences in the scale-wise control of NEE by vegetation with implications for model simplification and improvement.Item Multiscale model intercomparisons of CO2 and H2O exchange in a maturing southeastern U.S. pine forest(2006-07) Siqueira, Mario B. S.; Katul, Gabriel G.; Sampson, D. A.; Stoy, Paul C.; Juang, Jehn-Yih; McCarthy, Heather R.; Oren, RamWe compared four existing process‐based stand‐level models of varying complexity (physiological principles in predicting growth, photosynthesis and evapotranspiration, biogeochemical cycles, and stand to ecosystem carbon and evapotranspiration simulator) and a new nested model with 4 years of eddy‐covariance‐measured water vapor (LE) and CO2 (Fc) fluxes at a maturing loblolly pine forest. The nested model resolves the ‘fast’ CO2and H2O exchange processes using canopy turbulence theories and radiative transfer principles whereas slowly evolving processes were resolved using standard carbon allocation methods modified to improve leaf phenology. This model captured most of the intraannual variations in leaf area index (LAI), net ecosystem exchange (NEE), and LE for this stand in which maximum LAI was not at a steady state. The model comparisons suggest strong linkages between carbon production and LAI variability, especially at seasonal time scales. This linkage necessitates the use of multilayer models to reproduce the seasonal dynamics of LAI, NEE, and LE. However, our findings suggest that increasing model complexity, often justified for resolving faster processes, does not necessarily translate into improved predictive skills at all time scales. Additionally, none of the models tested here adequately captured drought effects on water and CO2 fluxes. Furthermore, the good performance of some models in capturing flux variability on interannual time scales appears to stem from erroneous LAI dynamics and from sensitivity to droughts that injects unrealistic flux variability at longer time scales.