Quantitative estimation of water storage and residence time in the epikarst with time‐lapse refraction seismic

The hydrodynamic characterization of the epikarst, the shallow part of the unsaturated zone in karstic systems, has always been challenging for geophysical methods. This work investigates the feasibility of coupling time‐lapse refraction seismic data with petrophysical and hydrologic models for the quantitative determination of water storage and residence time at shallow depth in carbonate rocks. The Biot–Gassmann fluid substitution model describing the seismic velocity variations with water saturation at low frequencies needs to be modified for this lithology. I propose to include a saturation‐dependent rock‐frame weakening to take into account water–rock interactions. A Bayesian inversion workflow is presented to estimate the water content from seismic velocities measured at variable saturations. The procedure is tested first with already published laboratory measurements on core samples, and the results show that it is possible to estimate the water content and its uncertainty. The validated procedure is then applied to a time‐lapse seismic study to locate and quantify seasonal water storage at shallow depth along a seismic profile. The residence time of the water in the shallow layers is estimated by coupling the time‐lapse seismic measurements with rainfall chronicles, simple flow equations, and the petrophysical model. The daily water input computed from the chronicles is used to constraint the inversion of seismic velocities for the daily saturation state and the hydrodynamic parameters of the flow model. The workflow is applied to a real monitoring case, and the results show that the average residence time of the water in the epikarst is generally around three months, but it is only 18 days near an infiltration pathway. During the winter season, the residence times are three times shorter in response to the increase in the effective rainfall.     

Impacts of human activities on recharge in a multilayered semiarid aquifer (Campo de Cartagena,SE Spain)

The development of intense agriculture in semiarid areas modifies intensity and spatial distribution of groundwater recharge by summing irrigation return flow to limited rainfall infiltration. Environmental tracers provide key information, but their interpretation is complicated by more complex groundwater flow patterns. In multilayered aquifers, the real origin of the groundwater samples is hard to assess because of local mixing processes occurring inside long‐screened boreholes. We use environmental tracers (¹⁴C, ¹³C, ²H, ¹⁸O, ³H) to investigate the long‐term evolution of recharge in the five‐layer Campo de Cartagena aquifer in South‐Eastern Spain, in addition to high‐resolution temperature loggings to identify the depth of origin of groundwater. Despite the complex background, this methodology allowed a reliable interpretation of the geochemistry and provided a better understanding of the groundwater flow patterns. The tritium method did not give good quantitative results because of the high variability of the recharge signal but remained an excellent indicator of recent recharge. Nonetheless, both pre‐anthropization and post‐anthropization recharge regime could be identified and quantified by radiocarbon. Before the development of agriculture, recharge varied from 17 mm.year^‐1 at the mountain ranges to 6 mm.year^‐1 in the plain, whereas the mean annual rainfall is about 300 mm. In response to the increase of agricultural activity, recharge fluxes to the plain were amplified and nowadays reach up to 210 mm.year^‐1 in irrigated areas. These values are strengthened by global water budget and local unsaturated zone studies. Copyright © 2013 John Wiley & Sons, Ltd.     

Sea Surface Salinity Observations from Space with the SMOS Satellite: A New Means to Monitor the Marine Branch of the Water Cycle

While it is well known that the ocean is one of the most important component of the climate system, with a heat capacity 1,100 times greater than the atmosphere, the ocean is also the primary reservoir for freshwater transport to the atmosphere and largest component of the global water cycle. Two new satellite sensors, the ESA Soil Moisture and Ocean Salinity (SMOS) and the NASA Aquarius SAC-D missions, are now providing the first space-borne measurements of the sea surface salinity (SSS). In this paper, we present examples demonstrating how SMOS-derived SSS data are being used to better characterize key land–ocean and atmosphere–ocean interaction processes that occur within the marine hydrological cycle. In particular, SMOS with its ocean mapping capability provides observations across the world’s largest tropical ocean fresh pool regions, and we discuss from intraseasonal to interannual precipitation impacts as well as large-scale river runoff from the Amazon–Orinoco and Congo rivers and its offshore advection. Synergistic multi-satellite analyses of these new surface salinity data sets combined with sea surface temperature, dynamical height and currents from altimetry, surface wind, ocean color, rainfall estimates, and in situ observations are shown to yield new freshwater budget insight. Finally, SSS observations from the SMOS and Aquarius/SAC-D sensors are combined to examine the response of the upper ocean to tropical cyclone passage including the potential role that a freshwater-induced upper ocean barrier layer may play in modulating surface cooling and enthalpy flux in tropical cyclone track regions.     

Edge Flow and Canopy Structure: A Large-Eddy Simulation Study

Sharp heterogeneities in forest structure, such as edges, are often responsible for wind damage. In order to better understand the behaviour of turbulent flow through canopy edges, large-eddy simulations (LES) have been performed at very fine scale (2 m) within and above heterogeneous vegetation canopies. A modified version of the Advanced Regional Prediction System (ARPS), previously validated in homogeneous conditions against field and wind-tunnel measurements, has been used for this purpose. Here it is validated in a simple forest-clearing-forest configuration. The model is shown to be able to reproduce accurately the main features observed in turbulent edge flow, especially the “enhanced gust zone” (EGZ) present around the canopy top at a few canopy heights downwind from the edge, and the turbulent region that develops further downstream. The EGZ is characterized by a peak in streamwise velocity skewness, which reflects the presence of intense intermittent wind gusts. A sensitivity study of the edge flow to the forest morphology shows that with increasing canopy density the flow adjusts faster and turbulent features such as the EGZ become more marked. When the canopy is characterized by a sparse trunk space the length of the adjustment region increases significantly due to the formation of a sub-canopy wind jet from the leading edge. It is shown that the position and magnitude of the EGZ are related to the mean upward motion formed around canopy top behind the leading edge, caused by the deceleration in the sub-canopy. Indeed, this mean upward motion advects low turbulence levels from the bottom of the canopy; this emphasises the passage of sudden strong wind gusts from the clearing, thereby increasing the skewness in streamwise velocity as compared with locations further downstream where ambient turbulence is stronger.     

Comparative behavior of Sr, Nd and Hf isotopic systems during fluid-related deformation at middle crust levels

We have carried out a comparative Rb-Sr, Sm-Nd and Lu-Hf isotopic study of a progressively deformed hercynian leucogranite from the French Massif Central, belonging to the La Marche ductile shear zone, in order to investigate the respective perturbation of these geochronometers with fluid induced deformation. The one-meter wide outcrop presents a strongly deformed and mylonitized zone at the center, and an asymmetric deformation pattern with a higher deformation gradient on the northern side of the zone. Ten samples have been carefully collected every 10 cm North and South away from the strongest deformed mylonitic zone. They have been analyzed for a complete major, trace element data set, oxygen isotopes, Rb-Sr, Sm-Nd and Lu-Hf isotopic systematics.We show that most of major and trace elements except SiO₂, alkaline elements (K₂O, Rb), and some metal transition elements (Cu), are progressively depleted with increasing deformation. This depletion includes REE + Y, but also HFS elements (Ti, Hf, Zr, Nb) which are commonly considered as immobile elements during upper level processes. Variations in elemental ratios with deformation, e.g. decrease in LREE/MREE- HREE, Nd/Hf, Th/Sr, increase in Rb/Sr, U/Th and constant Sr/Nd, lead to propose the following order of element mobility: U ? Th > Sr = Nd ? Hf + HREE. We conclude in agreement with previous tectonic and metallogenic studies that trace element patterns across the shear zone result from circulation of oxidizing F-rich hydrothermal fluids associated with deformation. A temperature of the fluid of 470-480 °C can be deduced from the δ¹⁸O equilibrium between quartz-muscovite pairs.Elemental fractionation induces perturbation of the Rb-Sr geochronometer. The well-defined ⁸⁷Rb/⁸⁶Sr-⁸⁷Sr/⁸⁶Sr correlation gives an apparent age of 294 ± 19 Ma, slightly younger than the 323 ± 4 Ma age of leucogranites in this area. This apparent age is interpreted as dating event of intense deformation and fluid circulation associated with mass transfer, and exhumation of the ductile crust shortly after the leucogranite emplacement. Sm-Nd and Lu-Hf isochron-type diagrams do not define any correlation, because of the low fractionated Sm/Nd and Lu/Hf ratios. Isotopic data demonstrate that only the Lu-Hf geochronometer system is not affected by fluid circulation and gives reliable T_DM age (1.29 ± 0.03 Ga) and ε_Hf signatures. By contrast, the Sm-Nd geochronometer system gives erroneous old T_DM ages of 2.84-4 Ga. There is no positive ε_Nd-ε_Hf correlation, because of decreasing ε_Nd values with deformation at constant ε_Hf values. However, ε_Nd-ε_Hf values remain in the broad ε_Nd-ε_Hf terrestrial array, which strongly indicates that fluid-induced fractionation can contribute to the width of the terrestrial array. The strong ε_Hf negative values of the leucogranite are similar to metasedimentary granulitic xenoliths from the French Massif Central and confirm the generation of the leucogranite by several episodes of reworking of the lower crust.     

Imprint of North-Atlantic abrupt climate changes on western European loess deposits as viewed in a dust emission model

Western European loess sequences of the last glaciation (100,000–15,000 years BP) exhibit strong, cyclic variations of the sedimentation rate, which are coeval to the Greenland stadial/interstadial cycles and the Heinrich events. These North-Atlantic rapid climate changes appear, thus, as a potential cause for the sedimentation variations, via changes in dust intensity cycle. Here we make a first step in testing this hypothesis, by modelling the impact of the North-Atlantic abrupt climate variations on dust emission. Our dust emission calculations use meteorological fields generated by the LMDZ atmospheric general circulation model at a resolution down to 60 km over Western Europe. Three numerical experiments are run, representing a Greenland stadial, an interstadial and a Heinrich event. Orbital parameters and ice-sheet configuration correspond to conditions from Marine Isotope Stage 3 (60,000–25,000 years BP), a period characterized by strong millennial-scale climate variability. The only differences we impose in the boundary conditions regard the North-Atlantic surface temperature and sea-ice cover in the latitudinal band 30°–63°N. The changes in wind, precipitation, soil moisture and snow cover from one simulated state to another result in small differences in dust emission intensity. In contrast, when the inhibition of the aeolian erosion by vegetation is taken into account, the dust fluxes for the cold climate states (Greenland stadial and Heinrich event) become generally more than twice higher than those for the relatively warmer Greenland interstadial, in agreement with the loess data. These results support the hypothesis that the North-Atlantic millennial-scale variability is imprinted in Western European loess profiles, and point to vegetation changes as the main factor responsible for millennial-scale sedimentation variations. An analysis for the English Channel and southern North Sea areas, major potential dust sources, shows that the seasonality of dust emission is not controlled by the wind speed, as in modern large deserts, but by the surface conditions. Consequently, the dusty season lasts from late winter to early summer, with maximum activity in April–May, and is shifted towards summer when the climate is colder.     

CESAM: a free code for stellar evolution calculations

The cesam code is a consistent set of programs and routines which perform calculations of 1D quasi-hydrostatic stellar evolution including microscopic diffusion of chemical species and diffusion of angular momentum. The solution of the quasi-static equilibrium is performed by a collocation method based on piecewise polynomials approximations projected on a B-spline basis; that allows stable and robust calculations, and the exact restitution of the solution, not only at grid points, even for the discontinuous variables. Other advantages are the monitoring by only one parameter of the accuracy and its improvement by super-convergence. An automatic mesh refinement has been designed for adjusting the localisations of grid points according to the changes of unknowns. For standard models, the evolution of the chemical composition is solved by stiffly stable schemes of orders up to four; in the convection zones mixing and evolution of chemical are simultaneous. The solution of the diffusion equation employs the Galerkin finite elements scheme; the mixing of chemicals is then performed by a strong turbulent diffusion. A precise restoration of the atmosphere is allowed for.     

长江三角洲近900年来气候变化对洪水事件的影响

Wide collection on the historic records of the climatic changes and flood events is performed in the Yangtze Delta. Man-Kendall (MK) method is applied to explore the changing trends of the time series of the flood discharge and the maximum high summer temperature. The research results indicate that the flood magnitudes increased during the transition from the medieval warm interval into the early Little Ice Age. Fluctuating climate changes of the Little Ice Age characterized by arid climate events followed by the humid and cold climate conditions give rise to the frequent flood hazards. Low-lying terrain made the study region prone to the flood hazards, storm tide and typhoon. MK analysis reveals that the jumping point of the time series of the flood discharge changes occurred in the mid-1960s, that of the maximum summer temperature changes in the mid-1990s, and the exact jump point in 1993. The flood discharge changes are on negative trend before the 1990s,they are on positive tendency after the 1990s; the maximum high summer temperature changes are on negative trend before the 1990s and on positive tendency after the 1990s. These results indicate that the trend of flood discharge matches that of the maximum high summer temperature in the Yangtze Delta. The occurrence probability of the maximum high summer temperature will be increasing under the climatic warming scenario and which will in turn increase the occurrence probability of the flood events. More active solar action epochs and the higher sea surface temperature index (SST index) of the south Pacific Ocean area lying between 4°N-4°S and 150°W-90°W correspond to increased annual precipitation, flood discharge and occurrence frequency of floods in the Yangtze Delta. This is partly because the intensified solar activities and the higher SST index give rise to accelerated hydrological circulation from ocean surface to the continent, resulting in increased precipitation on the continent.     

Separating wheat from chaff: Diatom taxon selection usingan artificial neural network pruning algorithm

This study addresses the question of what diatom taxa to includein a modern calibration set based on their relative contribution in apalaeolimnological calibration model. Using a pruning algorithm for ArtificialNeural Networks (ANNs) which determines the functionality of individual taxa interms of model performance, we pruned the Surface Water Acidification Project(SWAP) pH-diatom data-set until the predictive performance of thepruned set (as assessed by a jackknifing procedure) was statistically differentfrom the initial full-set. Our results, based on the validation at each5% data-set reduction, show that (i) 85% of the taxa canbe removed without any effect on the pH model calibration performance, and (ii)that the complexity and the dimensionality reduction of the model by theremoval of these non-essential or redundant taxa greatly improve therobustness of the calibration. A comparison between the commonly usedmarginal criteria for inclusion (species tolerance andHill's N2) and our functionality criterion shows that the importance ofeach taxon in an ANN palaeolimnological model calibration does not appear todepend on these marginal characteristics.