Modelling hydrological processes influenced by soil,rock and vegetation in a small karst basin of southwest China

Hydrological processes in karst basins are controlled by permeable multimedia, consisting of soil pores, epikarst fractures, and underground conduits. Distributed modelling of hydrological dynamics in such heterogeneous hydrogeological conditions is a challenging task. Basing on the multilayer structure of the distributed hydrology‐soil‐vegetation model (DHSVM), a distributed hydrological model for a karst basin was developed by integrating mathematical routings of porous Darcy flow, fissure flow and underground channel flow. Specifically, infiltration and saturated flow movement within epikarst fractures are expressed by the ‘cubic law’ equation which is associated with fractural width, direction, and spacing. A small karst basin located in Guizhou province of southwest China was selected for this hydrological simulation. The model parameters were determined on the basis of field measurement and calibrated against the observed soil moisture contents, vegetation interception, surface runoff, and underground flow discharges from the basin outlet. The results show that due to high permeability of the epikarst zone, a significant amount of surface runoff is only generated after heavy rainfall events during the wet season. Rock exposure and the epikarst zone significantly increase flood discharge and decrease evapotranspiration (ET) loss; the peak flood discharge is directly proportional to the size of the aperture. Distribution of soil moisture content (SMC) primarily depends on topographic variations just after a heavy rainfall, while SMC and actual ET are dominated by land cover after a period of consecutive non‐rainfall days. The new model was able to capture the sharp increase and decrease of the underground streamflow hydrograph, and as such can be used to investigate hydrological effects in such rock features and land covers. Copyright © 2011 John Wiley & Sons, Ltd.     

Are open fractures necessarily aligned with maximum horizontal stress?

Fluid flow in fractured rock is an increasingly central issue in recovering water and hydrocarbon supplies and geothermal energy, in predicting flow of pollutants underground, in engineering structures, and in understanding large-scale crustal behaviour. Conventional wisdom assumes that fluids prefer to flow along fractures oriented parallel or nearly parallel to modern-day maximum horizontal compressive stress, or S_Hmax. The reasoning is that these fractures have the lowest normal stresses across them and therefore provide the least resistance to flow. For example, this view governs how geophysicists design and interpret seismic experiments to probe fracture fluid pathways in the deep subsurface. Contrary to these widely held views, here we use core, stress measurement, and fluid flow data to show that S_Hmax does not necessarily coincide with the direction of open natural fractures in the subsurface (>3 km depth). Consequently, in situ stress direction cannot be considered to predict or control the direction of maximum permeability in rock. Where effective stress is compressive and fractures are expected to be closed, chemical alteration dictates location of open conduits, either preserving or destroying fracture flow pathways no matter their orientation.     

Simulation of steady-state flow in three-dimensional fracture networks using the boundary-element method

An efficient method for simulating steady-state flow in three-dimensional fracture networks is formulated with the use of the boundary-element method. The host rock is considered to be impervious, and the fractures can be of any orientation and areal extent. The fractures are treated as surfaces where fluid movement is essentially two-dimensional. Fracture intersections are regarded as one-dimensional fluid conduits. Hence, the three-dimensional geometric characteristics of the fracture geometry is retained in solutions of coupled sets of one- and two-dimentional equations. Use of the boundary-element method to evaluate the fluid responses in the fractures precludes the need to internally discretize the areal extent of the fractures.     

Classification of Thermal Patterns at Karst Springs and Cave Streams

Thermal patterns of karst springs and cave streams provide potentially useful information concerning aquifer geometry and recharge. Temperature monitoring at 25 springs and cave streams in southeastern Minnesota has shown four distinct thermal patterns. These patterns can be divided into two types: those produced by flow paths with ineffective heat exchange, such as conduits, and those produced by flow paths with effective heat exchange, such as small fractures and pore space. Thermally ineffective patterns result when water flows through the aquifer before it can equilibrate to the rock temperature. Thermally ineffective patterns can be either event‐scale, as produced by rainfall or snowmelt events, or seasonal scale, as produced by input from a perennial surface stream. Thermally effective patterns result when water equilibrates to rock temperature, and the patterns displayed depend on whether the aquifer temperature is changing over time. Shallow aquifers with seasonally varying temperatures display a phase‐shifted seasonal signal, whereas deeper aquifers with constant temperatures display a stable temperature pattern. An individual aquifer may display more than one of these patterns. Since karst aquifers typically contain both thermally effective and ineffective routes, we argue that the thermal response is strongly influenced by recharge mode.     

Conceptual models for carbonate aquifers. 1969

The very diverse types of ground‐water behavior in carbonate terrains can be classified by relating the flow type to a particular hydrogeologic environment each exhibiting a characteristic cave morphology. The ground water may move by diffuse flow, by retarded flow, or by free flow. Diffuse flow occurs in less soluble rocks such as extremely shaley limestones or crystalline dolomites. Integrated conduits are rare. Caves tend to be small, irregular, and often little more than solutionally widened joints. Retarded flows occur in artesian environments and in situations where unfavorable stratigraphy forces ground water to be confined to relatively thin beds. Network cave patterns are characteristic since hydrodynamic forces are damped by the external controls. Solution occurs along many available joints. Free flowing aquifers are those in which solution has developed a subsurface drainage system logically regarded as an underground extension of surface streams. These streams may have fully developed surface tributaries as well as recharge from sinkholes and general infiltration. Characteristic cave patterns are those of integrated conduit systems which are often truncated into linear, angulate, and branchwork caves. Free Flow aquifers may be further subdivided into Open aquifers lying beneath karst plains and Capped aquifers in which significant parts of the drainage net lie beneath an insoluble cap rock. Other geologic factors such as structure, detailed lithology, relief, and locations of major streams, control the details of cave morphology and orientation of the drainage network.     

Spatial characterization of salt accumulation in early stage limestone weathering using probe permeametry

This research characterizes the weathering of natural building stone using an unsteady‐state portable probe permeameter. Variations between the permeability properties of fresh rock and the same rocks after the early stages of a salt weathering simulation are used to examine the effects of salt accumulation on spatial variations in surface rock permeability properties in two limestones from Spain. The Fraga and Tudela limestones are from the Ebro basin and are of Miocene age. Both stone types figure largely in the architectural heritage of Spain and, in common with many other building limestones, they are prone to physical damage from salt crystallization in pore spaces. To examine feedbacks associated with salt accumulation during the early stages of this weathering process, samples of the two stone types were subjected to simulated salt weathering under laboratory conditions using magnesium sulphate and sodium chloride at concentrations of 5% and 15%. Permeability mapping and statistical analysis (aspatial statistics and spatial prediction) before and after salt accumulation are used to assess changes in the spatial variability of permeability and to correlate these changes with salt movement, porosity change, potential rock deterioration and textural characteristics. Statistical analyses of small‐scale permeability measurements are used to evaluate the drivers for decay and hence aid the prediction of the weathering behaviour of the two limestones. Copyright © 2010 John Wiley & Sons, Ltd.     

Effects of mining-induced seismic events on a deep underground mine

The performance of an underground excavation at the Lucky Friday Mine that was subjected to repetitive episodes of mining-induced seismic activities is assessed. Field measurements indicate that an underground excavation may respond to mining-induced seismicity with step changes in displacement of the rock mass around the excavation or step changes in excavation closures. However, step changes in displacement or closure were also observed without the presence of mine seismic events. This behavior may be explained using the concept of stick-slip on joints or bedding planes within the rock mass. The stick-slip mechanism is believed to be controlled by the state of stresses at a joint or bedding plane. Through this mechanism, joint displacement is accumulated in a progressive fashion leading to a weakening of rock mass around an excavation. As a result, the excavation becomes less stable and therefore is more susceptible to seismic impacts. The concept of using peak particle velocity to assess damage threshold for underground excavations may not be conservative when considering the effect of rock mass fatigue.     

Influence of Bedding and Mineral Composition on Mechanical Properties and Its Implication for Hydraulic Fracturing of Shale Oil Reservoirs

The premise of hydraulic fracturing is to have an accurate and detailed understanding of the rock mechanical properties and fracture propagation law of shale reservoirs. In this paper, a comprehensive evaluation of the mechanical properties of the shale oil reservoir in the south of Songliao Basin is carried out. Based on the experiments and the in-situ stress analysis, the fracture propagation law of three types of shale reservoirs is obtained, and the suggestions for fracturing are put forward. The results have shown that the fracture propagation of pure shale and low mature reservoir is easy to open along the bedding plane under compression loading, which is greatly influenced by the bedding. Sand-bearing shale is slightly better, the fractures of which are not easy to open along the bedding plane. The mechanical experimental results show that all the samples have the characteristics of low compressive strength, low Young''s modulus and strong anisotropy, indicating that the shale oil reservoir is certain plastic, which is related to its high clay mineral content and controlled by the bedding development. Compared with pure shale and low mature shale, the sand-bearing shale has less clay content and less developed bedding, which maybe the main reason for its slightly better brittleness. Overall, the expansion of hydraulic fracture is controlled by in-situ stress and bedding. Because of the development of bedding, it is easy to form horizontal fractures. Thus it is not suitable for horizontal well fracturing. Because of the high content of clay minerals, the applicability of conventional slick hydraulic fracturing fluid is poor. It is suggested to use vertical well or directional well to carry out volume fracturing. In this way, the effect of bedding can be effectively used to open and connect the bedding and form a larger fracture network.     

Karst spring responses examined by process-based modeling

Ground water in karst terrains is highly vulnerable to contamination due to the rapid transport of contaminants through the highly conductive conduit system. For contamination risk assessment purposes, information about hydraulic and geometric characteristics of the conduits and their hydraulic interaction with the fissured porous rock is an important prerequisite. The relationship between aquifer characteristics and short-term responses to recharge events of both spring discharge and physicochemical parameters of the discharged water was examined using a process-based flow and transport model. In the respective software, a pipe-network model, representing fast conduit flow, is coupled to MODFLOW, which simulates flow in the fissured porous rock. This hybrid flow model was extended to include modules simulating heat and reactive solute transport in conduits. The application of this modeling tool demonstrates that variations of physicochemical parameters, such as solute concentration and water temperature, depend to a large extent on the intensity and duration of recharge events and provide information about the structure and geometry of the conduit system as well as about the interaction between conduits and fissured porous rock. Moreover, the responses of solute concentration and temperature of spring discharge appear to reflect different processes, thus complementing each other in the aquifer characterization.     

Catchment‐scale conceptual modelling of water and solute transport in the dual flow system of the karst critical zone

Hydrological and hydrochemical processes in the critical zone of karst environments are controlled by the fracture‐conduit network. Modelling hydrological and hydrochemical dynamics in such heterogeneous hydrogeological settings remains a research challenge. In this study, water and solute transport in the dual flow system of the karst critical zone were investigated in a 73.5‐km² catchment in southwest China. We developed a dual reservoir conceptual run‐off model combined with an autoregressive and moving average model with algorithms to assess dissolution rates in the “fast flow” and “slow flow” systems. This model was applied to 3 catchments with typical karst critical zone architectures, to show how flow exchange between fracture and conduit networks changes in relation to catchment storage dynamics. The flux of bidirectional water and solute exchange between the fissure and conduit system increases from the headwaters to the outfall due to the large area of the developed conduits and low hydraulic gradient in the lower catchment. Rainfall amounts have a significant influence on partitioning the relative proportions of flow and solutes derived from different sources reaching the underground outlet. The effect of rainfall on catchment function is modulated by the structure of the karst critical zone (e.g., epikarst and sinkholes). Thin epikarst and well‐developed sinkholes in the headwaters divert more surface water (younger water) into the underground channel network, leading to a higher fraction of rainfall recharge into the fast flow system and total outflow. Also, the contribution of carbonate weathering to mass export is also higher in the headwaters due to the infiltration of younger water with low solute concentrations through sinkholes.     

Characterization of a hydraulically induced crystalline‐rock fracture

Hydraulic fracturing has become an important technique for enhancing the permeability of hydrocarbon source rocks and increasing aquifer transmissivity in many hard rock environments where natural fractures are limited, yet little is known about the nature or behaviour of these hydraulically induced fractures as conduits to flow and transport. We propose that these fractures tend to be smooth based on observed hydraulic and transport behaviour. In this investigation a multi‐faceted approach was used to quantify the properties and characteristics of an isolated hydraulically induced fracture in crystalline rocks. Packers were used to isolate the fracture that is penetrated by two separate observation wells located approximately 33 m apart. A series of aquifer tests and an induced gradient tracer test were performed to better understand the nature of this fracture. Aquifer test results indicate that full recovery is slow because of the overall low permeability of the crystalline rocks. Drawdown tests indicate that the fracture has a transmissivity of 1–2 m²/day and a specific storage on the order of 2–9 × 10^?7/m. Analysis of a potassium–bromide tracer test break through curve shows classic Fickian behaviour with minimal tailing analogous to parallel plate flow. Virtually all of the tracer was recovered, and the breakthrough curve dilution indicates that the swept area is only about 11% of a radial flow field and the estimated aperture is ≤0.5 mm, which implies a narrow linear flow region. These outcomes suggest that transport within these hydraulically induced ‘smooth’ fractures in crystalline rocks is rapid with minimal mixing, small local velocity fluctuations and no apparent diffusion into the host rock or secondary fractures. Copyright © 2016 John Wiley & Sons, Ltd.     

Quantitative measurement of channel-block hydraulic interactions by experimental saturation of a large, natural, fissured rock mass

The hydrodynamic behavior of fissured media relies on the relationships between a few very conductive fractures (channels) and the remaining low-conductivity fractures and matrix (blocks). We made a quantitative measurement of those relationships and their effect on water drainage and storage in a 19,000 m3 natural reservoir consisting of karstified limestones. This reservoir was artificially saturated with water by closing a water gate on the main outlet during a varying time (delta t) fixed by the operator. The water gate was completely or partly closed until the water pressure reached a particular specified value. If the water gate was left completely closed long enough, the water pressure was fixed by the elevation of temporary outlets at the site boundaries. The water elevation within the reservoir was monitored by means of pressure cells located on single fractures representative of the bedding plane and the two families of fractures of the massif network. The comparison of pressure variations with the network geometry allows us to identify a typical double permeability characterized by a few very conductive channels along 10 vertical faults. These channels limit blocks consisting of low-conductivity bedding planes and a rather impervious matrix. Depending on the closure interval, delta t, of the water gate, the total volume of water stored in the reservoir can vary from 0.8 m3 (delta t = 5 minutes) to 18.6 m3 (delta t = 2 days). Such a variance of storage versus closure time is explained by the reservoir's double permeability that is characterized by blocks that saturate much more slowly than channels. If plotted versus time, this injected volume fits a power relationship, according to the saturation state of the blocks. In less than 34 minutes, storage is close to zero in the blocks and to 1.6 to 2 m3 in the channels. For closing times shorter than 1 hour, only 1% of the volume that flows in the channels is stored into the blocks. Depending on the water pressure and for a given delta t = 3000 minutes, the volume of water stored is controlled by the geometry of the joint network and of the aquifer boundaries. Such an experiment shows that the flow is concentrated in about 10% of the fractured network (channels). The water storage that takes place in the 90% remaining fractures (blocks) depends mainly on time during which pressure remains high into the 10% channels. The accurate modeling of such typical double-permeability media needs a careful study of the geometry of the channels whose narrowings modify the flow and induce dam effects that maintain a high pressure gradient with surrounding blocks.     

The distribution rule and seepage effect of the fractures in the ultra-low permeability sandstone reservoir in east Gansu Province,Ordos Basin

To study the impact of the fractures on development in the ultra-low permeability sandstone reservoir of the Yangchang Formation of the Upper Triassic in the Ordos Basin,data on outcrops,cores,slices,well logging and experiments are utilized to analyze the cause of the formation of the fractures,their distribution rules and the control factors and discuss the seepage flow effect of the fractures. In the studied area developed chiefly high-angle tectonic fractures and horizontal bedding fractures,inter-granular fractures and grain boundary fractures as well. Grain boundary fractures and intragranular fractures serve as vital channels linking intragranular pores and intergranular solution pores in the reservoir matrix,thus providing a good connectivity between the pores in the ultra-low perme-ability sandstone reservoir. The formation of fractures and their distribution are influenced by such external and internal factors as the palaeo-tectonic stress field,the reservoir lithological character,the thickness of the rock layer and the anisotropy of a layer. The present-day stress field influences the preservative state of fractures and their seepage flow effect. Under the tec-tonic effect of both the Yanshan and Himalayan periods,in this region four sets of fractures are distributed,respectively assuming the NE-SW,NW-SE,nearly E-W and nearly S-N orientations,but,due to the effect of the rock anisotropy of the rock formation,in some part of it two groups of nearly orthogonal fractures are chiefly distributed. Under the effect of the present-day stress field,the fractures that assume the NE-SW direction have a good connectivity,big apertures,a high permeability and a minimum starting pressure,all of which are main advantages of the seepage fractures in this region. With the development of oilfields,the permeability of the fractures of dif-ferent directions will have a dynamic change.     

The geochemistry of sulphur in a mixed allogenic–autogenic karst catchment,Castleton, Derbyshire,UK

Analyses are presented of anion chemistry and sulphur isotopic compositions of sulphate in sinking streams and groundwaters in a mixed allogenic–autogenic karst catchment. Using the sulphur isotopic data, sources of sulphate from agriculture and the effects of sulphate reduction arising from slurry application can be distinguished from natural rock weathering sources. Within the aquifer, sulphate in known autogenic waters has isotopic compositions distinct from allogenic waters, the autogenic waters being dominated by sulphate from rainfall and rock weathering in these low agricultural intensity catchments. On this basis, water rising at low flow from Whirlpool Rising, Speedwell Cavern, has been identified as dominantly autogenic. Groundwater flow between the sinks and risings in Speedwell Cavern is believed to be along conduits following mineralized faults (rakes). During transit SO42−/Cl− in the water increases. Isotopic mass balance shows that this must be due to addition of sulphate from the oxidation of ore minerals by groundwater. Mass balance considerations show that the present rate of sulphide oxidation must be the result of enhancement by lead mining operations on the rakes. Copyright © 2000 John Wiley & Sons, Ltd.     

Can injection tests reveal the potential for fault movements?

The role of fluids in faulting mechanism and triggering earthquakes is widely accepted. The effective-stress law is the basis for the postulated theories. Using a generalized version of this law, applicable to both continuum and discontinuities, hydromechanical behaviour of a horizontal fracture in a hypothetical fluid-injection problem is investigated. In this problem the increasing intake flow rates, unpredictable by the traditional fluid-flow solutions, brings out another significant aspect of the role of the fluid pressure in rocks. By reducing the magnitude of the compressive effective stresses the fluid pressure causes elastic recovery in fractures. Simple rheologic models are used to demonstrate this fact. Such effects may lead to permeability increases in the rock mass, depending on the magnitude of the fluid pressure. Such variations in permeability, however, are governed by the path dependency of the fracture-deformation response. Therefore, a significant increase in permeability is an indication of comparability of the state of stress and the applied fluid pressure.This index may reveal the potential of hydroactivation of faults, as may arise in the regions of dam reservoirs, underground waste injections, and known faults, for certain ranges of working pressures relevant to each of the above-cited situations. Fluid-injection tests under constant working pressures are suggested as a means revealing the likelihood of movement on the faults.     

The significance of rainsplash in the surficial debris cascade of the colorado front range foothills

The surficial characteristics and precipitation regime of sparsely vegetated hillslopes in the montane zone of the Colorado Front Range suggest that rainsplash may be an important component of the surficial debris cascade. Site sediment flux data from two study periods reveal marked spatial and temporal variability. Comparison of these date with sediment movement data from splashcups suggests the following conclusions: (1) Detachment rates of the rainsplash process appear great enough to account for the sediment flux in the open Gerlach-type traps; (2) Areal extrapolation of rainsplash transport suggest that 88 per cent of the fine sediment flux in 1982 can be attributed to rainsplash; (3) Estimates of rainfall energy and changes in the potential energy of hillslopes by mass transport suggests a process efficiency of 0.05 per cent for rainsplash. If this procedure is applied to sediment flux values from open traps, the low precipitation-energy cascade of 1982 appears to be largely rainsplash-transported sediment. Extrapolation with the 1981 data suggests more aggressive overland flow erosion and transport.     

Experimental investigation of the strength and failure behavior of layered sandstone under uniaxial compression and Brazilian testing

As a typical inherently anisotropic rock, layered sandstones can differ from each other in several aspects, including grain size, type of material, type of cementation, and degree of compaction. An experimental study is essential to obtain and convictive evidence to characterize the mechanical behavior of such rock. In this paper, the mechanical behavior of a layered sandstone from Xuzhou, China, is investigated under uniaxial compression and Brazilian test conditions. The loading tests are conducted on 7 sets of bedding inclinations, which are defined as the angle between the bedding plane and horizontal direction. The uniaxial compression strength (UCS) and elastic modulus values show an undulatory variation when the bedding inclination increases. The overall trend of the UCS and elastic modulus values with bedding inclination is decreasing. The BTS value decreases with respect to the bedding inclination and the overall trend of it is approximating a linear variation. The 3D digital high-speed camera images reveal that the failure and fracture of a specimen are related to the surface deformation. Layered sandstone tested under uniaxial compression does not show a typical failure mode, although shear slip along the bedding plane occurs at high bedding inclinations. Strain gauge readings during the Brazilian tests indicate that the normal stress on the bedding plane transforms from compression to tension as the bedding inclination increases. The stress parallel to the bedding plane in a rock material transforms from tension to compression and agrees well with the fracture patterns; “central fractures” occur at bedding inclinations of 0°–75°, “layer activation” occurs at high bedding inclinations of 75°–90°, and a combination of the two occurs at 75°.     

Fluid flow through rough fractures in rocks. II: A new matching model for rough rock fractures

An understanding of fluid flow through natural fractures in rocks is important in many areas, such as in the hydrocarbon and water industries, and in the safe design of disposal sites for domestic, industrial and nuclear waste. It is often impractical to obtain this information by field or laboratory scale measurements, so numerical modelling of fluid flow must be carried out using synthetic fractures with rough fracture surfaces that are representative of the natural rock fractures. Clearly there are two practical requirements; (i) the development of a method for analysing natural rock fractures to obtain their characteristic parameters, and (ii) the development of techniques for creating high quality synthetic fractures using these parameters. We have implemented these practical requirements in two new software packages. The first, ParaFrac allows the analysis and parameterisation of fracture surfaces and apertures. The second, SynFrac, enables the numerical synthesis of fracture surfaces and apertures with basic prescribed parameters. Synthetic fractures are created using, (i) a new model, which takes full account of the complex matching properties of fracture surfaces using two new parameters, a minimum matching fraction and a transition length and (ii) an improved method of partially correlated random number generation. This model more closely captures the often complex matching properties of real rock fractures than previous more simplified models.     

Integrated cave drip monitoring for epikarst recharge estimation in a dry Mediterranean area,Sif Cave,Israel

Understanding recharge mechanisms and controls in karst regions is extremely important for managing water resources because of the dynamic nature of the system. The objective of this study was to evaluate water percolation through epikarst by monitoring water flow into a cave and conducting artificial irrigation and tracer experiments, at Sif Cave in Wadi Sussi, Israel from 2005 through 2007. The research is based on continuous high‐resolution direct measurements of both rainfall and water percolation in the cave chamber collected by three large PVC sheets which integrate drips from three different areas (17, 46, and 52 m²). Barrels equipped with pressure transducers record drip rate and volume for each of the three areas. The combined measured rainfall and cave data enables estimation of recharge into the epikarst and to better understand the relationship of rainfall‐recharge. Three distinct types of flow regimes were identified: (1) ‘Quick flow’ through preferential flow paths (large fractures and conduits); (2) ‘Intermediate flow’ through a secondary crack system; and (3) ‘Slow flow’ through the matrix. A threshold of ～100 mm of rain at the beginning of the rainy season is required to increase soil water content allowing later rainfall events to percolate deeper through the soil and to initiate dripping in the cave. During winter, as the soil water content rises, the lag time between a rain event and cave drip response decreases. Annual recharge (140–160 mm in different areas in the cave) measured represents 30–35% of annual rainfall (460 mm). Copyright © 2011 John Wiley & Sons, Ltd.     

Karst processes of the eastern upper Galilee,Northern Israel

Karst processes dominate most of the geomorphic activity in the Upper Galilee, consisting mainly of dolomites and limestones. Study of the chemical evolution of water passing through the karst hydrologic cycle clearly shows that the major portion of its carbonate solute is gained subaerially and in the upper part of the vadose zone. Most cave and spring water is already saturated with respect to aragonite and calcite.The karst depressions typical to surface morphology are mostly associated with fault-line traces. Their evolution is possible mainly in areas sloping initially less that 5°.The absence of evolved caves, representing well-developed karst of an earlier period, is attributed mainly to the marginal climate throughout the past combined with tectonic, and hence hydrologic, instability of the region.The discharge of the karst prings shows clearly dependence on annual precipitation, with a lag of about 2 years of the response to drought or more humid periods. Long-term fluctuations are larger in the smaller T'eo Spring than in the affluent 'Enan Springs.Most of the denuded material is extracted from the region as dissolved load via underground conduits and only small amounts as clastics. Mean long-term denudation is approximately 20 mm/1000 years, averaged for the surface area contributing to the springs.In spite of the above, most topographic forms are shaped by runoff erosion, active during medium to high intensity rainstorms. Solution processes prevail during low to medium rainfall intensities, while different parts of the region are denuded at similar rates. Even in karst depressions, erosion becomes dominant after their bottoms are covered by almost impervious terra-rossa mantle.