Extreme storm surge modelling in the North Sea

This study shows that storm surge model performance in the North Sea is mostly unaffected by the application of temporal variations of surface drag due to changes in sea state provided the choice of a suitable constant Charnock parameter in the sea-state-independent case. Including essential meteorological features on smaller scales and minimising interpolation errors by increasing forcing data resolution are shown to be more important for the improvement of model performance particularly at the high tail of the probability distribution. This is found in a modelling study using WAQUA/DCSMv5 by evaluating the influence of a realistic air-sea momentum transfer parameterization and comparing it to the influence of changes in the spatial and temporal resolution of the applied forcing fields in an effort to support the improvement of impact and climate analysis studies. Particular attention is given to the representation of extreme water levels over the past decades based on the example of the Netherlands. For this, WAQUA/DCSMv5 is forced with ERA-Interim reanalysis data. Model results are obtained from a set of different forcing fields, which either (i) include a wave-state-dependent Charnock parameter or (ii) apply a constant Charnock parameter (α _{C h} =?0.032) tuned for young sea states in the North Sea, but differ in their spatial and/or temporal resolution. Increasing forcing field resolution from roughly 79 to 12 km through dynamically downscaling can reduce the modelled low bias, depending on coastal station, by up to 0.25 m for the modelled extreme water levels with a 1-year return period and between 0.1 m and 0.5 m for extreme surge heights.     

Why and when it is useful to publish and share inconclusive results and failures: reply to “Reporting negative results to stimulate experimental hydrology”*

We thank van Emmerik et al. for their discussion of our opinion paper. We do not fully agree that publication of negative results will be very effective in promoting experimental hydrology. Instead, we think that experimental hydrology is considered to be more risky than modelling studies because of difficulties in publishing inconclusive results and the potential for failure.     

Clinker formation in basaltic and trachybasaltic lava flows

Clinker is a term used to describe massive or scoriaceous fragments commonly associated with ‘a‘ā lava flows. Clinker is generally considered to form by fragmentation of an upper vesiculated crust, due to an increase in apparent viscosity and/or to an increase in shear strain rate. Surface clinker is considered to be transported to the flow front and incorporated at the base by caterpillar motion. Clinker that we have observed on a variety of lava flows has very variable textures, which suggests several different mechanisms of formation. In order to study clinker formation, we examined several lava flows from the Chaîne des Puys Central France, where good sections, surface morphology and surface textures are widespread and clearly visible. We observed basal and surface ‘a‘ā clinker that has fragmentation textures similar to those observed in ash formed in eruptions under dry conditions. In two pāhoehoe flows we have observed basal clinker that formed in-situ. Two other flows display clinker features identical to those commonly observed in phreatomagmatic ash, such as adhering particles, blocky shapes, spherical glass and attached microphenocrysts. Another pāhoehoe flow has a flakey, angular basal breccia, with microfaulted and abraded clasts. These were probably formed at a cooled lava base by large amounts of simple shear and consequent intra-lava brittle faulting. Using these observations we propose three different ways of fragmentation. (1) Clinker can form at the surface and eventually produce roll-over basal breccia. (2) Water/lava interactions can form basal clinker by phreatomagmatic fragmentation. Water/lava ratio variations may produce different clinker structures, in a manner similar to observed textural changes in phreatomagmatic eruptions. (3) Clinker can be formed by brittle brecciation during basal simple shear. The different clinker can provide information about the mechanisms and environmental conditions during lava flow emplacement.     

Rapid assessment of field erosion and sediment transport pathways in cultivated catchments after heavy rainfall events

This paper explores a scale‐adapted erosion mapping method which aims at a rapid assessment of field erosion and sediment transport pathways in catchments up to several square kilometres and compares the results with the output of a well‐known erosion model (LISEM). The mapping method is based on an event‐defined classification scheme of erosion intensity (zero, weak, moderate and strong) that is applied to arable fields, in combination with incision measurements of erosion features for each erosion intensity class on a small sample of fields. Sediment deposition is classified on the basis of quantity indicators and abundance. In addition, relevant conditions and erosion factors are determined for each field. The method was applied to an agricultural catchment (4·2 km²) in the Sundgau (Alsace), after a short but violent thunderstorm in May 2001, to illustrate its potential use and its limitations. The rainfall event led to strong erosion on the arable fields and a muddy flow that caused significant damage in the built‐up area. On the basis of the analyses of the incision measurements in combination with the mapping of erosion intensity classes, total erosion for the catchment was estimated as 15 000 t (an average of about 36 t[sol ]ha). Sediment deposition was found to occur in three major locations: (1) in thalwegs at the interface between maize and downslope winter wheat fields, (2) in downslope headlands where the flow direction suddenly changed due to oriented tillage structures in the perpendicular direction, and (3) the lowest corners of fields which collect all the runoff from the field. Preliminary data analyses suggest that erosion intensity is related to field size and[sol ]or tillage direction and to slope morphology. Model output (LISEM) appeared to depend more strongly on slope gradient than the results obtained with the mapping method. The method yields a database, which can be used as a foundation for conservation strategies in small regions with similar land use and geomorphology. The mapping and modelling methods are compared, and their complementary aspects are highlighted. Copyright © 2005 John Wiley & Sons, Ltd.     

Preservation of the last aggradation phase in climate-driven terraces: Evidence from Late Quaternary reach-specific fluvial dynamics of the Allier River (France)

There is limited knowledge about the preservation of aggradation phases in Quaternary fluvial records. Previous numerical modelling of erosion and deposition in Late Quaternary Allier River (France) generated the prediction that this river has reach-specific fluvial dynamics related to climate-driven tributary sediment-flux dynamics. To test this prediction, new optically stimulated luminescence (OSL) samples were collected of the Late Quaternary Fx terrace at five locations along a stretch of 60 km. OSL dates of both quartz and feldspar sand grains indicate that all relatively basalt-poor sediments display significantly different ages for each reach (ranging from 36.3 ± 2.0 to 21.1 ± 2.3 ka). The more basalt-rich terrace body consistently yields ages in the range 21.1 ± 1.7 to 16.1 ± 1.5 ka, suggesting contemporaneous aggradation along the whole studied Allier reach during this interval. Our own new OSL date of a Tartaret eruption around 16.8 ± 2.5 ka also fits this time window, suggesting a direct link with volcanic activity. However, there are many more dated volcanic events that coincide with the older basalt-poor units, making it less likely that a direct link between terrace-sediment basalt content and volcanic activity exists. The timings of the dated depositional events in MIS 3 and 2 all match with simulated climate drivers and published landscape erosion rates. Counterintuitively, the volcanic Chaîne des Puys area supplied more sediment during the cold and dry Last Glacial Maximum. Basalt content in the Allier terrace sediments reflects climate-related sediment-flux dynamics upstream. The scarcity of older basalt-poor sediment bodies from MIS 4 and 3 in the Fx terrace suggests that less sediment was supplied and/or the intermittent erosional phases in the Allier were very effective at removing them. We hypothesize that this observation of predominant preservation of the last aggradation phase could be a common phenomenon in most climate-driven terraces. © 2020 The Authors. Earth Surface Processes and Landforms published by John Wiley & Sons Ltd     

The value of remotely sensed surface soil moisture for model calibration using SWAT

Remotely sensed (RS) data can add value to a hydrological model calibration. Among this, RS soil moisture (SM) data have mostly been assimilated into conceptual hydrological models using various transformed variable or indices. In this study, raw RS surface SM is used as a calibration variable in the Soil and Water Assessment Tool model. This means the SM values were not transformed into another variable (e.g., soil water index and root zone SM index). Using a nested catchment, calibration based only on RS SM and optimizing model parameters sensitive to SM using particle swarm optimization improved variations in streamflow predictions at some of the gauging stations compared to the uncalibrated model. This highlighted part of the catchments where the SM signal directly influenced the flow distribution. Additionally, highlighted high and low flow signals were mostly influenced. The seasonal breakdown indicates that the SM signal is more useful for calibrating in wetter seasons and in areas with higher variations in elevation. The results identified that calibration only on RS SM improved the general rainfall–runoff response simulation by introducing delays but cannot correct the overall routing effect. Furthermore, catchment characteristics (e.g., land use, elevation, soil types, and precipitation) regulating SM variation in different seasons highlighted by the model calibration are identified. This provides further opportunities to improve model parameterization.     

Effective stresses and shear failure pressure from in situ Biot's coefficient,Hejre Field,North Sea

We propose a combination of Biot's equations for effective stress and the expression for shear failure in a rock to obtain an expression for minimum pore pressure in a stable vertical well bore. We show that a Biot's coefficient calculated from logging data in the Hejre Field, North Sea, is significantly different from 1. The log‐derived Biot's coefficient is above 0.8 in the Shetland Chalk Group and in the Tyne Group, and 0.6–0.8 in the Heno Sandstone Formation. We show that the effective vertical and horizontal stresses obtained using the log‐derived Biot's coefficient result in a drilling window for a vertical well larger than if approximating Biot's coefficient by 1. The estimation of the Biot's coefficient is straightforward in formations with a stiff frame, whereas in formations such as shales, caution has to be taken. We discuss the consequence of assumptions made on the mineral composition of shales as unphysical results could be obtained when choosing inappropriate mineral moduli.