A comparison of hillslope drainage area estimation methods using high-resolution DEMs with implications for topographic studies of gullies

Topographic models provide a useful tool for understanding gully occurrence in the landscape but require reliable estimates of gully head drainage areas. Modern high-resolution topography data (collected using structure from motion photogrammetry or light detection and ranging) is increasingly used for topographic studies of gullies, but little work has been done to assess the variability of gully head drainage area estimates using different methods. This study evaluated alternative approaches to using high-resolution digital elevation models (DEMs) so that gully topographic models can be more readily applied to any area with suitably high-resolution data. Specifically, we investigated the impact of single- or multiple-direction flow routing algorithms, DEM hydrologic-enforcement procedures and spatial resolution on gully head drainage area estimation. We tested these methods on a 40 km² site centred on Weany Creek, a low-relief semi-arid landscape draining towards the Great Barrier Reef, Australia. Using a subroutine to separate gully heads into those with divergent or convergent flow patterns upslope, we found that divergent flow conditions occurred at half of 484 studied gullies. Drainage areas estimated by different flow routing algorithms were more variable in these divergent cases than for convergent cases. This variation caused a significant difference between topographic threshold parameters (slope b and intercept k) derived from single- or multiple-direction flow routing algorithms, respectively. Different methods of hydrologic enforcement (filling or breaching) also affected threshold analysis, resulting in estimates of the exponent b being ~188% higher if the DEM was filled than if breached. The testing of the methods to date indicates that a finer resolution (≤2 m) DEM and a multiple-direction flow routing algorithm achieve the most realistic drainage area estimates in low-relief landscapes. For Weany Creek we estimated threshold parameters k = 0.033 and b = 0.189, indicating that it is highly susceptible to gully erosion.     

Knowledge gaps in our perceptual model of Great Britain's hydrology

There is a no lack of significant open questions in the field of hydrology. How will hydrological connectivity between freshwater bodies be altered by future human alterations to the hydrological cycle? Where does water go when it rains? Or what is the future space–time variability of flood and drought events? However, the answers to these questions will vary with location due to the specific and often poorly understood local boundary conditions and system properties that control the functional behaviour of a catchment or any other hydrologic control volume. We suggest that an open, shared and evolving perceptual model of a region's hydrology is critical to tailor our science questions, as it would be for any other study domain from the plot to the continental scale. In this opinion piece, we begin to discuss the elements of and point out some knowledge gaps in the perceptual model of the terrestrial water cycle of Great Britain. We discuss six major knowledge gaps and propose four key ways to reduce them. While the specific knowledge gaps in our perceptual model do not necessarily transfer to other places, we believe that the development of such perceptual models should be at the core of the debate for all hydrologic communities, and we encourage others to have a similar debate for their hydrologic domain.     

Poro-viscoelasticity modelling based on upscaling quasistatic fluid-saturated solids

In this paper, quasistatic models are developed for the slow flow of compressible fluids through porous solids, where the solid exhibits fading memory viscoelasticity. Problems of this type are important in practical geomechanics contexts, for example, in the context of fluid flow through unconsolidated reservoir sands and of wellbore deformation behaviour in gas and oil shale reservoirs, all of which have been studied extensively. For slow viscous fluid flow in the poro-viscoelastic media we are able to neglect the dynamic effects related to inertia forces, as well as the dissipation associated with the viscous flows. This is in contrast to the vast body of work in the poro-elastic context, where much faster flow of the viscous fluids may give rise to memory effects associated with microflows in pores of the solid media. Such problems have been treated extensively in both the dynamic and quasistatic cases. We are taking a specific type of the porous medium subject to slow deformation processes possibly inducing moderate pressure gradients and flow rates characterised by negligible inertia effects. As the result of homogenisation of such a two-phase medium, we observe the fading memory behaviour in the Biot modulus which controls the pressure increase due to skeleton macroscopic deformation and pore fluid content. Although our derivation leads to a result which is consistent with the formal phenomenological approach proposed by Biot (J Appl Phys 23:1482–1498, 1962), we offer the reader more insight into the structure of the poro-viscoelastic constitutive relations obtained; in particular, we can show that the Biot compressibility evolves in time according to the creep function while the skeleton stiffness is driven by the relaxation function.     

Dating Micrometre‐Thin Rims Using a LA‐ICP‐MS Depth Profiling Technique on Zircon from an Archaean Metasediment: Comparison with the SIMS Depth Profiling Method

Laser ablation‐inductively coupled plasma‐mass spectrometry (LA‐ICP‐MS) was examined as a tool for measuring isotopic variation as a function of ablation depth in unpolished zircon from an Archaean metasediment specimen. This technique was able to identify micrometre‐thin (> 3 μm) isotopically distinct mineral domains characterised by ca. 100 Myr younger ²⁰⁷Pb/²⁰⁶Pb ages associated with 2s age uncertainties as low ~ 0.2%, as well as elevated U content relative to grain interiors (up to an order of magnitude). Our calculated drilling rate suggests that each laser pulse excavated depths of ~ 0.06 μm. Ages resolved through the LA‐ICP‐MS methods overlap the 2s uncertainties of ²⁰⁷Pb/²⁰⁶Pb ages measured using SIMS depth profiling on the same zircon population. The rims were further evinced by the detection of relative enrichment (> 3 orders of magnitude) in REE in the outermost micrometres of the same zircon, measured using a different and independent LA‐ICP‐MS depth profiling technique. We propose a LA‐ICP‐MS U–Pb technique capable of quickly identifying and quantifying rims, which are indication of late, yet geologically significant, fluid events that are otherwise undefined.     

Increased hydrocarbon recovery and CO2 storage in Neogene sandstones,a Croatian example: part II

During 2003–2006, a pilot project of alternating water and CO₂ injection was performed on a limited part of the Upper Miocene sandstone oil reservoir of the Ivani? Field. During the test period oil and gas recovery was significantly increased. Additionally 4,440 m³ of oil and 2.26 × 10⁶ m³ of gas were produced. It has initiated further modelling of sandstone reservoirs in the Ivani? Field in order to calculate volumes available for CO₂ injection for the purpose of increasing hydrocarbon production from depleted sandstone reservoirs in the entire Croatian part of the Pannonian Basin System. In the first phase, modelling was based on results of laboratory testing on the core samples. It considered applying analogies with world-known projects of CO₂ subsurface storage and its usage to enhance hydrocarbon production. In the second phase, reservoir variables were analysed by variograms and subsequently mapped in order to reach lithological heterogeneities and to determine reliable average values of reservoir volumes. Data on porosity, depth and reservoir thickness for the “Gamma 3” and the “Gamma 4” reservoirs, are mapped by the ordinary kriging technique. Calculated volume of CO₂ expressed at standard condition which can be injected in the main reservoirs of the Ivani? Field at near miscible conditions is above 15.5 billion m³.     

Geological interpretation of current subsidence and uplift in the London area,UK, as shown by high precision satellite-based surveying

Long term planning for flood risk management in coastal areas requires timely and reliable information on changes in land and sea levels. A high resolution map of current changes in land levels in the London and Thames estuary area has been generated by satellite-based persistent scatterer interferometry (PSI), aligned to absolute gravity (AG) and global positioning system (GPS) measurements. This map has been qualitatively validated by geological interpretation, which demonstrates a variety of controlling influences on the rates of land level change, ranging from near-surface to deep-seated mechanisms and from less than a decade to more than 100,000 years’ duration.     

Protomylonite evolution potentially revealed by the 3D depiction and fractal analysis of chemical data from a feldspar

Tourmalinization associated with peraluminous granitic intrusions in metapelitic host-rocks has been widely recorded in the Iberian Peninsula, given the importance of tourmaline as a tracer of granite magma evolution and potential indicator of Sn-W mineralizations. In the Penamacor-Monsanto granite pluton (Central Eastern Portugal, Central Iberian Zone), tourmaline occurs: (1) as accessory phase in two-mica granitic rocks, muscovite-granites and aplites, (2) in quartz (±mica)-tourmaline rocks (tourmalinites) in several exocontact locations, and (3) as a rare detrital phase in contact zone hornfels and metapelitic host-rocks. Electron microprobe and stable isotope (δ¹⁸O, δD, δ¹¹B) data provide clear distinctions between tourmaline populations from these different settings: (a) schorl–oxyschorl tourmalines from granitic rocks have variable foititic component (^X□ = 17–57 %) and Mg/(Mg + Fe) ratios (0.19–0.50 in two-mica granitic rocks, and 0.05–0.19 in the more differentiated muscovite-granite and aplites); granitic tourmalines have constant δ¹⁸O values (12.1 ± 0.1 ‰), with wider-ranging δD (?78.2 ± 4.7 ‰) and δ¹¹B (?10.7 to ?9.0 ‰) values; (b) vein/breccia oxyschorl [Mg/(Mg + Fe) = 0.31–0.44] results from late, B- and Fe-enriched magma-derived fluids and is characterized by δ¹⁸O = 12.4 ‰, δD = ?29.5 ‰, and δ¹¹B = ?9.3 ‰, while replacement tourmalines have more dravitic compositions [Mg/(Mg + Fe) = 0.26–0.64], close to that of detrital tourmaline in the surrounding metapelitic rocks, and yield relatively constant δ¹⁸O values (13.1–13.3 ‰), though wider-ranging δD (?58.5 to ?36.5 ‰) and δ¹¹B (?10.2 to ?8.8 ‰) values; and (c) detrital tourmaline in contact rocks and regional host metasediments is mainly dravite [Mg/(Mg + Fe) = 0.35–0.78] and oxydravite [Mg/(Mg + Fe) = 0.51–0.58], respectively. Boron contents of the granitic rocks are low (<650 ppm) compared to the minimum B contents normally required for tourmaline saturation in granitic melts, implying loss of B and other volatiles to the surrounding host-rocks during the late-magmatic stages. This process was responsible for tourmalinization at the exocontact of the Penamacor-Monsanto pluton, either as direct tourmaline precipitation in cavities and fractures crossing the pluton margin (vein/breccia tourmalinites), or as replacement of mafic minerals (chlorite or biotite) in the host-rocks (replacement tourmalinites) along the exocontact of the granite. Thermometry based on ¹⁸O equilibrium fractionation between tourmaline and fluid indicates that a late, B-enriched magmatic aqueous fluid (av. δ¹⁸O ~12.1 ‰, at ~600 °C) precipitated the vein/breccia tourmaline (δ¹⁸O ~12.4 ‰) at ~500–550 °C, and later interacted with the cooler surrounding host-rocks to produce tourmaline at lower temperatures (400–450 °C), and an average δ¹⁸O ~13.2 ‰, closer to the values for the host-rock. Although B-metasomatism associated with some granitic plutons in the Iberian Peninsula seems to be relatively confined in space, extending integrated studies such as this to a larger number of granitic plutons may afford us a better understanding of Variscan magmatism and related mineralizations.     

A relative sea-level history for Arviat,Nunavut, and implications for Laurentide Ice Sheet thickness west of Hudson Bay

Thirty-six new and previously published radiocarbon dates constrain the relative sea-level history of Arviat on the west coast of Hudson Bay. As a result of glacial isostatic adjustment (GIA) following deglaciation, sea level fell rapidly from a high-stand of nearly 170 m elevation just after 8000 cal yr BP to 60 m elevation by the mid Holocene (~ 5200 cal yr BP). The rate of sea-level fall decreased in the mid and late Holocene, with sea level falling 30 m since 3000 cal yr BP. Several late Holocene sea-level measurements are interpreted to originate from the upper end of the tidal range and place tight constraints on sea level. A preliminary measurement of present-day vertical land motion obtained by repeat Global Positioning System (GPS) occupations indicates ongoing crustal uplift at Arviat of 9.3 ± 1.5 mm/yr, in close agreement with the crustal uplift rate inferred from the inferred sea-level curve. Predictions of numerical GIA models indicate that the new sea-level curve is best fit by a Laurentide Ice Sheet reconstruction with a last glacial maximum peak thickness of ~ 3.4 km. This is a 30–35% thickness reduction of the ICE-5G ice-sheet history west of Hudson Bay.     

What causes the spread of model projections of ocean dynamic sea-level change in response to greenhouse gas forcing?

Sea levels of different atmosphere–ocean general circulation models (AOGCMs) respond to climate change forcing in different ways, representing a crucial uncertainty in climate change research. We isolate the role of the ocean dynamics in setting the spatial pattern of dynamic sea-level (ζ) change by forcing several AOGCMs with prescribed identical heat, momentum (wind) and freshwater flux perturbations. This method produces a ζ projection spread comparable in magnitude to the spread that results from greenhouse gas forcing, indicating that the differences in ocean model formulation are the cause, rather than diversity in surface flux change. The heat flux change drives most of the global pattern of ζ change, while the momentum and water flux changes cause locally confined features. North Atlantic heat uptake causes large temperature and salinity driven density changes, altering local ocean transport and ζ. The spread between AOGCMs here is caused largely by differences in their regional transport adjustment, which redistributes heat that was already in the ocean prior to perturbation. The geographic details of the ζ change in the North Atlantic are diverse across models, but the underlying dynamic change is similar. In contrast, the heat absorbed by the Southern Ocean does not strongly alter the vertically coherent circulation. The Arctic ζ change is dissimilar across models, owing to differences in passive heat uptake and circulation change. Only the Arctic is strongly affected by nonlinear interactions between the three air-sea flux changes, and these are model specific.