N and P constrain C in ecosystems under climate change: Role of nutrient redistribution, accumulation, and stoichiometry

We use the Multiple Element Limitation (MEL) model to examine responses of 12 ecosystems to elevated carbon dioxide (CO₂ ), warming, and 20% decreases or increases in precipitation. Ecosystems respond synergistically to elevated CO₂ , warming, and decreased precipitation combined because higher water-use efficiency with elevated CO₂ and higher fertility with warming compensate for responses to drought. Response to elevated CO₂ , warming, and increased precipitation combined is additive. We analyze changes in ecosystem carbon (C) based on four nitrogen (N) and four phosphorus (P) attribution factors: (1) changes in total ecosystem N and P, (2) changes in N and P distribution between vegetation and soil, (3) changes in vegetation C:N and C:P ratios, and (4) changes in soil C:N and C:P ratios. In the combined CO₂ and climate change simulations, all ecosystems gain C. The contributions of these four attribution factors to changes in ecosystem C storage varies among ecosystems because of differences in the initial distributions of N and P between vegetation and soil and the openness of the ecosystem N and P cycles. The net transfer of N and P from soil to vegetation dominates the C response of forests. For tundra and grasslands, the C gain is also associated with increased soil C:N and C:P. In ecosystems with symbiotic N fixation, C gains resulted from N accumulation. Because of differences in N versus P cycle openness and the distribution of organic matter between vegetation and soil, changes in the N and P attribution factors do not always parallel one another. Differences among ecosystems in C-nutrient interactions and the amount of woody biomass interact to shape ecosystem C sequestration under simulated global change. We suggest that future studies quantify the openness of the N and P cycles and changes in the distribution of C, N, and P among ecosystem components, which currently limit understanding of nutrient effects on C sequestration and responses to elevated CO₂ and climate change.

Local-scale Arctic tundra heterogeneity affects regional-scale carbon dynamics

In northern Alaska nearly 65% of the terrestrial surface is composed of polygonal ground, where geomorphic tundra landforms disproportionately influence carbon and nutrient cycling over fine spatial scales. Process-based biogeochemical models used for local to Pan-Arctic projections of ecological responses to climate change typically operate at coarse-scales (1km2-0.5°) at which fine-scale (<1km2) tundra heterogeneity is often aggregated to the dominant land cover unit. Here, we evaluate the importance of tundra heterogeneity for representing soil carbon dynamics at fine to coarse spatial scales. We leveraged the legacy of data collected near Utqiaġvik, Alaska between 1973 and 2016 for model initiation, parameterization, and validation. Simulation uncertainty increased with a reduced representation of tundra heterogeneity and coarsening of spatial scale. Hierarchical cluster analysis of an ensemble of 21st-century simulations reveals that a minimum of two tundra landforms (dry and wet) and a maximum of 4km2 spatial scale is necessary for minimizing uncertainties (<10%) in regional to Pan-Arctic modeling applications.

Mycobiont contribution to tundra plant acquisition of permafrost-derived nitrogen

As Arctic soils warm, thawed permafrost releases nitrogen (N) that could stimulate plant productivity and thus offset soil carbon losses from tundra ecosystems. Although mycorrhizal fungi could facilitate plant access to permafrost-derived N, their exploration capacity beyond host plant root systems into deep, cold active layer soils adjacent to the permafrost table is unknown. We characterized root-associated fungi (RAF) that colonized ericoid (ERM) and ectomycorrhizal (ECM) shrub roots and occurred below the maximum rooting depth in permafrost thaw-front soil in tussock and shrub tundra communities. We explored the relationships between root and thaw front fungal composition and plant uptake of a ¹⁵ N tracer applied at the permafrost boundary. We show that ERM and ECM shrubs associate with RAF at the thaw front providing evidence for potential mycelial connectivity between roots and the permafrost boundary. Among shrubs and tundra communities, RAF connectivity to the thaw boundary was ubiquitous. The occurrence of particular RAF in both roots and thaw front soil was positively correlated with ¹⁵ N recovered in shrub biomass Taxon-specific RAF associations could be a mechanism for the vertical redistribution of deep, permafrost-derived nutrients, which may alleviate N limitation and stimulate productivity in warming tundra.

Fuel-reduction management alters plant composition, carbon and nitrogen pools, and soil thaw in Alaskan boreal forest

Increasing wildfire activity in Alaska's boreal forests has led to greater fuel-reduction management. Management has been implemented to reduce wildfire spread, but the ecological impacts of these practices are poorly known. We quantified the effects of hand-thinning and shearblading on above- and belowground stand characteristics, plant species composition, carbon (C) and nitrogen (N) pools, and soil thaw across 19 sites dominated by black spruce (Picea mariana) in interior Alaska treated 2-12 years prior to sampling. The density of deciduous tree seedlings was significantly higher in shearbladed areas compared to unmanaged forest (6.4 vs. 0.1 stems/m² ), and unmanaged stands exhibited the highest mean density of conifer seedlings and layers (1.4 stems/m² ). Understory plant community composition was most similar between unmanaged and thinned stands. Shearblading resulted in a near complete loss of aboveground tree biomass C pools while thinning approximately halved the C pool size (1.2 kg C/m² compared to 3.1 kg C/m² in unmanaged forest). Significantly smaller soil organic layer (SOL) C and N pools were observed in shearbladed stands (3.2 kg C/m² and 116.8 g N/m² ) relative to thinned (6.0 kg C/m² and 192.2 g N/m² ) and unmanaged (5.9 kg C/m² and 178.7 g N/m² ) stands. No difference in C and N pool sizes in the uppermost 10 cm of mineral soil was observed among stand types. Total C stocks for measured pools was 2.6 kg C/m² smaller in thinned stands and 5.8 kg C/m² smaller in shearbladed stands when compared to unmanaged forest. Soil thaw depth averaged 13 cm deeper in thinned areas and 46 cm deeper in shearbladed areas relative to adjacent unmanaged stands, although variability was high across sites. Deeper soil thaw was linked to shallower SOL depth for unmanaged stands and both management types, however for any given SOL depth, thaw tended to be deeper in shearbladed areas compared to unmanaged forest. These findings indicate that fuel-reduction management alters plant community composition, C and N pools, and soil thaw depth, with consequences for ecosystem structure and function beyond those intended for fire management.

Polygonal tundra geomorphological change in response to warming alters future CO2 and CH4 flux on the Barrow Peninsula

The landscape of the Barrow Peninsula in northern Alaska is thought to have formed over centuries to millennia, and is now dominated by ice-wedge polygonal tundra that spans drained thaw-lake basins and interstitial tundra. In nearby tundra regions, studies have identified a rapid increase in thermokarst formation (i.e., pits) over recent decades in response to climate warming, facilitating changes in polygonal tundra geomorphology. We assessed the future impact of 100 years of tundra geomorphic change on peak growing season carbon exchange in response to: (i) landscape succession associated with the thaw-lake cycle; and (ii) low, moderate, and extreme scenarios of thermokarst pit formation (10%, 30%, and 50%) reported for Alaskan arctic tundra sites. We developed a 30 × 30 m resolution tundra geomorphology map (overall accuracy:75%; Kappa:0.69) for our ~1800 km² study area composed of ten classes; drained slope, high center polygon, flat-center polygon, low center polygon, coalescent low center polygon, polygon trough, meadow, ponds, rivers, and lakes, to determine their spatial distribution across the Barrow Peninsula. Land-atmosphere CO2 and CH4 flux data were collected for the summers of 2006-2010 at eighty-two sites near Barrow, across the mapped classes. The developed geomorphic map was used for the regional assessment of carbon flux. Results indicate (i) at present during peak growing season on the Barrow Peninsula, CO2 uptake occurs at -902.3 10(6) gC-CO2 day(-1) (uncertainty using 95% CI is between -438.3 and -1366 10(6) gC-CO2 day(-1)) and CH4 flux at 28.9 10(6) gC-CH4 day(-1) (uncertainty using 95% CI is between 12.9 and 44.9 10(6) gC-CH4 day(-1)), (ii) one century of future landscape change associated with the thaw-lake cycle only slightly alter CO2 and CH4 exchange, while (iii) moderate increases in thermokarst pits would strengthen both CO2 uptake (-166.9 10(6) gC-CO2 day(-1)) and CH4 flux (2.8 10(6) gC-CH4 day(-1)) with geomorphic change from low to high center polygons, cumulatively resulting in an estimated negative feedback to warming during peak growing season.

Age-related variation in carbon allocation at tree and stand scales in beech (Fagus sylvatica L.) and sessile oak (Quercus petraea (Matt.) Liebl.) using a chronosequence approach

Two types of physiological mechanisms can contribute to growth decline with age: (i) the mechanisms leading to the reduction of carbon assimilation (input) and (ii) those leading to modification of the resource economy. Surprisingly, the processes relating to carbon allocation have been little investigated as compared to research on the processes governing carbon assimilation. The objective of this paper was thus to test the hypothesis that growth decrease related to age is accompanied by changes in carbon allocation to the benefit of storage and reproductive functions in two contrasting broad-leaved species: beech (Fagus sylvatica L.) and sessile oak (Quercus petraea (Matt.) Liebl.). Age-related changes in carbon allocation were studied using a chronosequence approach. Chronosequences, each consisting of several even-aged stands ranging from 14 to 175 years old for beech and from 30 to 134 years old for sessile oak, were divided into five or six age classes. In this study, carbon allocations to growth, storage and reproduction were defined as the relative amount of carbon invested in biomass increment, carbohydrate increment and seed production, respectively. Tree-ring width and allometric relationships were used to assess biomass increment at the tree and stand scales. Below-ground biomass was assessed using a specific allometric relationship between root:shoot ratio and age, established from the literature review. Seasonal variations of carbohydrate concentrations were used to assess carbon allocation to storage. Reproduction effort was quantified for beech stands by collecting seed and cupule production. Age-related flagging of biomass productivity was assessed at the tree and stand scales, and carbohydrate quantities in trees increased with age for both species. Seed and cupule production increased with stand age in beech from 56 gC m(-)(2) year(-1) at 30 years old to 129 gC m(-2) year(-1) at 138 years old. In beech, carbon allocation to storage and reproductive functions increased with age to the detriment of carbon allocation to growth functions. In contrast, the carbon balance between growth and storage remained constant between age classes in sessile oak. The contrasting age-related changes in carbon allocation between beech and sessile oak are discussed with reference to the differences in growing environment, phenology and hydraulic properties of ring-porous and diffuse-porous species.