Search NASA⌕ Search

SEARCH · Search NASA

Results for “covariance”

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 523 records · Page 29
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Influence of Tundra Polygon Type and Climate Variability on CO 2 and CH 4 Fluxes Near Utqiagvik, Alaska

Arctic tundra has the potential to generate significant climate feedbacks, but spatial complexity makes it difficult to quantify the impacts of climate on ecosystem-atmosphere fluxes, particularly in polygonal tundra comprising wetter and drier polygon types on the scale of tens of meters. We measured CO 2 , CH 4 , and energy fluxes using eddy covariance for 7 yr (April to November, 2013–2019) in polygonal tundra near Utqiagvik, Alaska. This period saw the earliest snowmelt, latest snow accumulation, and hottest summer on record. To estimate fluxes by polygon type, we combined a polygon classification with a flux-footprint model. Methane fluxes were highest in the summer months but were also large during freeze-up and increased with the warming trend in August–November temperatures. While CO 2 respiration had a consistent, exponential relationship with temperature, net ecosystem exchange was more variable among years. CO 2 and CH 4 exchange (June–September) ranged between -0.83 (Standard error [SE] = 0.03) and -1.32 (SE = 0.04) μmol m -2 s-1 and 13.92 (SE = 0.26)—23.42 (SE = 0.45) nmol m -2 s -1 , respectively, and varied interannually (p ≤ 0.05). The maximum-influence method effectively attributed fluxes to polygon types. Areas dominated by low-centered polygons had higher CO 2 fluxes except in 2016–2017. Methane fluxes were highest in low-centered polygons 2013–2015 and in flat-centered polygons in subsequent years, possibly due to increasing temperature and precipitation. Sensible and latent heat fluxes also varied significantly among polygon types. Accurate characterization of Arctic fluxes and their climate dependencies requires spatial disaggregation and long term observations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Warming and Increased Respiration Have Transformed an Alpine Steppe Ecosystem on the Tibetan Plateau From a Carbon Dioxide Sink Into a Source

Cold region ecosystems store vast amounts of soil organic carbon (C), which upon warming and decomposition can affect the C balance and potentially change these ecosystems from C sinks to carbon dioxide (CO 2 ) sources. We quantified the decadal year-round CO 2 flux from an alpine steppe-ecosystem on the Tibetan Plateau using eddy covariance and automatic chamber approaches during a period of significant warming (0.13°C per 10 years; and 0.18°C in the non-growing season alone: 1st October to next 30th April). The results showed that ongoing climate change, mainly warming within the topsoil layers, is the main reason for the site’s change from a sink for to a source of CO 2 in the atmosphere. Non-growing-season ecosystem respiration accounted for 51% of the annual ecosystem respiration and has increased significantly. The growing seasons (1st May to 30th September) were consistent CO 2 sink periods without significant changes over the study period. Observations revealed high-emission events from the end of the non-growing season to early in the growing season (1st March to fifteenth May), which significantly (p < 0.01) increased at a rate of 22.6 g C m -2 decade -1 , ranging from 14.6 ± 10.7 g C m -2 yr -1 in 2012 to 35.3 ± 12.1 g C m -2 yr -1 in 2017. Further, structural equation modeling suggested that active layer warming was the key factor in explaining changes in ecosystem respiration, leading to significant changes in net ecosystem exchange over the period 2011–2020 and indicated that these changes have already transformed the ecosystem from a CO 2 sink into a source. These results can be used to improve our understanding of the sensitivity of ecosystem respiration to increased warming during the non-growing period.

58 GEOSCIENCES↗
DOE OSTI2022

Detecting Hot Spots of Methane Flux Using Footprint‐Weighted Flux Maps

Abstract In this study, we propose a new technique for mapping the spatial heterogeneity in gas exchange around flux towers using flux footprint modeling and focusing on detecting hot spots of methane (CH 4 ) flux. In the first part of the study, we used a CH 4 release experiment to evaluate three common flux footprint models: the Hsieh model (Hsieh et al., 2000), the Kljun model (Kljun et al., 2015), and the K & M model (Kormann and Meixner, 2001), finding that the K & M model was the most accurate under these conditions. In the second part of the study, we introduce the Footprint‐Weighted Flux Map, a new technique to map spatial heterogeneity in fluxes. Using artificial CH 4 release experiments, natural tracer approaches and flux chambers we mapped the spatial flux heterogeneity, and detected and validated a hot spot of CH 4 flux in a oligohaline restored marsh. Through chamber measurements during the months of April and May, we found that fluxes at the hot spot were on average as high as 6589 ± 7889 nmol m −2 s −1 whereas background flux from the open water were on average 15.2 ± 7.5 nmol m −2 s −1 . This study provides a novel tool to evaluate the spatial heterogeneity of fluxes around eddy‐covariance towers and creates important insights for the interpretation of hot spots of CH 4 flux, paving the way for future studies aiming to understand subsurface biogeochemical processes and the microbiological conditions that lead to the occurrence of hot spots and hot moments of CH 4 flux.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Drivers of Decadal Carbon Fluxes Across Temperate Ecosystems

Long-running eddy covariance flux towers provide insights into how the terrestrial carbon cycle operates over multiple timescales. Here, we evaluated variation in net ecosystem exchange (NEE) of carbon dioxide (CO 2 ) across the Chequamegon Ecosystem-Atmosphere Study AmeriFlux core site cluster in the upper Great Lakes region of the USA from 1997 to 2020. The tower network included two mature hardwood forests with differing management regimes (US-WCr and US-Syv), two fen wetlands with varying levels of canopy sheltering and vegetation (US-Los and US-ALQ), and a very tall (400 m) landscape-level tower (US-PFa). Together, they provided over 70 site-years of observations. The 19-tower Chequamegon Heterogenous Ecosystem Energy-balance Study Enabled by a High-density Extensive Array of Detectors 2019 campaign centered around US-PFa provided additional information on the spatial variation of NEE. Decadal variability was present in all long-term sites, but cross-site coherence in interannual NEE in the earlier part of the record became weaker with time as non-climatic factors such as local disturbances likely dominated flux time series. Average decadal NEE at the tall tower transitioned from carbon source to sink to near neutral over 24 years. Respiration had a greater effect than photosynthesis on driving variations in NEE at all sites. Declining snowfall offset potential increases in assimilation from warmer springs, as less-insulated soils delayed start of spring green-up. Higher CO 2 increased maximum net assimilation parameters but not total gross primary productivity. Stand-scale sites were larger net sinks than the landscape tower. Clustered, long-term carbon flux observations provide value for understanding the diverse links between carbon and climate and the challenges of upscaling these responses across space.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Atmospheric Turbulent Intermittency Over the Arctic Sea‐Ice Surface During the MOSAiC Expedition

Turbulent motions in the Arctic stable boundary layer are characterized by intermittency, but they are rarely investigated due to limited observations, in particular over the sea‐ice surface. In the present study, we explore the characteristics of turbulent intermittency over the Arctic sea‐ice surface using data collected during the Multidisciplinary drifting Observation for the Study of Arctic Climate expedition from October 2019 to September 2020. We first develop a new algorithm, which performs well in identifying the spectral gap over the Arctic sea‐ice surface. Then the characteristics of intermittency are investigated. It is found that the strength of intermittency increases under the conditions of light surface wind speed, small surface wind speed gradient, and strong surface air temperature gradient. The momentum flux, sensible heat flux, and latent heat flux calculated by raw eddy‐covariance fluctuations are overestimated by 3%, 10%, and 24%, respectively, because submesoscale motions are included. Furthermore, the characteristics of the atmospheric boundary layer structure under various intermittency conditions reveal that strong low‐level jets are favorable to surface turbulent motions that result in weak intermittency, while strong temperature inversions above the surface layer suppress surface turbulent motions and lead to strong intermittency.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Responses of Marginal and Intrinsic Water-Use Efficiency to Changing Aridity Using FLUXNET Observations

According to classic stomatal optimization theory, plant stomata are regulated to maximize carbon assimilation for a given water loss. A key component of stomatal optimization models is marginal water-use efficiency (mWUE), the ratio of the change of transpiration to the change in carbon assimilation. Although the mWUE is often assumed to be constant, variability of mWUE under changing hydrologic conditions has been reported. However, there has yet to be a consensus on the patterns of mWUE variabilities and their relations with atmospheric aridity. We investigate the dynamics of mWUE in response to vapor pressure deficit (VPD) and aridity index using carbon and water fluxes from 115 eddy covariance towers available from the global database FLUXNET. We demonstrate a non-linear mWUE-VPD relationship at a sub-daily scale in general; mWUE varies substantially at both low and high VPD levels. However, mWUE remains relatively constant within the mid-range of VPD. Despite the highly non-linear relationship between mWUE and VPD, the relationship can be informed by the strong linear relationship between ecosystem-level inherent water-use efficiency (IWUE) and mWUE using the slope, m *. We further identify site-specific m * and its variability with changing site-level aridity across six vegetation types. We suggest accurately representing the relationship between IWUE and VPD using Michaelis–Menten or quadratic functions to ensure precise estimation of mWUE variability for individual sites.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Is There a Scalar Atmospheric Surface Layer Within a Convective Boundary Layer? Implications for Flux Measurements

Top‐down entrainment shapes the vertical gradients of sensible heat, latent heat, and CO 2 fluxes, influencing the interpretation of eddy covariance (EC) measurements in the unstable atmospheric surface layer (ASL). Using large eddy simulations for convective boundary layer flows, we demonstrate that decreased temperature gradients across the entrainment zone increase entrainment fluxes by enhancing the entrainment velocity, amplifying the asymmetry between top‐down and bottom‐up flux contributions. These changes alter scalar flux profiles, causing flux divergence or convergence and leading to the breakdown of the constant flux layer assumption (CFLA) in the ASL. As a result, EC‐measured fluxes either underestimate or overestimate “true” surface fluxes during divergence or convergence phases, contributing to energy balance non‐closure. The varying degrees of the CFLA breakdown are a fundamental cause for the non‐closure issue. These findings highlight the underappreciated role of entrainment in interpreting EC fluxes, addressing non‐closure, and understanding site‐to‐site variability in flux measurements.

eddy covariance fluxes of scalars↗
DOE OSTI2022

The NEON Daily Isotopic Composition of Environmental Exchanges Dataset

The National Ecological Observatory Network (NEON) provides open-access measurements of stable isotope ratios in atmospheric water vapor (δ 2 H, δ 18 O) and carbon dioxide (δ 13 C) at different tower heights, as well as aggregated biweekly precipitation samples (δ 2 H, δ 18 O) across the United States. These measurements were used to create the NEON Daily Isotopic Composition of Environmental Exchanges (NEON-DICEE) dataset estimating precipitation (P; δ 2 H, δ 18 O), evapotranspiration (ET; δ 2 H, δ 18 O), and net ecosystem exchange (NEE; δ 13 C) isotope ratios. Statistically downscaled precipitation datasets were generated to be consistent with the estimated covariance between isotope ratios and precipitation amounts at daily time scales. Isotope ratios in ET and NEE fluxes were estimated using a mixing-model approach with calibrated NEON tower measurements. NEON-DICEE is publicly available on HydroShare and can be reproduced or modified to fit user specific applications or include additional NEON data records as they become available. The NEON-DICEE dataset can facilitate understanding of terrestrial ecosystem processes through their incorporation into environmental investigations that require daily δ 2 H, δ 18 O, and δ 13 C flux data.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Efficient Langevin dynamics for “noisy” forces

Efficient Boltzmann-sampling using first-principles methods is challenging for extended systems due to the steep scaling of electronic structure methods with the system size. Additionally, stochastic approaches provide a gentler system-size dependency at the cost of introducing “noisy” forces, which could limit the efficiency of the sampling. When the forces are deterministic, the first-order Langevin dynamics (FOLD) offers efficient sampling by combining a well-chosen preconditioning matrix S with a time-step-bias-mitigating propagator [G. Mazzola and S. Sorella, Phys. Rev. Lett. 118, 015703 (2017)]. However, when forces are noisy, S is set equal to the force-covariance matrix, a procedure that severely limits the efficiency and the stability of the sampling. Here, we develop a new, general, optimal, and stable sampling approach for FOLD under noisy forces. We apply it for silicon nanocrystals treated with stochastic density functional theory and show efficiency improvements by an order-of-magnitude.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Stability of Kuramoto networks subject to large and small fluctuations from heterogeneous and spatially correlated noise

Oscillatory networks subjected to noise are broadly used to model physical and technological systems. Due to their nonlinear coupling, such networks typically have multiple stable and unstable states that a network might visit due to noise. In this article, we focus on the assessment of fluctuations resulting from heterogeneous and spatially correlated noise inputs on Kuramoto model networks. We evaluate the typical, small fluctuations near synchronized states and connect the network variance to the overlap between stable modes of synchronization and the input noise covariance. Going beyond small to large fluctuations, we introduce the indicator mode approximation that projects the dynamics onto a single amplitude dimension. Such an approximation allows for estimating rates of fluctuations to saddle instabilities, resulting in phase slips between connected oscillators. Statistics for both regimes are quantified in terms of effective noise amplitudes that are compared and contrasted for several noise models. In conclusion, bridging the gap between small and large fluctuations, we show that a larger network variance does not necessarily lead to higher rates of large fluctuations.

42 ENGINEERING↗
DOE OSTI2023

Empirical probability and machine learning analysis of m , n = 2, 1 tearing mode onset parameter dependence in DIII-D H-mode scenarios

m, n = 2, 1 tearing mode onset empirical probability and machine learning analyses of a multiscenario DIII-D database of over 14 000 H-mode discharges show that the normalized plasma beta, the rotation profile, and the magnetic equilibrium shape have the strongest impact on the 2,1 tearing mode stability, in qualitative agreement with neoclassical tearing modes (m and n are the poloidal and toroidal mode numbers, respectively). In addition, 2,1 tearing modes are most likely to destabilize when n > 1 tearing modes are already present in the core plasma. The covariance matrix of tearing sensitive plasma parameters takes a nearly block-diagonal form, with the blocks incorporating thermodynamic, current and safety factor profile, separatrix shape, and plasma flow parameters, respectively. This suggests a number of paths to improved stability at fixed pressure and edge safety factor primarily by preserving a minimum of 1 kHz differential rotation, increasing the minimum safety factor above unity, using upper single null magnetic configuration, and reducing the core impurity radiation. In addition, lower triangularity, lower elongation, and lower pedestal pressure may also help to improve stability. Furthermore, the electron and ion temperature, collisionality, resistivity, internal inductance, and the parallel current gradient appear to only weakly correlate with the 2,1 tearing mode onsets in this database.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Simultaneous inference of equation of state parameters and unknown data errors with uncertainty quantification via hierarchical Bayesian posterior maximization

Equations of state (EOSs) are a key component in running hydrodynamic simulations as they relate the thermodynamic states for the material. The Davis reactants EOS is commonly used for modeling high explosives (HEs), and the EOS model parameters are calibrated using material specific data. The calibrations are often performed with uncertainty quantification via Bayesian inference to account for uncertainty in the data and generate ensembles of likely parameters. However, there are relatively few HE data sets to use for calibration and many are historical and lack error information. In this work, we simultaneously calibrate the Davis reactants EOS model parameters and unknown data error terms for the high explosive PBX 9501. To quantify the uncertainty in the models and the data, we use a Bayesian framework for the calibration and compute the hierarchical Bayesian posterior distribution with both a posteriori maximization approach and Markov Chain Monte Carlo. In general, we find that, given our assumptions, the two approaches result in similar calibrated parameters, posterior covariance matrices, and insights about the parameters but that the posterior maximization requires far less computational resources.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Nuclear Data–Induced Uncertainties in Criticality Safety Analyses for High-Burnup and Extended Enrichment Fuels

Criticality safety analyses are conducted to show compliance with regulatory standards and to demonstrate safe operational conditions during the storage and transportation of spent nuclear fuel. Given the increased interest in the industry in low-enriched uranium plus (LEU+) and higher-burnup fuel, it is important to study the impact of such fuels’ use on criticality safety analyses and the resulting nuclear data–induced uncertainties. Here, in this work, nominal pressurized water reactor assemblies with LEU+ fuel enrichments up to 8 wt% 235 U and high burnups up to 80 GWd/tonne U were studied. The assemblies were placed in a generic burnup credit cask GBC-32. As a result of the different covariance libraries, using the ENDF/B-VII.1 nuclear data library consistently resulted in lower nuclear data uncertainties than did the use of the ENDF/B-VIII.0 data library. The highest contribution in the nuclear data–induced uncertainties resulted from the major actinides, and their contribution increased with increasing burnup and enrichment.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Substantial carbon loss respired from a corn–soybean agroecosystem highlights the importance of careful management as we adapt to changing climate

Understanding agroecosystem carbon (C) cycle response to climate change and management is vital for maintaining their long-term C storage. We demonstrate this importance through an in-depth examination of a ten-year eddy covariance dataset from a corn–corn–soybean crop rotation grown in the Midwest United States. Ten-year average annual net ecosystem exchange (NEE) showed a net C sink of –0.39 Mg C ha –1 yr –1 . However, NEE in 2014 and 2015 from the corn ecosystem was 3.58 and 2.56 Mg C ha –1 yr –1 , respectively. Most C loss occurred during the growing season, when photosynthesis should dominate and C fluxes should reflect a net ecosystem gain. Partitioning NEE into gross primary productivity (GPP) and ecosystem respiration (ER) showed this C 'burp' was driven by higher ER, with a 51% (2014) and 57% (2015) increase from the ten-year average (15.84 Mg C ha –1 yr –1 ). GPP was also higher than average (16.24 Mg C ha –1 yr –1 ) by 25% (2014) and 37% (2015), but this was not enough to offset the C emitted from ER. This increased ER was likely driven by enhanced soil microbial respiration associated with ideal growing season climate, substrate availability, nutrient additions, and a potential legacy effect from drought.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Diel, seasonal, and inter-annual variation in carbon dioxide effluxes from lakes and reservoirs

Abstract Accounting for temporal changes in carbon dioxide (CO 2 ) effluxes from freshwaters remains a challenge for global and regional carbon budgets. Here, we synthesize 171 site-months of flux measurements of CO 2 based on the eddy covariance method from 13 lakes and reservoirs in the Northern Hemisphere, and quantify dynamics at multiple temporal scales. We found pronounced sub-annual variability in CO 2 flux at all sites. By accounting for diel variation, only 11% of site-months were net daily sinks of CO 2 . Annual CO 2 emissions had an average of 25% (range 3%–58%) interannual variation. Similar to studies on streams, nighttime emissions regularly exceeded daytime emissions. Biophysical regulations of CO 2 flux variability were delineated through mutual information analysis. Sample analysis of CO 2 fluxes indicate the importance of continuous measurements. Better characterization of short- and long-term variability is necessary to understand and improve detection of temporal changes of CO 2 fluxes in response to natural and anthropogenic drivers. Our results indicate that existing global lake carbon budgets relying primarily on daytime measurements yield underestimates of net emissions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Tractable Learning in Underexcited Power Grids

Estimating the structure of physical flow networks, such as power grids, is critical to secure delivery of energy. This article discusses statistical structure estimation in power grids in the “underexcited” regime, where a subset of internal nodes has zero injection fluctuations. Prior estimation algorithms based on nodal voltages fail for such grids as the voltage covariance matrix is not invertible. Here, we propose a novel topology learning algorithm for learning underexcited general networks. Our algorithm uses physics-informed conservation laws to first identify the zero-injection buses and their neighbors, and then estimates the remaining edges in the grid. We prove the asymptotic correctness of our algorithm for grids with nonadjacent internal zero-injection nodes. More important, we theoretically analyze our algorithm’s efficacy under noisy measurements, and determine bounds on maximum noise under which asymptotically correct recovery is guaranteed. Our approach is validated through simulations with voltage samples generated on test distribution grids with real injection data and nonlinear power flow models.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2020

Partitioning net carbon dioxide fluxes into photosynthesis and respiration using neural networks

The eddy covariance (EC) technique is used to measure the net ecosystem exchange (NEE) of CO 2 between ecosystems and the atmosphere, offering a unique opportunity to study ecosystem responses to climate change. NEE is the difference between the total CO 2 release due to all respiration processes (RECO), and the gross carbon uptake by photosynthesis (GPP). These two gross CO 2 fluxes are derived from EC measurements by applying partitioning methods that rely on physiologically based functional relationships with a limited number of environmental drivers. However, the partitioning methods applied in the global FLUXNET network of EC observations do not account for the multiple co-acting factors that modulate GPP and RECO flux dynamics. To overcome this limitation, we developed a hybrid data-driven approach based on combined neural networks (NN C-part ). NN C-part incorporates process knowledge by introducing a photosynthetic response based on the light-use efficiency (LUE) concept, and uses a comprehensive dataset of soil and micrometeorological variables as fluxes drivers. We applied the method to 36 sites from the FLUXNET2015 dataset and found a high consistency in the results with those derived from other standard partitioning methods for both GPP (R 2 > .94) and RECO (R 2 > .8). High consistency was also found for (a) the diurnal and seasonal patterns of fluxes and (b) the ecosystem functional responses. NNC-part performed more realistic than the traditional methods for predicting additional patterns of gross CO 2 fluxes, such as: (a) the GPP response to VPD, (b) direct effects of air temperature on GPP dynamics, (c) hysteresis in the diel cycle of gross CO 2 fluxes, (d) the sensitivity of LUE to the diffuse to direct radiation ratio, and (e) the post rain respiration pulse after a long dry period. In conclusion, NN C-part is a valid data-driven approach to provide GPP and RECO estimates and complementary to the existing partitioning methods.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Identifying dominant environmental predictors of freshwater wetland methane fluxes across diurnal to seasonal time scales

While wetlands are the largest natural source of methane (CH 4 ) to the atmosphere, they represent a large source of uncertainty in the global CH 4 budget due to the complex biogeochemical controls on CH 4 dynamics. Here we present, to our knowledge, the first multi-site synthesis of how predictors of CH 4 fluxes (FCH4) in freshwater wetlands vary across wetland types at diel, multiday (synoptic), and seasonal time scales. In this work we used several statistical approaches (correlation analysis, generalized additive modeling, mutual information, and random forests) in a wavelet-based multi-resolution framework to assess the importance of environmental predictors, nonlinearities and lags on FCH4 across 23 eddy covariance sites. Seasonally, soil and air temperature were dominant predictors of FCH4 at sites with smaller seasonal variation in water table depth (WTD). In contrast, WTD was the dominant predictor for wetlands with smaller variations in temperature (e.g., seasonal tropical/subtropical wetlands). Changes in seasonal FCH4 lagged fluctuations in WTD by ~17 ± 11 days, and lagged air and soil temperature by median values of 8 ± 16 and 5 ± 15 days, respectively. Temperature and WTD were also dominant predictors at the multiday scale. Atmospheric pressure (PA) was another important multiday scale predictor for peat-dominated sites, with drops in PA coinciding with synchronous releases of CH 4 . At the diel scale, synchronous relationships with latent heat flux and vapor pressure deficit suggest that physical processes controlling evaporation and boundary layer mixing exert similar controls on CH 4 volatilization, and suggest the influence of pressurized ventilation in aerenchymatous vegetation. In addition, 1- to 4-h lagged relationships with ecosystem photosynthesis indicate recent carbon substrates, such as root exudates, may also control FCH4. By addressing issues of scale, asynchrony, and nonlinearity, this work improves understanding of the predictors and timing of wetland FCH4 that can inform future studies and models, and help constrain wetland CH 4 emissions.

59 BASIC BIOLOGICAL SCIENCES↗