Stability and Mixing of a Vertical Axisymmetric Buoyant Jet in Shallow Water

The stability, mixing and effect of downstream control on axisymmetric turbulent buoyant jets discharging vertically into shallow stagnant water is studied using 3D Reynolds-averaged Navier–Stokes equations (RANS) combined with a buoyancy-extended k –ε model. The steady axisymmetric turbulent flow, temperature (or tracer concentration) and turbulence fields are computed using the finite volume method on a high resolution grid. The numerical predictions demonstrate two generic flow patterns for different turbulent heated jet discharges and environmental parameters (i) a stable buoyant discharge with the mixed fluid leaving the vertical jet region in a surface warm water layer; and (ii) an unstable buoyant discharge with flow recirculation and re-entrainment of heated water. A stratified counterflow region always appears in the far-field for both stable and unstable buoyant discharges. Provided that the domain radius L exceeds about 6H, the near field interaction and hence discharge stability is governed chiefly by the jet momentum length scale to depth ratio l_M/H, regardless of downstream control. The near field jet stability criterion is determined to be l_M/H = 3.5. A radial internal hydraulic jump always exists beyond the surface impingement region, with a 3- to 6-fold increase in dilution across the jump compared with vertical buoyant jet mixing. The predicted stability category, velocity and temperature/concentration fields are well-supported by experiments of all previous investigators.     

Economic impacts of climate change on agriculture: The importance of additional climatic variables other than temperature and precipitation

Climate change shifts the distributions of a set of climatic variables, including temperature, precipitation, humidity, wind speed, sunshine duration, and evaporation. This paper explores the importance of those additional climatic variables other than temperature and precipitation. Using the county-level agricultural data from 1980 to 2010 in China, we find that those additional climatic variables, especially humidity and wind speed, are critical for crop growth. Therefore, omitting those variables is likely to bias the predicted impacts of climate change on crop yields. In particular, omitting humidity tends to overpredict the cost of climate change on crop yields, while ignoring wind speed is likely to underpredict the effect. Our preferred specification indicates that climate change is likely to decrease the yields of rice, wheat, and corn in China by 36.25%, 18.26%, and 45.10%, respectively, by the end of this century.     

Combining data and simulated data for space–time fields: application to ozone

This paper presents a theory for modeling random environmental spatial-temporal fields that allows simulated data (numerical-physical model output) to be combined with measurements made at fixed monitoring sites. That theory involves Bayesian hierarchical models that provide temporal forecasts and spatial predictions along with appropriate credibility intervals. A by-product is a method for re-calibrating the simulated data to bring it into line with the measurements for certain applications. While the approach covers a broad domain of potential applications, this paper addresses a field of particular importance, ground level ozone concentrations over the eastern and central USA. A univariate model is developed and illustrated with hourly ozone fields. A multivariate alternative is also provided and illustrated with daily concentration fields. The forecasts and predictions they provide are compared with those from other approaches.     

Intercomparison of Two Models,ETA and RAMS,with TRACT Field Campaign Data

In this work a model intercomparison between RAMS and ETA models is carried out, with the aim of evaluating the quality and accuracy of these mesoscale models in reproducing the time evolution of the meteorology in real complex terrain. This is of great importance not only for meteorological forecast but also for air quality assessment. Numerical simulations are performed to reproduce the mean variables' fields and to compare them with measurements collected during the field campaign TRACT. The domain covers the Rhine valley and surrounding mountainous region and we consider a time period of two days. Results from simulations are compared to observations relative to ground stations and radiosoundings. A qualitative analysis is joined to a quantitative estimation of some reference statistical indexes. Both RAMS and ETA models performances are satisfactory when compared to the measured data and also their relative agreement is good. The mean variable fields are reproduced with a satisfactory degree of reliability, even if the simulated profiles are not able to describe the largest fluctuations of the variables. At the surface stations, the best agreement between predictions and observations is obtained for the wind velocity, while the quality of the results is lower for temperature and humidity.     

Preferential leaching of nitrate,chloride and phosphate in an Australian clay soil

This paper investigates leaching of water and nutrients (NO^‐ ₃, Cl^‐, PO^3‐ ₄) from the unsaturated layer in an Australian soil using a multisegment percolation system (MPS). Large undisturbed soil cores were collected from a clay‐based, basaltic plain, agricultural soil at Grassmere, 300 km west of Melbourne, Australia. Significant heterogeneity (or preferential flow) of effluent moisture and solutes was detected (one‐way ANO VA, p < 0.001). Fifty percent of the applied nitrate and chloride leached from the soil core within three days after initial application. Hundred percent of the applied nitrate and chloride leached from the soil core within 8 days after application. These results indicate little incorporation into the soil matrix, and possible denitrification or mineralisation. In contrast, after 18 days, less than 1 % of the total applied phosphates leached from the soil, indicating strong adsorption. Our experiments indicate considerable heterogeneity in water flow patterns and solute leaching on a small spatial scale. Very rapid transport of nitrate and chloride through the soil was evident, in comparison phosphate leaching was negligible. These results have important implications for the management of nutrient schedules in agricultural soils, particularly those located in the Western District of Victoria, Australia.     

A simplified computational analysis of turbulent plumes and jets

The application of computational fluid dynamics (CFD), particularly Large Eddy Simulation, for the modelling of buoyant turbulent plumes, has been demonstrated to be very accurate, but computationally expensive. Here a more basic, and therefore more generally practicable, approach is presented. Commercial CFD software is used to model such plumes using Reynolds-Averaged Navier-Stokes (RANS) turbulence models. A careful comparison is made between the numerical predictions and well-established results regarding the bulk properties of plumes. During this process, we are able to observe the well-known approximate Gaussian nature of the plume and achieve quantitative agreement with empirical plume spread coefficients. The use of numerical modelling allows for the investigation of the flow field and turbulence in those regions of the plume of most interest—the plume edge and near source regions. A comprehensive sensitivity study is conducted to identify the limits of applicability of this modelling approach. It is shown that the standard modelling approach of Morton, Taylor and Turner, which introduced the well-known entrainment assumption, pertains in a region well above the source region. At the plume edge, the levels of turbulence are contrasted with the value of the entrainment parameter. Finally, the effects of forcing the plumes with additional momentum at the source are considered, including the case of a pure jet. We show how these forced plumes eventually lose their momentum excess and tend to the behaviour of a pure, buoyant plume.     

Maximum feasible distance of windborne cross-pollination in Brassica napus: A ‘mass budget’ model

A theory of gene dispersal by wind pollination can make an important contribution to understanding the viability and evolution of important plant groups in the Earth's changing landscape and it can be applied to evaluate concerns about the spread of engineered genes from genetically modified (GM) crops into conventional varieties via windborne pollen. Here, we present a model of cross-pollination between plant populations due to the wind. We perform a ‘mass budget’ of pollen by accounting for the number of pollen grains from release in the source population, dispersal from the source to the sink population by the wind, and deposition on receptive surfaces in the sink population. Our model can be parameterised for any wind-pollinated species, but we apply it to Brassica napus (oilseed rape or canola) to investigate the threat posed by wind pollination to GM confinement in agriculture. Specifically, we calculate the maximum feasible distance at which a particular level of windborne gene dispersal could be attained. This is equivalent to the separation distance between populations or fields required to achieve a given threshold of gene dispersal or adventitious GM presence. As required, model predictions of the upper bounds on levels of wind-mediated gene dispersal exceed observations from a wide range of published studies. For a level of gene dispersal below 0.9%, which is the EU threshold for GM adventitious presence, we predict that the maximum feasible distance for agricultural fields of B. napus is 1000 m, regardless of field shape and direction of prevailing winds. For fields closer than 1000 m, our model results do not necessarily imply that the 0.9% threshold is likely to be breached, because in this instance we have conservatively set the values of parameters where current knowledge is limited. We also predict that gene dispersal is reduced by 50% when the lag in peak flowering between the source and sink populations is 13 days, and reduced by 90% when the lag is 24 days. We identify further measurements necessary to improve the accuracy of the model predictions.     

A Large-Eddy Simulation Study of the Convective Boundary Layer over Philadelphia during the 1999 Summer NE-OPS Campaign

This study presents a large-eddy simulation (LES) study of the convective boundary layer on August 1, 1999 over Philadelphia, PA during a summer ozone episode. The study is an evaluation of the Colorado State University's Regional Atmospheric Modeling System Version 4.3 (RAMS4.3) with the LES option using Northeast Oxidant and Particulate Study (NE-OPS) data. Simulations were performed with different imposed sensible heat fluxes at the ground surface. The model was initialized with the atmospheric sounding data collected at Philadelphia at 1230 UTC and model integrations continued till 2130 UTC. The resulting mean profiles of temperature and humidity obtained from the LES model were compared with atmospheric soundings, tethered balloon and aircraft data collected during the NE-OPS 1999 field campaign. Also the model-derived vertical profiles of virtual temperature were compared with NE-OPS Radio Acoustic Sounder System (RASS) data while the humidity profiles were compared with NE-OPS lidar data. The comparison of the radiosonde data with the LES model predictions suggests that the growth of the mixing layer is reasonably well simulated by the model. Overall, the agreement of temperature predictions of the LES model with the radiosonde observations is good. The model appears to underestimate humidity values for the case of higher imposed sensible heat flux. However, the humidity values in the mixing layer agree quite well with radiosonde observations for the case of lower imposed sensible heat flux. The model-predicted temperature and humidity profiles are in reasonable agreement with the tethered balloon data except for some small overestimation of temperature at lower layers and some underestimation of humidity values. However, the humidity profiles as simulated by the model agree quite well with the tethered balloon data for the case of lower imposed sensible heat flux. The model-predicted virtual temperature profile is also in better agreement with RASS data for the case of lower imposed sensible heat flux. The model-predicted temperature profile further agrees quite well with aircraft data for the case of lower imposed heat flux. However, the relative humidity values predicted by the model are lower compared with the aircraft data. The model-predicted humidity profiles are only in partial agreement with the lidar data. The results of this study suggest that the explicitly resolved energetic eddies seem to provide the correct forcing necessary to produce good agreement with observations for the case of an imposed sensible heat flux of 0.1 K m s^–1 at the surface.     

A σ-coordinate transport model coupled with rotational boussinesq-type equations

This paper describes a σ-coordinate scalar transport model coupled with a Boussinesq-type hydrodynamic model. The Boussinesq model has the ability to calculate both three-dimensional velocity distributions and the water surface motion. To capture ‘dispersion’ processes in open channel flow, horizontal vorticity effects induced by a bottom shear stress are included in the Boussinesq model. Thus, a reasonable representation of vertical flow structure can be captured in shallow and wavy flow fields. To solve the coupled Boussinesq and scalar transport system, a finite-volume method, based on a Godunov-type scheme with the HLL Riemann solver, is employed. Basic advection and advection–diffusion numerical tests in a non-rectangular domain were carried out and the computed results show good agreement with analytic solutions. With quantitative comparisons of dispersion experiments in an open channel, it is verified that the proposed coupled model is appropriate for both near and far field scalar transport predictions. From numerical simulations in the surf zone, physically reasonable results showing expected vertical variation are obtained.     

Optimization of ecosystem model parameters using spatio-temporal soil moisture information

Parameters in process-based terrestrial ecosystem models are often nonlinearly related to the water flux to the atmosphere, and they also change temporally and spatially. Therefore, for estimating soil moisture, process-based terrestrial ecosystem models inevitably need to specify spatially and temporally variant model parameters. This study presents a two-stage data assimilation scheme (TSDA) to spatially and temporally optimize some key parameters of an ecosystem model which are closely related to soil moisture. At the first stage, a simplified ecosystem model, namely the Boreal Ecosystem Productivity Simulator (BEPS), is used to obtain the prior estimation of daily soil moisture. After the spatial distribution of 0–10 cm surface soil moisture is derived from remote sensing, an Ensemble Kalman Filter is used to minimize the difference between the remote sensing model results, through optimizing some model parameters spatially. At the second stage, BEPS is reinitialized using the optimized parameters to provide the updated model predictions of daily soil moisture. TSDA has been applied to an arid and semi-arid area of northwest China, and the performance of the model for estimating daily 0–10 cm soil moisture after parameter optimization was validated using field measurements. Results indicate that the TSDA developed in this study is robust and efficient in both temporal and spatial model parameter optimization. After performing the optimization, the correlation (r²) between model-predicted 0–10 cm soil moisture and field measurement increased from 0.66 to 0.75. It is demonstrated that spatial and temporal optimization of ecosystem model parameters can not only improve the model prediction of daily soil moisture but also help to understand the spatial and temporal variation of some key parameters in an ecosystem model and the corresponding ecological mechanisms controlling the variation.     

Sensitivity of ecological models to their climate drivers: statistical ensembles for forcing.

Global and regional numerical models for terrestrial ecosystem dynamics require fine spatial resolution and temporally complete historical climate fields as input variables. However, because climate observations are unevenly spaced and have incomplete records, such fields need to be estimated. In addition, uncertainty in these fields associated with their estimation are rarely assessed. Ecological models are usually driven with a geostatistical model's mean estimate (kriging) of these fields without accounting for this uncertainty, much less evaluating such errors in terms of their propagation in ecological simulations. We introduce a Bayesian statistical framework to model climate observations to create spatially uniform and temporally complete fields, taking into account correlation in time and space, spatial heterogeneity, lack of normality, and uncertainty about all these factors. A key benefit of the Bayesian model is that it generates uncertainty measures for the generated fields. To demonstrate this method, we reconstruct historical monthly precipitation fields (a driver for ecological models) on a fine resolution grid for a climatically heterogeneous region in the western United States. The main goal of this work is to evaluate the sensitivity of ecological models to the uncertainty associated with prediction of their climate drivers. To assess their numerical sensitivity to predicted input variables, we generate a set of ecological model simulations run using an ensemble of different versions of the reconstructed fields. We construct such an ensemble by sampling from the posterior predictive distribution of the climate field. We demonstrate that the estimated prediction error of the climate field can be very high. We evaluate the importance of such errors in ecological model experiments using an ensemble of historical precipitation time series in simulations of grassland biogeochemical dynamics with an ecological numerical model, Century. We show how uncertainty in predicted precipitation fields is propagated into ecological model results and that this propagation had different modes. Depending on output variable, the response of model dynamics to uncertainty in inputs ranged from uncertainty in outputs that matched that of inputs to those that were muted or that were biased, as well as uncertainty that was persistent in time after input errors dropped.     

内蒙古河套灌区秋浇对不同类型农田土壤盐分淋失的影响

秋浇是河套灌区传统的秋后淋盐、春季保墒的一种特殊灌溉制度 ,秋浇对不同类型农田盐分损失影响的试验结果表明 ,秋浇前白菜地土壤含水率最高 ,向日葵地最低 ,秋浇后 5种农田土壤含水率均较高 ,且彼此间无显著性差异。秋浇后不同土壤剖面盐分的损失量按照表层 (0～ 40cm)、中层 (40～ 80cm)、深层 (80～ 1 2 0cm)依次递减(向日葵地除外 )。不同农田的盐分淋失量按照白菜地、向日葵地、玉米地、小麦 -玉米地、小麦地依次递减。非生长季灌溉下 ,农田土壤 0～ 1 0 0cm土层中盐分的淋失量与土壤储水量的增加量 (灌溉前土壤含水率 <2 50 g·kg- 1 ,r =0 .990 2 )密切相关 ,但与灌溉前土壤盐分初始含量无显著性关系     

A numerical model of wave- and current-driven nutrient uptake by coral reef communities

We developed a numerical model capable of simulating the spatial zonation of nutrient uptake in coral reef systems driven by hydrodynamic forcing (both from waves and currents). Relationships between nutrient uptake and bed stress derived from flume and field studies were added to a four-component biogeochemical model embedded within a three-dimensional (3-D) hydrodynamic ocean model coupled to a numerical wave model. The performance of the resulting coupled physical-biogeochemical model was first evaluated in an idealized one-dimensional (1-D) channel for both a pure current and a combined wave-current flow. Waves in the channel were represented by an oscillatory flow with constant amplitude and frequency. The simulated nutrient concentrations were in good agreement with the analytical solution for nutrient depletion along a uniform channel, as well as with existing observations of phosphate uptake across a real reef flat. We then applied this integrated model to investigate more complex two-dimensional (2-D) nutrient dynamics, firstly to an idealized coral reef-lagoon morphology, and secondly to a realistic section of Ningaloo Reef in Western Australia, where nutrients were advected into the domain via alongshore coastal currents. Both the idealized reef and Ningaloo Reef simulations showed similar patterns of maximum uptake rates on the shallow forereef and reef crest, and with nutrient concentration decreasing as water flowed over the reef flat. As a result of the cumulative outflow of nutrient-depleted water exiting the reef channels and then being advected down the coast by alongshore currents, both reef simulations exhibited substantial alongshore variation in nutrient concentrations. The coupled models successfully reproduced the observed spatial-variability in nitrate concentration across the Ningaloo Reef system.     

Stratification and Circulation in a Shallow Turbid Waterbody

Shallow waterbodies are often assumed to be well mixed in the vertical. However, when they are characterised by high turbidity levels, absoption of solar heating within a relatively thin surface layer can produce thermal stratification. Results from an intensive monitoring program have been combined with three-dimensional circulation modelling to examine the diurnal stratification cycles in a small turbid waterbody. The waterbody, known as Rushy Billabong, is located in southeastern Australia and its high turbidity coupled with forcing by wind and solar radiation resulted in regular diurnal cycles of stratification and overturning. Under conditions of light wind and high solar radiation, the model results were generally consistent with the observed temperature field. However, under stronger winds, preferential cooling and sinking of shallow water around the edge of the lake began to contribute significantly to the interior stratification. Model estimates then became more sensitive to the detailed bathymetry and the choice of turbulence parameterisation. The level of stratification is also shown to influence the circulation in the billabong by trapping the wind-driven flow near the surface. Insights provided by the observations and modelling may have broader implications for the management of small turbid systems such as settling ponds, aquaculture ponds, and some natural wetlands.     

Mating behaviour of the tiger blue butterfly (Tarucus theophrastus): competitive mate-searching when not all females are captured

Summary Females of the tiger blue butterfly (Tarucus theophrastus) fly upwind in search of the larval hostplant. Males perch or patrol the downwind edges, searching for incoming females or those already on the bush. A model of competitive mate-searching is developed for the case where not all receptive females are contacted by searching males: the model differs from the earlier ones of Parker (1970, 1974), particularly when few males are involved. Observed distributions of males upon bushes of different sizes agree better with the predictions of the game theory model than with a random distribution of males. The described model may be widely applicable to populations where females remate, but frequently evade detection by searching males.     

Advanced integral model for groups of interacting round turbulent buoyant jets

An integral model that combines all advantages of Superposition Method (SM), Entrainment Restriction Approach (ERA) and Second Order Approach (SOA) is proposed to predict the mean axial velocity and concentration fields of a group of N interacting vertical round turbulent buoyant jets. SM is successful in predicting the fields of mean axial velocity and mean concentration for a group of N interacting jets or plumes and ERA is advantageous in predicting the above fields for either two or large number (N → ∞) of interacting buoyant jets in the whole range of buoyancy. SOA takes into consideration in a dynamic way the turbulent contribution to the momentum and buoyancy fluxes and provides better accuracy than the usual procedures. A novelty of the proposed model is the production and utilisation of advanced profile distributions, convenient for the mean axial velocities and concentrations in a cross-section of the entire group of buoyant jets. These profiles are developed on the basis of flux conservation of momentum, buoyancy and kinetic energy for the mean motion. They enhance dynamic adaptation of the individual buoyant jet axes to the group centreline. Due to these profile distributions, the present model owns generality of application and better accuracy of predictions compared to usual integral models using simple Gaussian or top-hat profiles; thus it conferred the name Advanced Integral Model (AIM). AIM is herein applied to predict the mean flow properties of two different arrangement types of any number of buoyant jets: (a) linear diffusers and (b) rosette-type risers. Present results are compared to available experimental data and traditional solutions based on Gaussian profiles. Findings may be useful for design purposes and environmental impact assessment.     

The Biodiversity Integrity Index: An Illustration Using Ants in Western Australia

Although Western Australia is a relatively unpopulated region, considerable areas of native vegetation have been modified by agricultural clearing, rangeland grazing, urbanization, road construction, and mining. Ant diversity is reduced and community composition changed by each of these land uses. Road construction has the greatest long-term effect on the alpha diversity of ants, followed by agricultural clearing, mining, urbanization, and rangeland grazing. We present data on the extent of these various land uses in each major Western Australian vegetation association. Then, examples of ant diversity and community composition for each land use are coupled with geographic information system data on the extent of each land use in the various vegetation associations to calculate indices of "biodiversity integrity." The extent of biodiversity integrity in each region concurs with a subjective opinion of the condition of each unit. Agricultural clearing, followed by rangeland grazing, were found responsible for the greatest loss of ant biodiversity integrity. The findings relate to Australia in general and may serve as a framework for estimating losses of biodiversity integrity in other regions of the world in taxa other than ants.     

A comparative study of prognostic MM5 meteorological modeling with aircraft, wind profiler, lidar, tethered balloon and RASS data over Philadelphia during a 1999 summer episode

This study presents a comparative evaluation of the prognostic meteorological Fifth Generation NCAR Pennsylvania State University Mesoscale Model (MM5) using data from the Northeast Oxidant and Particle Study (NE-OPS) research program collected over Philadelphia, PA during a summer episode in 1999. A set of model simulations utilizing a nested grid of 36 km, 12 km and 4 km horizontal resolutions with 21 layers in the vertical direction was performed for a period of 101 h from July 15, 1999; 12 UTC to July 19, 1999; 17 UTC. The model predictions obtained with 4 km horizontal grid resolution were compared with the NE-OPS observations. Comparisons of model temperature with aircraft data revealed that the model exhibited slight underestimation as noted by previous investigators. Comparisons of model temperature with aircraft and tethered balloon data indicate that the mean absolute error varied up to 1.5 °C. The comparisons of model relative humidity with aircraft and tethered balloon indicate that the mean relative error varied from –11% to –22% for the tethered balloon and from –5% to –30% for the aircraft data. The mean relative error for water vapor mixing ratio with respect to the lidar data exhibited a negative bias consistent with the humidity bias corresponding to aircraft and tethered balloon data. The tendency of MM5 to produce estimates of very low wind speeds, especially in the early-mid afternoon hours, as noted by earlier investigators, is seen in this study also. It is indeed true that the initial fields as well as the fields utilized in the data assimilation also contribute to some of the differences between the model and observations. Studies such as these which compare the grid averaged mean state variables with observations have inherent difficulties. Despite the above limitations, the results of the present study broadly conform to the general traits of MM5 as noted by earlier investigators.     

水文因素影响稻田氮磷流失的研究进展

水文因素为农田氮磷元素的迁移提供了动力和载体,对稻田氮磷元素的流失具有重要影响。利用新型的稻田水分管理模式以代替传统的水分管理模式,对于有效控制面源污染具有重要意义。在综合调研国内外已有研究成果和最新进展的基础上,阐述了稻田氮磷流失特征、传统水分管理模式的弊端,并从新型稻田水分管理模式的种类、削减氮磷流失的效果与机理、与水平衡模型、营养负荷模型等结合应用等几个方面综述了国内外水文因素影响稻田氮磷流失的研究现状。进行稻田土壤吸附氮磷容量及人工调节机制的研究,在部分地区开展流域化水分管理系统研究,以及适用于我国稻田的水平衡模型和营养负荷模型的建立和深入研究应为今后我国关于水文因素影响稻田氮磷流失方面研究的一些方向。     

Greenhouse gas fluxes in southeastern U.S. coastal plain wetlands under contrasting land uses

Whether through sea level rise or wetland restoration, agricultural soils in coastal areas will be inundated at increasing rates, renewing connections to sensitive surface waters and raising critical questions about environmental trade-offs. Wetland restoration is often implemented in agricultural catchments to improve water quality through nutrient removal. Yet flooding of soils can also increase production of the greenhouse gases nitrous oxide and methane, representing a potential environmental trade-off. Our study aimed to quantify and compare greenhouse gas emissions from unmanaged and restored forested wetlands, as well as actively managed agricultural fields within the North Carolina coastal plain, USA. In sampling conducted once every two months over a two-year comparative study, we found that soil carbon dioxide flux (range: 8000-64 800 kg CO2 x ha(-1) x yr(-1)) comprised 66-100% of total greenhouse gas emissions from all sites and that methane emissions (range: -6.87 to 197 kg CH4 x ha(-1) x yr(-1)) were highest from permanently inundated sites, while nitrous oxide fluxes (range: -1.07 to 139 kg N2O x ha(-1) x yr(-1)) were highest in sites with lower water tables. Contrary to predictions, greenhouse gas fluxes (as CO2 equivalents) from the restored wetland were lower than from either agricultural fields or unmanaged forested wetlands. In these acidic coastal freshwater ecosystems, the conversion of agricultural fields to flooded young forested wetlands did not result in increases in greenhouse gas emissions.