Chemical weathering occurs faster in hot wet climates such as Rainforests. Once rocks are broken down, the particles are carried downhill under the force of gravity in a process called mass movement. In a mass movement , rocks roll, slide, or free-fall downhill under the force of gravity. What are the common forms of mass movement? It may only be noticed as curved tree trunks, bent telegram posts and broken retaining walls on roadsides.
A habitat for several plants and animals. A storage bank for water and nutrients. Foundation of the world's food chain. Provides an anchor for plants. Parent Material : the underlying bedrock over which. Soils forms from a variety of parent. Biological Activity Organisms : The activities of. Topography relief of the land : Various sections of. Time : Soil takes a long time to develop and mature.
Organic Matter Humus - decomposed materials :. Mineral Nutrients contained in rocks. Gas Content Oxygen, Nitrogen, C0 2 etc. Moisture Water. A horizontal layer of a soil is called a horizon. Sources of discrepancy can be that total soil thickness may include an organic layer O horizon , for example, which is deposited on top of the weathered soil, or that aeolian dust deposition or alluvial deposits, that are disconnected from the parent material, may be significant, neither of which is directly connected to soil production from weathering.
At many of the sites discussed in the present study the correlation has already been shown to be robust Hunt and Ghanbarian, ; Yu and Hunt, b , c. Furthermore, it has also been shown Hunt and Ghanbarian, that the time dependence of the soil formation rates of these particular studies corresponds closely to the time dependence of chemical weathering from completely different sources White and Brantley, , and over time scales from years to millions of years. Determination that chemical weathering and soil formation rates are proportional across a range of field experiments lends support to the assumption that such a proportionality holds more universally, but does not constitute a proof that extension to all sites under all conditions will be valid.
Our present study applies systematically this recent model of soil formation derived from the chemical weathering depth as limited mainly by solute transport percolation theory to two suites of soils; i.
What distinguishes this particular work is in its more systematic approach to the analysis, an improvement which is appropriate for the greater richness of the data analyzed. Such an approach also allows an increase in the corresponding richness of interpretation. The purpose of the study is fundamentally two-fold. The first goal is to check whether the data conform generally to the theory as published heretofore.
The second purpose, if the first should be verified, is to use any discrepancies between theory and observation to draw inferences regarding 1 in what ways the theory may be incomplete or inadequate, 2 what modifications of typical existing experimental procedure may be necessary to test the theoretical results properly, and 3 what parameters may have been estimated improperly.
The field data were drawn from alpine sites as well as from sites with Mediterranean climate. Data collected include what is required to calculate or estimate such relevant parameters as the infiltration rate, the erosion rate, and a characteristic particle size, d In order to incorporate any temporal variability in length and time scales into predictions, it will be necessary to make some straightforward extensions of the theoretical model and, possibly, to collect additional data as well.
Chemical weathering rates in the field decline according to a power-law by orders of magnitude over time White and Brantley, Such weathering rates are also demonstrated to be proportional to fluid flow rates Maher, Understanding this particular pair of results has posed problems for workers in this field.
In particular, progress in understanding reactive solute transport has been limited Hunt and Ewing, ; Hunt and Manzoni, ; Hunt et al. Although the ADE can be applied to predicting solute transport at the scale of a single pore e.
The only way the ADE can relate the observed reduction in reaction rates to a diminishing solute transport capability in time rather diffusion like is to abandon the observed proportionality to the flow rates.
In other words, one cannot have it both ways with the continuum approach. The solute velocity is obtained from a known scaling relationship between transit time and system length Lee et al. When Da I is larger than 1, the transport limitation for the chemical weathering rate is expected to be valid; for smaller values of Da I , reaction kinetics will dominate.
With increasing time of soil evolution soil depth usually increases. At shorter time scales, when Da I is small and the kinetics of the particular weathering reaction dominant at each specific site dominates, more structure is expected to be visible in the time-dependence of the weathering rates. It is necessary at the outset to be clear that the solute transport limitations discussed here arise from advective solute transport, not from diffusion, as has been argued by Bandopadhyay et al.
Arguments based on the Peclet number as calculated from characteristic instantaneous flow velocities at the scale of a single pore were used to justify this a priori Hunt and Manzoni, However, this may be criticized on the grounds that one should use a yearly average velocity, since diffusion may be relevant throughout the year, even in the vadose zone.
The structure of advective solute transport paths in heterogeneous porous media is fractal, though not, in general, related to any structure of the medium, and it is this fractality, which leads to the power-law decay in solute transport fluxes with time Hunt and Ewing, In applying percolation theory, it is assumed that the flow paths of least resistance can be calculated using the critical percolation probability.
As long as the optimal flow paths can be described using such critical paths, their fractal dimensionality is as given in percolation theory, explaining why the particular characteristics of the medium have reduced importance to flow path characteristics.
Percolation theory generates a suite of properties relevant to flow, diffusion, and dispersion, as well as to structure, although it is sometimes necessary to choose which percolation results are appropriate, or whether an alternate formulation, such as the effective medium approximation EMA , may be more suited to generating an accurate prediction in any specific case Hunt and Sahimi, Nevertheless, the general theory of percolation is best discussed in terms of its topological implications first.
Also, although its application need not be restricted to a regular grid, it is much easier to discuss under such restrictions, so only regular grids, also known as lattices, are considered here in the theory review.
The best review for understanding the basics of percolation theory is by Stauffer and Aharony , from which much of the following discussion ultimately comes. Consider a square lattice, on which sites may be occupied by either plastic or steel spheres of the same size. The choice, metal or plastic, at any given site is randomly generated, though the conclusions do not change fundamentally if spatial correlations are added.
Nearest neighbor spheres touch at one point. If sufficiently many of the spheres are metallic, a continuously connected infinitely long path through metallic spheres is produced.
This transition occurs, for any given medium, at a specific value of the fraction of spheres that is metallic, called p c. All other common systems have smaller values of p Hunt et al. Bond percolation values are smaller than site percolation values on the same lattice. In order for a bond percolation problem to be relevant to real systems, the hydraulic or electric conductance values connecting sites i. In a procedure called critical path analysis CPA the percolation probability can then be used in a way that generates the connected path of lowest total resistance to flow , i.
This is particularly relevant here, as Sahimi has emphasized how the topological structure of the critical paths will dominate solute transport whenever they describe such flow concentration as are observed.
The small values of p c in bond percolation in natural media imply that only a relatively small portion of the medium generates nearly all the fluid flow. It thus also means that percolation theory may generate connected flow paths that look like what hydrogeologists refer to as preferential flow.
It is quite generally believed that such preferential flow paths are important to solute transport National Research Council, Here, and in previous work, we have made this assumption regarding the relevance of such paths, as medium flow paths were assumed compatible with CPA. Then the percolation descriptions of, e. Once one has an infinitely large interconnected cluster of sites or bonds , it is possible to define the percolation backbone as that part of the infinite connected cluster which remains when all sites bonds that can be connected to it through only one point have been removed Stanley, All sites that connect through two points are retained.
Sites that connect only through one point are called dead-ends, since they do not support flow. The remaining structure has mass fractal dimensionality D b , a quantity which is relevant for solute transport Lee et al. Importantly, for a wide range of conditions, the value of D b remains the same, but it differs fundamentally depending on the dimensionality of the system studied Sheppard et al.
Flow in fractures as well as in strongly layered anisotropic soils may be fundamentally 2D, but, more generally, 3D conditions are observed. The value of D b does depend sensitively on the saturation condition, i. The latter two processes are known as invasion percolation Wilkinson and Willemsen, Most other medium variations, such as the particular lattice type chosen, or the particle size values, do not change the value of D b. However, in the present work, complications from such correlations are not addressed.
Moreover, evidence for the relevance of such correlations, though demonstrated for such frequently measured quantities as the electrical conductivity Hunt and Sahimi, , has not yet been found for problems of chemical weathering or soil formation. In the following, the scaling relationships relating time, distance, and velocity of solute transport are reproduced from known results for systems near the percolation threshold Sheppard et al.
What can make their applicability universal, however, is the tendency for water flow in disordered media to follow paths of least resistance, as defined in CPA. Then the solute transport is controlled by paths near the percolation threshold, as quantified by using percolation theory.
In Figure 1 , the concept of the percolation theory is schematically drawn and compared with a soil mass balance. Soil depth depends on mass input and output. Figure 1. Schematic overview of the applied concept. The weathering front e is intimately linked to these water fluxes. B Soil depth x as a function of the median granulometry x 0 as fundamental length scale of the soil scheme modified after Birkeland, Time constraints are given by the fundamental length scale, x 0 , and the corresponding time scale, t 0 , and x and the related time scale t.
Erosion E and chemical leaching W contribute to mass losses. Besides P Soil and the parameters E and W , soil depth is strongly related to the hydrologic water balance, granulometry of the medium and time or velocity. The specific role of percolation theory in describing solute transport is now discussed. When solute enters a medium at a point Lee et al. This proportionality can be represented as an equation, if appropriate values of constants representing a fundamental length scale, x 0 , and corresponding time scale, t 0 , can be identified,.
We refer to v 0 as a pore-scale flow rate. In the context of the following discussion, it will become clear that v 0 must be a value that characterizes the mean of the local vertical flow rate. Determination of one additional parameter, either x 0 or t 0 , completes the parameter determination in Equation 2. If the network representation is applied to a porous medium, essential suggestive choice for x 0 is a pore separation Hunt and Manzoni, , which should be more or less equal to a particle diameter, because this length scale defines the separations of the local connections between flow pathways.
In a highly disordered network, where particle sizes can vary widely, we have proposed Yu and Hunt, b , c that the best choice for x 0 is d 50 , the median particle size.
The larger instantaneous flow rates, which are limited primarily by the hydraulic conductivity, can often be on the order of 0. While the instantaneous flow rates must be large enough to neglect diffusion, the time-averaged flow rates must be small enough to limit reaction. Otherwise, reaction kinetics can be the limiting factor. Note that this definition corrects inconsistencies in Salehikhoo et al. The tendency for Da I to increase with column length is evident in Equation 3.
The pore volume is proportional to length, canceling the linear factor in L in the numerator, but the more rapid than linear increase in transport time with length leaves a second factor in length to the 0. In the previously analyzed case of MgCO 3 dissolution Yu and Hunt, a , for all field conditions the value of Da I was never smaller than tens of thousands.
However, reaction rates of silicate minerals, even under well-mixed conditions, can be orders of magnitude smaller Stumm and Morgan, than for MgCO 3 , as treated by Salehikhoo et al. The time derivative of x t yields the solute velocity, argued above to be a proxy for a chemical weathering rate, v. Equation 5 can be rewritten in a form that depends only on the distance, x.
Since the introduction of erosion can, in principle, make it impossible to define a unique time for a given transport distance, it is more useful to write Equation 5 in a form that eliminates time from the equation:. In the absence of erosion, it is a reasonable hypothesis that the total solute transport distance is equal to the soil depth.
Then, the temporal derivative of the soil depth, v , is the soil production function in units of depth divided by time , as also given by Equation 5 , making the soil production function proportional to the chemical weathering rate. The proportionality constant is equal to the ratio of the bulk density of the chemically weathered material to its molecular mass.
Equation 4 for the soil depth can hold only as long as erosion can be neglected, which is very rarely the case. Since the soil production rate Equation 5 declines with age or depth, Equation 6 , the period of time when erosion can be neglected is always limited.
When erosion, E , must be considered, one can construct an equation for the soil depth based on the concept of mass balance,. Here, R is given by Equation 6 and E is, for constant erosion rates, a parameter. In general, however, E is a function of time. The term E in fact is equaled to denudation that includes besides the output of solid material also chemical leaching of silicate particles. Compared to erosion, this type of leaching and, therefore, loss is in most cases fully subordinate Dixon and von Blanckenburg, In three dimensions with moisture conditions corresponding either to wetting or full saturation, D b has the value 1.
Owing to the fractional exponent of x introduced by R , Equation 7 does not have an analytical solution, but it may be readily solved numerically. One simply solves Equation 4 for some initial sufficiently small time and then calculates R from Equation 6.
Using the calculated value of R and the field value for E as well as an arbitrary time step one can calculate the change in soil depth and add it to the initial value.
Then one calculates R from Equation 6 using the new soil depth. This procedure is then simply followed until the total time elapsed is equal to the age of the soil, or until the depth no longer changes in time, at which point a steady-state soil depth has been generated.
As long as no significant changes in parameters occur, steady-state conditions will then prevail. At the opposite short time end of the time spectrum, an important complication can arise when chemical weathering is not solute transport-limited. In this case, a constant rate of weathering would ensue, reaction kinetics provide a limitation which is unchanging in time.
Since the limits imposed by kinetics are stronger in this case, at least at short time scales, than those due to solute transport, data would lie below the percolation predictions.
Further, a constant reaction rate would imply a linear increase in soil depth with increasing time up until the point that the predicted and observed depths were equal, at which time the transport-limited result would become valid.
Theoretical sensitivity of soil thickness to various parameters may be estimated from Equation 4 for short times when erosion might be neglected , or at long times from Equation 9. Owing to the power-law forms of these equations, the sensitivity relationships relate simply to the exponents. At intermediate times, sensitivities, like depths, must be obtained numerically. On account of the gradual evolution of the overall behavior from Equation 4 to Equation 9 through time, and the variable time period over which this change occurs due to variation in actual parameter values, actual sensitivities will exhibit somewhat more complex behavior as a function of time.
Data for soil depth as a function of age have been collected for a large number of sites in two distinct geographic environments: Alpine sites; Supplementary Table S1 , and Mediterranean 94 sites; Supplementary Table S2. Other data have been accessed from the literature.
All soils have developed from unconsolidated material. The sites are distributed on five continents, as shown in Figure 2. Some of these parameters are more or less easily accessible e. Accompanying data relevant to assess actual evapotranspiration, AET, and soil particle sizes, which are necessary to make concrete predictions, have also been collected for each site within these regions.
All data, together with their sources, are given in the Supplementary Tables S1 — S6. Soil depths are determined mostly through a process of excavation and measurement. We used datasets where information was available about soil horizons designation and thickness. In addition, data was collected where available about soil density, coarse fragment content and grain size distribution.
According to Sauer et al. Further, particle size distributions at the time of original exposure of the medium are assumed to be preserved in the C layer. Soil mass was determined using the thickness of the horizons, their corresponding bulk density and summed up over the entire profile. This mass was calculated with and without coarse fragments soil mass and mass of fine earth FE. For the stocks of FE we have:. Sites were selected where numerical indication about the surface and its soils was available e.
Erosion rates are sensitive to local vegetation, relief and slope angle, aspect, climate, and topographic curvature, as well as substrate. Not all of these influences can be quantified. In situ measurements of erosion or denudation were available only for a few sites e. Otherwise, present-day erosion rates had to be estimated from published maps e. Furthermore, related information was also available from specific investigations, e. Time-averaged over the entire soil evolution rates of soil erosion were measured for a few alpine sites some sites of the Wind River Range and European Alps.
For the other sites, erosion rates had to be estimated over the entire soil evolution. Soils on terraces many Mediterranean sites with a flat position have a very low to almost negligible water erosion rate Panagos et al. According to Raab et al.
The transition from the Pleistocene to the Holocene was accompanied by distinctly higher rates. Raab et al. These fluctuations, which occurred during the evolution of many of the soils considered, are shown in Figure 3 and had to be considered for sites where erosion values were derived from the previously mentioned maps mostly alpine sites. We assumed that the maps provide an average and relatively reliable erosion rate value for undisturbed sites for the entire Holocene.
The erosion rate of soils having an older age had to be corrected using the average trend given in Figure 3. Figure 3. A Climate oscillations over the last 2. Note, similarly to B that erosion increases at the transition from cold to warm periods. Using these data it was possible to determine directly d A mechanical disintegration of the rock material into small units facilitates chemical decay by increasing the total area of particle surfaces and surface reactive sites that are in contact with the solutions Stumm and Wollast, ; Lageat et al.
The d 50 value was then multiplied by 0. Tortuosity models for porous media Ghanbarian-Alavijeh et al. We therefore used a length factor of 0.
Data for many sites, however were not sufficient for the determination of d Previously published data, for example, often give only the percentages of the three fundamental size classes: sand, silt, and clay. Since it is necessary to know d 50 in order to make a concrete prediction of soil depths, we developed a regression routine over well-characterized soils for this purpose, using the percentages of sand, silt, and clay to calculate a mean diameter for the input, and the observed d 50 as an output.
Since we expected, at the very least, key differences in regression parameters between Alpine and Mediterranean sites, these regressions were performed separately. The results are given in Figures 4A,B. Once these relationships are established, we can use the appropriate regression relationship to generate automatically a reasonable median particle size for any soil in these two environments, as long as the sand, silt, and clay fractions are available.
Note that the particle size distribution is a function of height in the soil column. This result implies that the soil texture is changing over time. Specific results indicate that the studied soils either become finer over time, or do not change perceptibly Figure 5. These results are rather comparable in both Mediterranean and Alpine regions. From the theoretical sensitivity analysis, one should expect that those soils most severely impacted by diminishing particle sizes could be a factor 10 0.
Even use of a median grain size in such soils, rather than the final value, could still overpredict the depth by a factor roughly as high as the square root of 3. Figure 4. Transfer functions to calculate d 50 for sites having only indications about the three fundamental grain size classes sand, silt, clay based on detailed grain size data for A the alpine Supplementary Table S3 and B Mediterranean type of sites Supplementary Table S4.
Figure 5. Evolution of the median grain size d 50 of the soils alpine or Mediterranean environment over time. This detectable decrease is significant for both series of soils.
Infiltration is the parameter that is most difficult to obtain accurately, while it is also demonstrated below to be the parameter, to which the predicted soil depth responds most sensitively.
Strictly speaking, the infiltration rate required here, the fraction of precipitation actually relevant for soil weathering, is what penetrates to the bottom of the soil layer, and is the difference between precipitation and actual evapotranspiration plus whatever surface water runs on to the site less the amount that runs off.
Infiltration is seldom measured. It is furthermore difficult to estimate the local variability in this variable. These datasets, however, provide only an overview.
The infiltration rate, however, might have varied considerably over the period of soil evolution. Consequently, an estimate of the hydrologic mass balance had to be estimated for the duration of soil development. For this purpose, information about palaeoclimate had to be accessed. Basic data about climate variability were obtained from the following sources:.
Asia: Kigoshi et al. Southern Europe: Peyron , Blain et al. Andes: Graf , Schauwecker et al. In southern Europe and mid to southern USA, the climate was distinctly colder and precipitation rates varied from slightly lower to higher, depending on the area. The present-day hydrologic mass balance was used for soils that started to form during the Holocene.
For older soils, climatic conditions also of the Pleistocene had to be considered. Because the percolation theory approach uses the hydrologic mass balance, precipitation, overland flow, and evapotranspiration had to be estimated also for periods prior to the Holocene. Even though the climatic conditions were cold to very cold at the alpine sites, it does not mean that no weathering has occurred during this period. Zollinger et al. The hydrologic mass balance is relevant for the percolation theory approach—temperature is therefore at first sight less important.
The present-day and averaged precipitation, evapotranspiration, and infiltration rates over the entire soil formation period are given in the Supplementary Tables S1 , S2. Freshly created volcanic soils, called andisols , and clay-rich soils that hold nutrients and water are examples of productive soils. Water erosion is accentuated on sloped surfaces because fast-flowing water has higher eroding power than still water.
Raindrops can disaggregate exposed soil particles, putting the finer material e. Sheetwash , unchannelled flow across a surface carries suspended material away, and channels erode right through the soil layer, removing both fine and coarse material.
Wind erosion is exacerbated by the removal of trees that act as windbreaks and by agricultural practices that leave bare soil exposed. Tillage is also a factor in soil erosion, especially on slopes, because each time a cultivator lifts the soil, it is moved a few centimeters down the slope. Earth has two important carbon cycles. One is the biological one, wherein living organisms, mostly plants, consume carbon dioxide from the atmosphere to make their tissues, and then, after they die, that carbon is released back into the atmosphere when they decay over years or decades.
A small proportion of this biological-cycle carbon becomes buried in sedimentary rocks: during the slow formation of coal, as tiny fragments and molecules in organic-rich shale, and as the shells and other parts of marine organisms in limestone. The geological carbon cycle below shows the various steps in the process not necessarily in this order :. Under these conditions, the climate remains relatively stable.
This can happen during prolonged periods of higher than average volcanism. One example is the eruption of the Siberian Traps at around Ma, which appears to have led to strong climate warming over a few million years. A carbon imbalance is also associated with significant mountain-building events. For example, the Himalayan Range was formed between about 40 and 10 million years ago.
Over that period, and still today, the rate of weathering on Earth has been enhanced because those mountains are so high, and the range is so extensive. The weathering of these rocks, most importantly the hydrolysis of feldspar, has resulted in the consumption of atmospheric carbon dioxide and transfer of the carbon to the oceans and ocean-floor carbonate minerals.
The steady drop in carbon dioxide levels over the past 40 million years, which led to the Pleistocene glaciations, is partly attributable to the Himalayan Range formation. Another, non-geological form of carbon-cycle imbalance is happening today on a very rapid time scale.
Eolian soils are not very productive because they have very low water-holding capacity , are low in organic matter , and are nutrient deficient as compared to loess soils. Most are used for grass production or natural habitat. Geologic materials moved from the parent material by water are known as alluvium. Alluvial deposits are found in flood plain areas such as the Platte River and other stream valleys.
Since stream beds constantly change over time, alluvial parent materials are highly variable as are the soils that form them. Physical processes primarily result in the breakdown of rocks into smaller and smaller particles. As the particles become smaller, various living organisms begin to have a great impact on soil formation because they contribute organic matter. In addition, the smaller particles speed chemical processes which result in new chemical compounds. All of these processes are greatly influenced by climate , especially temperature and precipitation.
Precipitation, in particular, ranges from an average of 33 inches per year in southeastern Nebraska to 15 inches per year in western Nebraska Fig. The amount of water entering a soil influences the movement of calcium and other chemical compounds in the soil. Ultimately, if more chemicals are removed, the soils will be deeper and more developed. Precipitation influences vegetation and, therefore, greatly determines the organic matter content of soils.
Because of greater precipitation in eastern Nebraska, native vegetation included luxuriant growth of the tallgrass prairie. In western Nebraska where precipitation is about half that in the east, plants of the shortgrass prairies grow much less abundantly. Thus, soil organic matter content is greater in the east than in the west. Higher temperatures can speed the rate of organic matter decomposition.
Temperatures are typically higher in the southern portion of the state than in the northern portion Fig. Because of this trend, organic matter content decreases from north to south.
However, the change in organic matter content from north to south due to temperature is minuscule when compared to the change from east to west due to precipitation. Soils in eastern Nebraska commonly contain 3 percent organic matter as compared to about 1 to 2 percent in the west. The most abundant living organism in the soil is vegetation. Vegetation influences the kind of soil developed because plants differ in their root systems, size, above ground vegetative volume, nutrient content and life cycle.
Soils formed under trees are greatly different from soils formed under grass even though other soil-forming factors are similar. Trees and grass vary considerably in their search for food and water and in the amount of various chemicals taken up by roots and deposited in or on top of the soil when tree leaves and grass blades die. Soils formed under grass are much higher in organic matter than soils formed under forests because of their massive fibrous root structure and annual senescence of above ground vegetation.
Grassland soils tend to be darker, particularly to greater depths, and have a more stable structure than forest soils. Soils developed under grass are generally more fertile and best suited for crop production. Nebraska soils from any parent material are nearly all formed under grass and, with adequate water, can be very productive.
0コメント