Historical behaviour of Dome C and Talos Dome (East Antarctica) as investigated by snow accumulation and ice velocity measurements

Ice divide–dome behaviour is used for ice sheet mass balance studies and interpretation of ice core records. In order to characterize the historical behaviour (last 400 yr) of Dome C and Talos Dome (East Antarctica), ice velocities have been measured since 1996 using a GPS system, and the palaeo-spatial variability of snow accumulation has been surveyed using snow radar and firn cores. The snow accumulation distribution of both domes indicates distributions of accumulation that are non-symmetrical in relation to dome morphology. Changes in spatial distributions have been observed over the last few centuries, with a decrease in snow accumulation gradient along the wind direction at Talos Dome and a counter-clockwise rotation of accumulation distribution in the northern part of Dome C. Observations at Dome C reveal a significant increase in accumulation since the 1950s, which could correlate to altered snow accumulation patterns due to changes in snowfall trajectory. Snow accumulation mechanisms are different at the two domes: a wind-driven snow accumulation process operates at Talos Dome, whereas snowfall trajectory direction is the main factor at Dome C. Repeated GPS measurements made at Talos Dome have highlighted changes in ice velocity, with a deceleration in the NE portion, acceleration in the SW portion and migration of dome summit, which are apparently correlated with changes in accumulation distribution. The observed behaviour in accumulation and velocity indicates that even the most remote areas of East Antarctica have changed from a decadal to secular scale.     

The ichnology of carbonate drifts

Carbonate drifts have so far not been as intensely investigated as their siliciclastic equivalents, especially from an ichnological perspective. The aim of this work is therefore to provide an overview of the different bioturbation styles in carbonate drifts for ichnologists and sedimentologists working in such deposits. Different types of carbonate drifts from the Maldives were studied to address this objective. The cores recovered during International Ocean Discovery Program Expedition 359 were examined to provide the sedimentological and ichnological data for a detailed analysis of the ichnology of carbonate drifts. The ichnological characteristics of the Maldives drifts are compared to other carbonate drifts in order to discuss similarities and differences, and thus provide an overview of the general characteristics of carbonate drift ichnology. These drifts are located in the Santaren Channel which lines Great Bahama Bank, along the Marion Plateau in Australia, in the Limassol and Larnaca basins in Cyprus and in the Danish Basin in Denmark. The common characteristics of bioturbation in carbonate drifts are: (i) the complete bioturbation of the sediment with bioturbation indexes between four and six; (ii) the occurrence of distinctive trace fossils limited to facies contacts or condensed intervals; (iii) a typical ichnoassemblage consisting of Thalassinoides, Scolicia, Planolites, Zoophycos, Chondrites, Phycosiphon and Palaeophycus; (iv) the contiguous occurrence of ichnogenera from different tiers, with only Zoophycos and Chondrites as deep tiers; and (v) distinct infills of the traces including particulate organic matter, pyrite, silica and celestine. In addition, the main ichnofacies of carbonate drifts is the Zoophycos ichnofacies. Ichnofabrics grade from coarse‐grained and completely bioturbated to ichnofabrics with present to rare trace fossils and preserved sedimentary structures. The type and intensity of the bioturbation is controlled by the amount of organic matter and the oxygenation at the sea floor that is determined by the action of bottom currents and the sea‐level fluctuations affecting the carbonate factory in carbonate platforms bordering the basins where the carbonate drifts form. The study of the bioturbation in core and outcrop provides palaeoenvironmental information about carbonate‐drift deposits that complement the classical sedimentological data.     

Phosphate removal from water by naturally occurring shale,sandstone, and laterite: The role of iron oxides and of soluble species

Aqueous phosphate removal by three geomaterials from Ivory Coast was evaluated to determine their potential application as low-cost phosphate adsorbents in wastewater treatment. Batch experiments showed that phosphate uptake strongly depended on pH. Laterite and sandstone dissolution was less pronounced compared to shale. A correlation between concentrations of aqueous cation species released from shale and phosphate uptake was observed. The kinetics were well described using the pseudo-second-order model. Isotherms displayed a saturation level on shale, while phosphate uptake continuously increased for laterite and sandstone. The removal efficiency decreased in the following ranking order: laterite > sandstone > shale. Laterite was also the most efficient adsorbent in column experiments. The high phosphate removal efficiency of laterite (8.3 mg PO₄ g^?1) was attributed to the presence of superparamagnetic low grain sizes of goethite. Laterite is a particularly promising material for further investigation in wastewater treatment technology such as constructed wetlands.     

A concerted research effort to advance the hydrological understanding of tropical páramos

Páramos, a neotropical alpine grassland-peatland biome of the northern Andes and Central America, play an essential role in regional and global cycles of water, carbon, and nutrients. They act as water towers, delivering water and ecosystem services from the high mountains down to the Pacific, Caribbean, and Amazon regions. Páramos are also widely recognized as a biodiversity and climate change hot spots, yet they are threatened by anthropogenic activities and environmental changes. Despite their importance for water security and carbon storage, and their vulnerability to human activities, only three decades ago, páramos were severely understudied. Increasing awareness of the need for hydrological evidence to guide sustainable management of páramos prompted action for generating data and for filling long-standing knowledge gaps. This has led to a remarkably successful increase in scientific knowledge, induced by a strong interaction between the scientific, policy, and (local) management communities. A combination of well-established and innovative approaches has been applied to data collection, processing, and analysis. In this review, we provide a short overview of the historical development of research and state of knowledge of the hydrometeorology, flux dynamics, anthropogenic impacts, and the influence of extreme events in páramos. We then present emerging technologies for hydrology and water resources research and management applied to páramos. We discuss how converging science and policy efforts have leveraged traditional and new observational techniques to generate an evidence base that can support the sustainable management of páramos. We conclude that this co-evolution of science and policy was able to successfully cover different spatial and temporal scales. Lastly, we outline future research directions to showcase how sustainable long-term data collection can foster the responsible conservation of páramos water towers.     

Coupled groundwater flow and heat transport simulation for estimating transient aquifer–stream exchange at the lowland River Spree (Germany)

Subsurface flow and heat transport near Freienbrink, NE Germany, was simulated in order to study groundwater–surface water exchange between a floodplains aquifer and a section of the lowland River Spree and an adjacent oxbow. Groundwater exfiltration was the dominant process, and only fast surface water level rises resulted in temporary infiltration into the aquifer. The main groundwater flow paths are identified based on a 3D groundwater flow model. To estimate mass fluxes across the aquifer–surface water interfaces, a 2D flow and heat transport modelling approach along a transect of 12 piezometers was performed. Results of steady‐state and transient water level simulations show an overall high accuracy with a Spearman coefficient ρ = 0.9996 and root mean square error (RMSE) = 0.008 m. Based on small groundwater flow velocities of about 10^?7 to 10^?6 ms^?1, mean groundwater exfiltration rates of 233 l m^?2 d^?1 are calculated. Short periods of surface water infiltration into the aquifer do not exceed 10 days, and the infiltration rates are in the same range. The heat transport was modelled with slightly less accuracy (ρ = 0.8359 and RMSE = 0.34 °C). In contrast to the predominant groundwater exfiltration, surface water temperatures determine the calculated temperatures in the upper aquifer below both surface water bodies down to 10 m during the whole simulation period. These findings emphasize prevailing of heat conduction over advection in the upper aquifer zones, which seems to be typical for lowland streams with sandy aquifer materials and low hydraulic gradients. Moreover, this study shows the potential of coupled numerical flow and heat transport modelling to understand groundwater–surface water exchange processes in detail. Copyright © 2013 John Wiley & Sons, Ltd.     

Spatially distributed tracer‐aided modelling to explore water and isotope transport,storage and mixing in a pristine,humid tropical catchment

Rapidly transforming headwater catchments in the humid tropics provide important resources for drinking water, irrigation, hydropower, and ecosystem connectivity. However, such resources for downstream use remain unstudied. To improve understanding of the behaviour and influence of pristine rainforests on water and tracer fluxes, we adapted the relatively parsimonious, spatially distributed tracer‐aided rainfall–runoff (STARR) model using event‐based stable isotope data for the 3.2‐km² San Lorencito catchment in Costa Rica. STARR was used to simulate rainforest interception of water and stable isotopes, which showed a significant isotopic enrichment in throughfall compared with gross rainfall. Acceptable concurrent simulations of discharge (Kling–Gupta efficiency [KGE] ~0.8) and stable isotopes in stream water (KGE ~0.6) at high spatial (10 m) and temporal (hourly) resolution indicated a rapidly responding system. Around 90% of average annual streamflow (2,099 mm) was composed of quick, near‐surface runoff components, whereas only ~10% originated from groundwater in deeper layers. Simulated actual evapotranspiration (ET) from interception and soil storage were low (~420 mm/year) due to high relative humidity (average 96%) and cloud cover limiting radiation inputs. Modelling suggested a highly variable groundwater storage (~10 to 500 mm) in this steep, fractured volcanic catchment that sustains dry season baseflows. This groundwater is concentrated in riparian areas as an alluvial–colluvial aquifer connected to the stream. This was supported by rainfall–runoff isotope simulations, showing a “flashy” stream response to rainfall with only a moderate damping effect and a constant isotope signature from deeper groundwater (~400‐mm additional mixing volume) during baseflow. The work serves as a first attempt to apply a spatially distributed tracer‐aided model to a tropical rainforest environment exploring the hydrological functioning of a steep, fractured‐volcanic catchment. We also highlight limitations and propose a roadmap for future data collection and spatially distributed tracer‐aided model development in tropical headwater catchments.     

High-P granulites of the Songshugou area (Qinling Orogen,east-central China): Petrography,phase relations,and U/Pb zircon geochronology

High-pressure (HP) granulites provide telling records of mineral reactions at upper mantle to lower crustal levels and key information on the fate of material in subduction systems. The latter especially applies when they abut eclogite and mantle dunite because such rock associations are crucial for understanding the incompletely known processes at the interface of converging plates. A continental arc, active c. 520–395 Ma ago, formed an enigmatic example of such a rock association in the Songshugou area, Qinling Orogen. To unravel the juxtaposition of the distinct rocks, this study combines petrography, phase equilibria modelling, conventional thermobarometry, and zircon U–Th–Pb–Ti–REE analysis. Two mafic HP granulites, which contain the mineral assemblages garnet–clinopyroxene–plagioclase–rutile–mesoperthite–quartz and garnet–clinopyroxene–plagioclase–rutile, experienced peak metamorphic conditions of ≤1.4 GPa, 860°C and ~1.3 GPa, ≥910°C, respectively. During decompression and cooling, at 489 ± 4 Ma, amphibole lamellae unmixed from a clinopyroxene solid solution and orthopyroxene in part replaced garnet. A felsic HP granulite shows equilibration of garnet, perthite, antiperthite, kyanite, quartz, and rutile at 810–860°C, ~1.2 GPa, sillimanite growth during decompression, and upper amphibolite facies cooling at 510 ± 4 Ma. Though the thermobarometric data are just within the methodological errors, the U/Pb zircon ages imply the HP granulites did not evolve coherently. The HP granulites either represent foundered lower arc crust or originated from subduction erosion because their geochemistry is indistinguishable from that of the hanging-wall plate. Published and new pressure–temperature–time–deformation paths converge at ~710°C, ~0.9 GPa, and ≲470 Ma, implying exhumation tectonics juxtaposed the HP granulites with a mélange of eclogite and mantle dunite at lower crustal levels. This study highlights that lower arc crust can comprise material of diverse evolution.     

A rainfall‐runoff model parameterized from GIS and runoff data

A guiding principle in hydrological modelling should be to keep the number of calibration parameters to a minimum. A reduced number of parameters to be calibrated, while maintaining the accuracy and detail required by modern hydrological models, will reduce parameter and model structure uncertainty and improve model diagnostics. In this study, the dynamics of runoff are derived from the distribution of distances from points in the catchments to the nearest stream. This distribution is unique for each catchment and can be determined from a geographical information system. The distribution of distances, will, when a celerity of (subsurface) flow is introduced, provide a distribution of travel times, or a unit hydrograph (UH). For spatially varying levels of saturation deficit, we have different celerities and, hence, different UHs. Runoff is derived from the superposition of the different UHs. This study shows how celerities can be estimated if we assume that recession events represent the combined UHs for different levels of saturation deficit. A new soil moisture routine which estimates saturated and unsaturated volumes of subsurface water and with only one parameter to calibrate is included in the new model. The performance of the new model is compared with that of the Swedish HBV model and is found to perform equally well for eight Norwegian catchments although the number of parameters to be calibrated in the module concerning soil moisture and runoff dynamics is reduced from seven in the HBV model to one in the new model. It is also shown that the new model has a more realistic representation of the subsurface hydrology. Copyright © 2013 John Wiley & Sons, Ltd.     

On the equivalence of spherical splines with least-squares collocation and Stokes’s formula for regional geoid computation

This work is an investigation of three methods for regional geoid computation: Stokes’s formula, least-squares collocation (LSC), and spherical radial base functions (RBFs) using the spline kernel (SK). It is a first attempt to compare the three methods theoretically and numerically in a unified framework. While Stokes integration and LSC may be regarded as classic methods for regional geoid computation, RBFs may still be regarded as a modern approach. All methods are theoretically equal when applied globally, and we therefore expect them to give comparable results in regional applications. However, it has been shown by de Min (Bull Géod 69:223–232, 1995. doi: 10.1007/BF00806734) that the equivalence of Stokes’s formula and LSC does not hold in regional applications without modifying the cross-covariance function. In order to make all methods comparable in regional applications, the corresponding modification has been introduced also in the SK. Ultimately, we present numerical examples comparing Stokes’s formula, LSC, and SKs in a closed-loop environment using synthetic noise-free data, to verify their equivalence. All agree on the millimeter level.