Importance of geological information for assessing drain flow in a Danish till landscape

There has, in recent years, been an increasing interest in developing nutrient load mitigation measures focussing on tile drains. To plan the location of such tile drain measures, it is important to know where in the landscape drain flow is generated and to understand the key factors governing drain flow dynamics. In the present study, we test two approaches to assess spatial patterns in drain flow generation and thereby assess the importance of including geological information. The approaches are the widely used topographical wetness index (TWI), based solely on elevation data, and hydrological models that include the subsurface geology. We set‐up an ensemble of 20 hydrological models based on 20 stochastically generated geological models to predict drain flow dynamics in the clay till Norsminde catchment in Denmark and test the results against TWI. We find that the hydrological models predict observed daily drain flow reasonably well. High drain flow volumes were found in stream valleys and in wetlands and lower drain flow volumes in the more hilly parts of the catchment. In spite of the apparent connection to the landscape, there was no statistically significant correlation between TWI and drain flow at grid scale (100 × 100 m). TWI was therefore not found to be a sufficient index on its own to assess where drain flow is generated, especially in the highlands of the catchment. The geology below 3 m was found to have a large impact on the drain flow, and correlations between sand percentage in the subsurface geology and drain flow volume were found to be statistically significant. Geological uncertainty therefore give rise to uncertainty on simulated drain flow, and this uncertainty was found to be high at the model grid scale but decreasing with increasing scale.     

Catchment chemostasis revisited: Water quality responds differently to variations in weather and climate

Solute concentrations in streamflow typically vary systematically with stream discharge, and the resulting concentration–discharge relationships are important signatures of catchment biogeochemical processes. Solutes derived from mineral weathering often exhibit decreasing concentrations with increasing flows, suggesting dilution of a kinetically limited weathering flux by a variable flux of water. However, previous work showed that concentration–discharge relationships of weathering‐derived solutes in 59 headwater catchments were much weaker than this simple dilution model would predict. Instead, catchments behaved as chemostats, with rates of solute production and/or mobilization that were nearly proportional to water fluxes, on both event and interannual timescales. Here, we re‐examine these findings using data for a wider range of solutes from 2,186 catchments, ranging from ~10 to >1,000,000 km² in drainage area and spanning a wide range of lithologic and climatic settings. Concentration–discharge relationships among this much larger set of larger catchments are broadly consistent with the previously described chemostatic behaviour, at least on event and interannual timescales for weathering‐derived solutes. Among these same catchments, however, site‐to‐site variations in mean concentrations of weathering‐derived solutes exhibit strong negative correlations with long‐term average precipitation and discharge, reflecting strong climatic control on long‐term leaching of the critical zone. We use multiple regression of site characteristics including discharge to identify potential controls on long‐term mean concentrations and find that lithologic and land cover controls are significant predictors for many analytes. The picture that emerges is one in which, on event and interannual timescales, weathering‐derived stream solute concentrations are chemostatically buffered by groundwater storage and fast chemical reactions, but each catchment's chemostatic “set point” reflects site‐to‐site differences in climatically driven evolution of the critical zone. In contrast to these weathering products, some nutrients and particulates are often near‐chemostatic across all timescales, and their long‐term mean concentrations correlate more strongly with land use than climatic characteristics.     

Petrology, magnetostratigraphy and geochronology of the Miocene volcaniclastic Tepoztlán Formation: implications for the initiation of the Transmexican Volcanic Belt (Central Mexico)

The volcaniclastic Tepoztlán Formation (TF) represents an important rock record to unravel the early evolution of the Transmexican Volcanic Belt (TMVB). Here, a depositional model together with a chronostratigraphy of this Formation is presented, based on detailed field observations together with new geochronological, paleomagnetic, and petrological data. The TF consists predominantly of deposits from pyroclastic density currents and extensive epiclastic products such as tuffaceous sandstones, conglomerates and breccias, originating from fluvial and mass flow processes, respectively. Within these sediments fall deposits and lavas are sparsely intercalated. The clastic material is almost exclusively of volcanic origin, ranging in composition from andesite to rhyolite. Thick gravity-driven deposits and large-scale alluvial fan environments document the buildup of steep volcanic edifices. K-Ar and Ar-Ar dates, in addition to eight magnetostratigraphic sections and lithological correlations served to construct a chronostratigraphy for the entire Tepoztlán Formation. Correlation of the 577 m composite magnetostratigraphic section with the Cande and Kent (1995) Geomagnetic Polarity Time Scale (GPTS) suggests that this section represents the time intervall 22.8–18.8 Ma (6Bn.1n-5Er; Aquitanian-Burdigalian, Lower Miocene). This correlation implies a deposition of the TF predating the extensive effusive activity in the TMVB at 12 Ma and is therefore interpreted to represent its initial phase with predominantly explosive activity. Additionally, three subdivisions of the TF were established, according to the dominant mode of deposition: (1) the fluvial dominated Malinalco Member (22.8–22.2 Ma), (2) the volcanic dominated San Andrés Member (22.2–21.3 Ma) and (3) the mass flow dominated Tepozteco Member (21.3–18.8 Ma).     

On the role of domain size and resolution in the simulations with the HIRHAM region climate model

We investigate the simulated temperature and precipitation of the HIRHAM regional climate model using systematic variations in domain size, resolution and detailed location in a total of eight simulations. HIRHAM was forced by ERA-Interim boundary data and the simulations focused on higher resolutions in the range of 5.5–12 km. HIRHAM outputs of seasonal precipitation and temperature were assessed by calculating distributed model errors against a higher resolution data set covering Denmark and a 0.25° resolution data set covering Europe. Furthermore the simulations were statistically tested against the Danish data set using bootstrap statistics. The results from the distributed validation of precipitation showed lower errors for the winter (DJF) season compared to the spring (MAM), fall (SON) and, in particular, summer (JJA) seasons for both validation data sets. For temperature, the pattern was in the opposite direction, with the lowest errors occurring for the JJA season. These seasonal patterns between precipitation and temperature are seen in the bootstrap analysis. It also showed that using a 4,000 × 2,800 km simulation with an 11 km resolution produced the highest significance levels. Also, the temperature errors were more highly significant than precipitation. In similarly sized domains, 12 of 16 combinations of variables, observation validation data and seasons showed better results for the highest resolution domain, but generally the most significant improvements were seen when varying the domain size.     

Scaling Laws in Canopy Flows: A Wind-Tunnel Analysis

An analysis of velocity statistics and spectra measured above a wind-tunnel forest model is reported. Several measurement stations downstream of the forest edge have been investigated and it is observed that, while the mean velocity profile adjusts quickly to the new canopy boundary condition, the turbulence lags behind and shows a continuous penetration towards the free stream along the canopy model. The statistical profiles illustrate this growth and do not collapse when plotted as a function of the vertical coordinate. However, when the statistics are plotted as function of the local mean velocity (normalized with a characteristic velocity scale), they do collapse, independently of the streamwise position and freestream velocity. A new scaling for the spectra of all three velocity components is proposed based on the velocity variance and integral time scale. This normalization improves the collapse of the spectra compared to existing scalings adopted in atmospheric measurements, and allows the determination of a universal function that provides the velocity spectrum. Furthermore, a comparison of the proposed scaling laws for two different canopy densities is shown, demonstrating that the vertical velocity variance is the most sensible statistical quantity to the characteristics of the canopy roughness.     

Pufferungskurven der Gewässer in der Typologie

Ohne Zusammenfassung     

Effect of an oil spill in Spitzbergen in 1978

The oil content in the sediment and the marine life along the arctic shores of Van Mijenfjord, Spitzbergen, were investigated about two years after a spill from diesel storage tanks. High values of oil were recorded in the sediment along the shore near the tanks. The shore fauna is generally poor in these areas and the only biological effect detected was the disappearance of the amphipod Gammarus setosus from the surface layers.     

The Importance of Microbial Iron Sulfide Oxidation for Nitrate Depletion in Anoxic Danish Sediments

Nitrate (NO₃ ^?) reduction processes are important for depleting the NO₃ ^? load from agricultural source areas before the discharge water reaches surface waters or groundwater aquifers. In this study, we experimentally demonstrate the co-occurrence of microbial iron sulfide oxidation by NO₃ ^? (MISON) and other NO₃ ^?-depleting processes in a range of contrasting sediment types: sandy groundwater aquifer, non-managed minerotrophic freshwater peat and two brackish muddy sediments. Approximately 1/3 of the net NO₃ ^? reduction was caused by MISON in three of the four environments despite the presence of organic carbon in the sediment. An apparent salinity limitation to MISON was observed in the most brackish environment. Addition of high surface area synthetically precipitated iron sulfide (FeS _x ) to the aquifer sediment with the lowest natural FeS _x reactivity increased both the relative fraction of NO₃ ^? reduction linked to MISON from approximately 30–100 % and the absolute rates by a factor of 17, showing that the potential for MISON-related NO₃ ^? reduction is environmentally significant and rate limited by the availability of reactive FeS _x .     

Predicting riverine dissolved silica fluxes to coastal zones from a hyperactive region and analysis of their first-order controls

Silicate weathering and resulting transport of dissolved matter influence the global carbon cycle in two ways. First by the uptake of atmospheric/soil CO₂ and second by providing the oceanic ecosystems via the fluvial systems with the nutrient dissolved silica (DSi). Previous work suggests that regions dominated by volcanics are hyperactive or even “hot spots” concerning DSi-mobilization. Here, we present a new approach for predicting DSi-fluxes to coastal zones, emphasizing “first-order” controlling factors (lithology, runoff, relief, land cover and temperature). This approach is applied to the Japanese Archipelago, a region characterized by a high percentage of volcanics (29.1% of surface area). The presented DSi-flux model is based on data of 516 catchments, covering approximately 56.7% of the area of the Japanese Archipelago. The spatial distribution of lithology—one of the most important first order controls—is taken from a new high resolution map of Japan. Results show that the Japanese Archipelago is a hyperactive region with a DSi-yield 6.6 times higher than the world average of 3.3 t SiO₂ km⁻² a⁻¹, but with large regional variations. Approximately 10% of its area exceeds 10 times the world average DSi-yield. Slope constitutes another important controlling factor on DSi-fluxes besides lithology and runoff, and can exceed the influence of runoff on DSi-yields. Even though the monitored area on the Japanese Archipelago stretches from about 31° to 46°N, temperature is not identified as a significant first-order model variable. This may be due to the fact that slope, runoff and lithology are correlated with temperature due to regional settings of the Archipelago, and temperature information is substituted to a certain extent by these factors. Land cover data also do not improve the prediction model. This may partly be attributed to misinterpreted land cover information from satellite images. Implications of results for Earth System and global carbon cycle modeling are discussed.     

Petrographic,fluid inclusion and isotopic study of the Dikulushi Cu–Ag deposit,Katanga (D.R.C.): implications for exploration

The Dikulushi Cu–Ag vein-type deposit is located on the Kundelungu Plateau, in the southeastern part of the Democratic Republic of Congo (D.R.C.). The Kundelungu Plateau is situated to the north of the Lufilian Arc that hosts the world-class stratiform Cu–Co deposits of the Central African Copperbelt. A combined petrographic, fluid inclusion and stable isotope study revealed that the mineralisation at Dikulushi developed during two spatially and temporally distinct mineralising episodes. An early Cu–Pb–Zn–Fe mineralisation took place during the Lufilian Orogeny in a zone of crosscutting EW- and NE-oriented faults and consists of a sequence of sulphides that precipitated from moderate-temperature, saline H₂O–NaCl–CaCl₂-rich fluids. These fluids interacted extensively with the country rocks. Sulphur was probably derived from thermochemical reduction of Neoproterozoic seawater sulphate. Undeformed, post-orogenic Cu–Ag mineralisation remobilised the upper part of the Cu–Pb–Zn–Fe mineralisation in an oxidising environment along reactivated and newly formed NE-oriented faults in the eastern part of the deposit. This mineralisation is dominated by massive Ag-rich chalcocite that precipitated from low-temperature H₂O–NaCl–KCl fluids, generated by mixing of moderate- and low-saline fluids. The same evolution in mineralisation assemblages and types of mineralising fluids is observed in three other Cu deposits on the Kundelungu Plateau. Therefore, the recognition of two distinct types of (vein-type) mineralisation in the study area has a profound impact on the exploration in the Kundelungu Plateau region. The identification of a Cu–Ag type mineralisation at the surface could imply the presence of a Cu–Pb–Zn–Fe mineralisation at depth.