As humans continue to impact the Critical Zone, we need to project how our environment will evolve into the future. To model such change requires the ability to simulate interactions among the lithosphere, pedosphere, hydrosphere, biosphere, and atmosphere — including the activities of humans. Such projections, which some have called earthcasts, must be made with mechanistic models that capture the important phenomena, as well as scenarios of human behavior. As an example, we present earthcasts of future weathering in the mid-continent of the USA into the next century of projected warming. Rates of sequestration of CO2 from the atmosphere due to weathering will change in the future as carbonate and silicate minerals are dissolved or precipitated in soil. The downward or upward advance of the carbonate reaction front in the soil is an analogue of the oceanic lysocline. Like the movement of the oceanic lysocline in response to oceanic acidification, this terrestrial lysocline will likely move due to fluxes of CO2 driven by human activity. Understanding this and other responses to perturbations will best be achieved using multiple models for earthcasting.

Humans are changing Earth’s atmosphere, biota, water, soil, and sediments to the extent that a new geological epoch — the Anthropocene — has been proposed by Nobel laureate Paul Crutzen to emphasize the global and lasting impact of humans (Crutzen and Stoermer, 2000; Zalzsiewicz et al., 2008). Through processes such as deforestation, farming, urbanization, mining, and dam construction, humans transform water, vegetation, soils, and climate (Foster et al., 1998; Boyer et al., 2002; Pielke et al., 2002; McNeill and Winiwarter, 2004; Bond et al., 2005; Wilkinson, 2005; Brikowski, 2008; Walter and Merritts, 2008; van Mantgem et al., 2009). The pace of anthropogenic change threatens to overwhelm humankind’s ability to plan sustainable use of soils and waters (Vitousek et al., 1997; DeFries and Eschleman, 2004; McNeill and Winiwarter, 2004; Wilkinson, 2005; McNeill and Winiwarter, 2006; Holdren, 2008; Mann and Kump, 2008). Quantitative approaches are needed to project how these earth resources evolve naturally and under human impact.

Just as atmospheric scientists use general circulation models (GCMs) to compare multiple climate simulations, Earth scientists are developing and comparing models for forecasting —“Earthcasting”—the behavior of the zone from vegetation canopy to groundwater (Bachlet et al., 2003; Qu and Duffy, 2007; Minasny et al., 2008; Rasmussen et al., 2009). This realm, called the Critical Zone, represents the interface where rocks meet life (U.S. National Research Council Committee on Basic Research Opportunities in the Earth Sciences, 2001; Brantley et al., 2006; Brantley et al., 2008).

When this zone is observed over short timeframes, the response to perturbation is often linear, but over long timeframes, nonlinear responses emerge, because the system is driven by many coupled processes, each acting with its own characteristic response time (Chadwick and Chorover, 2001; Bachlet et al., 2003). To understand these responses, therefore, both short and long time records of changes to the Critical Zone must be investigated and model simulations must be refined. These models are extraordinarily difficult to develop and explore because they must include earth materials, fluids, and biota over extended lengths of time. Perhaps most daunting, the models must also incorporate scenarios of human behavior.

Soils and deeper weathering profiles play an important role in allowing scientists to decipher the behavior of the Critical Zone because they integrate reactions among rock, water, air, and biota over time (Figure 1) (Richter, 2007; Chesworth, 2008; Minasny et al., 2008). Indeed, it has long been known that the chemistries and textures of soils record the signature of environmental fluxes and changing conditions in response to tectonism, climate, and human activity (Jenny, 1980; Yaalon, 1983; Richter and Markewitz, 2001; Warkentin, 2006). Models can be used to interpret the stories written in soils by both natural and anthropogenic forcings to anticipate future changes to the Critical Zone (Minasny et al., 2008; Brantley and Lebedeva, 2011). These models should become increasingly useful as earthcasts now that the anthropogenic perturbation of the carbon fluxes from land to ocean have been clearly identified (Regnier et al., 2013).

Figure 1
Humans are acting as a geological force by perturbing flows of water, solutes, gases, and sediments.

The story of change of these flows is recorded in soils and sediment. When the complex surface earth system is perturbed, both linear and nonlinear “tipping point” responses result. Models, needed to project such responses, must first be tested against the “stories” written in the soil and sediment record. (Figure, inspired by Chesworth (2008), is reprinted here with permission from the original, Brantley and Lebedeva, 2011).

Figure 1
Humans are acting as a geological force by perturbing flows of water, solutes, gases, and sediments.

The story of change of these flows is recorded in soils and sediment. When the complex surface earth system is perturbed, both linear and nonlinear “tipping point” responses result. Models, needed to project such responses, must first be tested against the “stories” written in the soil and sediment record. (Figure, inspired by Chesworth (2008), is reprinted here with permission from the original, Brantley and Lebedeva, 2011).

Close modal

One example of an important process that must be earthcast is the process whereby minerals in soils weather and consume atmospheric CO2 (Banwart et al., 2009; Berner, 2006). During weathering, rock materials are catalyzed by bacteria to react with water and gas, equilibrating minerals to earth surface conditions. As minerals weather, CO2 from the atmosphere is solubilized as carbonic acid which reacts with minerals producing bicarbonate ions that are either precipitated as carbonates or transported in soil porewaters to rivers and ultimately the ocean. Over geological timescales, solubilization of CO2 during silicate weathering balances the production of CO2 to the atmosphere during volcanism and metamorphism, stabilizing the Earth climate within limits accessible for life (Walker et al., 1981; Kump et al., 2000; Berner, 2006).

Today, continental silicate and carbonate weathering together consume about 0.3 Gt of atmospheric CO2 each year (Gaillardet et al., 1999). This value is close to the preindustrial net carbon exchange between the atmosphere and the terrestrial biosphere (0.4 Gt C/year) (Solomon et al., 2007). The rate of continental weathering is sensitive to land use change (Raymond et al., 2008; Perrin et al., 2008) and ongoing climate change (Gislason et al., 2009; Beaulieu et al., 2012). As humans increase both the concentration of CO2 in the atmosphere as well as their impact on the continental ecosystems, the rate of weathering will respond. Quantifying how C fluxes will change in the future is a complex task for models because of the coupling between climate, hydrology, land use, and biogeochemical reactions, and because the models must be able to simulate a fast-changing system under transient conditions. In addition to its potential role in the global carbon cycle, continental weathering is also a key process in the global cycles of the major cations (Ca2+, Mg2+, K+, Na+) and thus contributes to the control of the future evolution of soil and riverine pH. Finally, at the centennial timescale, the continental biosphere derives nourishment from the elements released during weathering reactions. Despite its non-negligible magnitude, apparent sensitivity, and acknowledged importance in the biogeochemical cycle of many elements, continental weathering is generally neglected in studies dealing with the evolution of the global carbon cycle and climate over century timescales.

Considering the future of weathering is a new point of view for Earth scientists. Up to now, the science of weathering mostly focused on processes in the past. Indeed, weathering has been a key process during the geological evolution of the Earth system (Walker et al., 1981; Berner, 2006). We use earthcasting here to explore whether it might also be a key process for the short term future of our planet.

Parametric laws have been developed to predict the future evolution of weathering from mean annual temperature and runoff for several average lithologies (Hartmann et al., 2009). Such laws, generally derived from field observations, are chosen empirically and are used to assess integrated weathering fluxes for large continental surfaces without excessive computational time. As pointed out in Goddéris et al. (2009), this parametric method can be a reasonable snapshot of the current weathering system. However, parametric approaches are not necessarily adapted to describing decadal to centennial dynamics. In contrast, mechanistic models that incorporate physical descriptions of the below-ground hydrology and weathering reactions can be used to predict the response of weathering to increasing CO2 in the future.

Mechanistic models have been used and tested by focusing on natural experiments where climate and time are variable. For example the following phenomena have been modeled that are relevant to the goal of earthcasting: (1) mineralogical evolution of loess that was deposited approximately 13 ky ago along the Mississippi valley transect in the U.S.A (Williams et al., 2010; Goddéris et al., 2010), (2) soil development at the 226 ka marine terrace chronosequence near Santa Cruz, U.S.A. (Maher et al., 2009), and (3) soil development over a million years of weathering in the Merced chronosequence (Moore et al., 2012). These studies have shown that mechanistic weathering models can successfully reproduce mineralogical patterns observed in soils along a climate gradient (Williams et al., 2010; Goddéris et al., 2010; Goddéris et al., 2013) or along a gradient in time (Maher et al., 2009; Moore et al., 2012).

This approach documents that we can mechanistically model weathering as a function of time and changing climate. Given the current constraints of computational time, it is not possible to fully couple climate and reactive transport (weathering models). Instead, here we summarize generalizable lessons learned from those modeling efforts to highlight challenges for the future.

In our approach, we use a cascade of numerical models for climate (ARPEGE), vegetation (CARAIB) and weathering (WITCH), to explore the effect of an increase in CO2 into the future. Specifically, we simulate one of the carbon emission scenarios from the Intergovernmental Panel on Climate Change (Solomon et al 2007). The scenario we chose, A1B, assumes very rapid economic growth and a peak in global population mid-century, after which the population declines while new and efficient technologies are rapidly introduced to achieve balance across energy sources. With this scenario, CO2 increases from 315 to 700 ppmv from 1950 to 2100. Other models can be similarly used in a cascade to project weathering into the future (Figure 2).

Figure 2
To simulate rates of mineral weathering into the future requires a cascade of models as shown.

Climate models (GCM for General Circulation Model) in use today include GENESIS and ARPEGE, among others; vegetation models (DGVM for Dynamic Global Vegetation Model) include BIOME and CARAIB among others; weathering models include WITCH, SAFE, CRUNCH, and FLOTRAN among others. The solid arrows stand for processes that are currently included in numerical models. Dashed arrows represent processes and feedbacks that are not yet modeled within our current Earthcasting efforts.

Figure 2
To simulate rates of mineral weathering into the future requires a cascade of models as shown.

Climate models (GCM for General Circulation Model) in use today include GENESIS and ARPEGE, among others; vegetation models (DGVM for Dynamic Global Vegetation Model) include BIOME and CARAIB among others; weathering models include WITCH, SAFE, CRUNCH, and FLOTRAN among others. The solid arrows stand for processes that are currently included in numerical models. Dashed arrows represent processes and feedbacks that are not yet modeled within our current Earthcasting efforts.

Close modal

Factors that affect weathering

The reason that earthcasting is difficult is that there are so many processes and factors influencing the future of biogeochemical reactions. To illustrate this complexity, we focus here on two key forcing functions of the weathering reactions: the ambient temperature, which controls the kinetics and thermodynamics of the weathering reactions, and the soil drainage, which controls the transport of the weathering products and the availability of water required for weathering. The response of each of these forcing functions to the anthropogenic rise in CO2 can be contrasted.

  1. On average, increasing CO2 results in an increase in temperature globally, promoting silicate weathering and inhibiting carbonate dissolution.

  2. In contrast, increasing CO2 can cause either an increase or a decrease in porefluid drainage through soils, depending upon location at the Earth surface. Drainage can increase or decrease due to several competing effects, and these effects, which cannot be simulated by parametric models, require sophisticated climate and vegetation projections. For example, rainfall may either decrease or increase locally, depending upon how the atmospheric global circulation changes. In addition, however, CO2 fertilization will increase productivity of land plants which will in turn result in enhanced evapotranspiration, yielding less drainage. Evapotranspiration is also enhanced by temperature increases. Conversely, closure of leaf stomata at high CO2 levels may ultimately decrease evapotranspiration. It is impossible to know for any one location what the effect of all these factors will be without a numerical simulation.

Learning from our experience simulating weathering of loess along the Mississippi valley in the U.S.A. over the last 13 ky, we have run earthcasts to simulate weathering profiles along this transect out to the year 2100. We calculated that the drainage out of the profiles in the northern section of the transect in the state of Illinois (42°59.00′, 91°12.66′) could decrease by as much as 40% owing to the projected rise in temperature and evapotranspiration (Figure 3). This decreased drainage is predicted despite a moderate rise in projected rainfall. Conversely, the drainage in the southern end of the transect, in the state of Mississippi (32°24.46′, 90°49.30′), is calculated to increase (simulation not shown). Based on our observations of models for the future along the Mississippi valley, the hydrology and the temperature are the two variables that largely determine the earthcast of weathering (Figure 3).

Figure 3
Earthcasts showing how the sequestration of atmospheric CO2 will change into the future through the weathering of the northernmost pedon of the Mississippi transect.

The top plot shows contours of mean annual bicarbonate drainage flux, plotted on axes of depth (vertical) and time (horizontal), calculated by the ARPEGE-CARAIB-WITCH model cascade. The superimposed white curve is the calculated HCO3- drainage flux at the basis of the pedon, as a function of time (scale on the right). The green line is the linear fit of the model output, documenting the decrease in the bicarbonate flux in the future. The reaction front for dolomite, which is like a terrestrial lysocline, lies at about 3 meters depth, above which no dolomite is present. High fluxes of bicarbonate produced by silicate dissolution are depicted as light blue vertical pulses in the upper parts of the soil, above 3 meters depth. High vertical fluxes are predicted during wet periods. Although silicates dissolve slightly faster into the future, most of the bicarbonate flux out of the soil derives from dolomite dissolution below 3 meters depth, which slows into the future as temperature increases. The three bottom panels show simulations of the soil temperature, vertical drainage, and the CO2 level below the root zone (in Pre-industrial Atmospheric Level units: 1 PAL=280 ppmv).

Figure 3
Earthcasts showing how the sequestration of atmospheric CO2 will change into the future through the weathering of the northernmost pedon of the Mississippi transect.

The top plot shows contours of mean annual bicarbonate drainage flux, plotted on axes of depth (vertical) and time (horizontal), calculated by the ARPEGE-CARAIB-WITCH model cascade. The superimposed white curve is the calculated HCO3- drainage flux at the basis of the pedon, as a function of time (scale on the right). The green line is the linear fit of the model output, documenting the decrease in the bicarbonate flux in the future. The reaction front for dolomite, which is like a terrestrial lysocline, lies at about 3 meters depth, above which no dolomite is present. High fluxes of bicarbonate produced by silicate dissolution are depicted as light blue vertical pulses in the upper parts of the soil, above 3 meters depth. High vertical fluxes are predicted during wet periods. Although silicates dissolve slightly faster into the future, most of the bicarbonate flux out of the soil derives from dolomite dissolution below 3 meters depth, which slows into the future as temperature increases. The three bottom panels show simulations of the soil temperature, vertical drainage, and the CO2 level below the root zone (in Pre-industrial Atmospheric Level units: 1 PAL=280 ppmv).

Close modal

Of course, as the climate and CO2 in the atmosphere change into the future, the CO2 concentrations in soil solutions will also vary (Figure 3). Like the drainage, soil pCO2 is also difficult to project without numerical models. As temperature rises, the CO2 production rate by biota increases in soils, releasing more CO2 into the soil atmosphere. However, in response to global warming, the rate of diffusion back to the atmosphere also increases, as do rates of mineral weathering and formation of aqueous bicarbonate that drains to rivers. Thus, it is possible that the CO2 concentration in the soil atmosphere will either increase, remain roughly constant, or even decrease over the course of the next century depending upon location (Goddéris et al., 2013). For our model simulations in the northern pedon shown in Fig. 3, the soil pCO2 decreases into the future.

In addition to these more-or-less smooth changes, more pronounced modifications may happen when the global atmosphere attains higher CO2 levels in the future. For example, Beaulieu et al. (2010) showed that when atmospheric CO2 approaches 4 times the pre-industrial level, the changing climate may force rapid drifting of biomes. As biomes change, the below-ground hydrology or the vegetation productivity can change drastically. With mechanistic models, all of these can be modeled. In addition, we could make other earthcasts using different scenarios of human activity to further explore how land use change could affect weathering response.

Our simulations for the Mississippi transect predict that the temperature increase expected in the next 100 years will cause the rate constants and solubilities of silicate minerals to increase to the extent that weathering rates will increase regardless of whether drainage increases or decreases along the transect. As temperature increases during the period 1950–2100, the reaction front for the dissolution of silicates moves at increasingly faster rates of advance along the whole transect.

However, for the loess, the slow-dissolving silicates only sequester a minor amount of CO2 over the next 100 y. What matters in this geologically short timeframe in terms of carbon uptake by weathering is the dissolution of carbonate minerals. This is because the dissolution rates for carbonates are about 10 orders of magnitude faster than the dissolution rates of most common silicate minerals (Brantley et al., 2008). Of course, the net carbon consumption due to weathering of carbonates is equal to zero over 105 to 106 year timeframes since CO2 is returned to the atmosphere as carbonates precipitate at the ocean floor to release the consumed carbon back to the atmosphere (Berner, 2004). However, dissolution of continental carbonates nonetheless removes carbon from the atmosphere and stores it in the ocean for thousands of years before precipitation and degassing. In contrast to silicate minerals, the weathering rate of carbonate minerals such as calcite and dolomite decreases with increasing temperature because carbonate minerals have retrograde solubility. Retrograde solubility means that for a given pCO2, calcite solubility (for example) decreases with increasing temperature.

Although only tiny amounts of dolomite are found in the loess soils, the uptake of CO2 from the soil atmosphere as they weather dominates the production of bicarbonate in soil porewaters along the transect (Figure 3). As bicarbonate forms and CO2 is pulled out of the soil atmosphere, the dissolution of dolomite defines a reaction front in these slowly eroding soils that slowly moves downward with time. Given the relative importance of all of the competing factors, we project slower rates of advance of the carbonate reaction front into the subsurface for both the southernmost and northernmost pedons of the Mississippi transect. The slowing advance of the carbonate front is exemplified in a map of the concentrations of bicarbonate in loess soil porefluids for the northernmost pedon (Figure 3).

Terrestrial lysocline

The deceleration of the advance rate of the carbonate reaction front is caused by the slow response of soils to the perturbation of higher CO2 and changing climate. In fact, the carbonate reaction front can be likened to a terrestrial lysocline — it represents a depth interval over which carbonate dissolution rates on continents change rapidly. In the ocean, the lysocline represents the transition across which calcite reaction rates change drastically from a surface zone of calcite supersaturation to a deeper zone of calcite undersaturation. As the surface ocean responds to higher atmospheric CO2 and acidifies, the oceanic lysocline will shallow (Boudreau et al., 2010). In contrast, our earthcasts project a deeper lysocline (a deeper carbonate reaction front) in future soils. In the soils we simulated, calcite is undersaturated above the lysocline but supersaturated below. The advance of the terrestrial lysocline, like the acidification of the ocean, is the response to perturbation of increased CO2 concentration in the atmosphere.

On soils eroding at slow rates such as the Mississippi loess pedons, this lysocline deepens with time. The question is, how will deepening of the terrestrial lysocline evolve into the future? The slow deceleration of advance of the lysocline shown in Figure 4 does nothing to counteract rising CO2. An acceleration in deepening will be necessary for weathering systems to produce more alkalinity so as to oppose the rise in atmospheric CO2. In fact, faster rates are predicted in the central part of the Mississippi transect (Godderis et al., 2013). Only in that part of the transect does the model cascade project that drainage increases enough to counteract the temperature-induced inhibition of carbonate dissolution. Thus it is important to know how much of the soil worldwide will act like the central Mississippi Valley, where we calculate that the deepening of the soil lysocline will accelerate, versus how much will act like the Southern and Northern pedons where the deepening rate slows (Godderis et al., 2013). Deceleration is particularly striking in the North where the temperature rise and the drainage decrease are maximized. Figure 4 compares the average terrestrial lysocline deepening that we calculated over the entire Holocene due to natural climatic forcing (Goddéris et al., 2010) versus the calculated evolution under the next 100 years of anthropogenic forcing (Godderis et al., 2013). The figure shows that the rate of deepening of the lysocline is slowing below the “natural” rates, especially after 2040.

Figure 4
Rate of deepening of the terrestrial lysocline for the nothernmost pedon, calculated by a cascade of climate-biosphere-weathering models.

The grey area represents the calculated long term (kyr) average rate over the Holocene (Goddéris et al., 2010). The blue line is the future projection of the deepening rate. As the climate gets warmer and drier under the anthropogenic forcing, carbonate dissolution rate goes down, and the retreat rate of the terrestrial lysocline slows down together with the atmospheric CO2 consumption by weathering.

Figure 4
Rate of deepening of the terrestrial lysocline for the nothernmost pedon, calculated by a cascade of climate-biosphere-weathering models.

The grey area represents the calculated long term (kyr) average rate over the Holocene (Goddéris et al., 2010). The blue line is the future projection of the deepening rate. As the climate gets warmer and drier under the anthropogenic forcing, carbonate dissolution rate goes down, and the retreat rate of the terrestrial lysocline slows down together with the atmospheric CO2 consumption by weathering.

Close modal

The projections obtained for the whole transect illustrate the complex behavior of carbonate weathering in the face of short term global climate change. Predicting the global response of terrestrial weathering to increased atmospheric CO2 and temperature in the future will mostly depend upon our ability to make precise assessments of which areas of the globe increase or decrease in precipitation and soil drainage.

Here we have promoted the concept of a terrestrial lysocline for loess soils — essentially, the carbonate weathering front. This front deepens as carbonate dissolution proceeds in slowly eroding soils. Can we generalize this concept to a larger variety of continental soils? In drier areas, intense evaporation causes the formation of pedogenic carbonates (McFadden et al., 1991). In cases where dryness gets more pronounced in the future, precipitation of carbonate minerals will occur in increasingly shallower zones. The terrestrial lysocline then moves upward, with the result of a net release of CO2 to the atmosphere. In permafrost areas, the warming will increase the thaw depth, and consequently the depth of carbonate weathering. The warming may thus potentially impact the dynamics of the terrestrial lysocline in these regions (Keller et al., 2010). We must struggle to understand where on earth the terrestrial lysocline deepens and where it will shallow if we are to understand and project our environment into the next 100 y.

The model simulations also teach us that the weathering of carbonates responds quickly to climate change and will thus be important in the near future. In fact, the very high reactivity of carbonates means that their dissolution is mainly controlled by advective transport. For such reactions limited by solute transport, projections must rely on accurate models of water fluxes. For soils dominated by carbonate precipitation instead of dissolution (i.e., soils in drier climates), we will furthermore need more accurate knowledge of carbonate precipitation kinetics close to equilibrium (Schott et al., 2009).

Erosion

So far we have simulated soils that are not actively eroding at rates that are important within the 100-year timeframe. In a non-eroding soil, the terrestrial lysocline moves with time. The effect of a CO2 increase is to change the rate of advance of that front. In contrast, in an eroding soil, the depth of the terrestrial lysocline can remain constant if the rate of advance of the carbonate reaction front is equivalent to the rate of erosive loss of material. Parametric models have also been used to begin to simulate such soils as well (Brantley et al., 2013).

However, rates of soil erosion have increased by 30x above natural background due to anthropogenic processes (Wilkinson, 2005; Wilkinson and McElroy, 2007). Today, the rate at which sediment and soil particles are transported to the oceans is approximately 120 x 109 tons/yr. This is 24 times the natural background rate of 5 x 109 tons/yr. Like many of the other coupled processes in the water-soil-gas-biota system, a strong link couples physical erosion to the weathering of silicate minerals (Millot et al., 2002; West et al., 2005). Efficient dissolution of silicate mineral particles in seawater may be further evidence of the important link between physical erosion and chemical weathering (Jeandel et al., 2011; Jones et al., 2013). The significance of this link has been explored from a theoretical point of view within the framework of the geological evolution of the Earth surface (Gabet and Mudd, 2009; West, 2012). At shorter timescales, physical erosion might be important in accelerating chemical weathering of the slow-dissolving silicate minerals by increasing their surface area. In addition, intense erosion also strongly modifies the hydrologic behavior of the weathering profiles. If erosion globally outcompetes the advance of weathering and concomitant production of soil (Wilkinson, 2005), soils will become thinner and the CO2 and elemental depth profiles in soil will change. The impact on the global weathering rates remains to be explored.

In contrast to erosion that removes materials, the deposition and subsequent dissolution of calcareous dust transported over long distances can also modify the terrestrial lysocline advance rates in the future, slowing its downward shift. In the future, because large areas are expected to become drier, the pattern of dust deposition over the continents will be altered, with plausible consequences on the weathering system.

The consumption of CO2 as a function of depth into the subsurface is not the only gas reaction that is important in this context. Consumption of O2 at depth as sulfides and organic matter are oxidized also create reaction fronts in the subsurface that affect global mass balances (Kump, 2008; Brantley et al., 2013). Physical erosion may also accelerate the rate that reduced minerals come into contact with the oxidizing atmosphere and consume O2. In addition, for the specific case of the Mackenzie River (Canada), the oxidation of pyrite releases sulfuric acids into the environment, dissolving minerals such as carbonates without consuming atmospheric CO2 (Spence and Telmer, 2005; Calmels et al., 2007). The global warming of the Mackenzie watershed and melting of the permafrost may ultimately increase physical erosion, sulfuric acid release, and ultimately cause CO2 degassing from the ocean to the atmosphere (Calmels et al., 2007). Once again, only mechanistic models can be relied upon to explore the range of possibilities.

Today we rely upon forecasts of the weather to make daily decisions. Tomorrow we may be able to make longer range decisions by relying on earthcasts of geologic properties of our surface environment. As discussed here, we are already earthcasting the rates of mineral change in soils. Our earthcasts show that over the next 100 y, weathering will respond into the future mostly as a function of changes in temperature and drainage. Only to a much lesser extent will the soil atmospheric CO2 affect mineral reactions.

In general, as humans pump up the concentration of CO2 in the atmosphere, temperature will increase but drainage and soil CO2 could either decrease or increase in any given location. Carbonate weathering is highly sensitive to these future environmental changes over decadal to centennial timescales. Like the oceans where a depth zone of transition separates carbonate supersaturation from undersaturation, a terrestrial lysocline or carbonate reaction front separates zones of saturation in soils. We expect that the oceanic lysocline will shallow into the future while the terrestrial lysocline will continue to deepen in temperate systems where erosion is slow. However, the rate of advance of the terrestrial lysocline could either decelerate or accelerate in the future, mostly depending upon changes in drainage. Projections of alkalinity production due to movement of the terrestrial lysocline are thus highly dependent on changes in drainage, and hence on climate and vegetation dynamics which are inherently difficult to predict. Likewise, carbonates may be both dissolved and precipitated in soils depending on aridity. The relative importance of factors controlling pedogenic carbonates must be better understood.

Parametric models of weathering may work for the short term but are not mechanistically based and cannot yield accurate projections of weathering into the future. The best tool we have for exploring the future of weathering is a cascade of models. Such a cascade includes a climate model that feeds a dynamic global vegetation model, which in turn feeds a weathering model. But no single such cascade of models is expected to always yield accurate simulations. Just as we compare multiple models when we try to predict the weather, we need multiple cascades of earthcasting models that we can compare and improve to simulate future environmental systems. Ultimately, these model cascades will be replaced by truly coupled models. Coupling will be especially critical at the millenial timescale where continental vegetation is nourished by elements released during weathering. Models that quantify biotic processes and changes in soil organic matter are sorely needed. Incorporating the impact of soil erosion into future models should also be a priority as anthropogenic processes are accelerating the rates of erosion worldwide by 1 to 2 orders of magnitude. Furthermore, no model can predict the rate of production of erodible material — i.e. the production of mineral soil. The large range of possible responses of natural soil systems to human perturbations requires development and testing of increasingly sophisticated and mechanistically realistic models by multiple groups worldwide.

© 2013 Goddéris, Brantley. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

SLB acknowledges funding from U.S. Office of Basic Energy Sciences, Department of Energy grant DE-FG02-OSER15675 and from the U.S. National Science Foundation Critical Zone Observatory program grant EAR 12-39285. YG acknowledges funding from the French National Agency for Research, grant ANR-09-BLAN-0221.

Bachlet
D
,
Neilson
RP
,
Hickler
T
,
Drapek
RJ
,
Lenihan
JM
, et al
2003
.
Simulating past and future dynamics of natural ecosystems in the United States
.
Global Biogeochemical Cycles
17
:
1045
1053
.
Banwart
SA
,
Berg
A
,
Beerling
DJ
.
2009
.
Process-based modeling of silicate mineral weathering responses to increasing atmospheric CO2 and climate change
.
Global Biogeochemical Cycles
23
:
1
13
.
Beaulieu
E
,
Goddéris
Y
,
Donnadieu
Y
,
Labat
D
,
Roelandt
C
.
2012
.
High sensitivity of the continental weathering carbon dioxide to future climate change
.
Nature Climate Change
2
:
346
349
.
Beaulieu
E
,
Goddéris
Y
,
Labat
D
,
Roelandt
C
,
Oliva
P
, et al
2010
.
Impact of atmospheric CO2 levels on continental silicate weathering
.
Geochemistry, Geophysics, Geosystems
11
. doi: .
Berner
RA
.
2004
.
The Phanerozoic Carbon Cycle
.
New York
:
Oxford University Press
.
Berner
RA
.
2006
.
GEOCARBSULF: A combined model for Phanerozoic atmospheric O2 and CO2
.
Geochimica et Cosmochimica Acta
70
:
5653
5664
.
Bond
WJ
,
Woodward
FI
,
Midgley
GF
.
2005
.
The global distribution of ecosystems in a world without fire
.
New Phytologist
165
:
525
538
.
Boyer
EW
,
Goodale
CL
,
Jaworski
NA
,
Howarth
RW
.
2002
.
Anthropogenic nitrogen sources and relationships to riverine nitrogen export in the northeastern U.S.A
.
Biogeochemistry
57
:
137
169
.
Brantley
SL
,
Lebedeva
M
.
2011
.
Learning to read the chemistry of regolith to understand the Critical Zone
.
Annual Review Earth Planetary Science
39
:
387
416
.
Brantley
SL
,
Holleran
M
,
Jin
L
,
Bazilevskaya
E
.
2013
.
Probing deep weathering in the Shale Hills Critical Zone Observatory, Pennsylvania (U.S.A.): The hypothesis of nested chemical reaction fronts in the subsurface
.
Earth Surface Processes and Landforms
38
:
1280
1298
.
Brantley
SL
,
Kubicki
JD
,
White
AF
.
2008
.
Kinetics of Water-Rock Interaction
.
New York
:
Springer
.
Brantley
SL
,
White
TS
,
Ragnarsdottir
KV
.
2008
.
The critical zone: Where rock meets life
.
Elements
3
:
297
368
.
Brantley
SL
,
White
TS
,
White
AF
,
Sparks
D
,
Richter
D
, et al
2006
.
Frontiers in Exploration of the Critical Zone
,
An NSF-sponsored Workshop
Brantley
SL
,
Lebedeva
M
,
Bazilevskaya
E
.
2013
. Relating weathering fronts for acid neutralization and oxidation to pCO2 and pO2 . In
Farquhar
J
,
Kasting
J
, and
Candfield
D
eds,
Treatise in Geochemistry, The Atmosphere-History
.
Amsterdam, The Netherlands
:
Elsevier
: in press.
Brikowski
TH
.
2008
.
Doomed reservoirs in Kansas, USA? Climate change and groundwater mining on the Great Plains lead to unsustainable surface water storage
.
Journal of Hydrology
354
:
90
101
.
Calmels
D
,
Gaillardet
J
,
Brenot
A
,
France-Lanord
C
.
2007
.
Sustained sulfide oxidation by physical erosion processes in the Mackenzie River basin: climatic perspectives
.
Geology
35
:
1003
1006
.
Chadwick
OA
,
Chorover
J
.
2001
.
The chemistry of pedogenic thresholds
.
Geoderma
100
:
321
353
.
Chesworth
W
.
2008
. Biogeochemical Cycles. In
Chesworth
W
. ed,
Encyclopedia of Soil Science
. Vol , pp
56
60
.
Dordrecht, The Netherlands
:
Springer
.
Crutzen
PJ
,
Stoermer
EF
.
2000
.
The “Anthropocene”
.
Global Change Newsletter
41
:
17
18
.
DeFries
R
,
Eschleman
NK
.
2004
.
Land-use change and hydrologic processes: a major focus for the future
.
Hydrological Processes
18
:
2183
2186
.
Foster
DR
,
Motzkin
G
,
Slater
B
.
1998
.
Land-use history as long-term broad-scale disturbance: Regional forest dynamics in central New England
.
Ecosystems
1
:
96
119
.
Gabbet
EJ
,
Mudd
SM
.
2009
.
A theoretical model coupling chemical weathering rates with denudation rates
.
Geology
37
:
151
154
.
Gislason
SR
,
Oelkers
EH
,
Eiriksdottir
ES
,
Kardjilov
MI
,
Gisladottir
G
, et al
2008
.
Direct evidence of the feedback between climate and weathering
.
Earth and Planetary Science Letters
277
:
213
222
.
Goddéris
Y
,
Brantley
SL
,
François
LM
,
Schott
J
,
Pollard
D
, et al
2013
.
Rates of consumption of atmospheric CO2 through the weathering of loess during the next 100 yr of climate change
.
Biogeosciences
10
:
135
148
.
Goddéris
Y
,
Roelandt
C
,
Schott
J
,
Pierret
M-C
,
François
LM
.
2009
.
Towards an integrated model of weathering, climate, and biospheric processes
.
Reviews in Mineralogy and Geochemistry
70
:
411
434
.
Goddéris
Y
,
Williams
JZ
,
Schott
J
,
Pollard
D
,
Brantley
SL
.
2010
.
Time evolution of the mineralogical composition of the Mississippi Valley loess over the last 10 kyr: climate and geochemical modeling
.
Geochimica et Cosmochimica Acta
74
:
6357
6374
.
Hartmann
J
,
Jansen
N
,
Dürr
HH
,
Kempe
S
,
Köhler
P
.
2009
.
Global CO2-consumption by chemical weathering: what is the contribution of highly active weathering regions?
Global and Planetary Change
69
:
185
194
.
Holdren
JP
.
2008
.
AAAS Presidential Address: Science and technology for sustainable well-being
.
Science
319
:
424
434
.
Jeandel
C
,
Peucker-Ehrenbrink
B
,
Jones
MT
,
Pearce
CR
,
Oelkers
EH
, et al
2011
.
Ocean margins: the missing term in oceanic element budgets?
EOS Transactions AGU
95
:
217
.
Jenny
H
.
1980
.
The Soil Resource: Origin and Behavior
.
New York
:
Springer-Verlag
.
Jones
MT
,
Pearce
CR
,
Jeandel
C
,
Gislason
SR
,
Eiriksdottir
ES
, et al
2013
.
Riverine particulate material dissolution as a significant flux of strontium to the oceans
.
Earth and Planetary Science Letters
355–356
:
51
59
.
Keller
K
,
Blum
JD
,
Kling
GW
.
2010
.
Stream geochemistry as an indicator of increasing permafrost thaw depth in an arctic watershed
.
Chemical Geology
273
:
76
81
.
Kump
LR
.
2008
.
The rise of atmospheric oxygen
.
Nature
451
:
277
278
.
Kump
LR
,
Brantley
SL
,
Arthur
MA
.
2000
.
Chemical weathering, atmospheric CO2 and climate
.
Annual Reviews of Earth and Planetary Sciences
28
:
611
667
.
McFadden
LD
,
Amundson
RG
,
Chadwick
OA
.
1991
.
Numerical modeling, chemical and isotopic studies of carbonate accumulation in soils of arid regions
.
SSSA Special Publication
,
26
:
17
35
.
Maher
K
,
Steefel
CI
,
White
AF
,
Stonestrom
DA
.
2009
.
The role of reaction affinity and secondary minerals in regulating chemical weathering rates at the Santa Cruz Soil Chronosequence, California
.
Geochimica et Cosmochimica Acta
73
:
2804
2831
.
Mann
ME
,
Kump
LR
.
2008
.
Dire Predictions: Understanding Global Warming
.
London
:
Pearson Education
.
McNeill
JR
,
Winiwarter
V
.
2004
.
Breaking the sod: Humankind, history, and soil
.
Science
304
:
1627
1629
.
McNeill
JR
,
Winiwarter
V
.
2006
.
Soils and Societies Perspectives from Environmental History
.
Isle of Harris, UK
:
White Horse Press
.
Minasny
B
,
McBratney
AB
,
Salvador-Blanes
S
.
2008
.
Quantitative models for pedogenesis - A review
.
Geoderma
144
:
140
157
.
Moore
J
,
Lichtner
PC
,
White
AF
,
Brantley
SL
.
2012
.
Using a reactive transport model to elucidate differences between laboratory and field dissolution rates in regolith
.
Geochimica et Cosmochimica Acta
93
:
235
261
.
Perrin
A-S
,
Probst
A
,
Probst
J-L
.
2008
.
Impact of nitrogenous fertilizers on carbonate dissolution in small agricultural catchments: implications for weathering CO2 uptake at regional and global scales
.
Geochimica et Cosmochimica Acta
,
72
:
3105
3123
.
Pielke
RAS
,
Marland
G
,
Betts
RA
,
Chase
TN
,
Eastman
JL
, et al
2002
.
The influence of land-use change and landscape dynamics on the climate system: Relevance to climate-change policy beyond the radiative effect of greenhouse gases
.
Phil. Trans.: Mathematical, Physical, and Engineering Sciences
360
:
1705
1719
.
Qu
Y
,
Duffy
CJ
.
2007
.
A semidiscrete finite volume formulation for multiprocess watershed simulation
.
Water Resources Research
43
. doi: .
Rasmussen
C
,
Troch
PA
,
Chorover
J
,
Brooks
P
,
Pelletier
TE
, et al
2011
.
An open system framework for integrating Critical Zone structure and function
.
Biogeochemistry
102
:
15
29
.
Raymond
PA
,
Oh
N-H
,
Turner
RE
,
Broussard
W
.
2008
.
Anthropogenically enhanced fluxes of water and carbon from the Mississippi River
.
Nature
451
:
449
452
.
Regnier
P
,
Friedlingstein
P
,
Ciais
P
,
Mackenzie
FT
,
Gruber
N
, et al
2013
.
Anthropogenic perturbation of the carbon fluxes from land to ocean
.
Nature Geoscience
6
:
597
607
.
Richter
D
,
Markewitz
D
.
2001
.
Understanding soil change
.
U.K
.:
Cambridge University Press
.
Richter
DD
.
2007
.
Humanity’s transformation of Earth’s soil: a new pedological frontier
.
Soil Science
172
:
957
967
.
Schott
J
,
Pokrovsky
OS
,
Oelkers
EH
.
2009
. The Link Between Mineral Dissolution/Precipitation Kinetics and Solution Chemistry. In
Oelkers
EH
and
Schott
J
(ed)
Thermodynamics and Kinetics of Water-Rock Interaction. Reviews in Mineralogy and Geochemistry
70
:
207
258
.
Solomon
S
,
Qin
D
,
Manning
M
,
Chen
Z
,
Marquis
M
, et al
2007
.
Contribution of working group 1 to the fourth assessment report of the Intergovernmental Panel on Climate Change
,
Cambridge University Press
,
Cambridge, United Kingdom and New York, NY, USA
.
Spence
RE
,
Telmer
K
.
2005
.
The role of sulfur in chemical weathering and atmospheric CO2 fluxes: evidence from major ions, δ13CDIC, and δ34SSO4 in rivers of the Canadian Cordillera
.
Geochimica et Cosmochimica Acta
,
69
:
5441
5458
.
U.S. National Research Council Committee on Basic Research Opportunities in the Earth Sciences
,
2001
.
Basic Research Opportunities in Earth Science
.
Washington, D.C
.:
National Academy Press
.
vanMantgem
PJ
,
Stephenson
NL
,
Byrne
JC
,
Daniels
LD
,
Franklin
JF
, et al
2009
.
Widespread increase of tree mortality rates in the western United States
.
Science
323
:
521
524
.
Vitousek
PM
,
Mooney
HA
,
Lubchenco
J
,
Melillo
JM
.
1997
.
Human domination of Earth’s ecosystems
.
Science
277
:
494
499
.
Walker
JCG
,
Hays
PB
,
Kasting
JF
.
1981
.
A negative feedback mechanism for the long-term stabilization of Earth’s surface temperature
.
J. Geophys. Res.
,
86
:
9776
9782
.
Walter
RC
,
Merritts
DJ
.
2008
.
Natural streams and the legacy of water-powered mills
.
Science
18
:
299
304
.
West
AJ
.
2012
.
Thickness of the chemical weathering zone and implications for erosional and climatic drivers of weathering and for carbon-cycle feedbacks
.
Geology
,
40
:
811
814
.
West
AJ
,
Galy
A
,
Bickle
M
.
2005
.
Tectonic and climatic controls on silicate weathering
.
Earth and Planetary Science Letters
,
235
:
211
228
.
Wilkinson
BH
.
2005
.
Humans as geologic agents: A deep-time perspective
.
Geology
33
:
161
164
.
Wilkinson
BH
and
McElroy
BJ
.
2007
.
The impact of humans on continental erosion and sedimentation
.
Geological Society of American Bulletin
119
:
140
156
.
Williams
JZ
,
Bandstra
JZ
,
Pollard
D
,
Brantley
SL
.
2010
.
The temperature dependence of feldspar dissolution determined using a coupled weathering-climate model for Holocene-aged loess soils
.
Geoderma
,
156
:
11
19
.
Yaalon
DH
.
1983
. Climate, Time and Soil Development. In
Wilding
LP
,
Smeck
NE
, and
Hall
GF
ed,
Pedogenesis and Soil Taxonomy, part 1, Concepts and Interactions
. pp
233
251
.
Amsterdam
:
Elsevier, Developments in Soil Science
.
Zalzsiewicz
J
,
Williams
M
,
Smith
A
,
Barry
TL
,
Coe
AL
,
Bown
PR
,
Brenchley
P
,
Cantrill
D
,
Gale
A
,
Gibbard
P
,
Gregory
FJ
,
Hounslow
MW
,
Kerr
AC
,
Pearson
P
,
Know
R
,
Powell
J
,
Waters
C
,
Marshall
J
,
Oates
M
,
Rawson
P
, and
Stone
P
.
2008
.
Are we now living in the Anthropocene?
GSA Today
18
:
3
8
.

This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

Supplementary data