Hydrological response characteristics of Mediterranean catchments at different time scales: a meta-analysis

ABSTRACT

This work examines 140 hydrological studies conducted in the Mediterranean region. It identifies key characteristics of the hydrological responses of Mediterranean catchments at various time scales and compares different methods and modelling approaches used for individual-catchment studies. The study area is divided into the northwestern (NWM), eastern (EM) and southern (SM) Mediterranean. The analysis indicates regional discrepancies in which the NWM shows the most extreme rainfall regime. A tendency for reduced water resources driven by both anthropogenic and climatic pressures and a more extreme rainfall regime are also noticeable. Catchments show very heterogeneous responses over time and space, resulting in limitations in hydrological modelling and large uncertainties in predictions. However, few models have been developed to address these issues. Additional studies are necessary to improve the knowledge of Mediterranean hydrological features and to account for regional specificities.

Editor D. Koutsoyiannis Associate editor A. Efstratiadis     

Tectono‐sedimentary evolution of the western Corinth rift (Central Greece)

The Corinth rift (Greece) is one of the world's most active rifts. The early Plio‐Pleistocene rift is preserved in the northern Peloponnese peninsula, south of the active Corinth rift. Although chronostratigraphic resolution is limited, new structural, stratigraphic and sedimentological data for an area >400 km² record early rift evolution in three phases separated by distinct episodes of extension rate acceleration and northward fault migration associated with major erosion. Minimum total N–S extension is estimated at 6.4–7.7 km. The earliest asymmetrical, broad rift accommodated slow extension (0.6–1 mm a^?1) over >3 Myrs and closed to the west. North‐dipping faults with throws of 1000–2200 m defined narrow blocks (4–7 km) with little footwall relief. A N‐NE flowing antecedent river system infilled significant inherited relief (Lower group). In the earliest Pleistocene, significant fluvial incision coincided with a 15 km northward rift margin migration. Extension rates increased to 2–2.5 mm a^?1. The antecedent rivers then built giant Gilbert‐type fan deltas (Middle group) north into a deepening lacustrine/marine basin. N‐dipping, basin margin faults accommodated throws <1500 m. Delta architecture records initiation, growth and death of this fault system over ca. 800 ka. In the Middle Pleistocene, the rift margin again migrated 5 km north. Extension rate increased to 3.4–4.8 mm a^?1. This transition may correspond to an unconformity in offshore lithostratigraphy. Middle group deltas were uplifted and incised as new hangingwall deltas built into the Gulf (Upper group). A final increase to present‐day extension rates (11–16 mm a^?1) probably occurred in the Holocene. Fault and fault block dimensions did not change significantly with time suggesting control by crustal rheological layering. Extension rate acceleration may be due to strain softening or to regional tectonic factors.     

Direct dating of thick‐ and thin‐skin thrusts in the Peruvian Subandean zone through apatite (U–Th)/He and fission track thermochronometry

Although the structure of the central Peruvian Subandean zone is well defined, the timing of thrust‐related exhumation and Cenozoic sedimentation remain poorly constrained. In this study, we report new apatite (U–Th)/He (AHe) and fission track (AFT) ages from thrust‐belt and foreland strata along three published balanced cross sections. AHe data from the northern, thick‐skinned domain (i.e. Shira Mountain, Otishi Cordillera and Ucayali Basin) show young AHe ages (ranging from 2.6 ± 0.2 to 13.1 ± 0.8 Ma) compared with AFT ages (ranging from 101 ± 5 to 133 ± 11 Ma). In the southern Camisea Basin, where deformation is mainly thin‐skinned, AHe and AFT ages have been both reset and show young cooling ages (3.7 ± 0.8 Ma and 8 ± 2 Ma respectively). Using low‐temperature thermochronology data and the latest fission track annealing and He diffusion codes, the thermal history of the study area has been reconstructed using inverse modelling. This history includes two steps of erosion: Early Cretaceous and late Neogene, but only Neogene sedimentation and exhumation varies in the different sectors of the study area. From a methodological point of view, large AHe data dispersion point to the need for refinement of AHe damage and annealing models. The influence of grain chemistry on damage annealing, multiple age components and the possibility of fission tracks as traps for He need further consideration. For the central Peruvian Subandes, AHe and AFT ages combined with balanced cross sections emphasize the dominant control of Paleozoic inheritance rather than climate on Cenozoic infilling and exhumation histories. Finally, our data provide the first field example of how thick‐skinned thrust‐related deformation and exhumation in the Subandes can be directly dated through AHe thermochronology.     

Structure and nuclear density distribution in the cheralite—CaTh(PO4)2: studies of its behaviour under high pressure (36?GPa)

The crystal structure of the cheralite—CaTh(PO₄)₂—has been revisited by neutron diffraction and its behaviour under high pressure investigated by X-ray diffraction up to 36?GPa. The neutron diffraction data at ambient pressure gave a more accurate determination of the Ca/Th cation position than previous XRD data, taking advantage that the neutron scattering lengths of calcium and thorium are of same order of magnitude. The nuclear density distribution was also determined using the maximum entropy method (MEM) confirming that the two cations are not located at the same position in the unit cell but are slightly displaced from one another along a specific direction in order to minimize the electrostatic repulsion with the surrounding phosphorus atoms. At high pressure, the compound did not show any phase transition or amorphization. From the evolution of the unit-cell volume as a function of the pressure, the zero-pressure bulk modulus B₀ and its pressure derivative B₀ ^′ have been determined by fitting the experimental compressibility curve to the Birch–Murnaghan equation of state. The results are B₀?=?140(2) GPa and B ₀ ^′ ?=?4.4(4) GPa.     

Micro- and nano-textural evidence of Ti(–Ca–Fe) mobility during fluid–rock interactions in carbonaceous lawsonite-bearing rocks from New Zealand

Understanding the mobility of chemical elements during fluid–rock interactions is critical to assess the geochemical evolution of a rock undergoing burial and metamorphism and, more generally, to constrain the geochemical budget of the subduction factory. In particular, determining the behavior and mobility of Ti in aqueous fluids constitutes a great challenge that is still under scrutiny. Here, we study plant fossils preserved in blueschist metasedimentary rocks from the Marybank Formation (New Zealand). Using scanning and transmission electron microscopies (SEM and TEM), we show that the carbonaceous material (CM) composing the fossils contains abundant nano-inclusions of Ti- and Fe-oxides. These nanocrystals are mainly anatase, rutile, and Fe–Ti oxides. The mineral composition observed within the fossils is significantly different from that detected in the surrounding rock matrix. We propose that Ti and Fe might have been mobilized by the alteration of a detrital Ti–Fe-rich protolith during an early diagenetic event under acidic and reducing conditions. Aqueous fluids rich in organic ligands released by the degradation of organic matter may have been involved. Moreover, using mass balance and petrological observations, we show that the contrasted mineralogy between the rock matrix and the fossil CM might be the consequence of the chemical isolation of fossil CM during the prograde path of the rock. Such an isolation results from the early formation of quartz and Fe-rich phyllosilicate layers enclosing the fossil as characterized by SEM and TEM investigations. Overall, this study shows that investigating minerals associated with CM down to the nanometer scale in metamorphic rocks can provide a precious record of early prograde geochemical conditions.     

Links between NAO fluctuations and inter-annual variability of winter-months precipitation in the Seine River watershed (north-western France)

The inter-annual to multi-decadal winter variability (DJFM) of precipitation on the Seine River watershed (France) was analysed using continuous wavelet transform analysis and compared to the winter North Atlantic Oscillation Index (NAOI). Nine weather stations were used over the 1951 to 2004 period and confirmed the homogeneity of inter-annual fluctuations for all stations but one. Wavelet coherence between SLP over the Icelandic and Azores regions and precipitation highlighted coherence for different scales of variability according to the centre of action considered. Segmentation and wavelet analysis and coherence between precipitation and NAOI over a long period of time (1873–2004) showed: i) increasing variability across the last century at most time scales, especially for NAOI; ii) the existence of change points for the mean and variance of both signals; iii) overall discontinuity of the coherence whatever the scale considered, especially between ∼1910 and ∼1955 for inter-decennial to pluri-decennial scales.     

Improving the assessment of ICESat water altimetry accuracy accounting for autocorrelation

Given that water resources are scarce and are strained by competing demands, it has become crucial to develop and improve techniques to observe the temporal and spatial variations in the inland water volume. Due to the lack of data and the heterogeneity of water level stations, remote sensing, and especially altimetry from space, appear as complementary techniques for water level monitoring. In addition to spatial resolution and sampling rates in space or time, one of the most relevant criteria for satellite altimetry on inland water is the accuracy of the elevation data. Here, the accuracy of ICESat LIDAR altimetry product is assessed over the Great Lakes in North America. The accuracy assessment method used in this paper emphasizes on autocorrelation in high temporal frequency ICESat measurements. It also considers uncertainties resulting from both in situ lake level reference data. A probabilistic upscaling process was developed. This process is based on several successive ICESat shots averaged in a spatial transect accounting for autocorrelation between successive shots. The method also applies pre-processing of the ICESat data with saturation correction of ICESat waveforms, spatial filtering to avoid measurement disturbance from the land–water transition effects on waveform saturation and data selection to avoid trends in water elevations across space. Initially this paper analyzes 237 collected ICESat transects, consistent with the available hydrometric ground stations for four of the Great Lakes. By adapting a geostatistical framework, a high frequency autocorrelation between successive shot elevation values was observed and then modeled for 45% of the 237 transects. The modeled autocorrelation was therefore used to estimate water elevations at the transect scale and the resulting uncertainty for the 117 transects without trend. This uncertainty was 8 times greater than the usual computed uncertainty, when no temporal correlation is taken into account. This temporal correlation, corresponding to approximately 11 consecutive ICESat shots, could be linked to low transmitted ICESat GLAS energy and to poor weather conditions. Assuming Gaussian uncertainties for both reference data and ICESat data upscaled at the transect scale, we derived GLAS deviations statistics by averaging the results at station and lake scales. An overall bias of −4.6 cm (underestimation) and an overall standard deviation of 11.6 cm were computed for all lakes. Results demonstrated the relevance of taking autocorrelation into account in satellite data uncertainty assesment.     

Influence of the orographic roughness of glacier valleys across the Transantarctic Mountains in an atmospheric regional model

Glacier valleys across the Transantarctic Mountains are not properly taken into account in climate models, because of their coarse resolution. Nonetheless, glacier valleys control katabatic winds in this region, and the latter are thought to affect the climate of the Ross Sea sector, frsater formation to snow mass balance. The purpose of this paper is to investigate the role of the production of turbulent kinetic energy by the subgrid-scale orography in the Transantarctic Mountains using a 20-km atmospheric regional model. A classical orographic roughness length parametrization is modified to produce either smooth or rough valleys. A one-year simulation shows that katabatic winds in the Transantarctic Mountains are strongly improved using smooth valleys rather than rough valleys. Pressure and temperature fields are affected by the representation of the orographic roughness, specifically in the Transantarctic Mountains and over the Ross Ice Shelf. A smooth representation of escarpment regions shows better agreement with automatic weather station observations than a rough representation. This work stresses the need to improve the representation of subgrid-scale orography to simulate realistic katabatic flows. This paper also provides a way of improving surface winds in an atmospheric model without increasing its resolution.     

Influence of coupling on atmosphere, sea ice and ocean regional models in the Ross Sea sector, Antarctica

Air–sea ice–ocean interactions in the Ross Sea sector form dense waters that feed the global thermohaline circulation. In this paper, we develop the new limited-area ocean–sea ice–atmosphere coupled model TANGO to simulate the Ross Sea sector. TANGO is built up by coupling the atmospheric limited-area model MAR to a regional configuration of the ocean–sea ice model NEMO. A method is then developed to identify the mechanisms by which local coupling affects the simulations. TANGO is shown to simulate realistic sea ice properties and atmospheric surface temperatures. These skills are mostly related to the skills of the stand alone atmospheric and oceanic models used to build TANGO. Nonetheless, air temperatures over ocean and winter sea ice thickness are found to be slightly improved in coupled simulations as compared to standard stand alone ones. Local atmosphere ocean feedbacks over the open ocean are found to significantly influence ocean temperature and salinity. In a stand alone ocean configuration, the dry and cold air produces an ocean cooling through sensible and latent heat loss. In a coupled configuration, the atmosphere is in turn moistened and warmed by the ocean; sensible and latent heat loss is therefore reduced as compared to the stand alone simulations. The atmosphere is found to be less sensitive to local feedbacks than the ocean. Effects of local feedbacks are increased in the coastal area because of the presence of sea ice. It is suggested that slow heat conduction within sea ice could amplify the feedbacks. These local feedbacks result in less sea ice production in polynyas in coupled mode, with a subsequent reduction in deep water formation.     

Mo isotope composition in Mo-rich high- and low-T hydrothermal systems from the Swiss Alps

We analyzed the molybdenum (Mo) isotope compositions (IC) of 59 samples from two molybdenite mineralizations (Alpjahorn and Grimsel) and from a Mo-rich hydrothermal breccia (Grimsel) from the Aar Massif, Switzerland. The formation temperature of the Late Paleozoic Mo mineralizations (300-600 °C) is much higher than that of the Pliocene breccia (100-160 °C). The Mo IC of the molybdenites varies over 1.35‰. Even in a single hand specimen it spans 0.45‰, indicating that fractionation processes during molybdenite precipitation can vary on a cm scale. The Mo IC of most molybdenites analyzed here are significantly heavier than that of the host rock (δ^98/95Mo = (0.05 ± 0.1)‰) and show a bimodal distribution centered around δ^98/95Mo ≈ 1.1‰ and 0.2‰. This result rules out single stage Rayleigh fractionation as the relevant formation mechanism and instead, redox variations are suggested to be a main factor controlling the Mo IC of the studied high-temperature Mo deposits. The range of the Mo IC in one single deposit, the Alpjahorn, overlaps with the variation range of almost all other published values for Mo IC in Mo deposits. Compared to the molybdenites, the breccia shows an even wider variation of 3.0‰ (δ^98/95Mo between −1.6‰ and +1.4‰). In contrast to the high-T molybdenite deposits, here the Mo was transported via oxidized surface waters into the breccia system, where it was reduced and precipitated. This indicates that oxidation and reduction of Mo complexes may lead to highly variable Mo IC in hydrothermal systems.