USGS ScienceSearch

Geology topics

J. Harden

Publications and source records attributed to J. Harden.

7 recordsLinked to original sources

Simulating boreal forest carbon dynamics after stand-replacing fire disturbance: insights from a global process-based vegetation model

Stand-replacing fires are the dominant fire type in North American boreal forests. They leave a historical legacy of a mosaic landscape of different aged forest cohorts. This forest age dynamics must be included in vegetation models to accurately quantify the role of fire in the historical and current regional forest carbon balance. The present study adapted the global process-based vegetation model ORCHIDEE to simulate the CO 2 emissions from boreal forest fire and the subsequent recovery after a stand-replacing fire; the model represents postfire new cohort establishment, forest stand structure and the self-thinning process. Simulation results are evaluated against observations of three clusters of postfire forest chronosequences in Canada and Alaska. The variables evaluated include: fire carbon emissions, CO 2 fluxes (gross primary production, total ecosystem respiration and net ecosystem exchange), leaf area index, and biometric measurements (aboveground biomass carbon, forest floor carbon, woody debris carbon, stand individual density, stand basal area, and mean diameter at breast height). When forced by local climate and the atmospheric CO 2 history at each chronosequence site, the model simulations generally match the observed CO 2 fluxes and carbon stock data well, with model-measurement mean square root of deviation comparable with the measurement accuracy (for CO 2 flux ~100 g C m −2 yr −1 , for biomass carbon ~1000 g C m −2 and for soil carbon ~2000 g C m −2 ). We find that the current postfire forest carbon sink at the evaluation sites, as observed by chronosequence methods, is mainly due to a combination of historical CO 2 increase and forest succession. Climate change and variability during this period offsets some of these expected carbon gains. The negative impacts of climate were a likely consequence of increasing water stress caused by significant temperature increases that were not matched by concurrent increases in precipitation. Our simulation results demonstrate that a global vegetation model such as ORCHIDEE is able to capture the essential ecosystem processes in fire-disturbed boreal forests and produces satisfactory results in terms of both carbon fluxes and carbon-stock evolution after fire. This makes the model suitable for regional simulations in boreal regions where fire regimes play a key role in the ecosystem carbon balance.

Biogeosciences

Expert assessment of vulnerability of permafrost carbon to climate change

Approximately 1700 Pg of soil carbon (C) are stored in the northern circumpolar permafrost zone, more than twice as much C than in the atmosphere. The overall amount, rate, and form of C released to the atmosphere in a warmer world will influence the strength of the permafrost C feedback to climate change. We used a survey to quantify variability in the perception of the vulnerability of permafrost C to climate change. Experts were asked to provide quantitative estimates of permafrost change in response to four scenarios of warming. For the highest warming scenario (RCP 8.5), experts hypothesized that C release from permafrost zone soils could be 19–45 Pg C by 2040, 162–288 Pg C by 2100, and 381–616 Pg C by 2300 in CO 2 equivalent using 100-year CH 4 global warming potential (GWP). These values become 50 % larger using 20-year CH 4 GWP, with a third to a half of expected climate forcing coming from CH 4 even though CH 4 was only 2.3 % of the expected C release. Experts projected that two-thirds of this release could be avoided under the lowest warming scenario (RCP 2.6). These results highlight the potential risk from permafrost thaw and serve to frame a hypothesis about the magnitude of this feedback to climate change. However, the level of emissions proposed here are unlikely to overshadow the impact of fossil fuel burning, which will continue to be the main source of C emissions and climate forcing.

Climatic Change

The role of historical fire disturbance in the carbon dynamics of the pan-boreal region: A process-based analysis

[1] Wildfire is a common occurrence in ecosystems of northern high latitudes, and changes in the fire regime of this region have consequences for carbon feedbacks to the climate system. To improve our understanding of how wildfire influences carbon dynamics of this region, we used the process-based Terrestrial Ecosystem Model to simulate fire emissions and changes in carbon storage north of 45°N from the start of spatially explicit historically recorded fire records in the twentieth century through 2002, and evaluated the role of fire in the carbon dynamics of the region within the context of ecosystem responses to changes in atmospheric CO 2 concentration and climate. Our analysis indicates that fire plays an important role in interannual and decadal scale variation of source/sink relationships of northern terrestrial ecosystems and also suggests that atmospheric CO 2 may be important to consider in addition to changes in climate and fire disturbance. There are substantial uncertainties in the effects of fire on carbon storage in our simulations. These uncertainties are associated with sparse fire data for northern Eurasia, uncertainty in estimating carbon consumption, and difficulty in verifying assumptions about the representation of fires that occurred prior to the start of the historical fire record. To improve the ability to better predict how fire will influence carbon storage of this region in the future, new analyses of the retrospective role of fire in the carbon dynamics of northern high latitudes should address these uncertainties.

Journal of Geophysical Research G: Biogeosciences

Effect of tree line advance on carbon storage in NW Alaska

We investigated the size, distribution, and temporal dynamics of ecosystem carbon (C) pools in an area of recent tree line advance, northwest Alaska. Repeat aerial photographs show forest cover increased ???10% in our study area since 1949. We sampled C pools of four principal ecosystem types, tussock tundra, shrub tundra, woodland, and forest, all located on a 600-800 year old river terrace. Significant differences between ecosystem C pools, both above ground and below ground existed. Tundra sites store >22.2 kg C/m2, shrub tundra sites and woodland sites store 9.7 kg C/m2 and 14.3 kg C/m2, respectively, and forest sites store 14.4 kg C/m2. Landscape variation of total ecosystem C was primarily due to organic soil C and was secondarily due to C stored in trees. Soil C/N profiles of shrub tundra sites and woodland sites showed similarities with forest site soils at surface and tundra site soils at depth. We hypothesize that tundra systems transformed to forest systems in this area under a progression of permafrost degradation and enhanced drainage. On the basis of C pool estimates for the different ecosystem types, conversion of tundra sites to forest may have resulted in a net loss of > 7.8 kg C/m2, since aboveground C gains were more than offset by belowground C losses to decomposition in the tundra sites. Tree line advance therefore might not increase C storage in high-latitude ecosystems and thus might not, as previously suggested, act as a negative feedback to warming. Key to this hypothesis and to its projection to future climate response is the fate of soil carbon upon warming and permafrost drainage. Copyright 2006 by the American Geophysical Union.

Journal of Geophysical Research G: Biogeosciences

Soil N and 15N variation with time in a California annual grassland ecosystem

The %N and ??15N values of soils and plants were measured along a chronosequence spanning 3 to 3000 Ky in a California annual grassland. Total soil N decreased with increasing soil age (1.1 to 0.4 kg N m-2) while the mean ?? 15N values of the soil N increased by several ??? from the youngest to oldest sites (+3.5 to +6.2 ???). The ?? 15N values of plants varied along the gradient, reflecting changing soil N pools and differences in the form of N uptake. The decline in total N storage with time is hypothesized to be due to a shift from N to P limitation with increasing soil age. The general increase in ?? 15N values with time is interpreted using a N mass balance model, and appears to reflect a shift toward an increasing proportional losses of inorganic mineral forms of N (vs. organic forms) with increasing soil age. We develop a quantitative index of this trend (mineral vs. organic forms of N loss) using mass balance considerations and parameters. The %N and ?? 15N values along the California age gradient were compared to the published data for a comparably aged chronosequence in Hawaii. Most striking in this comparison is the observation that the California soil and plant ?? 15N values are several ??? greater than those on comparably aged Hawaiian sites. Multiple explanations are plausible, but assuming the sites have a similar range in ?? 15N values of atmospheric inputs, the isotopic differences suggest that N may be, at least seasonally, in greater excess in the strongly seasonal, semi-arid, California grassland. Copyright ?? 2001 Elsevier Science Ltd.

Geochimica et Cosmochimica Acta

Chemical weathering of a soil chronosequence on granitoid alluvium: II. Mineralogic and isotopic constraints on the behavior of strontium

The use of strontium isotopes to evaluate mineral weathering and identify sources of base cations in catchment waters requires an understanding of the behavior of Sr in the soil environment as a function of time. Our approach is to model the temporal evolution of 87 Sr/ 86 Sr of the cation exchange pool in a soil chronosequence developed on alluvium derived from central Sierra Nevada granitoids during the past 3 Ma. With increasing soil age, 87 Sr/ 86 Sr of ammonium-acetate extractable Sr initially decreases from values typical of K-feldspar to those of plagioclase and hornblende and then remains constant, even though plagioclase and hornblende are absent from the soils after approximately 1 Ma of weathering. The temporal variation of 87 Sr/ 86 Sr of exchangeable Sr is modeled by progressively equilibrating Sr derived from mineral weathering and atmospheric deposition with Sr on exchange sites as waters infiltrate a soil column. Observed decreases in quartz-normalized modal abundances of plagioclase, hornblende, and K-feldspar with time, and the distinct 87 Sr/ 86 Sr values of these minerals can be used to calculate Sr flux from weathering reactions. Hydrobiotites in the soils have nearly constant modal abundances, chemistry, and 87 Sr/ 86 Sr over the chronosequence and provide negligible Sr input to weathering solutions. The model requires time and soil horizon-dependent changes in the amount of exchangeable Sr and the efficiency of Sr exchange, as well as a biologic cycling term. The model predicts that exchangeable Sr initially has 87 Sr/ 86 Sr identical to that of K-feldspar, and thus could be dominated by Sr leached from K-feldspar following deposition of the alluvium. The maximum value of 87 Sr/ 86 Sr observed in dilute stream waters associated with granitoids of the Yosemite region is likewise similar to that of the K-feldspars, suggesting that K-feldspar and not biotite may be the dominant source of radiogenic Sr in the streams. This study reveals that, when attempting to use Strontium isotopes to identify sources of base cations in catchment waters and biomass, both preferential leaching of Sr from minerals during incipient soil development and changing Sr exchange efficiency must be considered along with chemical contributions due to mineral dissolution.

Geochimica et Cosmochimica Acta

10Be distribution in soils from Merced River terraces, California

The distribution and residence time of cosmogenic 10 Be in clay-rich soil horizons is fundamental to understanding and modelling the migration of 10 Be on terrestrial sediments and in groundwater solutions. We have analyzed seven profiles of clay-rich soils developed from terrace sediments of the Merced River, California. The terraces and soils of increasing age are used to compare the 10 Be inventory with a simple model of accumulation, decay and erosion. The data show that the distribution of 10 Be varies with soil horizon clay content, that the residence time of 10 Be in these horizons exceeds 10 5 years, and that to a rough approximation the inventory of 10 Be in a thoroughly sampled soil profile fits the equation: N = ( q − Em )(1 − e − λι )/ λ where q is delivery rate, E is erosion rate, m is the concentration of 10 Be in the eroding surface layer, λ is the decay constant, and t is the age of the depositional unit from which the soil has developed. The general applicability of this model is uncertain and warrants further testing in well-calibrated terrace sequences.

Geochimica et Cosmochimica Acta