Advancing isotope-enabled hydrological modelling for ungauged calibration of data-scarce humid tropical catchments

Realistic projections of the future climate and how this translates to water availability is crucial for sustainable water resource management. However, data availability constrains the capacity to simulate streamflow and corresponding hydrological processes. Developing more robust hydrological models and methods that can circumvent the need for large amounts of hydro-climatic data is crucial to support water-related decisions, particularly in developing countries. In this study, we use natural isotope tracers in addition to hydro-climate data within a newly developed version of the spatially-distributed J2000iso as an isotope-enabled rainfall-runoff model simulating both water and stable isotope (δ²H) fluxes. We pilot the model for the humid tropical San Carlos catchment (2500 km²) in northeastern Costa Rica, which has limited time series, but spatially distributed data. The added benefit of simulating stable isotopes was assessed by comparing different amounts of observation data using three model calibration strategies (i) three streamflow gauges, (ii) three gauges with stream isotopes and (iii) isotopes only. The J2000iso achieved a streamflow Kling–Gupta efficiency (KGE) of 0.55–0.70 across all the models and gauges, but differences in hydrological process simulations emerged when including stable water isotopes in the rainfall-runoff calibration. Hydrological process simulation varied between the standard J2000 rainfall-runoff model with a high simulated surface runoff proportion of 37% as opposed to the isotope version with 84%–89% simulated baseflow or interflow. The model solutions that used only isotope data for calibration exhibited differences in simulated interflow, baseflow and model performance but captured bulk water balances with a reasonable match between the simulated and observed hydrographs. We conclude that J2000iso has shown the potential to support water balance modelling for ungauged catchments using stable isotope, satellite and global reanalysis data sets.     

Evaluating morphological estimates of the aerodynamic roughness of debris covered glacier ice

Aerodynamic roughness length (z₀), the height above the ground surface at which the extrapolated horizontal wind velocity profile drops to zero, is one of the most poorly parameterised elements of the glacier surface energy balance equation. Microtopographic methods for estimating z₀ have become prominent in the literature in recent years, but are rarely validated against independent measures and are yet to be comprehensively analysed for scale or data resolution dependency. Here, we present the results of a field investigation conducted on the debris covered Khumbu Glacier during the post‐monsoon season of 2015. We focus on two sites. The first is characterised by gravels and cobbles supported by a fine sandy matrix. The second comprises cobbles and boulders separated by voids. Vertical profiles of wind speed recorded by a tower comprising five cup anemometers and deployed over both sites enable us to derive measurements of aerodynamic roughness that reflect their observed surface characteristics (0.0184 m and 0.0243 m, respectively). At the second site, z₀ also varied through time following snowfall (0.0055 m) and during its subsequent melt (0.0129 m), showing the importance of fine resolution topography for near‐surface airflow. To compare the wind profile data with microtopographic methods, we conducted structure from motion multi‐view stereo (SfM‐MVS) surveys across each patch and calculated z₀ using three previously published approaches. The fully three‐dimensional cloud‐based approach is shown to be most stable across different scales and these z₀ values are most correct in relative order when compared with the wind tower data. Popular profile‐based methods perform less well providing highly variable values across different scales and when using data of differing resolution. These findings hold relevance for all studies using microtopographic methods to estimate aerodynamic roughness lengths, including those in non‐glacial settings. Copyright © 2017 John Wiley & Sons, Ltd.     

Heterogeneous water storage and thermal regime of supraglacial ponds on debris‐covered glaciers

The water storage and energy transfer roles of supraglacial ponds are poorly constrained, yet they are thought to be important components of debris‐covered glacier ablation budgets. We used an unmanned surface vessel (USV) to collect sonar depth measurements for 24 ponds to derive the first empirical relationship between their area and volume applicable to the size distribution of ponds commonly encountered on debris‐covered glaciers. Additionally, we instrumented nine ponds with thermistors and three with pressure transducers, characterizing their thermal regime and capturing three pond drainage events. The deepest and most irregularly‐shaped ponds were those associated with ice cliffs, which were connected to the surface or englacial hydrology network (maximum depth = 45.6 m), whereas hydrologically‐isolated ponds without ice cliffs were both more circular and shallower (maximum depth = 9.9 m). The englacial drainage of three ponds had the potential to melt ~100 ± 20 × 10³ kg to ~470 ± 90 × 10³ kg of glacier ice owing to the large volumes of stored water. Our observations of seasonal pond growth and drainage with their associated calculations of stored thermal energy have implications for glacier ice flow, the progressive enlargement and sudden collapse of englacial conduits, and the location of glacier ablation hot‐spots where ponds and ice cliffs interact. Additionally, the evolutionary trajectory of these ponds controls large proglacial lake formation in deglaciating environments. Copyright © 2017 John Wiley & Sons, Ltd.     

Secondary Hydrogeologic Regions of the Conterminous United States

The U.S. Geological Survey (USGS) previously identified and mapped 62 Principal Aquifers (PAs) in the U.S., with 57 located in the conterminous states. Areas outside of PAs, which account for about 40% of the conterminous U.S., were collectively identified as “other rocks.” This paper, for the first time, subdivides this large area into internally-consistent features, defined here as Secondary Hydrogeologic Regions (SHRs). SHRs are areas of other rock within which the rocks or deposits are of comparable age, lithology, geologic or physiographic setting, and relationship to the presence or absence of underling PAs or overlying glacial deposits. A total of 69 SHRs were identified. The number and size of SHRs identified in this paper are comparable to the number and size of PAs previously identified by the USGS. From a two-dimensional perspective, SHRs are complementary to PAs, mapped only where the PAs were not identified on the USGS PA map and not mapped where the PAs were identified. SHRs generally consist of low permeability rocks or deposits, but can include locally productive aquifers. The two maps, taken together, provide a comprehensive, national-scale hydrogeologic framework for assessing and understanding groundwater systems.     

A multifrequency maximum-entropy joint analysis of COBE and Tenerife data

We use data from the Tenerife 10-, 15- and 33‐GHz beam-switching experiments along with the COBE 53- and 90‐GHz data to separate the cosmic microwave background (CMB) signal from the Galactic signal, and create two maps at high Galactic latitude. The new multi-MEM technique is used to obtain the best reconstruction of the two channels. The two maps are presented, and known features are identified within each. We find that the Galactic contribution to both the 15- and 33-GHz Tenerife data is small enough to be ignored when compared with the errors in the data and the magnitude of the CMB signal.     

Subsidies to high seas bottom trawl fleets and the sustainability of deep-sea demersal fish stocks

The life spans of demersal species of fishes occurring in deep-waters are much longer and their potential growth rates much lower than those of related shallow water species. As a result, deep-sea demersal fish species are more vulnerable to exploitation. This is because low growth rates relative to the available market discount rate for capital makes it desirable for fishing firms to mine, rather than sustainably exploit, these resources even in the absence of fisheries subsidies. However, it is common knowledge that governments around the world do provide subsidies to their fishing industries. The objective of this contribution is to estimate the global amount of subsidies paid to bottom trawl fleets operating in the high seas, i.e., outside of the Exclusive Economic Zones of maritime countries. Our study suggests that fisheries subsidies to these fleets stand at about US$152 million per year, which constitutes 25% of the total landed value of the fleet. Economic data for bottom trawlers suggest that the profit achieved by this vessel group is normally not more than 10% of landed value. The implication of this finding is that without subsidies, the bulk of the world's bottom trawl fleet operating in the high seas will be operating at a loss, and unable to fish, thereby reducing the current threat to deep-sea and high seas fish stocks.     

Diffusion in silicate melts: II. Multicomponent diffusion in CaOAl2O3SiO2 at 1500°C and 1 GPa

Chemical diffusion profiles in molten CaO---Al₂O₃---SiO₂ have been measured over a large range of compositions at 1500°C and l GPa. The diffusion profiles have been inverted for effective binary diffusion coefficients (EBDCs) and for the chemical diffusion matrix. The EBDCs are shown to depend strongly on both composition and direction of diffusion in composition space. The dependence of EBDCs on direction in composition space, which for the system studied here can be as large as a factor of seven, severely limits the applicability of EBDCs to interdiffusion in any direction other than the one used to derive the EBDCs.

The chemical diffusion matrix for molten CaO---Al₂O₃---SiO₂ was determined using diffusion profiles from two or three mutually orthogonal diffusion couples in the ternary composition space. All features of the diffusion profiles shown in this work can be reproduced by representing the chemical fluxes in the three-component system as a linear combination of concentration gradients via a 2 × 2 diffusion matrix. Chemical diffusion in molten CaO---Al₂O₃---SiO₂ shows clear evidence of strong diffusive coupling among the components. This can be seen in the uphill diffusion profiles of components that were initially uniform, in the fact that the apparent rate of diffusion of some components is a strong function of direction in composition space, and most quantitatively in the magnitude of off-diagonal elements of the diffusion matrix relative to the magnitude of the diagonal elements. SiO₂ for example, is found to be strongly coupled with CaO in relatively silicic melts, whereas Al₂O₃ is strongly coupled with CaO in less silicic melts. Furthermore, the coupling of CaO with either Al₂O₃ or Si02 reverses sign between more and less polymerized compositions. Interdiffusion profiles in natural melts have numerous features that suggest similar coupling between Al₂O₃ and CaO and between SiO₂ and CaO.     

Olivine/silicate melt partitioning of germanium: an example of a nearly constant partition coefficient

The partitioning of germanium between forsterite (Fo) and liquids in the diopside-anorthiteforsterite join was investigated by electron microprobe analysis of Ge-doped samples equilibrated at 1300°–1450°C. Germanium is somewhat incompatible in Fo relative to the haplobasaltic melts, with a grand mean for all simple partition coefficients (D_Fo-l^Ge) of 0.68 ± 0.06. For the melt composition range studied, D_Fo-l^Ge is virtually constant in isothermal series of experiments, and shows only minor overall temperature dependence. The exchange reaction partition coefficient $\text{K}{}_{\mathrm{D}}\text{= (}\text{Mg}{}_2\text{GeO}{}_4\text{)}{}_{\mathrm{Fo}}\text{(SiO}{}_2\text{)}{}_{\mathrm{l}}\text{(Mg}{}_2\text{SiO}{}_4\text{)}{}_{\mathrm{Fo}}\text{(GeO}{}_2\text{)}{}_{\mathrm{l}}$] is near unity in all cases, with a grand mean of 0.93 ± 0.11. One exploratory run at 20 kbar yielded a distinctly lower partition coefficient (D_Fo-l^Ge = 0.54 ± 0.04), which confirms the negative pressure dependence predicted by the thermodynamics of Ge ai Si exchange.These new data indicate that absolute Ge enrichment must occur in terrestrial magmas undergoing olivine fractionation, while $\text{Ge}\text{Si}$ remains nearly constant.     

A Back-arc Palaeotectonic Setting for the Augaro Neoproterozoic Magmatic Rocks of Western Eritrea

The Augaro volcano-sedimentary assemblages of western Eritrea are part of the Neoproterozoic, N-S trending belt of low-grade volcano-sedimentary and associated plutonic rocks. In contrast to the volcanic-dominated oceanic-arc assemblages in central Eritrea, the predominant rock types in the west are supracrustal sequences of sedimentary origin with subordinate volcanic rocks. These Augaro supracrustal rocks are overlain, unconformably, by a basin-fill metasedimentary succession known as the Gulgula Group. The Augaro metavolcanic rocks are tholeiitic and range in composition from basalt to basaltic andesite. Comparison of trace element characteristics and N-MORB-normalised spidergrams of these rocks with those of modern volcanic environments and age-comparable metavolcanic rocks of known tectonic association from the Arabian-Nubian Shield suggest that the volcanic assemblages from western Eritrea were generated in a back-arc tectonic setting.

Single zircon Pb-Pb evaporation and vapour-transfer U-Pb analyses of magmatic zircons from pre/syn-tectonic granites yield a mean ²⁰⁷Pb/²⁰⁶Pb age of 849±20 Ma and an upper concordia intercept age of 849±26 Ma. These ages are interpreted to represent the time of major magmatism in western Eritrea and are comparable to ages of early arc magmatism in central and northern Eritrea and in the southern Nubian Shield. Initial eNd values and initial Sr isotope ratios of whole-rock samples of magmatic rocks calculated for an age of 850 Ma range from +4.0 to +7.1 and 0.7026 to 0.7037, respectively. Single zircon ²⁰⁷Pb/²⁰⁶Pb ages, initial eNd value and Sr isotope ratio for a granitic clast in the Gulgula metaconglomerate suggest that the source area for the Gulgula metasedimentary rocks is similar to the surrounding Neoproterozoic rocks of western Eritrea.     

The periglacial origin of the drifts at Morfa-Bychan,near Aberystwyth

The superficial deposits at Morfa-bychan have previously been interpreted as deposited by Welsh ice. A reconsideration, including: stone orientation analysis, suggests however that they consist of screes and solifluction deposits with subordinate rain-wash gravels, and a thin loess—the whole laid down under periglacial conditions, when the Irish Sea ice no longer reached the area. There is no direct evidence for dating but their volume and field relations suggest they represent the deposits of the last glaciation.