Modelling weathering processes at the catchment scale: The WITCH numerical model

A numerical model of chemical weathering in soil horizons and underlying bedrock (WITCH) has been coupled to a numerical model of water and carbon cycles in forest ecosystems (ASPECTS) to simulate the concentration of major species within the soil horizons and the stream of the Strengbach granitic watershed, located in the Vosges Mountains (France). For the first time, simulations of solute concentrations in soil layers and in the catchment river have been performed on a seasonal basis. The model is able to reproduce the concentrations of most major species within the soil horizons, as well as catching the first-order seasonal fluctuations of aqueous calcium, magnesium and silica concentrations. However, the WITCH model underestimates concentrations of Mg²⁺ and silica at the spring of the catchment stream, and significantly underestimates Ca²⁺ concentration. The deficit in calculated calcium can be compensated for by dissolution of trace apatite disseminated in the bedrock. However, the resulting increased Ca²⁺ release yields important smectite precipitation in the deepest model layer (in contact with the bedrock) and subsequent removal of large amount of silica and magnesium from solution. In contrast, the model accurately accounts for the concentrations of major species (Ca, Mg and silica) measured in the catchment stream when precipitation of clay minerals is not allowed. The model underestimation of Mg²⁺ and H₄SiO₄ concentrations when precipitation of well crystallized smectites is allowed strongly suggests that precipitation of well crystallized clay minerals is overestimated and that more soluble poorly crystallized and amorphous materials may be forming. In agreement with observations on other watersheds draining granitic rocks, this study indicates that highly soluble trace calcic phases control the aqueous calcium budget in the Strengbach watershed.     

Multifractal behaviour of long‐term karstic discharge fluctuations

Karstic watersheds are highly complex hydrogeological systems that are characterized by a multiscale behaviour corresponding to the different pathways of water in these systems. The main issue of karstic spring discharge fluctuations consists in the presence and the identification of characteristic time scales in the discharge time series. To identify and characterize these dynamics, we acquired, for many years at the outlet of two karstic watersheds in South of France, discharge data at 3‐mn, 30‐mn and daily sampling rate. These hydrological records constitute to our knowledge the longest uninterrupted discharge time series available at these sampling rates. The analysis of the hydrological records at different levels of detail leads to a natural scale analysis of these time series in a multifractal framework. From a universal class of multifractal models based on cascade multiplicative processes, the time series first highlights two cut‐off scales around 1 and 16 h that correspond to distinct responses of the aquifer drainage system. Then we provide estimates of the multifractal parameters α and C₁ and the moment of divergence q_D corresponding to the behaviour of karstic systems. These results constitute the first estimates of the multifractal characteristics of karstic spingflows based on 10 years of high‐resolution discharge time series and should lead to several improvements in rainfall‐karstic springflow simulation models. Copyright © 2012 John Wiley & Sons, Ltd.     

Contrasting regional discharge evolutions in the Amazon basin (1974–2004)

Former hydrological studies in the Amazon Basin generally describe annual discharge variability on the main stem. However, the downstream Amazon River only represents the mean state of the Amazonian hydrological system. This study therefore uses a new data set including daily discharge in 18 sub-basins to analyze the variability of regional extremes in the Amazon basin, after recalling the diversity of the hydrological annual cycles within the Amazon basin. Several statistical tests are applied in order to detect trends and breaks in the time series. We show that during the 1974–2004 period, the stability of the mean discharge on the main stem in Óbidos is explained by opposite regional features that principally involve Andean rivers: a decrease in the low stage runoff, particularly important in the southern regions, and an increase in the high stage runoff in the northwestern region. Both features are observed from the beginning of the nineties. These features are also observed in smaller meridian sub-basins in Peru and Bolivia. Moreover we show that the changes in discharge extremes are related to the regional pluriannual rainfall variability and the associated atmospheric circulation as well as to tropical large-scale climatic indicators.     

High-resolution characterization and integrated study of a reservoir formation: the danian carbonate platform in the Aquitaine Basin (France)

This paper provides an example of an integrated multi-scale study of a carbonate reservoir. The Danian Lower R2 carbonate reservoir is located in the South of the Aquitaine Basin (France) and represents a potential underground gas storage site for Gaz de France. The Danian Lower R2 reservoir was deposited as a prograding carbonate platform bordered by a reef barrier. The effects of sedimentary and diagenetic events on the reservoir properties, particularly dolomitization, were evaluated. In this study, the reservoir quality has been assessed by seismic analyses at the basin scale, by log-analysis at the reservoir scale, by petrographic methods and by petrophysical tools at the pore-core scale.Two dolomitization stages, separated by a compaction event with associated fracturing and stylolites, have been identified. These diagenetic events have significantly improved the Lower R2 carbonate reservoir properties. It is demonstrated that the reservoir quality is mainly controlled by the pore-geometry, which is determined by various diagenetic processes. The permeability values of the reservoir range over 4 orders of magnitude, from 0.1 to 5600 mD and the porosity values range between 2 and 42%. Reservoir unit 4 (a karstic dolomite) shows the best reservoir properties with average porosity values ranging between 11.1% and 19.3% and an average permeability ranging between 379 and 766 mD. Reservoir unit 2 (a fine-grained limestone) shows the worst reservoir properties. The cementation factors range from 1.68 to 2.48. The dolomitic crystal carbonate texture (mainly units 3 and 4) shows the highest value of the cementation factor (1.98–2.48) and formation factor (9.54–36.97), which is due to its high degree of cementation. The saturation exponents vary between 1.2 and 3.4. Using these experimental electrical parameters and the resistivity laterolog tool we predicted the water saturation in the various reservoir units. The permeability was predicted by combining the formation factor with the micro-geometric characteristic length. The best fit is obtained with the Katz and Thompson's model and for a constant of 1/171.     

Assessment of climate change impacts on the hydrology of the Peruvian Amazon–Andes basin

In this article, we propose an investigation of the modifications of the hydrological response of two Peruvian Amazonas–Andes basins in relationship with the modifications of the precipitation and evapotranspiration rates inferred by the IPCC. These two basins integrate around 10% of the total area of the Amazonian basin. These estimations are based on the application of two monthly hydrological models, GR2M and MWB3, and the climatic projections come from BCM2, CSMK3 and MIHR models for A1B and B1 emission scenarios (SCE A1B and SCE B1). Projections are approximated by two simple scenarios (anomalies and horizon) and annual rainfall rates, evapotranspiration rates and discharge were estimated for the 2020s (2008–2040), 2050s (2041–2070) and 2080s (2071–2099). Annual discharge shows increasing trend over Requena basin (Ucayali river), Puerto Inca basin (Pachitea river), Tambo basin (Tambo river) and Mejorada basin (Mantaro river) while discharge shows decreasing trend over the Chazuta basin (Huallaga river), the Maldonadillo basin (Urubamba river) and the Pisac basin (Vilcanota river). Monthly discharge at the outlet of Puerto Inca, Tambo and Mejorada basins shows increasing trends for all seasons. Trends to decrease are estimated in autumn discharge over the Requena basin and spring discharge over Pisac basin as well as summer and autumn discharges over both the Chazuta and the Maldonadillo basins. Copyright © 2011 John Wiley & Sons, Ltd.     

Low‐frequency modulation and trend of the relationship between ENSO and precipitation along the northern to centre Peruvian Pacific coast

The relationship between El Niño Southern Oscillation (ENSO) and precipitation along the Peruvian Pacific coast is investigated over 1964–2011 on the basis of a variety of indices accounting for the different types of El Niño events and atmospheric and oceanographic manifestations of the interannual variability in the tropical Pacific. We show the existence of fluctuations in the ENSO/precipitation relationship at decadal timescales that are associated with the ENSO property changes over the recent decades. Several indices are considered in order to discriminate the influence of the two types of El Niño, namely, the eastern Pacific El Niño and the central Pacific El Niño, as well as the influence of large‐scale atmospheric variability associated to the Madden and Julian Oscillation, and of regional oceanic conditions. Three main periods are identified that correspond to the interleave periods between the main climatic transitions over 1964–2011, i.e. the shifts of the 1970s and the 2000s, over which ENSO experiences significant changes in its characteristics. We show that the relationship between ENSO and precipitation along the western coast of Peru has experienced significant decadal change. Whereas El Niño events before 2000 lead to increased precipitation, in the 2000s, ENSO is associated to drier conditions. This is due to the change in the main ENSO pattern after 2000 that is associated to cooler oceanic conditions off Peru during warm events (i.e. central Pacific El Niño). Our analysis also indicates that the two extreme El Niño events of 1982/1983 and 1997/1998 have overshadowed actual trends in the relationship between interannual variability in the tropical Pacific and precipitation along the coast of Peru. Overall, our study stresses on the complexity of the hydrological cycle on the western side of the Andes with regard to its relationship with the interannual to decadal variability in the tropical Pacific. Copyright © 2014 John Wiley & Sons, Ltd.     

Post‐stack stratigraphic inversion workflow applied to carbon dioxide storage: application to the saline aquifer of Sleipner field‡

In the Norwegian North Sea, the Sleipner field produces gas with a high CO₂ content. For environmental reasons, since 1996, more than 11 Mt of this carbon dioxide (CO₂) have been injected in the Utsira Sand saline aquifer located above the hydrocarbon reservoir. A series of seven 3D seismic surveys were recorded to monitor the CO₂ plume evolution. With this case study, time‐lapse seismics have been shown to be successful in mapping the spread of CO₂ over the past decade and to ensure the integrity of the overburden. Stratigraphic inversion of seismic data is currently used in the petroleum industry for quantitative reservoir characterization and enhanced oil recovery. Now it may also be used to evaluate the expansion of a CO₂ plume in an underground reservoir. The aim of this study is to estimate the P‐wave impedances via a Bayesian model‐based stratigraphic inversion. We have focused our study on the 1994 vintage before CO₂ injection and the 2006 vintage carried out after a CO₂ injection of 8.4 Mt. In spite of some difficulties due to the lack of time‐lapse well log data on the interest area, the full application of our inversion workflow allowed us to obtain, for the first time to our knowledge, 3D impedance cubes including the Utsira Sand. These results can be used to better characterize the spreading of CO₂ in a reservoir. With the post‐stack inversion workflow applied to CO₂ storage, we point out the importance of the a priori model and the issue to obtain coherent results between sequential inversions of different seismic vintages. The stacking velocity workflow that yields the migration model and the a priori model, specific to each vintage, can induce a slight inconsistency in the results.     

Basin-scale analysis of rainfall and runoff in Peru (1969–2004): Pacific,Titicaca and Amazonas drainages

Abstract

We present a procedure for estimating Q95 low flows in both gauged and ungauged catchments where Q95 is the flow that is exceeded 95% of the time. For each step of the estimation procedure, a number of alternative methods was tested on the Austrian data set by leave-one-out cross-validation, and the method that performed best was used in the final procedure. To maximise the accuracy of the estimates, we combined relevant sources of information including long streamflow records, short streamflow records, and catchment characteristics, according to data availability. Rather than deriving a single low flow estimate for each catchment, we estimated lower and upper confidence limits to allow local information to be incorporated in a practical application of the procedure. The components of the procedure consist of temporal (climate) adjustments for short records; grouping catchments into eight seasonality-based regions; regional regressions of low flows with catchment characteristics; spatial adjustments for exploiting local streamflow data; and uncertainty assessment. The results are maps of lower and upper confidence limits of low flow discharges for 21 000 sub-catchments in Austria.