ALBERT

Description: Biogeochemical cycles in the coastal ocean are changing and will continue to change in response to a changing climate. Effects on the oxygen and carbon cycles are particularly important, as either episodic or permanent shifts toward lower oxygen and/or higher inorganic carbon conditions can impact coastal ecosystems negatively. Here we study the sensitivity of these cycles to changes that may occur in the coastal ocean, focusing on a summer wind-driven upwelling region off southern Vancouver Island shelf. We use a quasi 2-D configuration of the Regional Ocean Modeling System (ROMS) to perform six sensitivity experiments. Results indicate that carbon and oxygen cycles in this region may be significantly affected by an altered upwelling season, a shallower offshore Oxygen Minimum Zone, and a carbon-enriched environment. Combinations of these scenarios suggest a potentially increasing risk for the development of coastal hypoxia and corrosive conditions in the region.

Description: Longitudinal and diel measurements of dual isotope composition (δ15N and δ18O) in nitrate (NO3-N) were made in the Ichetucknee River, a large (∼8 m3 s−1), entirely spring-fed river in North Florida, to determine whether isotopic variation can deconvolve assimilatory and dissimilatory removal. Comparing nitrate concentrations and isotope composition during the day and night we predicted (1) daytime declines in total fractionation due to low assimilatory fractionation and (2) diurnal variation in dual isotope coupling between 1:1 (assimilation) and 2:1 (denitrification). Five daytime longitudinal transects comprising 10 sampling stations showed consistent NO3-N removal (25–35% of inputs) and modest fractionation (15εtotal between −2 and −6‰, enriching the residual nitrate pool). Lower fractionation (by ∼1‰) during two nighttime transects, suggests higher fractionation due to assimilation than denitrification. Total fractionation was significantly negatively associated with discharge, input [NO3-N], N mass removal, and fractional water loss. Despite well-constrained mass balance estimates that denitrification dominated total N removal, isotope coupling was consistently 1:1, both for longitudinal and diel sampling. Hourly samples on two dates at the downstream location showed significant diel variation in concentration ([NO3-N] amplitude = 60 to 90 μg N L−1) and isotope composition (δ15N amplitude = −0.7‰ to −1.6‰). Total fractionation differed between day and night only on one date but estimated assimilatory fractionation assuming constant denitrification was highly variable and implausibly large (for N, 15ε = −2 to −25‰), suggesting that fractionation and removal due to denitrification is not diurnally constant. Pronounced counterclockwise hysteresis in the relationship between [NO3-N] and δ15N suggests diel variation in N isotope dynamics. Together, low fractionation, isotope versus concentration hysteresis, and consistent 1:1 isotope coupling suggests that denitrification is controlled by NO3− diffusion into the benthic sediments, the length of which is mediated by riverine oxygen dynamics. While using dual isotope behavior to deconvolve removal pathways was not possible, isotope measurements did yield valuable information about riverine N cycling and transformations.

Description: Thermokarst lakes are thought to have been an important source of methane (CH4) during the last deglaciation when atmospheric CH4 concentrations increased rapidly. Here we demonstrate that meltwater from permafrost ice serves as an H source to CH4 production in thermokarst lakes, allowing for region-specific reconstructions of δDCH4 emissions from Siberian and North American lakes. δDCH4 reflects regionally varying δD values of precipitation incorporated into ground ice at the time of its formation. Late Pleistocene-aged permafrost ground ice was the dominant H source to CH4 production in primary thermokarst lakes, whereas Holocene-aged permafrost ground ice contributed H to CH4 production in later generation lakes. We found that Alaskan thermokarst lake δDCH4 was higher (−334 ± 17‰) than Siberian lake δDCH4 (−381 ± 18‰). Weighted mean δDCH4 values for Beringian lakes ranged from −385‰ to −382‰ over the deglacial period. Bottom-up estimates suggest that Beringian thermokarst lakes contributed 15 ± 4 Tg CH4 yr−1 to the atmosphere during the Younger Dryas and 25 ± 5 Tg CH4 yr−1 during the Preboreal period. These estimates are supported by independent, top-down isotope mass balance calculations based on ice core δDCH4 and δ13CCH4 records. Both approaches suggest that thermokarst lakes and boreal wetlands together were important sources of deglacial CH4.

Description: A perturbation of the global carbon cycle has often been used for interpreting the Frasnian-Famennian (F-F) mass extinction. However, the changes of atmospheric CO2 level (pCO2) during this interval are much debatable. To illustrate the carbon cycle during F-F transition, paired inorganic (δ13Ccarb) and organic (δ13Corg) carbon isotope analyses were carried out on two late Devonian carbonate sequences (Dongcun and Yangdi) from south China. The larger amplitude shift of δ13Corg compared to δ13Ccarb and its resultant Δ13C (Δ13C = δ13Ccarb − δ13Corg) decrease indicate decreased atmospheric CO2 level around the F-F boundary. The onset of pCO2 level decrease predates that of marine regressions, which coincide with the beginning of conodont extinctions, suggesting that temperature decrease induced by decreased greenhouse effect of atmospheric CO2 might have contributed to the F-F mass extinction.

Description: A large increase in near-infrared (NIR) reflectance of Amazon forests during the light-rich dry season and a corresponding decrease during the light-poor wet season has been observed in satellite measurements. This increase has been variously interpreted as seasonal change in leaf area resulting from net leaf flushing in the dry season or net leaf abscission in the wet season, enhanced photosynthetic activity during the dry season from flushing new leaves and as change in leaf scattering and absorption properties between younger and older leaves covered with epiphylls. Reconciling these divergent views using theory and observations is the goal of this article. The observed changes in NIR reflectance of Amazon forests could be due to similar, but small, changes in NIR leaf albedo (reflectance plus transmittance) resulting from the exchange of older leaves for newer ones, but with the total leaf area unchanged. However, this argument ignores accumulating evidence from ground-based reports of higher leaf area in the dry season than the wet season, seasonal changes in litterfall and does not satisfactorily explain why NIR reflectance of these forests decreases in the wet season. More plausibly, the increase in NIR reflectance during the dry season and the decrease during the wet season would result from changes in both leaf area and leaf optical properties. Such change would be consistent with known phenological behavior of tropical forests, ground-based reports of seasonal changes in leaf area, litterfall, leaf optical properties and fluxes of evapotranspiration, and thus, would reconcile the various seemingly divergent views.

Description: In this study, tundra N2O and CH4 fluxes were measured from one seabird sanctuary (SBT) and two non-seabird colonies (NST-I and NST-II) in Ny-Ålesund (79°55′N, 11°56′E), Svalbard during the summers of 2008 and 2009. N2O and CH4 fluxes from SBT showed large temporal and spatial variations depending on the intensity of seabird activity. High seabird activity sites showed large N2O and CH4 emissions while low N2O and CH4 emissions, even CH4 uptake occurred at medium and low seabird activity sites. Overall the mean fluxes were 18.3 ± 3.6 μg N2O m−2 h−1 and 53.5 ± 20.3 μg CH4 m−2 h−1 from tundra SBT whereas tundra NST-I and NST-II represented a relatively weak N2O source (8.3 ± 13.2 μg N2O m−2 h−1) and strong CH4 sink (−82.8 ± 22.3 μg CH4 m−2 h−1). Seabird activity was the strongest control of N2O and CH4 fluxes compared with soil temperature and moisture, and high N2O and CH4 emissions were created by soil physical and chemical processes (the sufficient supply of nutrients NH4+–N, NO3−–N, total nitrogen, total phosphorus and total carbon from seabird guano, seabird tramp and appropriate water content) related to the seabird activity. Our work suggests that tundra ecosystems impacted by seabird activity are the potential “hotspots” for N2O and CH4 emissions although these sources have been largely neglected at present. Furthermore the combination of seabird activity and warming climate will likely further enhance N2O and CH4 emissions from the High Arctic tundra.

Description: Riparian groundwater can exhibit considerable patchiness in the concentration and reactivity of dissolved organic matter (DOM), which ultimately shapes subsurface biogeochemical transformations. Free and combined amino acids are bioavailable constituents of DOM, and their concentration and composition can provide valuable information about the diagenetic state of DOM. Based on riparian groundwater samples and relevant DOM end-member samples, we adapted the amino-acid-based marine DOM degradation index (DI) to groundwater. The groundwater DI was applied to evaluate the spatial and temporal variability in the bioavailability and diagenetic state of riparian DOM in a restored and a channelized section of the River Thur, Switzerland. Among different indicators for DOM diagenetic state (total hydrolysable amino acid concentrations, C-normalized yields, and the contribution of nonprotein amino acids), the groundwater DI correlated best with the activity of the enzyme leucine-aminopeptidase and bacterial secondary production in riparian groundwater. The “freshest” DOM was consistently found in the channel and during high-flow conditions in the groundwater of the restored riparian section and was spatially constrained to a zone inhabited by a dense willow population. The use of amino acid data and the newly developed DI for DOM in groundwater is a promising approach for characterizing the spatial and temporal dynamics of DOM reactivity and diagenesis within riparian groundwater.

Description: The responses of soil methane (CH4) net fluxes to nitrogen (N) addition in a N-fixing tree species (Acacia auriculiformis (AA)) and a non-N-fixing tree species (Eucalyptus citriodora (EU)) plantation were studied in southern China. Treatments were conducted at each plantation with three N levels (0, 50, and 100 kg N ha−1 yr−1 for control, medium-N, and high-N treatment, respectively, abbreviated as C, MN, and HN). From August 2010 to July 2011, CH4 flux was measured biweekly using a static chamber and gas chromatography technique. The soils of both sites acted as sink of atmospheric CH4. The CH4 uptake rate in control of the AA site (36.3 ± 3.2 μg CH4-C m−2 h−1) was greater than that of the EU plantation (29.9 ± 0.9 μg CH4-C m−2 h−1). In the AA plantation, the averaged rates of CH4 uptake for the MN (28.6 ± 2.3 μg CH4-C m−2 h−1) and HN treatment (23.8 ± 2.8 μg CH4-C m−2 h−1) were decreased by 21% and 35%, respectively, compared to the control. However, there was no change of soil CH4 uptake between N-treated plots and the controls in the EU site. Our results indicated that there might be large difference of inhibitive effect of N deposition on soil CH4 oxidation between the AA and EU plantations. The projected increase of N deposition would weaken the capability of N-fixing tree species plantations for atmospheric CH4 sink in tropical and subtropical regions.

Description: In several biomes, including croplands, wooded savannas, and tropical forests, many small fires occur each year that are well below the detection limit of the current generation of global burned area products derived from moderate resolution surface reflectance imagery. Although these fires often generate thermal anomalies that can be detected by satellites, their contributions to burned area and carbon fluxes have not been systematically quantified across different regions and continents. Here we developed a preliminary method for combining 1-km thermal anomalies (active fires) and 500 m burned area observations from the Moderate Resolution Imaging Spectroradiometer (MODIS) to estimate the influence of these fires. In our approach, we calculated the number of active fires inside and outside of 500 m burn scars derived from reflectance data. We estimated small fire burned area by computing the difference normalized burn ratio (dNBR) for these two sets of active fires and then combining these observations with other information. In a final step, we used the Global Fire Emissions Database version 3 (GFED3) biogeochemical model to estimate the impact of these fires on biomass burning emissions. We found that the spatial distribution of active fires and 500 m burned areas were in close agreement in ecosystems that experience large fires, including savannas across southern Africa and Australia and boreal forests in North America and Eurasia. In other areas, however, we observed many active fires outside of burned area perimeters. Fire radiative power was lower for this class of active fires. Small fires substantially increased burned area in several continental-scale regions, including Equatorial Asia (157%), Central America (143%), and Southeast Asia (90%) during 2001–2010. Globally, accounting for small fires increased total burned area by approximately by 35%, from 345 Mha/yr to 464 Mha/yr. A formal quantification of uncertainties was not possible, but sensitivity analyses of key model parameters caused estimates of global burned area increases from small fires to vary between 24% and 54%. Biomass burning carbon emissions increased by 35% at a global scale when small fires were included in GFED3, from 1.9 Pg C/yr to 2.5 Pg C/yr. The contribution of tropical forest fires to year-to-year variability in carbon fluxes increased because small fires amplified emissions from Central America, South America and Southeast Asia—regions where drought stress and burned area varied considerably from year to year in response to El Nino-Southern Oscillation and other climate modes.

Description: The Mackenzie River Delta is a lake-rich arctic floodplain that receives high inputs of dissolved organic matter (DOM) and suspended particulates from allochthonous and autochthonous sources, and may transfer carbon from dissolved to particulate phase via in situ formation of transparent exopolymer particles (TEP). TEP provides food for grazers, surfaces for bacteria, and increased potential for aggregation and sedimentation of organic matter. During open water 2006, we tracked TEP abundances in three Delta lakes representing gradients that include declining river-to-lake connection times, increasing levels of dissolved organic carbon (DOC), and declining chromophoric-DOM (CDOM). Unexpectedly, TEP abundances were highest immediately after the flood, when autochthonous autotrophic production was at a seasonal low and CDOM a seasonal high. Moreover, the lake with the strongest riverine influence and lowest levels of autochthonous autotrophic production had the highest mean TEP-carbon (TEP-C) concentrations among the lakes. The mean proportion of particulate organic carbon (POC) represented by TEP-C increased with increasing river connection time, and appears to represent a substantial proportion of POC in Mackenzie Delta Lakes. Unexpectedly, the TEP gradient was most strongly related to CDOM (river water source) rather than overall DOC. Variations in CDOM accounted for 53% of TEP-C variation among the lakes, indicating allochthonous matter was the most important source of TEP. DOC release from in situ macrophytes during periods of high photosynthesis may contribute to TEP formation in the lake with lowest riverine influence, but pH levels 〉9.5 driven by the high photosynthetic rates complicate the interpretation of results from this lake.

Description: At a semiarid steppe site located in the SE of Spain, relatively large CO2 emissions were measured that could not be attributed to the ecosystem activity alone. Since the study site was located in a tectonically active area, it was hypothesized that a part of the measured CO2 was of geologic origin. This investigation included a survey of soil CO2 efflux, together with carbon isotope analyses of the CO2 in the soil atmosphere, soil CO2 efflux (i.e., Keeling plots), groundwater and local thermal springs. These measurements confirmed the hypothesis of degassing from geologic sources. In areas with local faults and ancient volcanic structures, soil CO2 efflux rates were significantly higher (i.e., up to 6.3 and 1.4 μmol CO2 m−2 s−1) than measurements in a comparable site that was some distance from fault sites (means of 1.0 and 0.43 μmol CO2 m−2 s−1 in March and June, respectively). The CO2 concentration in the soil atmosphere at the eddy covariance site reached 0.14% v/v at 0.70 m soil depth with a 13C-enriched isotopic composition (δ13C from −10.2‰ to −16.6‰), consistent with the isotopic composition of the soil CO2 efflux estimated by Keeling plots (i.e., −16.6‰). 13C-enriched CO2 also occurred in local aquifers, and there was evidence of degassing from deep crust and mantle at regional scale by the helium isotopic ratio in spring waters located about 30 km (R/Ra: 0.12) and 200 km (R/Ra: 0.95) NW of the eddy covariance site. This study highlights the importance of considering CO2 sources of geologic origin when assessing the net ecosystem carbon balance of sites that may possibly be affected by circulation of such CO2-rich fluids.

Description: Rivers transport globally significant amounts of carbon (C) from terrestrial ecosystems to ocean margins. Understanding and quantifying the sources and respective contributions to riverine C has emerged as an important biogeochemical problem that can be approached through natural-abundance isotope mass balance. Traditionally, the sources of riverine C have been identified either qualitatively or quantitatively through application of static mixing models. However, both source signatures and contributions can vary significantly with time. Here we apply two time-based mixing models to a study of six rivers draining the northeast U.S. In the first model, a time-averaged mixing model (TAMM), we vary only the source isotopic (δ13C and Δ14C) signatures. In the second model, a time-varying mixing model (TVMM), we allow both isotopic signatures and contributions to vary with time. Based on results from the TVMM, drivers of variation in riverine particulate organic C (POC), dissolved organic C (DOC), and dissolved inorganic C (DIC) include stream discharge, stream discharge and water temperature, and water temperature and vegetation phenology, respectively. Major sources include C3 plant material, algal material and slow-turnover soil OC (“slow SOC”), which together account for 50%–100% (95% CI) of riverine POC; C3 plant material and slow SOC, which together typically account for 60%–100% (95% CI) of DOC; and atmospheric exchange which alone typically accounts for 40%–60% (95% CI) of DIC. Seasonal change in relative contributions from algal material, slow SOC, and photosynthesis (in response to the identified drivers) dominates the observed variation in POC, DOC and DIC, respectively. The TVMM is a novel tool to identify component contributions under more realistic non-static conditions, and with potential application to a broad range of biogeochemical studies.

Description: Hyporheic flow in streams has typically been studied separately from geomorphic processes. We investigated interactions between bed mobility and dynamic hyporheic storage of solutes and fine particles in a sand-bed stream before, during, and after a flood. A conservatively transported solute tracer (bromide) and a fine particles tracer (5 μm latex particles), a surrogate for fine particulate organic matter, were co-injected during base flow. The tracers were differentially stored, with fine particles penetrating more shallowly in hyporheic flow and retained more efficiently due to the high rate of particle filtration in bed sediment compared to solute. Tracer injections lasted 3.5 h after which we released a small flood from an upstream dam one hour later. Due to shallower storage in the bed, fine particles were rapidly entrained during the rising limb of the flood hydrograph. Rather than being flushed by the flood, we observed that solutes were stored longer due to expansion of hyporheic flow paths beneath the temporarily enlarged bedforms. Three important timescales determined the fate of solutes and fine particles: (1) flood duration, (2) relaxation time of flood-enlarged bedforms back to base flow dimensions, and (3) resulting adjustments and lag times of hyporheic flow. Recurrent transitions between these timescales explain why we observed a peak accumulation of natural particulate organic matter between 2 and 4 cm deep in the bed, i.e., below the scour layer of mobile bedforms but above the maximum depth of particle filtration in hyporheic flow paths. Thus, physical interactions between bed mobility and hyporheic transport influence how organic matter is stored in the bed and how long it is retained, which affects decomposition rate and metabolism of this southeastern Coastal Plain stream. In summary we found that dynamic interactions between hyporheic flow, bed mobility, and flow variation had strong but differential influences on base flow retention and flood mobilization of solutes and fine particulates. These hydrogeomorphic relationships have implications for microbial respiration of organic matter, carbon and nutrient cycling, and fate of contaminants in streams.

Description: The forest area in the western United States that burns annually is increasing with warmer temperatures, more frequent droughts, and higher fuel densities. Studies that examine fire effects for regional carbon balances have tended to either focus on individual fires as examples or adopt generalizations without considering how forest type, fire severity, and regional climate influence carbon legacies. This study provides a more detailed characterization of fire effects and quantifies the full carbon impacts in relation to direct emissions, slow release of fire-killed biomass, and net carbon uptake from forest regrowth. We find important variations in fire-induced mortality and combustion across carbon pools (leaf, live wood, dead wood, litter, and duff) and across low- to high-severity classes. This corresponds to fire-induced direct emissions from 1984 to 2008 averaging 4 TgC yr−1 and biomass killed averaging 10.5 TgC yr−1, with average burn area of 2723 km2 yr−1 across the western United States. These direct emission and biomass killed rates were 1.4 and 3.7 times higher, respectively, for high-severity fires than those for low-severity fires. The results show that forest regrowth varies greatly by forest type and with severity and that these factors impose a sustained carbon uptake legacy. The western U.S. fires between 1984 and 2008 imposed a net source of 12.3 TgC yr−1 in 2008, accounting for both direct fire emissions (9.5 TgC yr−1) and heterotrophic decomposition of fire-killed biomass (6.1 TgC yr−1) as well as contemporary regrowth sinks (3.3 TgC yr−1). A sizeable trend exists toward increasing emissions as a larger area burns annually.

Description: From September 2009 to August 2010, intensive monthly sampling of nutrients was conducted at two stations at the mouth of the Changjiang (Yangtze River). Particulate organic carbon (POC), particulate nitrogen (PN), and their stable isotope values (δ13C and δ15N) were also measured in selected samples of all months. Most nutrients (nitrate, phosphate, ammonia, and nitrite) as well as POC, PN, and δ13C displayed peak values when the highest or lowest Changjiang monthly discharges occurred, suggesting the Changjiang discharges strongly influence the seasonal variations of these chemicals. The sharply increases in concentrations of ammonia and nitrite in winter probably suggest nitrification was greatly depressed during this cold period. Using five interpolation methods, the annual discharge fluxes of nutrients, POC, and PN from the Changjiang to the East China Sea shelf were calculated. Combining this nutrient data with data from previous studies, the seasonal Mann-Kendall test, in which the influence of seasonal variation was considered, suggests concentrations of nitrate and phosphate in the Changjiang have significantly increased during recent decades at rates of 2.2 μM yr−1 and 0.03 μM yr−1, respectively; no significant trend for silicate was noted. Decreased POC annual fluxes along with sharply decreased suspended particulate matter yields were also seen in recent years (1993–2010). However, no distinct changes of δ13C, δ15N, and the POC/PN ratio, which describe the particulate organic matter properties, were observed during this period.

Description: Hydrological exchange between a river and its floodplain plays a critical role in maintaining key ecosystem services like habitat formation, nutrient transformation, and flood attenuation. We studied the spatial and temporal patterns of river-floodplain exchange in the Kafue Flats, a 6500-km2, dam-impacted floodplain ecosystem in Zambia. In addition, we characterized the effects of floodplain runoff on river biogeochemistry and assessed dam-related changes in the hydrological regime. The basic flood pulse concept poorly describes conditions in the Kafue Flats. Instead, high resolution measurements of discharge and tracers (specific conductivity, δ18O-H2O) along 410 km of river revealed substantial spatial variations in both the magnitude and direction of river-floodplain exchange. During peak discharge, a river channel constriction, 230 km into the floodplain, diverted as much as 80% of the river's ∼700 m3 s−1 discharge into the floodplain. As a net result, 〉80% of the water exiting the Kafue Flats via the river, either passed through the floodplain or originated from precipitation on the floodplain. This floodplain-derived water had a strong impact on river water quality, resulting in a seasonally recurring sharp decline in dissolved oxygen levels to

Description: River channel dynamics over many decades provide a physical control on the age structure of floodplain vegetation as a river occupies and abandons locations. Floodplain reoccupation by a river, in particular, determines the interval of time during which vegetation can establish and mature. A general framework for analyzing floodplain reoccupation and a time series model are developed and applied to five alluvial rivers in the United States. Channel dynamics in these rivers demonstrate time-scale dependence with short-term oscillations in active channel area in response to floods and subsequent vegetation growth and progressive lateral movement that accounts for much of the cumulative area occupied by the rivers over decades. Rivers preferentially reoccupy locations recently abandoned causing a decreasing probability of reoccupation with time since abandonment. For a typical case, a river is 10 times more likely to reoccupy an area it abandoned in the past decade than it is to reoccupy an area it abandoned 30 yrs ago. The decreasing probability of reoccupation over time is consistent with observations of persistent stands of late seral stage floodplain forest. A power function provides a robust approach for estimating the cumulative area occupied by a river and the age structure of riparian forests resulting from a specific historical sequence of streamflow in comparison to either linear or exponential alternatives.

Description: Concentration versus runoff relationships can reveal how land use and watershed hydrology interactively regulate solute inputs to streams and downstream aquatic ecosystems. In six mountainous southern California coastal watersheds, consistent nitrate-runoff patterns within three broad land use classes exist: dilution in agricultural watersheds, invariance in urban watersheds, and enrichment in an undeveloped watershed. Locally weighted scatterplot smoothing (LOWESS) revealed these patterns in nitrate-runoff relationships. A hyperbolic equation reproduced these relationships, also identifying the most common nitrate concentration (i.e., nitrate mode) observed in stream water during periods of low runoff (i.e., base flow) and high runoff (i.e., stormflow). LOWESS and the hyperbolic equation were also used to reveal and reproduce electrical conductance-runoff relationships, which contrary to nitrate-runoff relationships, demonstrated uniform behavior (significant dilution; p 〈 0.0001) in all watersheds regardless of land use. Nitrate-electrical conductance plots revealed seasonal shifts in stormflow nitrate modes, indicating nitrate flushing behavior at the beginning of the winter wet season for all watersheds except the two with the highest agricultural land usage. Despite variation in land uses between watersheds, we found a consistent reduction in the variability of stormflow nitrate (∼12%) and electrical conductance (∼21%) modes relative to base flow modes indicating a common water and nitrate source during periods of high runoff. We propose the undeveloped, mountainous upland regions of the watersheds as this source, and suggest that this region plays an important role in determining watershed stream nitrate concentrations and nitrate flux to the Santa Barbara Channel (Pacific Ocean).

Description: Freshwater lakes can emit significant quantities of methane to the atmosphere by bubbling. The high spatial and temporal heterogeneity of ebullition, combined with a lack of high-resolution field measurements, has made it difficult to accurately estimate methane fluxes or determine the underlying mechanisms for bubble release. We use a high-temporal resolution data set of ebullitive fluxes from the eutrophic Upper Mystic Lake, Massachusetts to understand the triggers that lead to bubbling from submerged sediments. A wavelet approach is introduced to detect ebullition events for multiple time-scales, and is complemented with traditional statistical methods for data analyses. We show that bubble release from lake sediments occurred synchronously at several sites, and was closely associated with small, aperiodic drops in total hydrostatic pressure. Such results are essential to constrain mechanistic models and to design future measurement schemes, particularly with respect to the temporal scales that are needed to accurately observe and quantify ebullition in aquatic ecosystems.

Description: Thermokarst lakes and peat-accumulating drained lake basins cover a substantial portion of Arctic lowland landscapes, yet the role of thermokarst lake drainage and ensuing peat formation in landscape-scale carbon (C) budgets remains understudied. Here we use measurements of terrestrial peat thickness, bulk density, organic matter content, and basal radiocarbon age from permafrost cores, soil pits, and exposures in vegetated, drained lake basins to characterize regional lake drainage chronology, C accumulation rates, and the role of thermokarst-lake cycling in carbon dynamics throughout the Holocene on the northern Seward Peninsula, Alaska. Most detectable lake drainage events occurred within the last 4,000 years with the highest drainage frequency during the medieval climate anomaly. Peat accumulation rates were highest in young (50–500 years) drained lake basins (35.2 g C m−2 yr−1) and decreased exponentially with time since drainage to 9 g C m−2 yr−1 in the oldest basins. Spatial analyses of terrestrial peat depth, basal peat radiocarbon ages, basin geomorphology, and satellite-derived land surface properties (Normalized Difference Vegetation Index (NDVI); Minimum Noise Fraction (MNF)) from Landsat satellite data revealed significant relationships between peat thickness and mean basin NDVI or MNF. By upscaling observed relationships, we infer that drained thermokarst lake basins, covering 391 km2 (76%) of the 515 km2 study region, store 6.4–6.6 Tg organic C in drained lake basin terrestrial peat. Peat accumulation in drained lake basins likely serves to offset greenhouse gas release from thermokarst-impacted landscapes and should be incorporated in landscape-scale C budgets.

Description: Land use/land cover change often leads to increased nutrient loading to streams; however, its influence on stream ecosystem nutrient transport remains poorly understood. Given the deleterious impacts elevated nutrient loading can have on aquatic ecosystems, it is imperative to improve understanding of nutrient retention capacities across stream scales and watershed development gradients. We performed 17 nutrient addition experiments on six streams across the West Fork Gallatin Watershed, Montana, USA, to quantify nitrogen uptake kinetics and retention dynamics across stream sizes (first to fourth order) and along a watershed development gradient. We observed that stream nitrogen (N) uptake kinetics and spiraling parameters varied across streams of different development intensity and scale. In more developed watersheds we observed a fertilization affect. This fertilization affect was evident as increased ash-free dry mass, chlorophyll a, and ambient and maximum uptake rates in developed as compared to undeveloped streams. Ash-free dry mass, chlorophyll a, and the number of structures in a subwatershed were significantly correlated to nutrient spiraling and kinetic parameters, while ambient and average annual N concentrations were not. Additionally, increased maximum uptake capacities in developed streams contributed to low in-stream nutrient concentrations during the growing season, and helped maintain watershed export at low levels during base flow. Our results indicate that land use/land cover change can enhance in-stream uptake of limiting nutrients and highlight the need for improved understanding of the watershed dynamics that control nutrient export across scales and development intensities for mitigation and protection of aquatic ecosystems.

Description: Tropical and subtropical wetlands are considered to be globally important sources of greenhouse gases, but their capacity to store carbon is presumably limited by warm soil temperatures and high rates of decomposition. Unfortunately, these assumptions can be difficult to test across long timescales because the chronology, cumulative mass, and completeness of a sedimentary profile are often difficult to establish. We therefore made a detailed analysis of a core from the principal drainage outlet of the Everglades of South Florida in order to assess these problems and determine the factors that could govern carbon accumulation in this large subtropical wetland. Accelerator mass spectroscopy dating provided direct evidence for both hard-water and open-system sources of dating errors, whereas cumulative mass varied depending upon the type of method used. Radiocarbon dates of gastropod shells, nevertheless, seemed to provide a reliable chronology for this core once the hard-water error was quantified and subtracted. Long-term accumulation rates were then calculated to be 12.1 g m−2 yr−1 for carbon, which is less than half the average rate reported for northern and tropical peatlands. Moreover, accumulation rates remained slow and relatively steady for both organic and inorganic strata, and the slow rate of sediment accretion (0.2 mm yr−1) tracked the correspondingly slow rise in sea level (0.35 mm yr−1) reported for South Florida over the past 4000 years. These results suggest that sea level and the local geologic setting may impose long-term constraints on rates of sediment and carbon accumulation in the Everglades and other wetlands.

Description: Climate change in semi-arid, midlatitude mountain environments is expected to shift the spatial patterns of temperature, water availability, and vegetation upslope. Vegetation growing near its low-elevation range limit may prove especially vulnerable to mortality and decline. We investigated the altitudinal pattern of conifer mortality that occurred from 2002 to 2004 in Southern California's San Jacinto Mountains. We found that conifer mortality was focused in the lower portion of the midmontane conifer range, which drove the midmontane conifer distribution upslope. We investigated past reports of conifer mortality in Southern California by searching historical newspaper accounts. We found evidence of previous episodes of conifer mortality that coincided with past droughts, and which may have caused vegetation redistribution in the past. We interpret the early 2000s mortality and associated vegetation redistribution as a response to natural decadal to centennial climate variability. Moreover, we hypothesize this response mode will dominate the early impact of global climate change on semi-arid forest, which, in turn, may complicate efforts to distinguish between ecological changes attributable to natural climate variability and those attributable to global climate change.

Description: Most natural freshwater lakes are net greenhouse gas (GHG) emitters. Compared to natural systems, human perturbations such as watershed wood harvesting and long-term reservoir impoundment lead to profound alterations of biogeochemical processes involved in the aquatic cycle of carbon (C). We exploited these anthropogenic alterations to describe the C dynamics in five lakes and two reservoirs from the boreal forest through the analysis of dissolved carbon dioxide (CO2), methane (CH4), oxygen (O2), and organic carbon (DOC), as well as total nitrogen and phosphorus. Dissolved and particulate organic matter, forest soil/litter and leachates, as well as dissolved inorganic carbon were analyzed for elemental and stable isotopic compositions (atomic C:N ratios, δ13Corg, δ13Cinorg and δ15Ntot). We found links between the export of terrestrial organic matter (OM) to these systems and the dissolved CO2 and O2 concentrations in the water column, as well as CO2 fluxes to the atmosphere. All systems were GHG emitters, with greater emissions measured for systems with larger inputs of terrestrial OM. The differences in CO2 concentrations and fluxes appear controlled by bacterial activity in the water column and the sediment. Although we clearly observed differences in the aquatic C cycle between natural and perturbed systems, more work on a larger number of water bodies and encompassing all four seasons should be undertaken to better understand the controls, rates, and spatial as well as temporal variability of GHG emissions, and to make quantitatively meaningful comparisons of GHG emissions (and other key variables) from natural and perturbed systems.

Description: The effect of air pollution on vegetation and the consequent changes in atmospheric chemistry are largely under-investigated; a new generation of chemical transport models fully coupled with complex land surface models is needed to better represent the feedbacks between vegetation and atmospheric chemistry. In this context, we coupled at high spatial resolution (30 km) the chemistry transport model CHIMERE with the land surface model ORCHIDEE to study the regional impact of tropospheric ozone on Euro-Mediterranean vegetation and the consequent changes in biogenic emission and ozone dry deposition owing to modifications in canopy conductance and LAI due to the ozone stress on vegetation. Results for the year 2002 show that the effect of tropospheric ozone on vegetation leads to a significant reduction of about 23% in the annual gross primary production, followed by a reduction in leaf area index. In addition, results show that CHIMERE does not correctly reproduce the activity of evergreen forests, grassland and crops during winter and fall, and consequently the dry deposition velocity is affected by this wrong pattern. On the other hand, in the coupled model, we have a better representation of vegetation activity during cold months, and the general performance of the model is improved compared to local site observations.

Description: In terrestrial systems limited by water availability the spatial distribution of vegetation can self-organize into a mosaic of vegetated patches and bare soil. Spatially extensive competition for water and short-range facilitation underpin many models that describe the process of vegetation pattern formation. Earlier studies investigating this self-organized patchiness have largely considered smooth landscapes. However, topographic variations can significantly alter the redistribution of surface water flow and therefore the pattern-forming process. Here, we consider how microtopographic variations, at the scale of individual plants, alters self-organized vegetation patterns with the use of a simple ecohydrological model. We show that increasing microtopography can induce a change from banded vegetation, oriented across the slope, to irregular drainage patterns, oriented in the downslope direction. The mechanism responsible is shown to be a change in the spatial redistribution of infiltration around plants and plant patches. Only small increases in microtopography are required to cause banded systems with weak facilitation to change to downslope-oriented patterns. When non-periodic boundary conditions were considered, band orientation tended to become oblique to the topographic contour and in some circumstances their migration upslope ceased. These results suggest that diffusive sediment transport processes may be essential for the maintenance of regular periodic vegetation patterns, which implies that erosion may be critical for understanding the susceptibility of these ecosystems to catastrophic shifts.

Description: The use of smart tracers to study hydrologic systems is becoming more widespread. Smart tracers are compounds that irreversibly react in the presence of a process or condition under investigation. Resazurin (Raz) is a smart tracer that undergoes an irreversible reduction to resorufin (Rru) in the presence of cellular metabolic activity. We quantified the relationship between the transformation of Raz and aerobic bacterial respiration in pure culture experiments using two obligate aerobes and two facultative anaerobes, and in colonized surface and shallow (

Description: Methane (CH4) is an important anthropogenic greenhouse gas, up to 15% of which is consumed by terrestrial soils. In this field study of the CH4 cycle of a pine forest, 18 plots were established at each of two sites, located 40 m apart. The upper site was well-drained and the lower site was poorly drained, but they shared similar overstory tree composition. Nitrogen was added as NH4NO3 incrementally across the 2009 growing season in a high (67 kg NH4NO3 ha−1 yr−1) and a low (5 kg NH4NO3 ha−1 yr−1) concentration. The sites were monitored for soil and air temperature, soil moisture, precipitation, air pressure, and NH4 and NO2+NO3 concentrations throughout the growing season. Across all treatments for the duration of the field season, average CH4 flux showed consumption of −0.84 kg CH4 ha−1 yr−1, but CH4 flux differed between the upper and lower sites. Across all treatments, upper site CH4 flux averaged −5.38 kg CH4 ha−1 yr−1, while lower site flux averaged 3.72 kg CH4 ha−1 yr−1, with greater variance than was observed at the upper site. High N treatments caused greater CH4 release than the control in the lower, but not the upper, site. The main correlated variable with CH4 flux was soil moisture; however, it accounted for

Description: The role of aquatic ecosystems in regional and global carbon cycles is becoming increasingly apparent, and lakes and reservoirs may be particularly important to the retention and processing of organic carbon. If this is the case, then lakes and reservoirs may act as control points that decrease OC concentrations and fluxes in downstream aquatic ecosystems. We tested this hypothesis at a regional scale by comparing dissolved organic carbon (DOC) concentrations and fluxes in 52 randomly selected streams and rivers with and without upstream lakes in the water-rich Northern Highlands Lake District (NHLD), Wisconsin, USA. DOC concentrations were significantly higher (p 〈 0.01) in drainage networks that did not contain lakes (25.02 mg/L) than they were in networks with upstream lakes (10.38 mg/L). However, when accounting for differences in wetland extent between watersheds, we were unable to detect a lake effect on downstream DOC concentrations (p 〉 0.49). Likewise, there were no significant differences in DOC:DON or DOC:DOP ratios, or in yields from watersheds with and without upstream lakes after compensating for wetland influences. We suggest that lake OC storage or processing may be limited by high hydrologic flushing in lakes with stream outlets and overwhelmed by larger scale influences of landscape composition in the NHLD. Consequently, drainage lakes in carbon-rich regions like the NHLD may have limited influence on terrigenous carbon exports to the ocean.

Description: Although land-water carbon (C) transport represents a critical link in the global C cycle, rare attempts have been made to compare hydrologic controls over storm pulses of dissolved organic C (DOC) and particulate organic C (POC) in mountainous watersheds. An immersible UV/Vis spectrophotometer was used to comparatively investigate the rapid storm responses of stream water DOC and POC in a small mountainous forested watershed in South Korea. High-frequency measurements at 5-min intervals during 42 hydrologic events, including monsoon storms and winter snowmelts, showed consistent patterns: POC concentrations were lower than DOC concentrations during base flow and small storm events but exceeded them during the peak flow periods of intense storm events. Although both the DOC and POC concentrations had hysteretic relationships with discharge, the POC concentrations showed larger increases and variations after crossing a threshold discharge on the rising limb of the storm hydrograph. Stronger responses to intense storms resulted in a disproportionately large export of POC at high flow, whereas a large portion of the total DOC flux was exported under prevailing low-flow conditions. The results demonstrate the potential of in situ optical measurements for investigating fine-resolution dynamics of the DOC and POC export during storm events. Stronger storm responses of the POC export compared to the limited response range of the DOC export suggest that erosion-induced POC export will become more important as a major pathway for the hydrologic soil C loss from mountainous watersheds in response to an increasing occurrence of extreme storm events.

Description: Dryland vegetation frequently shows self-organized spatial patterns as mosaic-like structures of sources (bare areas) and sinks (vegetation patches) of water runoff and sediments with variable interconnection. Good examples are banded landscapes displayed by Mulga in semiarid Australia, where the spatial organization of vegetation optimizes the redistribution and use of water (and other scarce resources) at the landscape scale. Disturbances can disrupt the spatial distribution of vegetation causing a substantial loss of water by increasing landscape hydrological connectivity and consequently, affecting ecosystem function (e.g., decreasing the rainfall-use efficiency of the landscape). We analyze (i) connectivity trends obtained from coupled analysis of remotely sensed vegetation patterns and terrain elevations in several Mulga landscapes subjected to different levels of disturbance, and (ii) the rainfall-use efficiency of these landscapes, exploring the relationship between rainfall and remotely sensed Normalized Difference Vegetation Index. Our analyses indicate that small reductions in the fractional cover of vegetation near a particular threshold can cause abrupt changes in ecosystem function, driven by large nonlinear increases in the length of the connected flowpaths. In addition, simulations with simple vegetation-thinning algorithms show that these nonlinear changes are especially sensitive to the type of disturbance, suggesting that the amount of alterations that an ecosystem can absorb and still remain functional largely depends on disturbance type. In fact, selective thinning of the vegetation patches from their edges can cause a higher impact on the landscape hydrological connectivity than spatially random disturbances. These results highlight surface connectivity patterns as practical indicators for monitoring landscape health.

Description: Accurately simulating gross primary productivity (GPP) in terrestrial ecosystem models is critical because errors in simulated GPP propagate through the model to introduce additional errors in simulated biomass and other fluxes. We evaluated simulated, daily average GPP from 26 models against estimated GPP at 39 eddy covariance flux tower sites across the United States and Canada. None of the models in this study match estimated GPP within observed uncertainty. On average, models overestimate GPP in winter, spring, and fall, and underestimate GPP in summer. Models overpredicted GPP under dry conditions and for temperatures below 0°C. Improvements in simulated soil moisture and ecosystem response to drought or humidity stress will improve simulated GPP under dry conditions. Adding a low-temperature response to shut down GPP for temperatures below 0°C will reduce the positive bias in winter, spring, and fall and improve simulated phenology. The negative bias in summer and poor overall performance resulted from mismatches between simulated and observed light use efficiency (LUE). Improving simulated GPP requires better leaf-to-canopy scaling and better values of model parameters that control the maximum potential GPP, such as εmax (LUE), Vcmax (unstressed Rubisco catalytic capacity) or Jmax (the maximum electron transport rate).

Description: Changes of vegetation phenology in response to climate change in the temperate forests have been well documented recently and have important implications on the regional and global carbon and water cycles. Predicting the impact of changing phenology on terrestrial ecosystems requires an accurate phenology model. Although species-level phenology models have been tested using a small number of vegetation species, they are rarely examined at the regional level. In this study, we used remotely sensed phenology and meteorological data to parameterize the species-level phenology models. We used a remotely sensed vegetation index (Two-band Enhanced Vegetation Index, EVI2) derived from the Moderate Resolution Spectroradiometer (MODIS) 8-day reflectance product from 2000 to 2010 of New England, United States to calculate remotely sensed vegetation phenology (start/end of season, or SOS/EOS). The SOS/EOS and the daily mean air temperature data from weather stations were used to parameterize three budburst models and one senescence model. We compared the relative strengths of the models to predict vegetation phenology and selected the best model to reconstruct the “landscape phenology” in New England from year 1960 to 2010. Of the three budburst models tested, the spring warming model showed the best performance with an averaged Root Mean Square Deviation (RMSD) of 4.59 days. The Akaike Information Criterion supported the spring warming model in all the weather stations. For senescence modeling, the Delpierre model was better than a null model (the averaged phenology of each weather station, averaged model efficiency = 0.33) and has a RMSD of 8.05 days. A retrospective analysis using the spring warming model suggests a statistically significant advance of SOS in New England from 1960 to 2010 averaged as 0.143 days per year (p = 0.015). EOS calculated using the Delpierre model and growing season length showed no statistically significant advance or delay between 1960 and 2010 in this region. These results suggest the applicability of species-level phenology models at the regional level (and potentially terrestrial biosphere models) and the feasibility of using these models in reconstructing and predicting vegetation phenology.

Description: Using ground-penetrating radar (GPR) to map subsurface patterns in peat physical properties, we investigated the developmental history of meso-scale surface patterning of microforms within a raised bog. Common offset GPR measurements were obtained along a 45-m transect, at frequencies ranging from 100 to 900 MHz. We found that low-frequency (central frequency 〈 240 MHz) GPR could not adequately represent the subsurface structures of the peatland because individual peat layers were too thin. However, more detailed high-frequency measurements (central frequency ≥ 240 MHz) showed a striking pattern of subsurface reflections that dip consistently in a northerly direction. The angle of these dipping reflectors is calculated using a semblance algorithm and was shown to average 3.9° between a depth of 1.0 and 2.5 m. These dipping reflectors may indicate downslope migration of surface microforms during the development of the peatland. Based on the estimated angle and the rate of peat accumulation, the average rate of downslope propagation of these surface microforms is calculated at 9.8 mm per year. Further survey work is required to establish whether the downslope migration is common across the peatland.

Description: The Orinoco River is the fourth largest in the world in terms of water discharge and organic carbon export to the ocean. River export of organic carbon is a key component of the carbon cycle and the global carbon budget. Here, we examined the seasonal transport of organic carbon by the Orinoco River into the eastern Caribbean using the conservative relationship of colored dissolved organic matter (CDOM) and dissolved organic carbon (DOC) in low salinity coastal waters influenced by river plumes. In situ measurements of CDOM absorption, DOC, and salinity were used to develop an empirical model for DOC concentration at the Orinoco River Plume. Satellite remote sensing reflectances were used with empirical models to determine DOC and Particulate organic carbon (POC) river transport. Our estimates of CDOM and DOC significantly correlated with in situ measurements and were within the expected ranges for the river. Total organic carbon transport by the Orinoco River during the period of 1998 to 2010 was 7.10 ×1012 g C y−1, from 5.29 × 1012 g C y−1 of DOC and 1.81 × 1012 g C y−1 of POC, representing ∼6% increase to previous published estimates. The variability in organic carbon transport responded to the seasonality in river flow more than to changes in organic carbon concentration in the river. Our results corroborate that is possible to estimate organic carbon transport using ocean color data at global scales. This is needed to reduce the uncertainties of land–ocean carbon fluxes.

Description: During November–December 2009 community rates of gross photosynthesis (Pg), respiration (R) and net calcification (Gnet) were estimated from low-tide slack water measurements of dissolved oxygen, dissolved inorganic carbon and total alkalinity at the historical station DK13 One Tree Island reef, Great Barrier Reef, Australia. Compared to measurements made during the 1960s–1970s at DK13 in the same season, Pg increased from 833 to 914 mmol O2·m−2·d−1 and Pg:R increased from 1.14 to 1.30, indicating that the reef has become more autotrophic. In contrast, Gnet decreased from 133 mmol C·m−2·d−1 to 74 ± 24 mmol C·m−2·d−1. This decrease stems primarily from the threefold increase in nighttime CaCO3 dissolution from −2.5 mmol·m−2·h−1 to −7.5 mmol·m−2·h−1. Comparison of the benthic community survey results from DK13 and its vicinity conducted during this study and in studies from the 1970s, 1980s and 1990s suggest that there have been no significant changes in the live coral coverage during the past 40 years. The reduced Gnet most likely reflects the almost threefold increase in dissolution rates, possibly resulting from increased bioerosion due to changes in the biota (e.g., sea cucumbers, boring organisms) and/or from greater chemical dissolution produced by changing abiotic conditions over the past 40 years associated with climate change, such as increased temperatures and ocean acidification. However, at this stage of research on One Tree Island the effects of these changes are not entirely understood.

Description: The roots of many trees in temperate and boreal forests are sheathed with ectomycorrhizal fungi (EMF) that extend into the soil, forming intimate contact with soil minerals, from which they absorb nutrient elements required by the plants and, in return, are supported by the organic carbon photosynthesized by the trees. While EMF are strongly implicated in mineral weathering, their effects on mineral surfaces at the nanoscale are less documented. In the present study, we investigated the effects of symbiotic EMF on the topography of a chlorite mineral using atomic force microscopy. A cleaning protocol was successfully applied to remove fungal hyphae without altering the underlying mineral structure and topography. Examination of the exposed chlorite surface showed the presence of primary channels, of the order of a micron in width and up to 50 nm in depth, the morphology of which strongly indicates a fungal-induced origin. Smaller secondary channels were observed extending from the primary channels and would appear to be involved in their enlargement. The presence of channels is the first nanoscale demonstration of the effects of fungal interaction, fuelled by plant photosynthate, on the topography of a chlorite mineral, and it provides clear evidence of the ability of EMF to enhance mineral dissolution.

Description: In the North Pacific Subtropical Gyre (NPSG), the regular occurrence of summer phytoplankton blooms contributes to marine ecosystem productivity and the annual carbon export. The mechanisms underlying the formation, maintenance, and decay of these blooms remain largely unknown; nitrogen fixation, episodic vertical mixing of nutrients, and meso- (

Description: In a suburban neighborhood of Minneapolis–Saint Paul, Minnesota, USA, we simultaneously measured net CO2 exchange of trees using sap flow and leaf gas exchange measurements, net CO2 exchange of a turfgrass lawn using eddy covariance from a portable tower, and total surface-atmosphere CO2 fluxes (FC) using an eddy covariance system on a tall tower. Two years of continuous measurements showed that net CO2 exchange varied among vegetation types, with the largest growing-season (Apr–Nov) net CO2 uptake on a per cover area basis from evergreen needleleaf trees (−603 g C m−2), followed by deciduous broadleaf trees (−216 g C m−2), irrigated turfgrass (−211 g C m−2), and non-irrigated turfgrass (−115 g C m−2). Vegetation types showed seasonal patterns of CO2 exchange similar to those observed in natural ecosystems. Scaled-up net CO2 exchange from vegetation and soils (FC(VegSoil)) agreed closely with landscape FC measurements from the tall tower at times when fossil fuel emissions were at a minimum. Although FC(VegSoil) did not offset fossil fuel emissions on an annual basis, the temporal pattern of FC(VegSoil) did significantly alter the seasonality of FC. Total growing season FC(VegSoil) in recreational land-use areas averaged −165 g C m−2 and was dominated by turfgrass CO2 exchange (representing 77% of the total), whereas FC(VegSoil) in residential areas averaged −124 g C m−2 and was dominated by trees (representing 78% of the total). Our results suggest urban vegetation types can capture much of the variability required to predict seasonal patterns and differences in FC(VegSoil) that could result from changes in land use or vegetation composition in temperate cities.

Description: Stream ecotones, specifically the lateral floodplain and subsurface hyporheic zone, can be important sites for nitrogen (N) removal via denitrification, but their role in streams with constructed floodplains has not been examined. We studied denitrification in the hyporheic zone and floodplains of an agriculturally influenced headwater stream in Indiana, USA, that had floodplains added as part of a “two-stage ditch” restoration project. To examine the potential for N removal in the hyporheic zone, we seasonally measured denitrification rates and nitrate concentrations by depth into the stream sediments. We found that nitrate concentration and denitrification rates declined with depth into the hyporheic zone, but denitrification was still measureable to a depth of at least 20 cm. We also measured denitrification rates on the restored floodplains over the course of a flood (pre, during, and post-inundation), and also compared denitrification rates between vegetated and non-vegetated areas of the floodplain. We found that floodplain denitrification rates increased over the course of a floodplain inundation event, and that the presence of surface water increased denitrification rates when vegetation was present. Stream ecotones in midwestern, agriculturally influenced streams have substantial potential for N removal via denitrification, particularly when they are hydrologically connected with high-nitrate surface water.

Description: We aimed to investigate the temporal variation of in-stream net dissolved inorganic nitrogen (DIN) areal uptake rates (UDIN, in μg N m−2 min−1) and its implications on regulating catchment N export, under base flow conditions. To do so, we estimated UDIN from longitudinal profiles of ambient DIN concentration (nitrate + ammonium) in two streams on a monthly basis during two hydrological years (n = 45). We found that in-stream DIN uptake and release did not offset each other (UDIN ≠ 0) in half of the dates, and that UDIN 〉 0 occurred mostly in autumn. Based on these reach-scale uptake rates, we performed empirical calculations and model simulations to assess the potential of stream network DIN retention to regulate DIN export from catchments on an annual scale. The empirical approach consisted in up-scaling UDIN by means of a dynamic stream network analysis that considered temporal and spatial variation of UDIN. The modeling approach consisted in applying different scenarios with the INCA model based on the natural range of empirical UDIN values. Our results showed that the contribution of stream network DIN retention to catchment DIN export increased when calculations accounted for the temporal variation of UDIN. Both approaches suggested that stream network DIN retention can significantly reduce DIN export from headwater catchments under base flow conditions (from 4% to 38%).

Description: Many estimates of nitrogen removal in streams and watersheds do not include or account for nitrate removal in deep sediments, particularly in gaining streams. We developed and tested a conceptual model for nitrate removal in deep sediments in a nitrogen-rich river network. The model predicts that oxic, nitrate-rich groundwater will become depleted in nitrate as groundwater upwelling through sediments encounters a zone that contains buried particulate organic carbon, which promotes redox conditions favorable for nitrate removal. We tested the model at eight sites in upwelling reaches of lotic ecosystems in the Waupaca River Watershed that varied by three orders of magnitude in groundwater nitrate concentration. We measured denitrification potential in sediment core sections to 30 cm and developed vertical nitrate profiles to a depth of about 1 m with peepers and piezometer nests. Denitrification potential was higher on average in shallower core sections. However, core sections deeper than 5 cm accounted for 70% on average of the depth-integrated denitrification potential. Denitrification potential increased linearly with groundwater nitrate concentration up to 2 mg NO3-N/L, but the relationship broke down at higher concentrations (〉5 mg NO3-N/L), a pattern that suggests nitrate saturation. At most sites groundwater nitrate declined from high concentrations at depth to much lower concentrations prior to discharge into the surface water. The profiles suggested that nitrate removal occurred at sediment depths between 20 and 40 cm. Dissolved oxygen concentrations were much higher in deep sediments than in pore water at 5 cm sediment depth at most locations. The substantial denitrification potential in deep sediments coupled with the declines in nitrate and dissolved oxygen concentrations in upwelling groundwater suggest that our conceptual model for nitrate removal in deep sediments is applicable to this river network. Our results suggest that nitrate removal rates can be high in deep sediments of upwelling stream reaches, which may have implications for efforts to understand and quantify nitrogen transport and removal at larger scales.

Description: Woody plant encroachment, a worldwide phenomenon, is a major driver of land degradation in desert grasslands. Woody plant encroachment by shrub functional types ultimately leads to the formation of a patchy landscape with fertile shrub patches interspaced with nutrient-depleted bare soil patches. This is considered to be an irreversible process of land and soil degradation. Recent studies have indicated that in the early stages of shrub encroachment, when there is sufficient herbaceous connectivity, fires (prescribed or natural) might provide some reversibility to the shrub encroachment process by negatively affecting shrub demography and homogenizing soil resources across patches within weeks to months after burning. A comprehensive understanding of longer term changes in microtopography and spatial patterning of soil properties following fire in shrub-encroached grasslands is desirable. Here, we investigate the changes in microtopography with LiDAR (light detection and ranging), vegetation recovery, and spatial pattering of soil properties in replicated burned, clipped, and control areas in a shrub-grass transition zone in the northern Chihuahuan Desert four years after prescribed fire or clipping. Results indicate a greater homogeneity in soil, microtopography, and vegetation patterning on burned relative to clipped and control treatments. Findings provide further evidence that disturbance by prescribed fire may allow for reversal of the shrub encroachment process, if the event occurs in the early stages of the vegetation shift. Improved understanding of longer-term effects of fire and associated changes in soil patterning can inform the use and role of fire in the context of changing disturbance regimes and climate.

Description: The Community Land Model version 4 overestimates gross primary production (GPP) compared with estimates from FLUXNET eddy covariance towers. The revised model of Bonan et al. (2011) is consistent with FLUXNET, but values for the leaf-level photosynthetic parameter Vcmax that yield realistic GPP at the canopy-scale are lower than observed in the global synthesis of Kattge et al. (2009), except for tropical broadleaf evergreen trees. We investigate this discrepancy between Vcmax and canopy fluxes. A multilayer model with explicit calculation of light absorption and photosynthesis for sunlit and shaded leaves at depths in the canopy gives insight to the scale mismatch between leaf and canopy. We evaluate the model with light-response curves at individual FLUXNET towers and with empirically upscaled annual GPP. Biases in the multilayer canopy with observed Vcmax are similar, or improved, compared with the standard two-leaf canopy and its low Vcmax, though the Amazon is an exception. The difference relates to light absorption by shaded leaves in the two-leaf canopy, and resulting higher photosynthesis when the canopy scaling parameter Kn is low, but observationally constrained. Larger Kn decreases shaded leaf photosynthesis and reduces the difference between the two-leaf and multilayer canopies. The low model Vcmax is diagnosed from nitrogen reduction of GPP in simulations with carbon-nitrogen biogeochemistry. Our results show that the imposed nitrogen reduction compensates for deficiency in the two-leaf canopy that produces high GPP. Leaf trait databases (Vcmax), within-canopy profiles of photosynthetic capacity (Kn), tower fluxes, and empirically upscaled fields provide important complementary information for model evaluation.

Description: Recent climate variability (increasing temperature, droughts) and atmospheric composition changes (nitrogen deposition, rising CO2 concentration) along with harvesting, wildfires, and insect infestations have had significant effects on U.S. forest carbon (C) uptake. In this study, we attribute C changes in the conterminous U.S. forests to disturbance and non-disturbance factors with the help of forest inventory data, a continental stand age map, and an updated Integrated Terrestrial Ecosystem Carbon Cycle model (InTEC). We grouped factors into disturbances (harvesting, fire, insect infestation) and non-disturbances (CO2 concentration, N deposition, and climate variability) and estimated their subsequent impacts on forest regrowth patterns. Results showed that on average, the C sink in the conterminous U.S. forests from 1950 to 2010 was 206 Tg C yr−1 with 87% (180 Tg C yr−1) of the sink in living biomass. Compared with the simulation of all factors combined, the estimated C sink would be reduced by 95 Tg C yr−1 if disturbance factors were omitted, and reduced by 50 Tg C yr−1 if non-disturbance factors were omitted. Our study also showed diverse regional patterns of C sinks related to the importance of driving factors. During 1980–2010, disturbance effects dominated the C changes in the South and Rocky Mountain regions, were almost equal to non-disturbance effects in the North region, and had minor effects compared with non-disturbance effects in the West Coast region.

Description: We demonstrate the effect of an ecosystem differentiated insulation by snow on the soil thermal regime and on the terrestrial soil carbon distribution in the pan-Arctic area. This is done by means of a sensitivity study performed with the land surface model ORCHIDEE, which furthermore provides a first quantification of this effect. Based on field campaigns reporting higher thermal conductivities and densities for the tundra snowpack than for taiga snow, two distributions of near-equilibrium soil carbon stocks are computed, one relying on uniform snow thermal properties and the other using ecosystem-differentiated snow thermal properties. Those modeled distributions strongly depend on soil temperature through decomposition processes. Considering higher insulation by snow in taiga areas induces warmer soil temperatures by up to 12 K in winter at 50 cm depth. This warmer soil signal persists over summer with a temperature difference of up to 4 K at 50 cm depth, especially in areas exhibiting a thick, enduring snow cover. These thermal changes have implications on the modeled soil carbon stocks, which are reduced by 8% in the pan-Arctic continental area when the vegetation-induced variations of snow thermal properties are accounted for. This is the result of diverse and spatially heterogeneous ecosystem processes: where higher soil temperatures lift nitrogen limitation on plant productivity, tree plant functional types thrive whereas light limitation and enhanced water stress are the new constrains on lower vegetation, resulting in a reduced net productivity at the pan-Arctic scale. Concomitantly, higher soil temperatures yield increased respiration rates (+22% over the study area) and result in reduced permafrost extent and deeper active layers which expose greater volumes of soil to microbial decomposition. The three effects combine to produce lower soil carbon stocks in the pan-Arctic terrestrial area. Our study highlights the role of snow in combination with vegetation in shaping the distribution of soil carbon and permafrost at high latitudes.

Description: Humified soil organic matter storage in boreal forests is large, and its responses to warming over relatively long timescales is critical for predicting soil feedbacks to climate change. To derive information relevant across decades to centuries from manipulative short-term experiments, we conducted incubations of soils from two forested sites along the Newfoundland-Labrador Boreal Ecosystem Latitude Transect in eastern Canada and assessed linkages between incubation data and these sites' profile characteristics. The sites differ in mean annual temperature by 3.4°C, but vegetation and soil types are similar. Organic soils (Oe + Oa) were incubated for 120 days at 15°C and 20°C, with and without a replaced Oi subhorizon possessing a distinct δ13C signature. Laboratory warming induced significantly greater mineralization and leaching of humified SOM relative to replaced Oi, congruent with greater warming-induced increases in phenol oxidase activity relative to enzymes associated with labile C acquisition (percent increases of 101% versus 50%, respectively). These data suggest that warming can influence microbial communities and their enzymatic dynamics such that relative losses of humified SOM are disproportionately enhanced. This is consistent with stable isotopic, C:N, and radiocarbon profile differences between the two sites, which suggest a greater degree of microbial processing and greater relative losses of older SOC over the preceding decades at the warmer site, given our knowledge of organic inputs in these soils. This study is a first step toward linking the divergent timescales represented by soil profiles and laboratory manipulations, an important goal for biogeochemists assessing climate change impacts on SOM dynamics.

Description: Much progress has been made in methane modeling for the Arctic. However, there is still large uncertainty in emissions estimates due to the spatial variability in water table depth resulting from complex topographic gradients, and due to variations in methane production and oxidation due to complex freezing and thawing processes. Here we extended an extant methane emission module within a biogeochemistry model, the Terrestrial Ecosystem Model (TEM), to include a large-scale hydrology model, the variable infiltration capacity (VIC) model. The VIC model provides the required inputs, including freezing and thawing fronts, soil temperature and moisture, to the methane module. The effect of topography on the redistribution of soil moisture and water table depth was explicitly modeled using the TOPMODEL approach. The coupled modeling framework was applied to the Yukon River basin at a spatial resolution of 1 km from 1986 to 2005. The simulations show that the average annual net emissions of CH4 from the region are 4.01 Tg CH4 yr−1. El Niño phenomena usually lead to positive emission anomalies, while decreases in net CH4 emissions may be associated with strong La Niña events. Precipitation was found to be more closely related to CH4 dynamics than to soil temperature and active layer depth during the study period. This study suggests that the effects of soil freezing and thawing processes and the effects of microtopography on hydrology should be considered in the quantification of CH4 emissions.

Description: We measured the δ13C of assimilated carbon (foliage organic matter (δCOM), soluble carbohydrates (δCSC), and waxes (δCW)) and respiratory carbon (foliage (δCFR), soil (δCSR) and ecosystem 13CO2 (δCER)) for two years at adjacent ecosystems in the southeastern U.S.: a regenerated 32 m tall mature Pinus palustris forest, and a mid-rotation 13 m tall Pinus elliottii stand. Carbon pools and foliage respiration in P. palustris were isotopically enriched by 2‰ relative to P. elliottii. Despite this enrichment, mean δCER values of the two sites were nearly identical. No temporal trends were apparent in δCSC, δCFR, δCSR and δCER. In contrast, δCOM and δCW at both sites declined by approximately 2‰ over the study. This appears to reflect the adjustment in the δ13C of carbon storage reserves used for biosynthesis as the trees recovered from a severe drought prior to our study. Unexpectedly, the rate of δ13C decrease in the secondary C32–36 n-alkanoic acid wax molecular cluster was twice that observed for δCOM and the predominant C22–26 compound cluster, and provides new evidence for parallel but separate wax chain elongation systems utilizing different carbon precursor pools in these species. δCFR and δCER were consistently enriched relative to assimilated carbon but, in contrast to previous studies, showed limited variations in response to changes in vapor pressure deficit (D). This limited variability in respiratory fluxes and δCSC may be due to the shallow water table as well as the deep taproots of pines, which limit fluctuations in photosynthetic discrimination arising from changes in D.

Description: Soil organic matter (SOM) processes in dynamic landscapes are strongly influenced by soil erosion and sedimentation. We determined the contribution of physical isolation of organic matter (OM) inside aggregates, chemical interaction of OM with soil minerals, and molecular structure of SOM in controlling storage and persistence of SOM in different types of eroding and depositional landform positions. By combining density fractionation with elemental and spectroscopic analyses, we showed that SOM in depositional settings is less transformed and better preserved than SOM in eroding landform positions. However, which environmental factors exert primary control on storage and persistence of SOM depended on the nature of the landform position considered. In an annual grassland watershed, protection of SOM by physical isolation inside aggregates and chemical association of organic matter (complexation) with soil minerals, as assessed by correlation with radiocarbon concentration, were more effective in the poorly drained, lowest-lying depositional landform positions, compared to well-drained landform positions in the upper parts of the watershed. Results of this study demonstrated that processes of soil erosion and deposition are important mechanisms of long-term OM stabilization.

Description: In this study, we analyzed the influence of soil mineral characteristics (e.g., clay concentration and mineralogical composition, iron and aluminum oxide concentration and crystallinity, specific surface area, and exchangeable cation concentration) on (i) organic carbon (OC) content (kg m−2) and (ii) the concentration (g kg−1), composition, and stability of the mineral-associated organic matter (OM) of arable and forest topsoils. We selected seven soil types with different mineral characteristics for this study. For each soil type, samples were taken from topsoils of a deciduous forest and an adjacent arable site. The arable and forest sites have been used continuously for more than 100 years. Na-pyrophosphate soluble OM fractions (OM(PY)), representing mineral-associated OM, were extracted, analyzed for OC and 14C concentrations, and characterized by FTIR spectroscopy. For the forest and arable topsoils, a linear relationship was found between the OC content and exchangeable Ca. For the arable topsoils (pH 6.7–7.5), correlation analyses indicated that the OCPY concentration increased with an increase in oxalate soluble Fe and Al, exchangeable Ca, and Na-pyrophosphate soluble Mg and Fe concentrations. The stability of OM(PY) determined by the 14C measurements of the near-neutral arable topsoils was shown to increase with the specific surface area and the concentration of exchangeable Ca. For the acidic forest topsoils (pH

Description: Disturbance processes of various types substantially modify ecosystem carbon dynamics both temporally and spatially, and constitute a fundamental part of larger landscape-level dynamics. Forests typically lose carbon for several years to several decades following severe disturbance, but our understanding of the duration and dynamics of post-disturbance forest carbon fluxes remains limited. Here we capitalize on a recent North American Carbon Program disturbance synthesis to discuss techniques and future work needed to better understand carbon dynamics after forest disturbance. Specifically, this paper addresses three topics: (1) the history, spatial distribution, and characteristics of different types of disturbance (in particular fire, insects, and harvest) in North America; (2) the integrated measurements and experimental designs required to quantify forest carbon dynamics in the years and decades after disturbance, as presented in a series of case studies; and (3) a synthesis of the greatest uncertainties spanning these studies, as well as the utility of multiple types of observations (independent but mutually constraining data) in understanding their dynamics. The case studies—in the southeast U.S., central boreal Canada, U.S. Rocky Mountains, and Pacific Northwest—explore how different measurements can be used to constrain and understand carbon dynamics in regrowing forests, with the most important measurements summarized for each disturbance type. We identify disturbance severity and history as key but highly uncertain factors driving post-disturbance carbon source-sink dynamics across all disturbance types. We suggest that imaginative, integrative analyses using multiple lines of evidence, increased measurement capabilities, shared models and online data sets, and innovative numerical algorithms hold promise for improved understanding and prediction of carbon dynamics in disturbance-prone forests.

Description: The benthic, mat-forming diatom Didymosphenia geminata has the unique ability to produce large amounts of algal biomass under oligotrophic conditions in cold, fast flowing streams and rivers. This presents an ecological paradox that challenges our current understanding of stream ecosystem dynamics. Our understanding of the drivers of D. geminata ecology is still limited. Here we present a conceptual model for the blooming behavior and persistence of this species to advance scientific understanding of strategies for life in fast flowing oligotrophic waters and support the design of future research and mitigation measures for nuisance algal blooms. The conceptual model is based on a synthesis of data and ideas from a range of disciplines including hydrology, geomorphology, biogeochemistry, and ecology. The conceptual model highlights the role of water chemistry, river morphology, and flow thresholds in defining the habitat window for D. geminata. We propose that bed disturbance is a primary control on accumulation and persistence of D. geminata and that the removal threshold can be determined by synthesizing site-specific information on hydrology and geomorphology. Further, we propose that a key to understanding the didymo paradox is the separation of cellular reproduction and mat morphology with specific controls acting in respect of the different processes.

Description: Transitions between aquatic and terrestrial ecosystems represent zones where soil moisture is a dominant factor influencing vegetation composition. Niche models based on hydrological and vegetation observations can be powerful tools for guiding management of these zones, especially when they are linked with physically based hydrological models. Floodplain restoration represents a unique opportunity to utilize such a predictive vegetation tool when a site's hydrology is altered to create a wetter environment. A variably saturated groundwater flow model was developed and used to simulate the soil moisture regime across a floodplain in Wisconsin where post-settlement alluvium was removed with the intent of increasing regionally threatened wetland plant species. Hydrological niche models based on simultaneous observations of vegetation composition and surface effective saturation were used to predict probability of presence for two plant species (Carex vulpinoidea (fox sedge) and Elymus canadensis (Canada wildrye)) and wetland indicator score (a composite indicator of relative frequency of species in five habitat categories) based on simulated surface effective saturation. The vegetation predictions following restoration are more wetland-species dominant overall. However, zones of the study site where a confining layer is present that decouples groundwater from the near-surface soil zone tend to be drier following restoration due to restricted upward groundwater flow and less soil water storage above the confining layer. As reflected by an increase in the interquartile range in the predicted wetland indicator score, this restoration technique may increase the site-scale spatial diversity of plant community types while simultaneously accomplishing the goal of increasing wetland plant species occurrence.

Description: The emission of biogenic gases from large rivers can be an important component of regional greenhouse gas budgets. However, emission rate estimates are often poorly constrained due to uncertainties in the air-water gas exchange rate. We used the floating chamber method to estimate the gas transfer velocity (k) of carbon dioxide (CO2), methane (CH4), and nitrous oxide (N2O) in the Markland Pool of the Ohio River, a large tributary of the Mississippi River (U.S.A). We measured k every two weeks for a year at one site and at 15 additional sites distributed across the length of the pool during two summer surveys. We found that k was positively related to both water currents and wind speeds, with 46% of the gas transfer attributable to water currents at low wind speeds (e.g., 0.5 m s−1) and 11% at higher wind speeds (e.g., 〉2.0 m s−1). Gas transfer velocity was highly sensitive to wind, possibly because the direction of river flow was often directly opposed to the wind direction. Gas transfer velocity values derived for CH4 were consistently greater than those derived for CO2 when standardized to a Schmidt number of 600 (k600), possibly because the transfer of CH4, a poorly soluble gas, was enhanced by surfacing microbubbles. Additional research to determine the conditions that support microbubble enhanced gas transfer is merited.

Description: We used the carbon isotope composition (14C and δ13C) to measure the source and age of DOC, POC, dissolved CO2 and CH4 (δ13C only) released from three natural peat pipes and the downstream catchment outlet of a small peatland in northern England. Sampling under different hydrological extremes (high flows associated with storm events and low flows before or after storms) was used to explore variability in C sources as flow paths change over short periods of time. The δ13C composition of organic C differed (δ13C-DOC −28.6‰ to −27.6‰; δ13C-POC −28.1‰ to −26.1‰) from that of the dissolved gases (δ13C-CO2 −20.5‰ to +1.1‰; δ13C-CH4 −67.7‰ to −42.0‰) and showed that C leaving the catchment was a mixture of shallow/deep pipe and non-pipe sources. The isotopic composition of the dissolved gases was more variable than DOC and POC, with individual pipes either showing 13C enrichment or depletion during a storm event. The 14C age of DOC was consistently modern at all sites; POC varied from modern to 653 years BP and evasion CO2 from modern to 996 years BP. Differences in the isotopic composition of evasion CO2 at pipe outlets do not explain the variability in δ13C and 14C at the catchment outlet and suggest that overland flow is likely to be an important source of CO2. Our results also show that the sources of CO2 and CH4 are significantly more variable and dynamic than DOC and POC and that natural pipes vent old, deep peat CO2 and POC (but not DOC) to the atmosphere.

Description: The 15N-enrichment of plants and soils is believed to indicate characteristics of the open nitrogen (N) cycle in terrestrial ecosystems because N lost from an ecosystem is presumably 15N-depleted through isotopic fractionation. However, because of a lack of an appropriate analytical methodology to confirm that supposition, the δ15N value for total dissolved nitrogen (TDN, the sum of ammonium, nitrate, and dissolved organic N) in stream water from forests has been measured only rarely. This report describes the δ15N values for TDN, ammonium, and nitrate in precipitation and stream water, together with those for soil-emitted nitrous oxide (N2O; measured once) in an N-saturated subtropical forest in southern China. Concentration-weighted δ15N values of TDN were −0.7‰ in precipitation and +1.2‰ in stream water. The difference in δ15N between soil (+3.9‰) and TDN in the stream water was 2.7‰. In contrast, soil-emitted N2O was strongly 15N-depleted (−14.3‰): 18‰ lower than that of the soil. Our results demonstrate that the discharged N loss is 15N-depleted only slightly compared with soil N, and gaseous N losses can be a strong driver for raising the terrestrial ecosystem δ15N. Our findings suggest that the relation between ecosystem δ15N and the open N cycle can be interpreted better by considering the net discrimination against 15N determined by the balance between gaseous and discharge N losses. Steady state 15N budget calculations proposed by Houlton and Bai (2009) can provide important information about the gaseous N fluxes, which are difficult to measure directly. The steady state calculation for the relationships among gaseous N loss, apparent isotopic fractionation during gaseous N loss, and isotopic signature of N inputs suggests that precise measurements of unmeasured components (e.g., dry deposition, NO and N2 emission) are quite important for better estimation of gaseous N losses from the ecosystem.

Description: This study aims to compare and validate two soil-vegetation-atmosphere-transfer (SVAT) schemes: TERRA-ML and the Community Land Model (CLM). Both SVAT schemes are run in standalone mode (decoupled from an atmospheric model) and forced with meteorological in-situ measurements obtained at several tropical African sites. Model performance is quantified by comparing simulated sensible and latent heat fluxes with eddy-covariance measurements. Our analysis indicates that the Community Land Model corresponds more closely to the micrometeorological observations, reflecting the advantages of the higher model complexity and physical realism. Deficiencies in TERRA-ML are addressed and its performance is improved: (1) adjusting input data (root depth) to region-specific values (tropical evergreen forest) resolves dry-season underestimation of evapotranspiration; (2) adjusting the leaf area index and albedo (depending on hard-coded model constants) resolves overestimations of both latent and sensible heat fluxes; and (3) an unrealistic flux partitioning caused by overestimated superficial water contents is reduced by adjusting the hydraulic conductivity parameterization. CLM is by default more versatile in its global application on different vegetation types and climates. On the other hand, with its lower degree of complexity, TERRA-ML is much less computationally demanding, which leads to faster calculation times in a coupled climate simulation.

Description: Wildfire represents the single largest disturbance to the ecohydrological function of northern peatlands. Alterations to peatland thermal behavior as a result of wildfire will modify the carbon balance of these important long-term global carbon stores and regulate post-fire ecosystem recovery. We simulate the 3-D thermal behavior of a peatland that has been disturbed by wildfire to identify how changes in peat temperatures emerge from changes to the surface energy balance and peat thermal properties. Peat temperatures are simulated within two adjacent peatlands, one area having burned 4 years previously, the second which has been wildfire-free for 75 years. We demonstrate that there is only a small alteration to the thermal response in Sphagnum fuscum hummocks that are not severely burnt within the wildfire. In contrast, wildfire produces important changes to the energy balance of Sphagnum hollows. A large reduction in the latent heat flux post-fire increases surface temperatures by up to 30°C during daytime summer conditions. However, temperatures through the peat profile are insensitive to these increases in surface temperature. The low surface moisture content of near-surface peat insulates the profile from these higher temperatures and, at depths below 0.015 m, only small differences are identifiable between burned and unburned hollow temperatures. Nevertheless, we argue that these alterations to near-surface temperatures and evaporation rates likely substantially influence the thermal and hydrological conditions post-wildfire, impacting the peatland recovery.

Description: The conversion of grasslands to shrublands by woody plant encroachment is a common occurrence in arid and semi-arid regions that can have significant effects on ecosystem hydrology and biogeochemistry, including changes in soil organic carbon stocks. The mechanisms determining the direction and magnitude of soil carbon change due to encroachment are unknown but appear to depend on factors that vary along a gradient of mean annual precipitation. We present a model of coupled steady-state soil moisture and carbon dynamics that accounts for the effects of structural heterogeneity at the scale of vegetation patches by representing the local effects of root systems and canopies. We applied the model to paired grasslands and woody encroachment-produced shrublands spanning a climate gradient in the American Southwest, and model results were consistent with measured values. Modeled transpiration, which was used to determine productivity, was 5%–57% greater in shrublands, with the greatest difference occurring where the relative difference in leaf area index between grass and shrub patches was small and vegetation density was low. Patterns of soil carbon abundance were mainly driven by patterns of productivity, but decomposition rates were also affected by vegetation structure. Woody encroachment increases heterogeneity in soil carbon decomposition rates and has the net effect of increasing soil carbon residence times by 3.6–4.9 years. Our model represents an important step toward a mechanistic and quantitative understanding of the effect of vegetation structure on soil moisture and carbon dynamics and highlights the need for a better description of belowground vegetation structure in ecosystem models.

Description: We determined microbial decomposition of dissolved organic carbon (DOC) over 3.7 year long dark bioassays of six Swedish lake waters. The overall lost DOC fraction was similar in clearwater lakes (34.8 ± 2.4%) and in brownwater lakes (37.8 ± 1.9%). Reactivity continuum modeling revealed that the most labile DOC fraction, degrading at rates 〉0.01 d−1, was larger in the clearwater lakes (11.1 ± 1.2%) than in the brownwater lakes (0.8 ± 0.1%). The initial apparent first-order decay coefficient k was fivefold larger in the clearwater lakes (0.0043 ± 0.0012 d−1) than in the brownwater lakes (0.0009 ± 0.0003 d−1). Over time, k decreased more steeply in the clearwater lakes than in the brownwater lakes, reaching the k of the brownwater lakes within 5 months. Finally, k averaged 0.0001 d−1 in both lake categories. In the brownwater lakes, colored dissolved organic matter (CDOM) absorption decayed with an initial k twice as large (0.0018 ± 0.0008 d−1) as that of DOC. The initial k was inversely correlated with initial specific UV absorption and CDOM absorption and positively correlated with initial tryptophan-like fluorescence as proxy for autochthonous DOC. Exposure to simulated sunlight at the end of the incubations caused loss of color in the clearwater lakes and loss of DOC in the brownwater lakes, where subsequent mineralization was also stimulated. The DOC lost in the absence of photochemical processes fell below previously reported watershed-scale losses in Sweden by 25% at most. This suggests that a major part of the in situ DOC loss could potentially be attributed to dark reactions alone.

Description: Despite the wide acceptance of the “big-leaf” upscaling strategy in evapotranspiration modeling (e.g., the Penman-Monteith model), its usefulness in simulating canopy photosynthesis may be limited by the underlying assumption of homogeneous response of carbon assimilation light-response kinetics through the canopy. While previous studies have shown that the separation of the canopy into sunlit and shaded parts (i.e., two-leaf model) is typically more effective than big-leaf models for upscaling photosynthesis from leaf to canopy, a systematic comparison between these two upscaling strategies among multiple ecosystems has not been presented. In this study, gross primary productivity was modeled using two-leaf and big-leaf upscaling approaches in the Boreal Ecosystem Productivity Simulator for shrublands, broadleaf, and conifer forest types. When given the same leaf-level photosynthetic parameters, the big-leaf approach significantly underestimated canopy-level GPP while the two-leaf approach more closely predicted both the magnitude and day-to-day variability in eddy covariance measurements. The underestimation by the big-leaf approach is mostly caused by its exclusion of the photosynthetic contributions of shaded leaves. Tests of the model sensitivity to a foliage clumping index revealed that the contribution of shaded leaves to the total simulated productivity can be as high as 70% for highly clumped stands and seldom decreases below ∼40% for less-clumped canopies. Our results indicate that accurate upscaling of photosynthesis across a broad array of ecosystems requires an accurate description of canopy structure in ecosystem models.

Description: Watershed export of phosphorus (P) from anthropogenic sources has contributed to eutrophication in freshwater and coastal ecosystems. We explore impacts of watershed urbanization on the magnitude and export flow distribution of P along an urban-rural gradient in eight watersheds monitored as part of the Baltimore Ecosystem Study Long-Term Ecological Research site. Exports of soluble reactive phosphorus (SRP) and total P (TP) were lowest in small watersheds with forest and low-density residential land use (2.8–3.1 kg−1 km−2 yr−1). In contrast, SRP and TP exports increased with watershed impervious surface coverage and reached highest values in a small urban watershed (24.5–83.7 kg−1 km−2 yr−1). Along the Gwynns Falls, a larger watershed with mixed land use, the greatest proportion of SRP (68%) and TP (75%) was contributed from the lower watershed, where urban areas were the dominant land use. Load duration curve analysis showed that increasing urbanization in watersheds was associated with shifts in P export to high-flow conditions (〉2 mm d−1). SRP concentrations during low-flow conditions at urban headwater sites were highest during summer and lowest during winter. This seasonal pattern was consistent with sediment incubation experiments showing that SRP release from sediments was temperature dependent. Our results suggest that shifts in streamflow and alterations in water temperatures owing to urbanization and climate can influence stream water P concentrations and P export from urban watersheds.

Description: In large parts of the Southern Ocean, primary production is limited due to shortage of iron (Fe). We measured vertical Fe profiles in the western Weddell Sea, Weddell-Scotia Confluence, and Antarctic Circumpolar Current (ACC), showing that Fe is derived from benthic Fe diffusion and sediment resuspension in areas characterized by high turbulence due to rugged bottom topography. Our data together with literature data reveal an exponential decrease of dissolved Fe (DFe) concentrations with increasing distance from the continental shelves of the Antarctic Peninsula and the western Weddell Sea. This decrease can be observed 3500 km eastward of the Antarctic Peninsula area, downstream the ACC. We estimated DFe summer fluxes into the upper mixed layer of the Atlantic sector of the Southern Ocean and found that horizontal advection dominates DFe supply, representing 54 ± 15% of the total flux, with significant vertical advection second most important at 29 ± 13%. Horizontal and vertical diffusion are weak with 1 ± 2% and 1 ± 1%, respectively. The atmospheric contribution is insignificant close to the Antarctic continent but increases to 15 ± 10% in the remotest waters (〉1500 km offshore) of the ACC. Translating Southern Ocean carbon fixation by primary producers into biogenic Fe fixation shows a twofold excess of new DFe input close to the Antarctic continent and a one-third shortage in the open ocean. Fe recycling, with an estimated “fe” ratio of 0.59, is the likely pathway to balance new DFe supply and Fe fixation.

Description: In temperate regions, the budburst date of deciduous trees is mainly regulated by temperature variation, but the exact nature of the temperature dependence has been a matter of debate. One hypothesis is that budburst date depends purely on the accumulation of warm temperature; a competing hypothesis states that exposure to cold temperatures is also important for budburst. In this study, variability in budburst is evaluated using 15 years of budburst data for 17 tree species at Harvard Forest. We compare two budburst hypotheses through reversible jump Markov chain Monte Carlo. We then investigate how uncertainties in budburst date mapped onto uncertainties in ecosystem carbon using the Geophysical Fluid Dynamics Laboratory's LM3 land model. For 15 of 17 species, we find that more complicated budburst models that account for a chilling period are favored over simpler models that do not include such dependence. LM3 simulations show that the choice of budburst model induces differences in the timing of carbon uptake commencement of ∼11 days, in the magnitude of April–May carbon uptake of ∼1.03 g C m−2 day−1, and in total ecosystem carbon stocks of ∼2 kg C m−2. While the choice of whether to include a chilling period in the budburst model strongly contributes to this variability, another important factor is how the species-dependent field data gets mapped onto LM3's single deciduous plant functional type (PFT). We conclude budburst timing has a strong impact on simulated CO2 fluxes, and uncertainty in the fluxes can be substantially reduced by improving the model's representation of PFT diversity.

Description: The future carbon balance of high-latitude ecosystems is dependent on the sensitivity of biological processes (photosynthesis and respiration) to the physical changes occurring with permafrost thaw. Predicting C exchange in these ecosystems is difficult because the thawing of permafrost is a heterogeneous process that creates a complex landscape. We measured net ecosystem exchange of C using eddy covariance (EC) in a tundra landscape visibly undergoing thaw during two 6 month campaigns in 2008 and 2009. We developed a spatially explicit quantitative metric of permafrost thaw based on variation in microtopography and incorporated it into an EC carbon flux estimate using a generalized additive model (GAM). This model allowed us to make predictions about C exchange for the landscape as a whole and for specific landscape patches throughout the continuum of permafrost thaw and ground subsidence. During June through November 2008, the GAM predicted that the landscape on average took up 337.1 g C m−2 via photosynthesis and released 289.5 g C m−2 via respiration, resulting in a net C gain of 47.5 g C m−2 by the tundra ecosystem. During April through October 2009, the landscape on average took up 498.7 g C m−2 and released 410.3 g C m−2, resulting in a net C gain of 87.8 g C m−2. On average, between the years, areas with the highest permafrost thaw and ground subsidence photosynthesized 17.7% more and respired 3.3% more C than the average landscape. Areas with the least thaw and subsidence photosynthesized 15% less and respired 5.1% less than the landscape on average. By incorporating spatial variation into the EC C estimate, we were able to estimate the C balance of a heterogeneous landscape and determine the collective effect of permafrost thaw on the plant and soil processes that drive ecosystem C flux. In these study years, permafrost thaw appeared to increase the amplitude of the C cycle by stimulating both C release and sequestration, while the ecosystem remained a C sink at the landscape scale.

Description: High-latitude northern rivers export globally significant quantities of dissolved organic carbon (DOC) to the Arctic Ocean. Climate change, and its associated impacts on hydrology and potential mobilization of ancient organic matter from permafrost, is likely to modify the flux, composition, and thus biogeochemical cycling and fate of exported DOC in the Arctic. This study examined DOC concentration and the composition of dissolved organic matter (DOM) across the hydrograph in Siberia's Kolyma River, with a particular focus on the spring freshet period when the majority of the annual DOC load is exported. The composition of DOM within the Kolyma basin was characterized using absorbance-derived measurements (absorbance coefficient a330, specific UV absorbance (SUVA254), and spectral slope ratio SR) and fluorescence spectroscopy (fluorescence index and excitation-emission matrices (EEMs)), including parallel factor analyses of EEMs. Increased surface runoff during the spring freshet led to DOM optical properties indicative of terrestrial soil inputs with high humic-like fluorescence, SUVA254, and low SR and fluorescence index (FI). Under-ice waters, in contrast, displayed opposing trends in optical properties representing less aromatic, lower molecular weight DOM. We demonstrate that substantial losses of DOC can occur via biological (∼30% over 28 days) and photochemical pathways (〉29% over 14 days), particularly in samples collected during the spring freshet. The emerging view is therefore that of a more dynamic and labile carbon pool than previously thought, where DOM composition plays a fundamental role in controlling the fate and removal of DOC at a pan-Arctic scale.

Description: Fluvial transport of dissolved organic carbon (DOC) is an important link in the global carbon cycle. Previous studies largely increased our knowledge of fluvial exports of carbon to the marine system, but considerable uncertainty remains about in-stream/in-river losses of organic carbon. This study presents an empirical method to assess the nonconservative behavior of fluvial DOC at continental scale. An empirical DOC flux model was trained on two different subsets of training catchments, one with catchments smaller than 2,000 km2 (n = 246, avg. 494 km2) and one with catchments larger than 2,000 km2 (n = 207, avg. 26,525 km2). A variety of potential predictors and controlling factors of fluvial DOC fluxes is discussed. The predictors retained for the final DOC flux models are runoff, slope gradient, land cover, and areal proportions of wetlands. According to the spatially explicit extrapolation of the models, in North America south of 60°N, the total fluvial DOC flux from small catchments (25.8 Mt C a−1, std. err.: 12%) is higher than that from large catchments (19.9 Mt C a−1, std. err.: 10%), giving a total DOC loss of 5.9 Mt C a−1 (std. err.: 78%). As DOC losses in headwaters are not represented in this budget, the estimated DOC loss is rather a minimum value for the total DOC loss within the fluvial network.

Description: Severity of burning can influence multiple aspects of forest composition, carbon cycling, and climate forcing. We quantified how burn severity affected vegetation recovery and albedo change during early succession in Canadian boreal regions by combining satellite observations from the Moderate Resolution Imaging Spectroradiometer (MODIS) and the Canadian Large Fire Database. We used the MODIS-derived difference Normalized Burn Ratio (dNBR) and initial changes in spring albedo as measures of burn severity. We found that the most severe burns had the greatest reduction in summer MODIS Enhanced Vegetation Index (EVI) in the first year after fire, indicating greater loss of vegetation cover. By 5–8 years after fire, summer EVI for all severity classes had recovered to within 90%–108% of prefire levels. Spring and summer albedo progressively increased during the first 7 years after fire, with more severely burned areas showing considerably larger postfire albedo increases during spring and more rapid increases during summer as compared with moderate- and low-severity burns. After 5–7 years, increases in spring albedo above prefire levels were considerably larger in high-severity burns (0.20 ± 0.06; defined by dNBR percentiles greater than 75%) as compared to changes observed in moderate- (0.16 ± 0.06; for dNBR percentiles between 45% and 75%) or low-severity burns (0.13 ± 0.06; for dNBR percentiles between 20% and 45%). The sensitivity of spring albedo to dNBR was similar in all ecozones and for all vegetation types along gradients of burn severity. These results suggest carbon losses associated with increases in burn severity observed in some areas of boreal forests may be at least partly offset, in terms of climate impacts, by increases in negative forcing associated with changes in surface albedo.

Description: Desertification impacts a large proportion of drylands and can be driven by a variety of climate and land use factors. Most conceptual models of desertification include the underlying assumption that when herbaceous cover is reduced, increased erosion from bare patches is redistributed to shrub canopy patches, resulting in self-reinforcing “islands of fertility.” Notably, however, this underlying assumption has not been explicitly tested with direct field measurements. Here we provide direct measurements of horizontal sediment flux moving into and out of bare-, herbaceous-, and shrub-dominated patch types in a semiarid ecosystem for both simulated and natural dust events, as well as in response to simulated disturbance. Horizontal sediment flux out of the bare patches was ∼20% greater than the herbaceous patches and ∼50% greater than sediment flux out of the shrub-dominated patches. Differences among vegetation patch types indicate that shrub patches capture more sediment than herbaceous patches and, importantly, that bare patches serve as amplified sediment sources following disturbance. Our results provide explicit support for the pervasive but untested desertification redistribution assumption, highlighting that loss of grass cover is a compounding problem that not only increases dust emissions but also precludes capture, and may have global relevance for coupled human-environmental systems at risk due to current or potential desertification.

Description: The Mekong River ranks within the top ten rivers of the world in terms of water discharge and sediment load to the ocean, yet its organic matter (OM) composition remains unstudied. This river is experiencing anthropogenically forced changes due to land use and impoundment, and these changes are expected to intensify in the future. Accordingly, we monitored the composition (including vascular-plant signatures) of Mekong River fine particulate organic matter (FPOM) over a one-year period. Autochthonous production comprises a greater proportion of FPOM during the dry season than in the rainy season, as demonstrated by higher percent organic carbon values (7.9 ± 2.4 versus 2.2 ± 0.4%), lower yields of lignin normalized to carbon (0.40 ± 0.05 versus 1.1 ± 0.3 mg (100 mg OC)−1), and an increase in N:C ratios toward phytoplankton values during the dry season (from 0.06 to 0.12). Changes in the lignin-phenol composition of FPOM suggest that gymnosperms contribute more toward FPOM composition during the dry season, with angiosperms dominating in the wet season. This is supported by calculations of the lignin phenol vegetation index of riverine FPOM, which increases between the dry to wet seasons (dry: 29.4 ± 15.0 versus wet: 74.6 ± 17.3). These changes likely reflect seasonal differences in the proportion of flow that is coming from the Upper and Lower Basin, corresponding to compositional differences between the vegetation of these regions. Therefore, this work provides a baseline understanding of FPOM variability that can be used to assess how future change will affect this river.

Description: Approaches are needed to better predict spatial variation in riverine Hg concentrations across heterogeneous landscapes that include mountains, wetlands, and open waters. We applied multivariate linear regression to determine the landscape factors and chemical variables that best account for the spatial variation of total Hg (THg) and methyl Hg (MeHg) concentrations in 27 sub-basins across the 493 km2 upper Hudson River basin in the Adirondack Mountains of New York. THg concentrations varied by sixfold, and those of MeHg by 40-fold in synoptic samples collected at low-to-moderate flow, during spring and summer of 2006 and 2008. Bivariate linear regression relations of THg and MeHg concentrations with either percent wetland area or DOC concentrations were significant but could account for only about 1/3 of the variation in these Hg forms in summer. In contrast, multivariate linear regression relations that included metrics of (1) hydrogeomorphology, (2) riparian/wetland area, and (3) open water, explained about 66% to 〉90% of spatial variation in each Hg form in spring and summer samples. These metrics reflect the influence of basin morphometry and riparian soils on Hg source and transport, and the role of open water as a Hg sink. Multivariate models based solely on these landscape metrics generally accounted for as much or more of the variation in Hg concentrations than models based on chemical and physical metrics, and show great promise for identifying waters with expected high Hg concentrations in the Adirondack region and similar glaciated riverine ecosystems.

Description: Pulsed hydrologic inputs interact with antecedent moisture conditions to shape biogeochemical dynamics in many ecosystems, but the outcomes of these interactions remain difficult to predict. Hydrologic pulses may influence biogeochemical activity through several mechanisms: by providing water as a resource, providing limiting nutrients or substrates that fuel particular biogeochemical pathways, or determining redox conditions. Antecedent moisture conditions may modify the relative importance of each of these potential mechanisms, by influencing accumulation of labile carbon and nutrients, the severity of water limitation to biological processes, and longer-term effects on abiotic conditions, including redox. We experimentally applied hydrologic pulses of different sizes (1-cm and 20-cm events) to soils of desert floodplains and assessed responses of trace gases (CO2, CH4, NO, and N2O) in dry and monsoon seasons to test these mechanisms. Size of the hydrologic pulse strongly interacted with antecedent soil-moisture conditions to determine emissions of some trace gases. Following dry antecedent conditions, water addition stimulated emissions of CO2, CH4, and NO, but not N2O, and larger experimental pulses resulted in larger fluxes. In the monsoon season, responses to water addition were muted and size of the hydrologic pulse had no effect, except for CH4 emission, which increased in response to the 20-cm event. Seasonal contrasts indicated that antecedent moisture conditions constrain the effects of hydrologic pulses on biogeochemical processes, whereas contrasts among responses of different trace gases demonstrated that mechanisms controlling emissions of particular gases are water limitation (CO2), in situ production of nitrogen substrates (NO), or redox conditions (CH4). Strong and predictable interactive effects of water inputs and antecedent conditions indicate that extended droughts may cause elevated emissions of gaseous C and NO following the return of precipitation, whereas larger floods or longer wet seasons are expected to dampen gaseous fluxes, which may contribute to conserving soil C and nutrients within floodplains.

Description: Thermokarst lakes alter landscape topography and hydrology in widespread permafrost regions and mobilize significant permafrost carbon pools, including releasing methane (CH4) to the atmosphere. Despite this, the dynamics of lake evolution, permafrost thawing, and carbon mobilization are not well known. We present a 3-D numerical model of thermokarst lakes on organic-rich yedoma permafrost terrains with surface water flow and pooling naturally defining lakes that deepen, expand laterally, and drain due to talik formation, bank retreat, and both gradual and catastrophic drainage. We predict the 3-D pattern of microbial methane production within the talik over time. As a first model test and calibration, beginning with small protolakes, we simulated 10,000 years of evolution of Pear and Claudi lakes, two neighboring thermokarst features on the northern Seward Peninsula, Alaska. Simulated lakes approximated observed bathymetry, but results are sensitive to initial topography and soil ice content. Local topography caused markedly different dynamics for the two lakes. Pear expanded rapidly across low-relief topography, fully drained multiple times, and released little methane in later stages due to Pleistocene carbon depletion by the first and largest lake generation. Claudi grew slowly and continuously across high-relief topography, forming high subaerial banks; partial drainages left remnant horseshoe lakes that continued to expand into virgin yedoma, mobilizing carbon at roughly the same rate irrespective of lake drainage. The ∼2× discrepancy between simulated CH4 production and observed emission rates in Claudi likely results from misestimation of hot spot ebullition, labile carbon content, CH4:CO2 production ratio, or microbial CH4 oxidation.

Description: A leaf stomatal conductance model was combined with a hydrological tree and soil water flow model and a spatially explicit three-dimensional canopy light model. The model was applied to single, old-growth Fagus sylvatica L. trees, and the measured daily values of stem sap flux could be reproduced with a normalized root mean square error of 0.10 for an observation period of 32 days in the summer of 2009. The high temporal resolution of the model also makes it possible to simulate the diurnal dynamics of transpiration, stem sap flux, and root water uptake. We applied new data-processing algorithms to information from terrestrial laser scans to represent the canopies of the functional-structural model. The high spatial resolution of the root and branch geometry and connectivity makes the detailed modeling of the water usage of single trees possible and allows for the analysis of the interaction between single trees and the influence of the canopy light regime on the water flow inside the xylem. In addition to the laser scans of the observed trees, the model needs tree-species-specific physiological input parameters, which are easy to obtain. The model can be applied at various sites and to different tree species, allowing the up-scaling of the water usage of single trees to the total transpiration of mixed stands.

Description: Terrestrially derived dissolved organic matter (DOM) can impact the fate of persistent organic pollutants (POPs) that are transported to the Arctic via global distillation. Interactions between DOM and POPs through hydrophobic binding processes may influence their photofate in arctic surface waters. We examined the DOM-mediated photodegradation of γ-hexachlorocyclohexane (lindane) and hexachlorobenzene (HCB) in arctic surface waters. These two halogenated organic compounds are commonly detected in the Arctic. We examined how different sources of DOM affect the indirect photolysis of these compounds. We conducted our study using DOM from arctic streams and lakes near the Toolik Lake Long-term Ecological Research Site. HCB or lindane was irradiated in the presence and absence of DOM from these sources, both at the surface of an arctic lake and at 10 cm depth, to investigate the indirect phototransformation of these two compounds and the depth dependence of the observed chemistry. In both artificial and natural sunlight, two of four DOM sources studied stimulated the photodegradation of HCB but not of lindane, suggesting that the indirect phototransformation is a selective process depending on the interactions between DOM and POPs. Through solubility studies, we found that HCB readily partitions to isolated Toolik Lake DOM, while lindane shows no affinity for DOM, findings that corroborate results previously reported in the literature. We demonstrate for the first time DOM's role as a sensitizer for photodegradation of some POPs under field conditions, thus confirming that this process may be an important control on the fate of POPs exhibiting a strong affinity for DOM.

Description: Eddy covariance flux towers provide continuous measurements of ecosystem-level net exchange of carbon, water, energy, and other trace gases between land surface and the atmosphere. The upscaling of flux observations from towers to broad regions provides a new and independent approach for quantifying these fluxes over regions, continents, or the globe. The seven contributions of this special section reflect the most recent advances in the upscaling of fluxes from towers to these broad regions. The section mainly stems from presentations at the recent North American Carbon Program (NACP), FLUXNET, and AGU meetings. These studies focus on different aspects of upscaling: (1) assessing the representativeness of flux networks; (2) upscaling fluxes from towers to broad spatial scales; (3) examining the magnitude, distribution, and interannual variability of fluxes over regions, continents, or the globe; and (4) evaluating the impacts of spatial heterogeneity and parameter variability on flux estimates. Collectively, this special issue provides a timely update on upscaling science and also generates gridded flux data that can be used for model evaluations. Future upscaling studies are expected to advance toward incorporating the impacts of disturbance on ecosystem carbon dynamics, quantifying uncertainties associated with gridded flux estimates, and comparing various upscaling methods and the resulting gridded flux fields.

Description: Mercury (Hg) is one of the leading water quality concerns in surface waters of the United States. Although watershed-scale Hg cycling research has increased in the past two decades, advances in modeling watershed Hg processes in diverse physiographic regions, spatial scales, and land cover types are needed. The goal of this study was to assess Hg cycling in a Coastal Plain system using concentrations and fluxes estimated by multiple watershed-scale models with distinct mathematical frameworks reflecting different system dynamics. We simulated total mercury (HgT, the sum of filtered and particulate forms) concentrations and fluxes from a Coastal Plain watershed (McTier Creek) using three watershed Hg models and an empirical load model. Model output was compared with observed in-stream HgT. We found that shallow subsurface flow is a potentially important transport mechanism of particulate HgT during periods when connectivity between the uplands and surface waters is maximized. Other processes (e.g., stream bank erosion, sediment re-suspension) may increase particulate HgT in the water column. Simulations and data suggest that variable source area (VSA) flow and lack of rainfall interactions with surface soil horizons result in increased dissolved HgT concentrations unrelated to DOC mobilization following precipitation events. Although flushing of DOC-HgT complexes from surface soils can also occur during this period, DOC-complexed HgT becomes more important during base flow conditions. TOPLOAD simulations highlight saturated subsurface flow as a primary driver of daily HgT loadings, but shallow subsurface flow is important for HgT loads during high-flow events. Results suggest limited seasonal trends in HgT dynamics.

Description: Shallow coastal bays provide habitat for diverse fish and invertebrate populations and are an important source of sediment for surrounding marshes. The sediment dynamics of these bays are strongly affected by seagrass meadows, which limit sediment resuspension, thereby providing a more favorable light environment for their own survival and growth. Due to this positive feedback between seagrass and light conditions, it has been suggested that bare sediment and seagrass meadows are potential alternate stable states of the benthos in shallow coastal bays. To investigate the stability and resilience of seagrass meadows subjected to variation in environmental conditions (e.g., light, temperature), a coupled model of vegetation–sediment–water flow interactions and vegetation growth was developed. The model was used to examine the effect of dynamically varying seasonal and interannual seagrass density on sediment resuspension, water column turbidity, and the subsequent light environment on hourly time steps and then run over decadal time scales. A daily growth model was designed to capture both belowground biomass and the growth and senescence of aboveground biomass structural components (e.g., leaves and stems). This allowed us to investigate how the annual and seasonal variability in shoot and leaf density within a meadow affects the strength of positive feedbacks between seagrass and their light environment. The model demonstrates both the emergence of bistable behavior from 1.6 to 1.8 m mean sea level due to the strength of the positive feedback, as well as the limited resilience of seagrass meadows within this bistable range.

Description: A greater abundance of shrubs in semiarid grasslands affects the spatial patterns of soil temperature, moisture, and litter, resulting in fertile islands with potentially enhanced soil metabolic activity. The goal of this study was to quantify the microsite specificity of soil respiration in a semiarid riparian ecosystem experiencing shrub encroachment. We quantified the response of soil respiration to different microsite conditions created by big mesquite shrubs (near the trunk and the canopy edge), medium-sized mesquite, sacaton bunchgrasses, and open spaces. We hypothesized that soil respiration would be more temperature sensitive and less moisture sensitive and have a greater magnitude in shrub microsites compared with grass and open microsites. Field and incubation soil respiration data were simultaneously analyzed in a Bayesian framework to quantify the microsite-specific temperature and moisture sensitivities and magnitude of respiration. The analysis showed that shrub expansion increases the heterogeneity of respiration. Respiration has greater temperature sensitivity near the shrub canopy edge, and respiration rates are higher overall under big mesquite compared with those of the other microsites. Respiration in the microsites beneath medium-sized mesquites does not behave like a downscaled version of big mesquite microsites. The grass microsites show more similarity to big mesquite microsites than medium-sized shrubs. This study shows there can be a great deal of fine-scale spatial heterogeneity that accompanies shifts in vegetation structure. Such complexity presents a challenge in scaling soil respiration fluxes to the landscape for systems experiencing shrub encroachment, but quantifying this complexity is significantly important in determining overall ecosystem metabolic behavior.

Description: We studied the photochemical transformation of dissolved organic matter (DOM) to dissolved inorganic carbon (DIC), ammonium (NH4+), and labile organic substrates supporting bacterial carbon biomass along a salinity gradient throughout the Baltic Sea during summer, autumn, and spring. The photoproduced DIC, NH4+, and labile DOM supporting bacterial biomass were related to the number of photons absorbed during the irradiations of biologically recalcitrant DOM to determine apparent quantum yields. The apparent quantum yields for the photoproduction of DIC and NH4+ lacked seasonal variation, but behaved differently along the salinity gradient; the photoproduction of DIC decreased, while photoammonification increased with increasing salinity. The apparent quantum yield for the photoproduction of labile DOM supporting bacterial biomass was highest in summer and unaffected by salinity. The annual photoammonification rate over the entire Baltic Sea ranged from 0.038 to 0.049 Tg N, equivalent to 13%–23% of the annual atmospheric deposition of inorganic N. The annual phototransformation of dissolved organic carbon (DOC), including the direct photomineralization and indirect bacterial mineralization of photoproduced labile DOM (total of 2.71–3.94 Tg C), exceeded the annual river loading of photoreactive DOC, assuming that half of the total river DOC input to the Baltic Sea is photoreactive. As the annual photomineralization of DOC exceeded the annual terrestrial input of photoreactive DOC to the Baltic Sea, the photochemical transformation is a major sink for terrestrial DOC in such coastal systems.

Description: Loss of permafrost can modify the export and composition of dissolved organic carbon (DOC) from subarctic peatlands by changing the hydrological regime and altering ecosystem structure and function. In Stordalen peatland complex (68.20°N, 19.03°E) recent permafrost thaw has caused a conversion of the palsa parts (an ombrotrophic, permafrost affected peatland type) into both bog and flow-through fen peatland types. Within the Stordalen peatland complex we estimated the DOC mass balance and assessed DOC composition for one palsa catchment, one bog catchment and two fens in order to assess the possible impacts of permafrost thaw on peatland complex DOC export. The fens were found to have higher net DOC export rates at 8.1 and 7.0 g C m−2 yr−1 than either the palsa or bog catchments, at 3.2 and 3.5 g C m−2 yr−1, respectively. The snowmelt period was more important for the annual DOC export from the palsa and bog catchments than for the fens, representing 65–100% of the palsa and bog catchment exports while 35–60% of the net fen exports. DOC exported from the palsa and bog catchments were characterized by high aromaticity, molecular weight, C/N ratios, and contained DOC of primarily terrestrial origin. The fens exhibited a shift in DOC composition between inflows and outflows that suggested that fens act as catchment locations for degradation and transformation of DOC. Permafrost thaw can thus alter the magnitude, timing, and composition of subarctic peatland DOC exports due to interactions among peatland type, permafrost conditions, and hydrological setting.

Description: For some predominately phosphorus (P)-limited ecosystems, vegetation can be sustained under steady state conditions by atmospheric P inputs. The structure of the canopy influences deposition via the ability of the canopy to trap airborne P. This dependence suggests that a positive feedback may exist, which would have important impacts on the process of forest regeneration. One source of disturbance in the tropics is shifting cultivation. Over multiple cycles, studies have shown that shifting cultivation can lead to a large net loss of soil P, with a concomitant decline in forest biomass. In this study a model was developed to assess how shifting cultivation affects vegetation and phosphorus dynamics. This model is applied as a case study to a dry tropical forest system in the Southern Yucatan that is primarily P limited and has experienced shifting cultivation over several decades. Results show that following the second cycle, recovery would only be achieved after 100 years or longer in comparison to ∼30 years after one cycle. Examining the stable states of the system suggests that two stable states exist and that state changes as brought about by repeated disturbance can cause a shift to the other “low vegetation” stable state. Thus, for predominately P-limited ecosystems undergoing repeated disturbance, the depletion of soil P can significantly affect the long-term ability of the forest vegetation to recover. Results from this study have widespread implications, as 400 million ha of forest are affected by shifting cultivation.

Description: The purpose of this study was to quantify the effects of clear-cutting and site preparation on dissolved organic carbon (DOC) concentrations and export in four boreal headwater streams in northern Sweden. The data set included intensive stream water monitoring from 2 years of pretreatment conditions (2004–2005), a 2 year post-clear-cut period (2006–2007), and a 2 year period after site preparation (2008–2009). To investigate differences in [DOC], an analysis of variance on ranks was performed on the data sets. Clear-cutting increased the median DOC concentrations significantly from 15.9 to 20.4 mg L−1, which represents a net increase (treatment versus control) of 3.0 mg L−1 in the 2006–2007 period. Site preparation had an even more profound effect on DOC levels; an increase from 20.4 to 27.6 mg L−1 was found in the site-prepared catchments, whereas the control sites increased slightly from 17.4 to 21.4 mg L−1 during the wetter years of 2008–2009. Riverine C fluxes increased significantly by 100% after clear-cutting and by 79% after site preparation (92% and 195%, respectively, if compared to pretreatment conditions). When comparing these yearly C fluxes (183 kg C ha−1 yr−1 after clear-cutting; 280 kg C ha−1 yr−1 after site preparation) to the net ecosystem exchange (NEE) of a forest in the region, the DOC flux represented 10% of NEE before harvest, increased to 18% after the clear-cut, and increased to 28% after site preparation. These results underline the large impact of forestry operations on stream water quality as well as DOC exports leaving managed boreal forests.

Description: Post-fire storage of carbon (C) in organic-soil horizons was measured in one Canadian and three Alaskan chronosequences in black spruce forests, together spanning stand ages of nearly 200 yrs. We used a simple mass balance model to derive estimates of inputs, losses, and accumulation rates of C on timescales of years to centuries. The model performed well for the surface and total organic soil layers and presented questions for resolving the dynamics of deeper organic soils. C accumulation in all study areas is on the order of 20–40 gC/m2/yr for stand ages up to ∼200 yrs. Much larger fluxes, both positive and negative, are detected using incremental changes in soil C stocks and by other studies using eddy covariance methods for CO2. This difference suggests that over the course of stand replacement, about 80% of all net primary production (NPP) is returned to the atmosphere within a fire cycle, while about 20% of NPP enters the organic soil layers and becomes available for stabilization or loss via decomposition, leaching, or combustion. Shifts toward more frequent and more severe burning and degradation of deep organic horizons would likely result in an acceleration of the carbon cycle, with greater CO2 emissions from these systems overall.

Description: Wind erosion and large dust plumes are an increasingly important attribute in cold-desert rangelands, particularly as wildfire increases. Fire reduces vegetation, which increases erosivity. Whether sediment supply increases after fire has not been determined in this environment. We asked how sediment supply varied among sites burned 2-months to 5-years previously, in comparison to unburned sagebrush steppe, across 500 km of southern Idaho, USA. We measured potential dust emissions (PM10, particles

Description: A dynamic model is constructed for interactive silicon, nitrogen, sulfur processing in and below Arctic sea ice, by ecosystems residing in the lower few centimeters of the distributed pack. A biogeochemically active bottom layer supporting sources/sinks for the pennate diatoms is appended to thickness categories of a global sea ice code. Nutrients transfer from the ocean mixed layer to drive algal growth, while sulfur metabolites are reinjected from the ice interface. Freeze, flux, flush and melt processes are linked to multielement geocycling for the entire high-latitude regime. Major element kinetics are optimized initially to reproduce chlorophyll observations, which extend across the seasons. Principal influences on biomass are solute exchange velocity at the solid interface, optical averaging in active ice and cell retention against ablation. The sulfur mechanism encompasses open water features such as accumulation of particulate dimethyl sulfoniopropionate, grazing and other disruptive releases, plus bacterial/enzymatic conversion to volatile dimethyl sulfide. For baseline settings, the mixed layer trace gas distribution matches sparging measurements where they are available. However, concentrations rise to well over 10 nM in remote, unsampled locations. Peak contributions are supported by ice grazing, mortality and fractional melting. The model bottom layer adds substantially to a ring maximum of reduced sulfur chemistry that may be dominant across the marginal Arctic environment. Sensitivity tests on this scenario include variation of cell sulfur composition and remineralization, routings/chemical time scales, and the physical dimension of water layers. An alternate possibility that peripheral additions are small cannot be excluded from the outcomes. It is concluded that seagoing dimethyl sulfide data are far too sparse at the present time to distinguish sulfur-ice production levels.

Description: This study presents observations of atmospheric boundary layer CO2 mole fraction from a nine-tower regional network deployed during the North American Carbon Program's Mid-Continent Intensive (MCI) during 2007–2009. The MCI region is largely agricultural, with well-documented carbon exchange available via agricultural inventories. By combining vegetation maps and tower footprints, we show the fractional influence of corn, soy, grass, and forest biomes varies widely across the MCI. Differences in the magnitude of CO2 flux from each of these biomes lead to large spatial gradients in the monthly averaged CO2 mole fraction observed in the MCI. In other words, the monthly averaged gradients are tied to regional patterns in net ecosystem exchange (NEE). The daily scale gradients are more weakly connected to regional NEE, instead being governed by local weather and large-scale weather patterns. With this network of tower-based mole fraction measurements, we detect climate-driven interannual changes in crop growth that are confirmed by satellite and inventory methods. These observations show that regional-scale CO2 mole fraction networks yield large, coherent signals governed largely by regional sources and sinks of CO2.

Description: Dissolved organic matter (DOM) fuels the majority of in-stream microbial processes, including the removal of nitrate via denitrification. However, little is known about how the chemical composition of DOM influences denitrification rates. Water and sediment samples were collected across an ecosystem gradient, spanning the alpine to plains, in central Colorado to determine whether the chemical composition of DOM was related to denitrification rates. Laboratory bioassays measured denitrification potentials using the acetylene block technique and carbon mineralization via aerobic bioassays, while organic matter characteristics were evaluated using spectroscopic and fractionation methods. Denitrification potentials under ambient and elevated nitrate concentrations were strongly correlated with aerobic respiration rates and the percent mineralized carbon, suggesting that information about the aerobic metabolism of a system can provide valuable insight regarding the ability of the system to additionally reduce nitrate. Multiple linear regressions (MLR) revealed that under elevated nitrate concentrations denitrification potentials were positively related to the presence of protein-like fluorophores and negatively related to more aromatic and oxidized fractions of the DOM pool. Using MLR, the chemical composition of DOM, carbon, and nitrate concentrations explained 70% and 78% of the observed variability in denitrification potential under elevated and ambient nitrate conditions, respectively. Thus, it seems likely that DOM optical properties could help to improve predictions of nitrate removal in the environment. Finally, fluorescence measurements revealed that bacteria used both protein and humic-like organic molecules during denitrification providing further evidence that larger, more aromatic molecules are not necessarily recalcitrant in the environment.

Description: Wetland soil oxygen (O2) is rarely measured, which limits our understanding of a key regulator of nitrogen loss through denitrification. We asked: (1) How does soil [O2] vary in riparian wetlands? (2) How does this [O2] variation affect denitrification rates and end products? and (3) How does [O2] variation and previous exposure to O2 affect trace gas fluxes? We collected a continuous seven-month record of [O2] dynamics in a “wet” and “dry” riparian zone. In April 2009, soil [O2] ranged from 0 to 13% and consistently increased with increasing distance from the stream. [O2] gradually declined in all sensors until all sensors went anoxic in early September 2009. In mid-fall, a dropping water table increased soil [O2] to 15–20% within a 2–3 day period. We measured denitrification using the Nitrogen-Free Air Recirculation Method (N-FARM), a direct measurement of N2 production against a helium background. Denitrification rates were significantly higher in the wetter areas, which correlated to lower O2 conditions. Denitrification rates in the drier areas correlated with [O2] in the early spring and summer, but significantly decreased in late summer despite decreasing O2 concentrations. Increasing [O2] significantly increased core N2O production, and therefore may be an important control on nitrous oxide yield. Field N2O fluxes, however, were highly variable, ranging from 0 to 800 ug N m−2 hr−1 with no differences between the wet and dry sites. Future research should focus on understanding the biotic and abiotic controls on O2 dynamics, and O2 dynamics should be included in models of soil N cycling and trace gas fluxes.

Description: Nitrogen export from small forested watersheds is known to be affected by N deposition but with high regional variability. We studied 10 headwater catchments in the northeastern United States across a gradient of N deposition (5.4 − 9.4 kg ha−1 yr−1) to determine if soil nitrification rates could explain differences in stream water NO3− export. Average annual export of two years (October 2002 through September 2004) varied from 0.1 kg NO3−-N ha−1 yr−1 at Cone Pond watershed in New Hampshire to 5.1 kg ha−1 yr−1 at Buck Creek South in the western Adirondack Mountains of New York. Potential net nitrification rates and relative nitrification (fraction of inorganic N as NO3−) were measured in Oa or A soil horizons at 21–130 sampling points throughout each watershed. Stream NO3− export was positively related to nitrification rates (r2 = 0.34, p = 0.04) and the relative nitrification (r2 = 0.37, p = 0.04). These relationships were much improved by restricting consideration to the 6 watersheds with a higher number of rate measurements (59–130) taken in transects parallel to the streams (r2 of 0.84 and 0.70 for the nitrification rate and relative nitrification, respectively). Potential nitrification rates were also a better predictor of NO3− export when data were limited to either the 6 sampling points closest to the watershed outlet (r2 = 0.75) or sampling points

Description: There is current debate about whether the balance of photosynthesis and respiration has any impact on the net accumulation of organic matter on glacier surfaces. This study assesses controls on rates of net ecosystem production (NEP), respiration, and photosynthesis in cryoconite holes during the main melt season (June–August 2009) on three valley glaciers in Svalbard. Cryoconite thickness and organic matter content explained 87% of the total variation in rates of respiration (in units of volume), and organic matter (but not sediment depth) was a significant (p 〈 0.05) control on photosynthesis (by volume). The average rates of respiration and gross photosynthesis within the cryoconite holes were overall closely balanced, ranging from net autotrophic to heterotrophic. Sediment depth explained over half the variation of NEP, with net autotrophic rates typical only in sediment

Description: Increasing soil temperature has the potential to alter the activity of the extracellular enzymes that mobilize nitrogen (N) from soil organic matter (SOM) and ultimately the availability of N for primary production. Proteolytic enzymes depolymerize N from proteinaceous components of SOM into amino acids, and their activity is a principal driver of the within-system cycle of soil N. The objectives of this study were to investigate whether the soils of temperate forest tree species differ in the temperature sensitivity of proteolytic enzyme activity over the growing season and the role of substrate limitation in regulating temperature sensitivity. Across species and sampling dates, proteolytic enzyme activity had relatively low sensitivity to temperature with a mean activation energy (Ea) of 33.5 kJ mol−1. Ea declined in white ash, American beech, and eastern hemlock soils across the growing season as soils warmed. By contrast, Ea in sugar maple soil increased across the growing season. We used these data to develop a species-specific empirical model of proteolytic enzyme activity for the 2009 calendar year and studied the interactive effects of soil temperature (ambient or +5°C) and substrate limitation (ambient or elevated protein) on enzyme activity. Declines in substrate limitation had a larger single-factor effect on proteolytic enzyme activity than temperature, particularly in the spring. There was, however, a large synergistic effect of increasing temperature and substrate supply on proteolytic enzyme activity. Our results suggest limited increases in N availability with climate warming unless there is a parallel increase in the availability of protein substrates.

Description: Wetlands provide important ecohydrological services by regulating fluxes of nutrients and pollutants to receiving waters, which can in turn mitigate adverse effects on water quality. Turnover of redox-sensitive solutes in wetlands has been shown to take place in distinct spatial and temporal patterns, commonly referred to as hot spots and hot moments. Despite the importance of such patterns for solute fluxes the mechanistic understanding of their formation is still weak and their existence is often explained by variations in soil properties and diffusive transport only. Here we show that surface micro-topography in wetlands can cause the formation of biogeochemical hot spots solely by the advective redistribution of infiltrating water as a result of complex subsurface flow patterns. Surface and subsurface flows are simulated for an idealized section of a riparian wetland using a fully integrated numerical code for coupled surface-subsurface systems. Biogeochemical processes and transport along advective subsurface flow paths are simulated kinetically using the biogeochemical code PHREEQC. Distinct patterns of biogeochemical activity (expressed as reaction rates) develop in response to micro-topography induced subsurface flow patterns. Simulated vertical pore water profiles for various redox-sensitive species resemble profiles observed in the field. This mechanistic explanation of hot spot formation complements the more static explanations that relate hot spots solely to spatial variability in soil characteristics and can account for spatial as well as temporal variability of biogeochemical activity, which is needed to assess future changes in the biogeochemical turnover of wetland systems.