Soil cracking effects on hydrological and erosive processes: a study case in Mediterranean cultivated vertisols

Shrink–swell soils, such as those in a Mediterranean climate regime, can cause changes in terms of hydrological and erosive responses because of the changing soil water storage conditions. Only a limited number of long‐term studies have focused on the impacts on both hydrological and erosive responses and their interactions in an agricultural environment. In this context, this study aims to document the dynamics of cracks, runoff and soil erosion within a small Mediterranean cultivated catchment and to quantify the influence of crack processes on the water and sediment supplied to a reservoir located at the catchment outlet using water and sediment measurements at a cultivated field outlet as baseline. Detailed monitoring of the presence of topsoil cracks was conducted within the Kamech catchment (ORE OMERE, Tunisia), and runoff and suspended sediment loads were continuously measured over a long period of time (2005–2012) at the outlets of a field (1.32 ha) and a catchment (263 ha). Analysis of the data showed that topsoil cracks were open approximately half of the year and that the rainfall regime and water table level conditions locally control the seasonal cracking dynamics. Topsoil cracks appeared to seriously affect the generation of runoff and sediment concentrations and, consequently, sediment yields, with similar dynamics observed at the field and catchment outlets. A similar time lag in the seasonality between water and sediment delivery was observed at these two scales: although the runoff rates were globally low during the presence of topsoil cracks, most sediment transport occurred during this period associated with very high sediment concentrations. This study underlines the importance of a good prediction of runoff during the presence of cracks for reservoir siltation considerations. In this context, the prediction of cracking effects on runoff and soil erosion is a key factor for the development of effective soil and water management strategies and downstream reservoir preservation. Copyright © 2016 John Wiley & Sons, Ltd.     

Subsurface DC resistivity mapping: approximate 1-D interpretation

Resistivity prospecting is the main tool used to investigate the shallow structure of the ground. A series of new techniques for determining the 2-D and 3-D geometry of the ground is now finding increasing use, but the light and simple Wenner prospecting technique remains a practical and efficient tool for rapidly mapping lateral variations in resistivity. When the resistivity changes are smooth, 1-D modelling can be used to interpret the data, and the criteria governing this approximation can be defined from synthetic data generated by a 3-D slab-model. For a Wenner array, two quadripole configurations can be used, Normal and Dipole-Dipole. For these two configurations the width of the transition zone, the apparent anisotropy effect and the precision of the resistivity values recovered from 1-D inversion differ. However the simultaneous inversion of both sets of data gives better results than for either configuration by itself. Two examples illustrate that in geological contexts where the thickness of the weathered zone causes the changes in the apparent resistivity value, this parameter can be recovered from 1-D inversion.     

Internal erosion of chemically reinforced granular materials: a mathematical modeling approach

Internal erosion (IE) affects the stability of natural and reinforced materials by causing instability within their granular structure. The dislodgement and transport of eroded particles affect both the particulate concentration of eroding fluid and the pore network of eroded material. In this study, we examined these modifications using a transport model with a finite element code. First, IE tests on chemically reinforced sand columns were performed to obtain information about eroded material loss of mass, particulate concentration of effluent, porosity and permeability modifications, and existing IE stages. Second, based on experimental results, a mathematical one‐dimensional model has been formulated to monitor the evolution and spatial distribution of erodible solids, fluidized particles, porosity, permeability, and seepage stresses. The model consists of a set of coupled nonlinear differential equations solved in sequence. It provides valuable information about the extent and the dynamics of structural changes, which can be used to estimate an IE time for the hydraulic work to reach failure. Copyright © 2011 John Wiley & Sons, Ltd.     

An automated salt‐tracing gauge for flow‐velocity measurement

This article introduces the SVG (salt‐velocity gauge), a novel automated technique for measuring flow velocity by means of salt tracing. SVG allows a high measuring rate (up to one every 2 seconds), short control section length (down to 10 cm), high accuracy (+[sol ]?1·5 cm s^?1), and unbiased calculation of the mean velocity in experimental conditions with turbulent, supercritical flow. A few cubic centimetres of saturated salt solution (NaCl) are injected into the flow at regular time intervals using a programmable solenoid valve. The tracer successively passes two conductivity probes placed a short distance downstream. The transformation of the signal between the two probes is modelled as a one‐dimensional diffusion wave equation. Model calibration gives an estimation of the mean velocity and the diffusion for each salt plume. Two implementations of the SVG technique are described. The first was an outdoors simulated rainfall experiment in Senegal (conductivity probes at 40 cm apart, 8 Hz measurement rate, salt injections at 10 second intervals). Mean velocity was estimated to range between 0·1 and 0·3 m s^?1. The second was a laboratory‐based flume experiment (conductivity probes at 10 cm apart, 32 Hz, salt injections at 2 second intervals). Another SVG with probes at 34 cm apart was used for comparison. An acoustic Doppler velocimeter (ADV) was also used to give an independent assessment of velocity. Using the 10 cm salt gauge, estimated mean velocity ranged from 0·6 to 0·9 m s^?1 with a standard deviation of 1·5 cm s^?1. Comparisons between ADV, 10 cm SVG and 34 cm SVG were consistent and demonstrated that the salt‐tracing results were unbiased and independent of distance between probes. Most peaks were modelled with r² > 90 per cent. The SVG technology offers an alternative to the dye‐tracing technique, which has been severely criticized in the literature because of the wide interval of recommended values for the correction factor α to be applied to the timings. This article demonstrates that a fixed value of α is inappropriate, since the correction factor varies with velocity, diffusion and the length of the control section. Copyright © 2005 John Wiley & Sons, Ltd.     

Displacements and transformations of nitrate-rich and nitrate-poor water masses in the tropical Pacific during the 1997 El Niño

A Lagrangian analysis was applied to the outputs of a coupled physical-biogeochemical model to describe the redistribution of nitrate-rich and nitrate-poor surface water masses in the tropical Pacific throughout the major 1997 El Niño. The same tool was used to analyze the causes of nitrate changes along trajectories and to investigate the consequences of the slow nitrate uptake in the high nutrient low chlorophyll (HNLC) region during the growth phase of the event. Three patterns were identified during the drift of water masses. The first mechanism is well known along the equator: oligotrophic waters from the western Pacific are advected eastward and retain their oligotrophic properties along their drift. The second concerns the persistent upwelling in the eastern basin. Water parcels have complex trajectories within this retention zone and remain mesotrophic. This study draws attention to the third process which is very specific to the HNLC region and to the El Niño period. During the 1997 El Niño, horizontal and vertical inputs of nitrate decreased so dramatically that nitrate uptake by phytoplankton became the only mechanism driving nitrate changes along pathways. The study shows that because of the slow nitrate uptake characteristic of the tropical Pacific HNLC system, nitrate in the pre-El Niño photic layer can support biological production for a period of several months. As a consequence, the slow nitrate uptake delays the gradual onset of oligotrophic conditions over nearly all the area usually occupied by upwelled waters. Owing to this process, mesotrophic conditions persist in the tropical Pacific during El Niño events.     

Extreme value modelling of daily areal rainfall over Mediterranean catchments in a changing climate

Heavy rainfall events during the fall season are causing extended damages in Mediterranean catchments. A peaks‐over‐threshold model is developed for the extreme daily areal rainfall occurrence and magnitude in fall over six catchments in Southern France. The main driver of the heavy rainfall events observed in this region is the humidity flux (FHUM) from the Mediterranean Sea. Reanalysis data are used to compute the daily FHUM during the period 1958–2008, to be included as a covariate in the model parameters. Results indicate that the introduction of FHUM as a covariate can improve the modelling of extreme areal precipitation. The seasonal average of FHUM can improve the modelling of the seasonal occurrences of heavy rainfall events, whereas daily FHUM values can improve the modelling of the events magnitudes. In addition, an ensemble of simulations produced by five different general circulation models are considered to compute FHUM in future climate with the emission scenario A1B and hence to evaluate the effect of climate change on the heavy rainfall distribution in the selected catchments. This ensemble of climate models allows the evaluation of the uncertainties in climate projections. By comparison to the reference period 1960–1990, all models project an amplification of the mean seasonal FHUM from the Mediterranean Sea for the projection period 2070–2099, on average by +22%. This increase in FHUM leads to an increase in the number of heavy rainfall events, from an average of 2.55 events during the fall season in present climate to 3.57 events projected for the period 2070–2099. However, the projected changes have limited effects on the magnitude of extreme events, with only a 5% increase in the median of the 100‐year quantiles. Copyright © 2011 John Wiley & Sons, Ltd.     

An enhanced constitutive model for crushable granular materials

Studies in the past have tried to reproduce the mechanical behaviour of granular materials by proposing constitutive relations based on a common assumption that model parameters and parameters describing the properties, including gradation of individual grains are inevitably linked. However successful these models have proved to be, they cannot account for the changes in granular assembly behaviour if the grains start to break during mechanical loading. This paper proposes to analyse the relation between grading change and the mechanical behaviour of granular assembly. A way to model the influence of grain breakage is to use a critical state‐based model. The influence of the amount of grain breakage during loading, depending on the individual grain strength and size distribution, can be introduced into constitutive relations by means of a new parameter that controls the evolution of critical state with changes in grain size distribution. Experimental data from a calcareous sand, a quartz sand, and a rockfill material were compared with numerical results and good‐quality simulations were obtained. The main consequences of grain breakage are increased compressibility and a gradual dilatancy disappearance in the granular material. The critical state concept is also enriched by considering its overall relation to the evolution of the granular material. Copyright © 2009 John Wiley & Sons, Ltd.     

Fault zone geometry of a mature active normal fault: A potential high permeability channel (Pirgaki fault, Corinth rift, Greece)

We present results from petrophysical analysis of a normal fault zone with the aim of defining the flow pathways and their behavior during seismic and interseismic periods. Data are obtained on porosity geometry, strain structure and mineralogy of different domains of a normal fault zone in the Corinth rift. Data point out a close relationship between mineralogy of the clayey minerals, porosity network and strain structures and allow definition of a macroscopic anisotropy of the flow parameters with a strong control by microscopic ultracataclasite structures. The Pirgaki fault zone, developed within pelagic limestone, has a sharp asymmetric porosity profile, with a high porosity volume in the fault core and in the damage zone of the hanging wall. From porosity volumes and threshold measurements, a matrix permeability variation of 6 orders of magnitude could be expected between the protolith and the fault core. Modifications of this pathway during seismic and interseismic phases are depicted. Healing of cracks formed during seismic slip events occurred in the fault core zone and the porous network in the damage zone is sealed in a second step. The lens geometry of the fault core zone is associated with dissolution surfaces and open conduits where dissolved matter could move out of the fault core zone. These elementary processes are developed in particular along Riedel's structures and depend on the orientation of the strain surfaces relative to the local stress and depend also on the roughness of each surface type. P-surfaces are smooth and control shearing process. R-surfaces are rough and present two wavelengths of roughness. The long one controls localization of dissolution surfaces and conduits; the short one is characteristic of dissolution surfaces. The dissolved matter can precipitate in the open structures of the hanging wall damage zone, decreasing the connectivity of the macroscopic conduit developed within this part of the fault zone.     

Snow cover monitoring in the Northern Patagonia Icefield using MODIS satellite images (2000–2006)

The snow cover of the Northern Patagonia Icefield (NPI) was monitored after applying the Normalized Difference Snow Index (NDSI) and the Red/NIR band ratio to 134 Moderate Resolution Imaging Spectroradiometer (MODIS) images captured between 2000 and 2006. The final results show that the snow cover extent of the NPI fluctuates a lot in winter, in addition to its seasonal behaviour. The minimum snow cover extent of the period (3600 km²) was observed in March 2000 and the maximum (11,623 km²) in August 2001. We found that temperature accounts for approximately 76% of the variation of the snow cover extent over the entire icefield. We also show two different regimes of winter snow cover fluctuations corresponding to the eastern and the western sides of the icefield. The seasonality of the snow cover on the western side was determined by temperature rather than precipitation, while on the east side the seasonality of the snow cover was influenced by the seasonal behaviour of both temperature and precipitation. This difference can be explained by the two distinct climates: coastal and continental. The fluctuations in the winter snow cover extent were more pronounced and less controlled by temperature on the western side than on the eastern side of the icefield. Snow cover extent was correlated with temperature R² = 0.75 and R² = 0.74 for the western and eastern sides, respectively. Since limited meteorological data are available in this region, our investigation confirmed that the change in snow cover is an interesting climatic indicator over the NPI providing important insights in mass balance comprehension. Since snow and ice were distinguished snow cover fluctuations can be associated to fluctuations in the snow accumulation area of the NPI. In addition, days with minimum snow covers of summer season can be associated to the period in which Equilibrium Line Altitude (ELA) is the highest.