Integration of geospatial technologies with RUSLE for analysis of land use/cover change impact on soil erosion: case study in Rib watershed,north-western highland Ethiopia

In recent times, soil erosion interlocked with land use and land cover (LULC) changes has become one of the most important environmental issues in developing countries. Evaluation of this complex interaction between LULC change and soil erosion is indispensable in land use planning and conservation works. This paper analysed the impact of LULC change on soil erosion in the north-western highland Ethiopia over the period 1986–2016. Rib watershed, the area with dynamic LULC change and severe soil erosion problem, was selected as a case study site. Integrated approach that combined geospatial technologies with revised universal soil loss equation model was utilized to evaluate the spatio-temporal dynamics of soil loss over the study period. Pixel-based overlay of soil erosion intensity maps with LULC maps was carried out to understand the change in soil loss due to LULC change. Results showed that the annual soil loss in the study area varied from 0 to 236.5 t ha^?1 year^?1 (tons per hectare per year) in 1986 and 0–807 t ha^?1 year^?1 in 2016. The average annual soil loss for the entire watershed was estimated about 40 t ha^?1 year^?1 in 1986 comparing with 68 t ha^?1 year^?1 in 2016, a formidable increase. Soil erosion potential that was estimated to exceed the average soil loss tolerance level increased from 34.5% in 1986 to 66.8% in 2016. Expansion of agricultural land at the expense of grassland and shrubland was the most detrimental factor for severe soil erosion in the watershed. The most noticeable change in soil erosion intensity was observed from cropland with mean annual soil loss amount increased to 41.38 t ha^?1 year^?1 in 2016 from 26.60 in 1986. Moreover, the most successive erosion problems were detected in eastern, south-eastern and northern parts of the watershed. Therefore, the results of this study can help identify the soil erosion hot spots and conservation priority areas at local and regional levels.     

Geospatial assessment of soil erosion intensity and sediment yield: a case study of Potohar Region,Pakistan

Estimation of spatial extent of soil erosion, one of the most serious forms of land degradation, is critical because soil erosion has serious implications on soil fertility, water ecosystem, crop productivity and landscape beauty. The primary objective of the current study was to assess and map the soil erosion intensity and sedimentation yield of Potohar region of Pakistan. Potohar is the rainfed region with truncated and complex topography lying at the top of the Indus Basin, the world’s largest irrigation networks of canals and barrages. Spatially explicit Revised Universal Soil Loss Equation (RUSLE) Model integrated with Remote Sensing-GIS techniques was used for detecting/mapping of erosion prone areas and quantification of soil losses. The results show that the Potohar region is highly susceptible to soil erosion with an average annual soil loss of 19 tons ha^?1 year^?1 of which the maximum erosion (70–208 tons ha^?1 year^?1) was near the river channels and hilly areas. The sediment yield due to the erosion is as high as 148 tons ha^?1 year^?1 with an average of 4.3 tons ha^?1 year^?1. It was found that 2.06% of the total area falls under severe soil erosion, 13.34% under high erosion, 15.35% under moderate soil erosion while 69.25% of the area lies in the low (tolerable) soil erosion. Chakwal and Jhelum districts of the region are seriously affected by erosion owing to their topography and soil properties. The information generated in this study is a step forward towards proper planning and implementation of strategies to control the erosion and for protection of natural resources. It is, hence, necessary that suitable water harvesting structures be made to control water to prevent soil erosion and provision of water in the lean season in this region. Tree plantation and other erosion control practices such as strip cropping can also minimize soil erosion in this region.     

Multi-temporal annual soil loss risk mapping employing Revised Universal Soil Loss Equation (RUSLE) model in Nun Nadi Watershed, Uttrakhand (India)

The present comparative study is multi-temporal in nature. The Revised Universal Soil Loss Equation (RUSLE), remote sensing, and GIS were used to model the soil loss estimation for soil conservation and vegetation rehabilitation in Nun Nadi watershed for the years 2000 and 2009. The estimated mean soil loss for the year 2000 and 2009 is 3,283.11 and 1,419.39 Mg?ha^?1 year^?1, respectively. The study finds that about 80 % area has low or least risk of erosion and about 7 % is exposed to high or very high risk which indicates the improvement in terms of soil loss if we compare the data of both the time periods. The findings show that the rainfall, LULC change, and elevation are the main responsible factors for the soil loss in Nun Nadi watershed. Conservation measures have been adopted; however, the problem still remains serious and demands urgent attention.     

Unprecedented rates of landslide and surface erosion along a newly constructed road in Yunnan, China

Field measurements conducted 4 years after the construction of a new portion of the Weixi?CShangri-La road in Yunnan, China, reveal that unprecedented rates of mass wasting occurred along the road with much of this sediment directly impacting the headwaters of the Mekong River. Landslide erosion (including dry ravel) exceeded 33,000 t ha^?1 year^?1 along the most severely eroded sections of the road and averaged more than 9,600 t ha^?1 year^?1 along the surveyed 23.5 km of road; these values are the highest ever reported for road-related landslides. While surface erosion was only about 7% of the total erosion from the road, it is still more than an order of magnitude higher than typical surface erosion rates from disturbed lands in Southeast Asia. Combined landslide and surface erosion from this road delivered an estimated 19 times more sediment to the river than the remaining 99.6% of the contributing catchment. These sediment inputs are aggrading local channels, promoting downstream sediment transport, degrading aquatic habitat, and creating the possibility for a future debris flood or hyperconcentrated flow.     

Integration of GIS and remote sensing for estimation of soil loss and prioritization of critical sub-catchments: a case study of Tapacurá catchment

Mapping of erosion risk areas is an important tool for the planning of natural resources management, allowing researchers to propose the modification of land use properly and implement more sustainable long-term management strategies. The objective of this study was to assess and identify critical sub-catchments for soil conservation management using the USLE, GIS, and remote sensing techniques. The Tapacurá catchment is one of the planning units for water resource management of the Recife Metropolitan Region. Maps of the erosivity (R), erodibility (K), slope (LS), cover-management (C), and support practice (P) factors were derived from the climate database, digital elevation model, and soil and land-use maps. In order to validate the simulation process, total sediment delivery ratio was estimated. The results showed a mean sediment delivery ratio (SDR) of around 11.5?% and a calculated mean sediment yield of 0.108?t?ha^?1?year^?1, which is close to the observed one, 0.169?t?ha^?1?year^?1. The obtained soil loss map could be considered as a useful tool for environmental monitoring and water resources management. The methodology applied showed acceptable precision and allowed the identification of the most susceptible areas to soil erosion by water, constituting an important predictive tool for soil and environmental management in this region, which is highly relevant for the prediction of varying development scenarios for Tapacurá catchment. This approach can be applied to other areas for simple and reliable identification of critical areas of soil erosion in catchments.     

Erosion rates evaluated by the <Superscript>137</Superscript>Cs technique in the high altitude area of the Qinghai–Tibet plateau of China

Studying spatial and temporal variation of soil loss is of great importance because of global environmental concerns. Understanding the spatial distribution of soil erosion and deposition in the high-cold steppe is important for designing soil and water conservation measures. Measured ¹³⁷Cs losses (Bq m⁻²) from long-term high altitude (4,000 m above sea level) watershed plots on the Qinghai–Tibet plateau and derived soil erosion estimates (Mg ha⁻¹ year⁻¹) were significantly correlated to directly measured soil losses from the same plots, over the same period (1963–2005). The local reference inventory was estimated to be 2,468 Bq m⁻². The result of analyzing ¹³⁷Cs distribution and its intensity in the soil profiles in this area shows similarities to ¹³⁷Cs distribution in other areas. ¹³⁷Cs is basically distributed in the topsoil layer of 0–0.3 m. Soil erosions vary greatly in the entire sampled area, ranging from 5.5 to 23 Mg ha⁻¹ year⁻¹, with an average of 16.5 Mg ha⁻¹ year⁻¹ which is a moderate rate of erosion.     

Organic carbon stocks and erosion in the soils of Guangdong,South China

Soil organic carbon (SOC) storage and erosion in South China at the regional scale in the past decades remains far from being understood. This paper calculated the SOC density, storage and erosion in 14 soil classes in Guangdong Province, South China, based on statistical data from the soil survey and soil erosion survey of Guangdong, which was performed in the 1990s. The purpose of this study is to understand the relationships between soil classes and SOC erosion at the regional scale. The results indicated that the SOC density in the soils of Guangdong varied from 12.7 to 144.9 Mg ha^?1 over the entire profile and from 12.6 to 68.4 Mg ha^?1 in the top 20-cm soil layer. The average area-weighted SOC density in the topsoil (0–20 cm) and the entire profile was 32 ± 3 and 86 ± 4 Mg ha^?1, respectively. The total SOC storage was 1.27 ± 0.06 Pg, with 35.6 % (0.46 ± 0.04 Pg) located in the topsoil. The average area-weighted strength of the SOC erosion in the 1990s was 20.6 ± 0.8 Mg km^?2 year^?1. The results indicated that SOC erosion was strongly related to soil class.     

Estimation of soil erosion using RUSLE in Caijiamiao watershed,China

Soil erosion is a serious environmental and production problem in China. In particular, natural conditions and human impact have made the Chinese Loess Plateau particularly prone to intense soil erosion area. To decrease the risk on environmental impacts, there is an increasing demand for sound, and readily applicable techniques for soil conservation planning in this area. This work aims at the assessment of soil erosion and its spatial distribution in hilly Loess Plateau watershed (northwestern China) with a surface area of approximately 416.31 km². This study was conducted at the Caijiamiao watershed to determine the erosion hazard in the area and target locations for appropriate initiation of conservation measures using the revised universal soil loss equation (RUSLE). The erosion factors of RUSLE were collected and processed through a geographic information system (GIS)-based approach. The soil erosion parameters were evaluated in different ways: The R-factor map was developed from the rainfall data, the K-factor map was obtained from the soil map, the C-factor map was generated based on Landsat-5 Thematic Mapper image and spectral mixture analysis, and a digital elevation model with a spatial resolution of 25 m was derived from topographic map at the scale of 1:50,000 to develop the LS-factor map. Support practice P factor was from terraces that exist on slopes where crops are grown. By integrating the six-factor maps in GIS through pixel-based computing, the spatial distribution of soil loss in the study area was obtained by the RUSLE model. The results showed that spatial average soil erosion at the watershed was 78.78 ton ha^?1 year^?1 in 2002 and 70.58 ton ha^?1 year^?1 in 2010, while the estimated sediment yield was found to be 327.96 × 10⁴ and 293.85 × 10⁴ ton, respectively. Soil erosion is serious, respectively, from 15 to 35 of slope degree, elevation area from 1,126 to 1,395 m, in the particular area of soil and water loss prevention. As far as land use is concerned, soil losses are highest in barren land and those in waste grassland areas are second. The results of the study provide useful information for decision maker and planners to take appropriate land management measures in the area. It thus indicates the RUSLE–GIS model is a useful tool for evaluating and mapping soil erosion quantitatively and spatially at a river watershed scale on a cell basis in Chinese Loess Plateau and for planning of conservation practices.     

Using beryllium-7 to evaluate soil erosion processes in agricultural lands in semiarid regions of Central Argentina

This work presents the results of a soil erosion study using the ⁷Be technique. This technique estimates the water erosion/deposition from the comparison between ⁷Be soil content of a reference site and an eroded or sedimented site. The soil samples were collected from an agricultural area of the semiarid region of Argentina near San Luis City, which has a marked rainfall season. The area has been used for crop cultivation, being subjected to plowing practices. The ⁷Be in the Reference Site was in the first centimeter of soil, showing the typical exponential decreasing of ⁷Be soil content with depth, with the ⁷Be inventories value being 340?±?50 Bq m^?2 for the dry season and 571?±?48 Bq m^?2 for the rainy season. The ⁷Be technique was applied to a potential eroded site subjected to traditional tillage practices (plowing). A net soil erosion value of 13.5 t ha^?1 (1.35 kg m^?2) was obtained. From the assumptions of the applied technique, we can draw the inference that this erosion was caused by rains produced in the month prior to the date of soil sampling.     

Effectiveness of Living Shorelines as an Erosion Control Method in North Carolina

Living shorelines are a shoreline stabilization strategy encompassing a range of vegetative to structural materials and serve as an alternative approach to the use of structures like bulkheads, which are known to aggravate erosion. Living shorelines are often installed with little to no long-term monitoring for effectiveness; specifically, there is a lack of quantitative data regarding their performance as a shoreline stabilization strategy. This study sought to assess the performance of living shorelines with sills, with respect to shoreline protection, by determining shoreline change rates (SCR) using geospatial analysis. Shoreline surveys were conducted using a real-time kinematic (RTK)-GPS unit at a total of 17 living shoreline projects and nine control segments at 12 sites along the coast of North Carolina. Current shoreline position was compared to historic (pre-installation) shoreline positions obtained from aerial imagery, dating to 1993. The average SCR among northern sites before installation was ??0.45?±?0.49 m year^?1, and in southern sites, it was ??0.21?±?0.52 m year^?1. After installation, average SCR was significantly less erosive at northern and southern sites with living shorelines, 0.17?±?0.47 and ??0.01?±?0.51 m year^?1, respectively. Of the 17 living shoreline project segments, 12 exhibited a reduction in the rate of erosion; of those 12, six were observed to be accreting. This study supports the convention that living shorelines can reduce the rate of erosion and potentially restore lost shore zone habitat.     

Carbon Sequestration and Sediment Accretion in San Francisco Bay Tidal Wetlands

Tidal wetlands play an important role with respect to climate change because of both their sensitivity to sea-level rise and their ability to sequester carbon dioxide from the atmosphere. Policy-based interest in carbon sequestration has increased recently, and wetland restoration projects have potential for carbon credits through soil carbon sequestration. We measured sediment accretion, mineral and organic matter accumulation, and carbon sequestration rates using ¹³⁷Cs and ²¹⁰Pb downcore distributions at six natural tidal wetlands in the San Francisco Bay Estuary. The accretion rates were, in general, 0.2?C0.5?cm?year^?1, indicating that local wetlands are keeping pace with recent rates of sea-level rise. Mineral accumulation rates were higher in salt marshes and at low-marsh stations within individual sites. The average carbon sequestration rate based on ²¹⁰Pb dating was 79?g?C?m^?2?year^?1, with slightly higher rates based on ¹³⁷Cs dating. There was little difference in the sequestration rates among sites or across stations within sites, indicating that a single carbon sequestration rate could be used for crediting tidal wetland restoration projects within the Estuary.     

Integrated mass balance of some heavy metals fluxes in Yaakob village,south Sohag,Egypt

The pollution of soil with heavy metals has direct or indirect adverse effect on human health. The present work was conducted to identify all the expected sources and sinks for heavy metals by applying mass balance model to identify the retention rate of metals by soils in Yaakob village, south Sohag Governorate, Egypt. The studied inputs (sources) include P-fertilizers, irrigation water and dustfall, while the main outputs (sinks) are drainage water and harvested plants. The measurements indicate that soil, clover, dustfall and P-fertilizers contain considerable concentration of Cd, Cr, Co, Cu and Pb. The mass balance measurements indicate that the accumulation rate of Cd, Cr and Co in soil was 5.4, 54.6 and 16.3 g ha^?1 year^?1, respectively. However, depletion trend of Pb and Cu was about 1.4 and 5.2 g ha^?1 year^?1, respectively. The main source of Cd, Pb, Cr and Co in the study area is P-fertilizers with input flux 14.9, 89.9, 198.6 and 18.5 g ha^?1 year^?1, while Cu source was dustfall with 19.33 g ha^?1 year^?1. The index of geoaccumulation calculations indicates different degrees of contamination with Cd, Cr, Co and Cu. On the other hand, the main sink for the studied heavy metals was the Egyptian clover (Trifolium alexandrinum) which can be considered a good bioaccumulator of heavy metals.     

Assessment of soil erosion risk using SWAT model

Soil erosion is one of the most serious land degradation problems and the primary environmental issue in Mediterranean regions. Estimation of soil erosion loss in these regions is often difficult due to the complex interplay of many factors such as climate, land uses, topography, and human activities. The purpose of this study is to apply the Soil and Water Assessment Tool (SWAT) model to predict surface runoff generation patterns and soil erosion hazard and to prioritize most degraded sub-catchment in order to adopt the appropriate management intervention. The study area is the Sarrath river catchment (1,491 km²), north of Tunisia. Based on the estimated soil loss rates, the catchment was divided into four priority categories for conservation intervention. Results showed that a larger part of the watershed (90 %) fell under low and moderate soil erosion risk and only 10 % of the watershed was vulnerable to soil erosion with an estimated sediment loss exceeding 10 t?ha^?1?year^?1. Results indicated that spatial differences in erosion rates within the Sarrath catchment are mainly caused by differences in land cover type and gradient slope. Application of the SWAT model demonstrated that the model provides a useful tool to predict surface runoff and soil erosion hazard and can successfully be used for prioritization of vulnerable areas over semi-arid catchments.     

Relating soil erosion and sediment yield to geomorphic features and erosion processes at the catchment scale in the Spanish Pre-Pyrenees

Erosion and sediment redistribution are important processes in landscape changes in the short and long term. In this study, the RMMF model of soil erosion and the SEDD model of sediment delivery were used to estimate annual soil loss and sediment yield in an ungauged catchment of the Spanish Pre-Pyrenees and results were interpreted in the context of the geomorphic features. The Estaña Catchment is divided into 15 endorheic sub-catchments and there are 17 dolines. Gullies and slopes were the main erosive geomorphic elements, whereas the colluvial, alluvial, valley floor, and doline deposits were depositional elements. Spatially distributed maps of gross soil erosion, sediment delivery ratio (SDR), and sediment yield (SY) were generated in a GIS. Severe erosion rates (>100 Mg ha^?1 year^?1) were found in gullies, whereas mean and maximum erosion rates were very high on slopes developed on Keüper Facies and high in soils on Muschelkalk Facies. Where crops are grown, the depositional-type geoforms were predicted by the models to have an erosive dynamic. Those results were consistent with the rates of erosion quantified by ¹³⁷Cs which reflects the significant role of human activities in triggering soil erosion. Catchment area was positively correlated with erosion rate, but negatively correlated with SDR and SY. The latter were negatively correlated with the proportion of the surface catchment covered with forests and scrublands. The topography of the area influenced the high SDR and SY in the dolines and valley floors near the sinks. Intra-basin stored sediment was 59.2% of the total annual eroded soil in the catchment. The combination of the RMMF and SEDD models was an appropriate means of assessing the effects of land uses on soil erosion and obtaining a better understanding of the processes that underlie the geomorphic changes occurring in mountainous environments of the Mediterranean region.     

Soil erosion determinations by using 137Cs technique in the agricultural regions of Gediz Basin, Western Turkey

Gediz Basin is one of the regions where intense agricultural activities take place in Western Turkey. Erosion and soil degradation have long been causing serious problems to cultivated fields in the basin. This work describes the application of two different ¹³⁷Cs models for estimating soil erosion rates in cultivated sites of the region. Soil samples were collected from five distinct cultivated regions subject to soil erosion. The variations of ¹³⁷Cs concentrations with depth in soil profiles were investigated. Soil loss rates were calculated from ¹³⁷Cs inventories of the samples using both proportional model (PM) and simplified mass balance model (SMBM). When PM was used, erosion and deposition rates varied from −15 to −28 t ha⁻¹ year⁻¹ and from +5 to +41 t ha⁻¹ year⁻¹, respectively; they varied from −16 to −33 t ha⁻¹ year⁻¹ and from +5 to +55 t ha⁻¹ year⁻¹ with SMBM. A good agreement was observed between the results of two models up to 30 t ha⁻¹ year⁻¹ soil loss and gain in the study area. Ulukent, a small representative agricultural field, was selected to compare the present data of ¹³⁷Cs techniques with the results obtained by universal soil loss equation (USLE) applied in the area before.     

Short-Term and Interannual Variability in Primary Production in the Low-Salinity Zone of the San Francisco Estuary

We measured primary production during spring?Csummer 2006?C2007 to determine the carbon supply to the low-salinity pelagic food web of the San Francisco Estuary (SFE). Weekly or biweekly samples were taken at three stations of fixed salinity for size-fractionated primary production and biomass, both as chlorophyll and from biovolume based on counts. Error variance in productivity estimates arose mainly from the depth integration of ¹⁴C uptake, showing the importance of productivity measurements at high light levels for estimates of depth-integrated production. Temporal and spatial variability in production were surprisingly small. Combining data from this study with long-term monitoring data, productivity and biomass were variable in time and salinity but without persistent patterns and with infrequent blooms. Production within the low-salinity zone was unresponsive to variation in freshwater flow, in contrast to findings in other estuaries where nutrient loading drives variability in production and other regions of the SFE where production responds to residence time or to stratification. Estimated annual primary production was only 25 and 31?g?C?m^?2?year^?1 during 2006 and 2007, only half of it in cells >5???m. These results imply that phytoplankton provided poor food web support for higher trophic levels, probably contributing to the long-term decline in fish abundance in the brackish to freshwater region of the estuary.     

Rainfall erosivity and erosivity density in Eastern Ghats Highland of east India

The rainfall erosivity (R-factor in USLE) is the long-term average of the sum of the product of rainfall kinetic energy and its maximum 30-min intensity. Therefore, at most 30-min time intervals pluviograph records are required to calculate R-factor. But, such high-resolution data are scarce in many parts of the world and require lengthy processing period. In this study, R-factor was correlated with daily, monthly and annual rainfall, and its spatial variability in Eastern Ghats Highland of east India was mapped. The result showed that power regression models predicted satisfactorily the daily, monthly and annual R-factor, of which annual R-factor model performed best (model efficiency 0.93). Mean monsoon season R-factor was 15.6 and 10.0 times higher than the pre- and post-monsoon season R-factor, and thus remained highly critical with respect to erosion. Annual R-factor values ranged from 3040 to 10,127 MJ mm ha^?1 h^?1 year^?1, with standard deviation of 1981 MJ mm ha^?1 h^?1 year^?1. Rainfall intensity was positively correlated with erosivity density, and numerical value of rainfall intensity was almost double of the erosivity density value. The combination of rainfall and erosivity density was used to identify flood, erosion and landslide-prone areas. The developed iso-erosivity, erosivity density and risk maps can be opted as a tool for policy makers to take suitable measures against natural hazards in Eastern Ghats Highland of east India and elsewhere with similar rainfall characteristics.

     

GIS-based lake sediment budget estimation taking into consideration land use change in an urbanizing catchment area

The objective of this study was to assess the lake sediment budget of land use changes using the Universal Soil Loss Equation (USLE), sediment delivery ratio (SDR), and trap efficiency (TE). The geographic information system was combined with the USLE to estimate the soil erosion of the Lake Asan watershed. Spatial data for each of the USLE factors were obtained from the land use, soil, and 1/25,000 scale digital contour maps. Landsat-5 TM images were selected for analyzing soil erosion changes due to land use changes. The sediment yield to Lake Asan was estimated using the SDR and TE. The estimated sediment budget was compared with observed data from the Lake Asan watershed between 1974 and 2003. The total estimated annual mean sediment budgets from Lake Asan in 1986, 1992, and 2000 were 0.267, 0.301, and 0.339 × 10⁶ ton, respectively, with an average of 0.302 × 10⁶ ton. The average measured sediment budget was 3.15 × 10⁶ ton year^?1. The average estimated value shows reasonable agreement with the observed sediment balance. The average estimated and measured sediment budgets contain uncertainties due to both the methods and the approach used by the observers. The simulated results indicated that soil erosion in the Lake Asan watershed increased at a rate of approximately 2 % per year from 1986 to 2000 due to land use change. This study may be useful for managers to identify reservoir rehabilitation management methods for stable irrigation water supply.     

Rates of Biotite Weathering,and Clay Mineral Transformation and Neoformation,Determined from Watershed Geochemical Mass-Balance Methods for the Coweeta Hydrologic Laboratory,Southern Blue Ridge Mountains,North Carolina,USA

Biotite is a common constituent of silicate bedrock. Its weathering releases plant nutrients and consumes atmospheric CO₂. Because of its stoichiometric relationship with its transformational weathering product and sensitivity to botanical activity, calculating biotite weathering rates using watershed mass-balance methods has proven challenging. At Coweeta Hydrologic Laboratory the coupling of biotite to its transformational weathering product is only valid if the stoichiometric relationship for the two phases is known; this relationship is unlikely layer-for-layer. Rates of biotite weathering and transformation of its secondary weathering product at the Coweeta Hydrological Laboratory are comparable with other Appalachian watersheds. The magnitude and sign of the difference between field- and laboratory-determined biotite weathering rates are similar to those of other silicate minerals. The influence of major-cation proportions in biomass on the rates of biotite weathering and transformational weathering product is greatest for watersheds with high biomass aggradation rates. The watershed with the lowest bedrock reactivity and highest flushing rate yielded the highest gibbsite formation rate of ~500 mol ha^?1 year^?1 and lowest kaolin-group mineral formation rates of 4–78 mol ha^?1 year^?1. The kaolin-group mineral formation rate increases as bedrock reactivity increases and flushing rate decreases to a maximum of ~300 mol ha^?1 year^?1, with a similar minimum gibbsite formation rate. The relative differences in bedrock reactivity and flux of water through Coweeta Hydrological Laboratory watersheds studied appear to be invariant over geologic timescales.     

Fluid production rate during the regional metamorphism of a pelitic schist

Phase equilibria modeling of the pressure–temperature (P–T) path of regional metamorphism and associated fluid expulsion, combined with constraints on the timescale of garnet growth by Sm–Nd geochronology, elucidates the fluid production rate and fluid flux during Barrovian metamorphism of pelitic rocks from Townshend Dam, VT, USA. This modeling builds on a published companion study that utilized Sm–Nd geochronology of concentric growth zones in multiple garnet grains, to constrain the duration of garnet growth in a large sample of schist at Townshend Dam to 3.8?±?2.2 million years (Gatewood et al., Chem Geol 401:151–168, 2015). P–T pseudosections combined with observed mineral compositions constrain garnet growth conditions, and were utilized to construct P–T path-dependent thermodynamic forward models. These models determine that garnet growth was initiated at ~?0.6 GPa and ~?525 °C, with a roughly linear loading and heating P–T trajectory to >?0.8 GPa and ~?610 °C. Loading and heating rates of 2.4 km·Myear^?1 (with a range of 1.6 to 5.8 km·million year^?1) and 23 °C·million year^?1 (with a range of 14 to 54 °C·million year^?1), respectively, are consistent with model estimates and chronologic constraints for tectono-metamorphic rates during orogenesis. Phase equilibria modeling also constrains the amount of water release during garnet growth to be ~?0.7 wt% (or >?2 vol%), largely resulting from the complete consumption of chlorite. Coupling this estimate with calculated garnet growth durations provides a fluid production rate of 5.2 kg·m^?3·million year^?1 (with a range of 3.2 to 12.2 kg·m^?3·million year^?1) and when integrated over the overlying crustal column, a regional-scale fluid flux of 0.07–0.37 kg·m^?2·million year^?1. This range of values is consistent with those derived by numerical models and theory for regional-scale, pervasive fluid flow. This study signifies the first derivation of a fluid production rate and fluid flux in regional metamorphism using a direct chronology of water-producing (garnet-forming) reactions and can provide a framework for future studies on elucidating the nature and timescales of fluid release.