Estimation of vertical water fluxes from temperature time series by the inverse numerical computer program FLUX‐BOT

The application of heat as a hydrological tracer has become a standard method for quantifying water fluxes between groundwater and surface water. The typical application is to estimate vertical water fluxes in the shallow subsurface beneath streams or lakes. For this purpose, time series of temperatures in the surface water and in the sediment are measured and evaluated by a vertical 1D representation of heat transport by advection and conduction. Several analytical solutions exist to calculate the vertical water flux from the measured temperatures. Although analytical solutions can be easily implemented, they are restricted to specific boundary conditions such as a sinusoidal upper temperature boundary. Numerical solutions offer higher flexibility in the selection of the boundary conditions. This, in turn, reduces the effort of data preprocessing, such as the extraction of the diurnal temperature variation from the raw data. Here, we present software to estimate water fluxes based on temperatures—FLUX‐BOT. FLUX‐BOT is a numerical code written in MATLAB that calculates vertical water fluxes in saturated sediments based on the inversion of measured temperature time series observed at multiple depths. FLUX‐BOT applies a centred Crank–Nicolson implicit finite difference scheme to solve the one‐dimensional heat advection–conduction equation. FLUX‐BOT includes functions for the inverse numerical routines, functions for visualizing the results, and a function for performing uncertainty analysis. We present applications of FLUX‐BOT to synthetic and to real temperature data to demonstrate its performance.     

Proposed methods for potential evapotranspiration calculation of the Red River basin (North Vietnam)

Potential evapotranspiration (ETP) is an important part of a climatic water balance and a crucial variable in many kinds of models in computing actual evapotranspiration. The objective of this study was to find a reasonable approach of ETP calculation for a height‐differentiated landscape in subtropical climate. From the pool of diverse approaches, six common methods [Hamon, Priestley–Taylor, Thornthwaite, Blaney–Criddle, Turc and Food and Agricultural Organization Irrigation and Drainage Paper No. 56 (FAO‐56)] were selected. With the meteorological data for long‐term period (1964–2008), the calculation of ETP values was performed for 12 different meteorological stations in the Red River basin. Among the applied ETP calculation methods, the Turc and the FAO‐56 methods agreed well at most stations and represent best the expectations for the ETP values of the Thao and Da subbasins. The findings of our investigations indicate that in highly structured (land use and elevation) regions, not all methods provide satisfying results. Copyright © 2011 John Wiley & Sons, Ltd.     

Indo‐Antarctic derived detritus on the northern margin of Gondwana: evidence for continental‐scale sediment transport

Provenance studies from Cambro‐Ordovician sediments of the North Gondwana passive margin typically ascribe a North African source, a conclusion that cannot be reconciled with all observations. We present new U‐Pb ages from detrital rutile and zircon from Late Ordovician sediments from Saxo‐Thuringia, Germany. Detrital zircons yield age populations of 500–800 Ma, 900–1050 Ma and 1800–2600 Ma. The detrital rutile age spectra are unimodal with ages between 500 and 650 Ma and likely represent, together with the 500–800 Ma and 1800–2600 Ma zircon populations, detritus sourced predominantly from North Africa. In contrast, the c. 950 Ma zircons, which are persistently found in Cambro‐Ordovician sediments of North Gondwana, have no obvious African source. We propose that these zircons are sourced from the Rayner Complex–Eastern Ghats regions of Antarctica and India. An Indo‐Antarctic source indicates either continental‐scale sedimentary transport from central Gondwana to its peripheries or multiple cycles of sediment reworking and redeposition.     

Runoff evolution in the Swiss Alps: projections for selected high‐alpine catchments based on ENSEMBLES scenarios

The Alps are often referred to as the ‘water tower of Europe’. In Switzerland, many branches of the economy, especially the hydropower industry, are closely linked to and dependent on the availability of water. Assessing the impact of climate change on streamflow runoff is, thus, of great interest. Major efforts have already been made in this respect, but the analyses often focus on individual catchments and are difficult to intercompare. In this article, we analysed nine high‐alpine catchments spread over the Swiss Alps, selected for their relevance to a wide range of morphological characteristics. Runoff projections were carried out until the end of the current century by applying the Glacier Evolution Runoff Model (GERM) and climate scenarios generated in the framework of the ENSEMBLES project. We focused on assessing the uncertainty induced by the unknown climate evolution and provided general, statistically based statements, which should be useful as a ‘rule of thumb’ for analyses addressing questions related to water management. Catchments with a high degree of glacierization will undergo the largest changes. General statements about absolute variations in discharge are unreliable, but an overall pattern, with an initial phase of increased annual discharge, followed by a phase with decreasing discharge, is recognizable for all catchments with a significant degree of glacierization. In these catchments, a transition from glacial and glacio‐nival regime types to nival will occur. The timing of maximal annual runoff is projected to occur before 2050 in all basins. The time of year with maximal daily discharges is expected to occur earlier at a rate of 4·4 ± 1·7 days per decade. Compared to its present level, the contribution of snow‐ and icemelt to annual discharge is projected to drop by 15 to 25% until the year 2100. Copyright © 2011 John Wiley & Sons, Ltd.     

Late Holocene Baltic Sea outflow changes reconstructed using C37:4 content from marine cores

The Baltic Sea is an intra‐continental brackish water body. Low saline surface water, the so‐called Baltic outflow current, exits the Baltic Sea through the Kattegat into the Skagerrak. Ingressions of saline oxygen‐rich bottom water enter the Baltic Sea basins via the narrow and shallow Kattegat and are of great importance for the ecological and ventilation state of the Baltic Sea. Over recent decades, progress has been made in studying Holocene changes in saline water inflow. However, reconstructions of past variations in Baltic Sea outflow changes are sparse and hampered because of the lack of suitable proxies. Here, we used the relative proportion of tetra‐unsaturated C₃₇ ketones (C_37:4 %) in long‐chain alkenones produced by coccolithophorids as a proxy for outflowing Baltic Sea water in the Skagerrak. To evaluate the applicability of the proxy, we compared the biomarker results with grain‐size records from the Kattegat and Mecklenburg Bay in addition to previously published salinity reconstructions from the Kattegat over the last 5000 years. All Skagerrak records showed an increase in C_37:4 % that is accompanied by enhanced bottom water currents in the Kattegat and western Baltic Sea over the past 3500 cal. a BP, indicating an increase in Baltic Sea outflow. This probably reflects higher precipitation in the Baltic Sea catchment area owing to a re‐organization of North Atlantic atmospheric circulation with an increased influence of wintertime Westerlies over the Baltic catchment from the mid‐ to the late Holocene.     

In situ U–Pb rutile dating by LA-ICP-MS: <Superscript>208</Superscript>Pb correction and prospects for geological applications

Rutile is a common accessory mineral that occurs in a wide spectrum of metamorphic rocks, such as in blueschists, eclogites, and granulites and as one of the most stable detrital heavy minerals in sedimentary rocks. The advent of rutile trace element thermometry has generated increased interest in a better understanding of rutile formation. This study documents important analytical advances in in situ LA-ICP-MS U/Pb geochronology of rutile: (1) Matrix matching, necessary for robust in situ dating is fulfilled by calibrating and testing several rutile standards (R10, R19, WH-1), including the presentation of new TIMS ages for the rutile standard R19 (489.5 ± 0.9 Ma; errors always stated as 2 s). (2) Initial common lead correction is routinely applied via ²⁰⁸Pb, which is possible due to extremely low Th/U ratios (usually <0.003) in most rutiles. Employing a 213 nm Nd:YAG laser coupled to a quadrupole ICP-MS and using R10 as a primary standard, rutile U/Pb concordia ages for the two other rutile standards (493 ± 10 Ma for R19; 2640 ± 50 Ma for WH-1) and four rutile-bearing metamorphic rocks (181 ± 4 Ma for Ivrea metapelitic granulite; 339 ± 7 Ma for Saidenbach coesite eclogite; 386 ± 8 Ma for Fjortoft UHP metapelite; 606 ± 12 Ma for Andrelandia metepelitic granulite) always agree within 2% with the reported TIMS ages and other dating studies from the same localities. The power of in situ U/Pb rutile dating is illustrated by comparing ages of detrital rutile and zircon from a recent sediment from the Christie Domain of the Gawler Craton, Australia. While the U/Pb age spectrum from zircons show several pronounced peaks that are correlated with magmatic episodes, rutile U/Pb ages are marked by only one pronounced peak (at ca 1,675 Ma) interpreted to represent cooling ages of this part of the craton. Rutile thermometry of the same detrital grains indicates former granulite-facies conditions. The methods outlined in this paper should find wide application in studies that require age information of single spots, e.g., provenance studies, single-crystal zoning and texturally controlled dating.     

Imaging a shallow salt diapir using ambient seismic vibrations beneath the densely built-up city area of Hamburg,Northern Germany

Salt diapirs are common features of sedimentary basins. If close to the surface, they can bear a significant hazard due to possible dissolution sinkholes, karst formation and collapse dolines or their influence on ground water chemistry. We investigate the potential of ambient vibration techniques to map the 3-D roof morphology of shallow salt diapirs. Horizontal-to-vertical (H/V) spectral peaks are derived at more than 900 positions above a shallow diapir beneath the city area of Hamburg, Germany, and are used to infer the depth of the first strong impedance contrast. In addition, 15 small-scale array measurements are conducted at different positions in order to compute frequency-dependent phase velocities of Rayleigh waves between 0.5 and 25 Hz. The dispersion curves are inverted together with the H/V peak frequency to obtain shear-wave velocity profiles. Additionally, we compare the morphology derived from H/V and array measurements to borehole lithology and a gravity-based 3-D model of the salt diapir. Both methods give consistent results in agreement with major features indicated by the independent data. An important result is that H/V and array measurements are better suited to identify weathered gypsum caprocks or gypsum floaters, while gravity-derived models better sample the interface between sediments and homogeneous salt. We further investigate qualitatively the influence of the 3-D subsurface topography of the salt diapir on the validity of local 1-D inversion results from ambient vibration dispersion curve inversion.     

Consistent scale-dependency of future increases in hourly extreme precipitation in two convection-permitting climate models

Climate Dynamics - Convection-permitting models (CPMs) have been proven successful in simulating extreme precipitation statistics. However, when such models are used to study climate change,...     

A submarine landslide complex affecting the Jan Mayen Ridge, Norwegian–Greenland Sea: slide-scar morphology and processes of sediment evacuation

Swath-bathymetry data and 2D multichannel seismics reveal for the first time an up to ~60 km wide amphitheater-shaped slide scar on the eastern flank of the Jan Mayen Ridge, a micro-continent in the Norwegian–Greenland Sea. The scar opens southeastward where it continues as a narrower, topographically controlled translational area. It includes secondary scars, as well as channels and escarpments. Based on the identification of secondary scars, the slide is classified as a slide complex and the total volume of missing sediments was estimated at ~60 km³. From the overall shape of the scar, the upslope widening from a bottleneck- or channel-like bypass-area, the failure is inferred to have had a retrogressive development. The absence of ridges, slabs and sediment blocks indicates that the failed sediments have been evacuated entirely. The smaller channels were formed from single or repetitive smaller flows post-dating the large failure events.