Refine
Year of publication
- 2019 (79) (remove)
Document Type
- Article (55)
- Doctoral Thesis (7)
- Postprint (7)
- Other (6)
- Part of a Book (1)
- Habilitation Thesis (1)
- Report (1)
- Review (1)
Is part of the Bibliography
- yes (79) (remove)
Keywords
- climate change (3)
- hydrology (3)
- Abandonment (2)
- Biodiversity (2)
- Determinants (2)
- Drivers (2)
- Ecosystem services (2)
- GLOF (2)
- Intensity (2)
- National-Parks (2)
Institute
- Institut für Umweltwissenschaften und Geographie (79) (remove)
Bank filtration (BF) is an established indirect water-treatment technology. The quality of water gained via BF depends on the subsurface capture zone, the mixing ratio (river water versus ambient groundwater), spatial and temporal distribution of subsurface travel times, and subsurface temperature patterns. Surface-water infiltration into the adjacent aquifer is determined by the local hydraulic gradient and riverbed permeability, which could be altered by natural clogging, scouring and artificial decolmation processes. The seasonal behaviour of a BF system in Germany, and its development during and about 6 months after decolmation (canal reconstruction), was observed with a long-term monitoring programme. To quantify the spatial and temporal variation in the BF system, a transient flow and heat transport model was implemented and two model scenarios, 'with' and 'without' canal reconstruction, were generated. Overall, the simulated water heads and temperatures matched those observed. Increased hydraulic connection between the canal and aquifer caused by the canal reconstruction led to an increase of similar to 23% in the already high share of BF water abstracted by the nearby waterworks. Subsurface travel-time distribution substantially shifted towards shorter travel times. Flow paths with travel times <200 days increased by similar to 10% and those with <300 days by 15%. Generally, the periodic temperature signal, and the summer and winter temperature extrema, increased and penetrated deeper into the aquifer. The joint hydrological and thermal effects caused by the canal reconstruction might increase the potential of biodegradable compounds to further penetrate into the aquifer, also by potentially affecting the redox zonation in the aquifer.
The ability to reflect is considered an essential element of Education for Sustainable Development (ESD) and a key competence for learners and educators in ESD (UNECE Strategy for ESD, 2012). In contrast to its high importance, little is known about how reflective thinking can be identified, influenced or increased in the classroom. Therefore, the objective of this study is to address this need by developing an empirical multi-stage model designed to help educators diagnose different levels of reflective thinking and to identify factors that influence students’ reflective thinking about sustainability. Based on a 4–8-week project with grade 10 and 11 students studying sustainability, reflective thinking performance using weblogs as reflective journals was analysed. In addition, qualitative semi-structured interviews were conducted with the teachers to comprehend the learning environment and the personal value they assigned to ESD in their geography class. To determine the levels of reflective thinking achieved by the students, the study built on the work of Dewey (1933) and pre-existing multi-stage models of reflective thinking (Bain, Ballantyne, & Packer, 1999; Chen, Wei, Wu, & Uden, 2009). Using a qualitative, iterative data analysis, the study adapted the stage models to be applicable in ESD and found great differences in the students’ reflection levels. Furthermore, the study identified eight factors that influence students’ reflective thinking about sustainability. The outcomes of this study may be valuable for educators in high school and higher education, who seek to diagnose their students’ reflective thinking performance and facilitate reflection about sustainability.
Sustained glacier melt in the Himalayas has gradually spawned more than 5,000 glacier lakes that are dammed by potentially unstable moraines. When such dams break, glacier lake outburst floods (GLOFs) can cause catastrophic societal and geomorphic impacts. We present a robust probabilistic estimate of average GLOFs return periods in the Himalayan region, drawing on 5.4 billion simulations. We find that the 100-y outburst flood has an average volume of 33.5(+3.7)/(-3.7) x 10(6) m(3) (posterior mean and 95% highest density interval [HDI]) with a peak discharge of 15,600(+2.000)/(-1,800) m(3).S-1. Our estimated GLOF hazard is tied to the rate of historic lake outbursts and the number of present lakes, which both are highest in the Eastern Himalayas. There, the estimated 100-y GLOF discharge (similar to 14,500 m(3).s(-1)) is more than 3 times that of the adjacent Nyainqentanglha Mountains, and at least an order of magnitude higher than in the Hindu Kush, Karakoram, and Western Himalayas. The GLOF hazard may increase in these regions that currently have large glaciers, but few lakes, if future projected ice loss generates more unstable moraine-dammed lakes than we recognize today. Flood peaks from GLOFs mostly attenuate within Himalayan headwaters, but can rival monsoon-fed discharges in major rivers hundreds to thousands of kilometers downstream. Projections of future hazard from meteorological floods need to account for the extreme runoffs during lake outbursts, given the increasing trends in population, infrastructure, and hydropower projects in Himalayan headwaters.
A digital filter is introduced which treats the problem of predictability versus time averaging in a continuous, seamless manner. This seamless filter (SF) is characterized by a unique smoothing rule that determines the strength of smoothing in dependence on lead time. The rule needs to be specified beforehand, either by expert knowledge or by user demand. As a result, skill curves are obtained that allow a predictability assessment across a whole range of time-scales, from daily to seasonal, in a uniform manner. The SF is applied to downscaled SEAS5 ensemble forecasts for two focus regions in or near the tropical belt, the river basins of the Karun in Iran and the Sao Francisco in Brazil. Both are characterized by strong seasonality and semi-aridity, so that predictability across various time-scales is in high demand. Among other things, it is found that from the start of the water year (autumn), areal precipitation is predictable with good skill for the Karun basin two and a half months ahead; for the Sao Francisco it is only one month, longer-term prediction skill is just above the critical level.
A digital filter is introduced which treats the problem of predictability versus time averaging in a continuous, seamless manner. This seamless filter (SF) is characterized by a unique smoothing rule that determines the strength of smoothing in dependence on lead time. The rule needs to be specified beforehand, either by expert knowledge or by user demand. As a result, skill curves are obtained that allow a predictability assessment across a whole range of time-scales, from daily to seasonal, in a uniform manner. The SF is applied to downscaled SEAS5 ensemble forecasts for two focus regions in or near the tropical belt, the river basins of the Karun in Iran and the Sao Francisco in Brazil. Both are characterized by strong seasonality and semi-aridity, so that predictability across various time-scales is in high demand. Among other things, it is found that from the start of the water year (autumn), areal precipitation is predictable with good skill for the Karun basin two and a half months ahead; for the Sao Francisco it is only one month, longer-term prediction skill is just above the critical level.
During the earliest Triassic microbial mats flourished in the photic zones of marginal seas, generating widespread microbialites. It has been suggested that anoxic conditions in shallow marine environments, linked to the end-Permian mass extinction, limited mat-inhibiting metazoans allowing for this microbialite expansion. The presence of a diverse suite of proxies indicating oxygenated shallow sea-water conditions (metazoan fossils, biomarkers and redox proxies) from microbialite successions have, however, challenged the inference of anoxic conditions. Here, the distribution and faunal composition of Griesbachian microbialites from China, Iran, Turkey, Armenia, Slovenia and Hungary are investigated to determine the factors that allowed microbialite-forming microbial mats to flourish following the end-Permian crisis. The results presented here show that Neotethyan microbial buildups record a unique faunal association due to the presence of keratose sponges, while the Palaeotethyan buildups have a higher proportion of molluscs and the foraminifera Earlandia. The distribution of the faunal components within the microbial fabrics suggests that, except for the keratose sponges and some microconchids, most of the metazoans were transported into the microbial framework via wave currents. The presence of both microbialites and metazoan associations were limited to oxygenated settings, suggesting that a factor other than anoxia resulted in a relaxation of ecological constraints following the mass extinction event. It is inferred that the end-Permian mass extinction event decreased the diversity and abundance of metazoans to the point of significantly reducing competition, allowing photosynthesis-based microbial mats to flourish in shallow water settings and resulting in the formation of widespread microbialites.
Fire and grazing shape biodiversity in savannah landscapes. In land use management, knowing the effects of fire and grazing on biodiversity are important in order to ensure environmental sustainability. Beetles specifically are indicators of the biodiversity response to fire and grazing. A grazing exclusion and burning experiment in a split-plot design was used in order to investigate the interacting effects of fire and wildlife grazing on biomass, diversity, and species composition of darkling beetles (Coleoptera, Tenebrionidae) over time after fire. Darkling beetle species richness and diversity were responding in a three-way-interaction to fire, grazing, and time after fire, whereby biomass of darkling beetles remained unaffected and species compositional changes were attributed to seasonal changes of time only. Fire on ungrazed plots had a negative effect on species diversity and richness 2 weeks and 6 months post fire, whereas fire on grazed plots had no impact on species diversity and richness. Grazing only lowered species diversity and richness 6 months after fire treatments. Results suggest that grazing overrides the effects of fire and that the similar effects caused by fire and grazing are due to niche and assemblage simplification of the habitat.
In crop modeling and yield predictions, the heterogeneity of agricultural landscapes is usually not accounted for. This heterogeneity often arises from landscape elements like forests, hedges, or single trees and shrubs that cast shadows. Shading from forested areas or shrubs has effects on transpiration, temperature, and soil moisture, all of which affect the crop yield in the adjacent arable land. Transitional gradients of solar irradiance can be described as a function of the distance to the zero line (edge), the cardinal direction, and the height of trees. The magnitude of yield reduction in transition zones is highly influenced by solar irradiance-a factor that is not yet implemented in crop growth models on a landscape level. We present a spatially explicit model for shading caused by forested areas, in agricultural landscapes. With increasing distance to forest, solar irradiance and yield increase. Our model predicts that the shading effect from the forested areas occurs up to 15 m from the forest edge, for the simulated wheat yields, and up to 30 m, for simulated maize. Moreover, we estimated the spatial extent of transition zones, to calculate the regional yield reduction caused by shading of the forest edges, which amounted to 5% to 8% in an exemplary region.
Influence of the Main Border Faults on the 3D Hydraulic Field of the Central Upper Rhine Graben
(2019)
The Upper Rhine Graben (URG) is an active rift with a high geothermal potential. Despite being a well-studied area, the three-dimensional interaction of the main controlling factors of the thermal and hydraulic regime is still not fully understood. Therefore, we have used a data-based 3D structural model of the lithological configuration of the central URG for some conceptual numerical experiments of 3D coupled simulations of fluid and heat transport. To assess the influence of the main faults bordering the graben on the hydraulic and the deep thermal field, we carried out a sensitivity analysis on fault width and permeability. Depending on the assigned width and permeability of the main border faults, fluid velocity and temperatures are affected only in the direct proximity of the respective border faults. Hence, the hydraulic characteristics of these major faults do not significantly influence the graben-wide groundwater flow patterns. Instead, the different scenarios tested provide a consistent image of the main characteristics of fluid and heat transport as they have in common: (1) a topography-driven basin-wide fluid flow perpendicular to the rift axis from the graben shoulders to the rift center, (2) a N/NE-directed flow parallel to the rift axis in the center of the rift and, (3) a pronounced upflow of hot fluids along the rift central axis, where the streams from both sides of the rift merge. This upflow axis is predicted to occur predominantly in the center of the URG (northern and southern model area) and shifted towards the eastern boundary fault (central model area).
Above and underground hydrological processes depend on soil moisture (SM) variability, driven by different environmental factors that seldom are well-monitored, leading to a misunderstanding of soil water temporal patterns. This study investigated the stability of the SM temporal dynamics to different monitoring temporal resolutions around the border between two soil types in a tropical watershed. Four locations were instrumented in a small-scale watershed (5.84 km(2)) within the tropical coast of Northeast Brazil, encompassing different soil types (Espodossolo Humiluvico or Carbic Podzol, and Argissolo Vermelho-Amarelo or Haplic Acrisol), land covers (Atlantic Forest, bush vegetation, and grassland) and topographies (flat and moderate slope). The SM was monitored at a temporal resolution of one hour along the 2013-2014 hydrological year and then resampled a resolutions of 6 h, 12 h, 1 day, 2 days, 4 days, 7 days, and 15 days. Descriptive statistics, temporal variability, time-stability ranking, and hierarchical clustering revealed uneven associations among SM time components. The results show that the time-invariant component ruled SM temporal variability over the time-varying parcel, either at high or low temporal resolutions. Time-steps longer than 2 days affected the mean statistical metrics of the SM time-variant parcel. Additionally, SM at downstream and upstream sites behaved differently, suggesting that the temporal mean was regulated by steady soil properties (slope, restrictive layer, and soil texture), whereas their temporal anomalies were driven by climate (rainfall) and hydrogeological (groundwater level) factors. Therefore, it is concluded that around the border between tropical soil types, the distinct behaviour of time-variant and time-invariant components of SM time series reflects different combinations of their soil properties.
Groundwater travel time distributions (TTDs) provide a robust description of the subsurface mixing behavior and hydrological response of a subsurface system. Lagrangian particle tracking is often used to derive the groundwater TTDs. The reliability of this approach is subjected to the uncertainty of external forcings, internal hydraulic properties, and the interplay between them. Here, we evaluate the uncertainty of catchment groundwater TTDs in an agricultural catchment using a 3-D groundwater model with an overall focus on revealing the relationship between external forcing, internal hydraulic properties, and TTD predictions. Eight recharge realizations are sampled from a high-resolution dataset of land surface fluxes and states. Calibration-constrained hydraulic conductivity fields (Ks fields) are stochastically generated using the null-space Monte Carlo (NSMC) method for each recharge realization. The random walk particle tracking (RWPT) method is used to track the pathways of particles and compute travel times. Moreover, an analytical model under the random sampling (RS) assumption is fit against the numerical solutions, serving as a reference for the mixing behavior of the model domain. The StorAge Selection (SAS) function is used to interpret the results in terms of quantifying the systematic preference for discharging young/old water. The simulation results reveal the primary effect of recharge on the predicted mean travel time (MTT). The different realizations of calibration-constrained Ks fields moderately magnify or attenuate the predicted MTTs. The analytical model does not properly replicate the numerical solution, and it underestimates the mean travel time. Simulated SAS functions indicate an overall preference for young water for all realizations. The spatial pattern of recharge controls the shape and breadth of simulated TTDs and SAS functions by changing the spatial distribution of particles' pathways. In conclusion, overlooking the spatial nonuniformity and uncertainty of input (forcing) will result in biased travel time predictions. We also highlight the worth of reliable observations in reducing predictive uncertainty and the good interpretability of SAS functions in terms of understanding catchment transport processes.
Arctic lowlands are characterized by large numbers of small waterbodies, which are known to affect surface energy budgets and the global carbon cycle. Statistical analysis of their size distributions has been hindered by the shortage of observations at sufficiently high spatial resolutions. This situation has now changed with the high-resolution (<5 m) circum-Arctic Permafrost Region Pond and Lake (PeRL) database recently becoming available. We have used this database to make the first consistent, high-resolution estimation of Arctic waterbody size distributions, with surface areas ranging from 0.0001 km(2) (100 m(2)) to 1 km(2). We found that the size distributions varied greatly across the thirty study regions investigated and that there was no single universal size distribution function (including power-law distribution functions) appropriate across all of the study regions. We did, however, find close relationships between the statistical moments (mean, variance, and skewness) of the waterbody size distributions from different study regions. Specifically, we found that the spatial variance increased linearly with mean waterbody size (R-2 = 0.97, p < 2.2e-16) and that the skewness decreased approximately hyperbolically. We have demonstrated that these relationships (1) hold across the 30 Arctic study regions covering a variety of (bio)climatic and permafrost zones, (2) hold over time in two of these study regions for which multi-decadal satellite imagery is available, and (3) can be reproduced by simulating rising water levels in a high-resolution digital elevation model. The consistent spatial and temporal relationships between the statistical moments of the waterbody size distributions underscore the dominance of topographic controls in lowland permafrost areas. These results provide motivation for further analyses of the factors involved in waterbody development and spatial distribution and for investigations into the possibility of using statistical moments to predict future hydrologic dynamics in the Arctic.
The improvement of process representations in hydrological models is often only driven by the modelers' knowledge and data availability. We present a comprehensive comparison between two hydrological models of different complexity that is developed to support (1) the understanding of the differences between model structures and (2) the identification of the observations needed for model assessment and improvement. The comparison is conducted on both space and time and by aggregating the outputs at different spatiotemporal scales. In the present study, mHM, a process‐based hydrological model, and ParFlow‐CLM, an integrated subsurface‐surface hydrological model, are used. The models are applied in a mesoscale catchment in Germany. Both models agree in the simulated river discharge at the outlet and the surface soil moisture dynamics, lending their supports for some model applications (drought monitoring). Different model sensitivities are, however, found when comparing evapotranspiration and soil moisture at different soil depths. The analysis supports the need of observations within the catchment for model assessment, but it indicates that different strategies should be considered for the different variables. Evapotranspiration measurements are needed at daily resolution across several locations, while highly resolved spatially distributed observations with lower temporal frequency are required for soil moisture. Finally, the results show the impact of the shallow groundwater system simulated by ParFlow‐CLM and the need to account for the related soil moisture redistribution. Our comparison strategy can be applied to other models types and environmental conditions to strengthen the dialog between modelers and experimentalists for improving process representations in Earth system models.
The Value of Empirical Data for Estimating the Parameters of a Sociohydrological Flood Risk Model
(2019)
In this paper, empirical data are used to estimate the parameters of a sociohydrological flood risk model. The proposed model, which describes the interactions between floods, settlement density, awareness, preparedness, and flood loss, is based on the literature. Data for the case study of Dresden, Germany, over a period of 200years, are used to estimate the model parameters through Bayesian inference. The credibility bounds of their estimates are small, even though the data are rather uncertain. A sensitivity analysis is performed to examine the value of the different data sources in estimating the model parameters. In general, the estimated parameters are less biased when using data at the end of the modeled period. Data about flood awareness are the most important to correctly estimate the parameters of this model and to correctly model the system dynamics. Using more data for other variables cannot compensate for the absence of awareness data. More generally, the absence of data mostly affects the estimation of the parameters that are directly related to the variable for which data are missing. This paper demonstrates that combining sociohydrological modeling and empirical data gives additional insights into the sociohydrological system, such as quantifying the forgetfulness of the society, which would otherwise not be easily achieved by sociohydrological models without data or by standard statistical analysis of empirical data.
Remote Sensing technologies allow to map biophysical, biochemical, and earth surface parameters of the land surface. Of especial interest for various applications in environmental and urban sciences is the combination of spectral and 3D elevation information. However, those two data streams are provided separately by different instruments, namely airborne laser scanner (ALS) for elevation and a hyperspectral imager (HSI) for high spectral resolution data. The fusion of ALS and HSI data can thus lead to a single data entity consistently featuring rich structural and spectral information. In this study, we present the application of fusing the first pulse return information from ALS data at a sub-decimeter spatial resolution with the lower-spatial resolution hyperspectral information available from the HSI into a hyperspectral point cloud (HSPC). During the processing, a plausible hyperspectral spectrum is assigned to every first-return ALS point. We show that the complementary implementation of spectral and 3D information at the point-cloud scale improves object-based classification and information extraction schemes. This improvements have great potential for numerous land cover mapping and environmental applications.
Private precaution is an important component in contemporary flood risk management and climate adaptation. However, quantitative knowledge about vulnerability reduction via private precautionary measures is scarce and their effects are hardly considered in loss modeling and risk assessments. However, this is a prerequisite to enable temporally dynamic flood damage and risk modeling, and thus the evaluation of risk management and adaptation strategies. To quantify the average reduction in vulnerability of residential buildings via private precaution empirical vulnerability data (n = 948) is used. Households with and without precautionary measures undertaken before the flood event are classified into treatment and nontreatment groups and matched. Postmatching regression is used to quantify the treatment effect. Additionally, we test state-of-the-art flood loss models regarding their capability to capture this difference in vulnerability. The estimated average treatment effect of implementing private precaution is between 11 and 15 thousand EUR per household, confirming the significant effectiveness of private precautionary measures in reducing flood vulnerability. From all tested flood loss models, the expert Bayesian network-based model BN-FLEMOps and the rule-based loss model FLEMOps perform best in capturing the difference in vulnerability due to private precaution. Thus, the use of such loss models is suggested for flood risk assessments to effectively support evaluations and decision making for adaptable flood risk management.
Ecosystem services inherently involve people, whose values help define the benefits of nature's services. It is thus important for researchers to involve stakeholders in ecosystem services research. However, a simple and practicable framework to guide such engagement, and in particular to help researchers anticipate and consider key issues and challenges, has not been well explored. Here, we use experience from the 12 case studies in the European Operational Potential of Ecosystem Research Applications (OPERAs) project to propose a stakeholder engagement framework comprising three key elements: creating space, aligning motivations, and building trust. We argue that involving stakeholders in research demands thoughtful reflection from the researchers about what kind of space they want to create, including if and how they want to bring different interests together, how much space they want to allow for critical discussion, and whether there is a role for particular stakeholders to serve as conduits between others. In addition, understanding their own motivations—including values, knowledge, goals, and desired benefits—will help researchers decide when and how to involve stakeholders, identify areas of common ground and potential disagreement, frame the project appropriately, set expectations, and ensure each party is able to see benefits of engaging with each other. Finally, building relationships with stakeholders can be difficult but considering the roles of existing relationships, time, approach, reputation, and belonging can help build mutual trust. Although the three key elements and the paths between them can play out differently depending on the particular research project, we suggest that a research design that considers how to create the space in which researchers and stakeholders will meet, align motivations between researchers and stakeholders, and build mutual trust will help foster productive researcher–stakeholder relationships.
The link between streamflow extremes and climatology has been widely studied in recent decades. However, a study investigating the effect of large-scale circulation variations on the distribution of seasonal discharge extremes at the European level is missing. Here we fit a climate-informed generalized extreme value (GEV) distribution to about 600 streamflow records in Europe for each of the standard seasons, i.e., to winter, spring, summer and autumn maxima, and compare it with the classical GEV distribution with parameters invariant in time. The study adopts a Bayesian framework and covers the period 1950 to 2016. Five indices with proven influence on the European climate are examined independently as covariates, namely the North Atlantic Oscillation (NAO), the east Atlantic pattern (EA), the east Atlantic-western Russian pattern (EA/WR), the Scandinavia pattern (SCA) and the polar-Eurasian pattern (POL). It is found that for a high percentage of stations the climate-informed model is preferred to the classical model. Particularly for NAO during winter, a strong influence on streamflow extremes is detected for large parts of Europe (preferred to the classical GEV distribution for 46% of the stations). Climate-informed fits are characterized by spatial coherence and form patterns that resemble relations between the climate indices and seasonal precipitation, suggesting a prominent role of the considered circulation modes for flood generation. For certain regions, such as northwestern Scandinavia and the British Isles, yearly variations of the mean seasonal climate indices result in considerably different extreme value distributions and thus in highly different flood estimates for individual years that can also persist for longer time periods.
Continental rift systems form by propagation of isolated rift segments that interact, and eventually evolve into continuous zones of deformation. This process impacts many aspects of rifting including rift morphology at breakup, and eventual ocean-ridge segmentation. Yet, rift segment growth and interaction remain enigmatic. Here we present geological data from the poorly documented Ririba rift (South Ethiopia) that reveals how two major sectors of the East African rift, the Kenyan and Ethiopian rifts, interact. We show that the Ririba rift formed from the southward propagation of the Ethiopian rift during the Pliocene but this propagation was short-lived and aborted close to the Pliocene-Pleistocene boundary. Seismicity data support the abandonment of laterally offset, overlapping tips of the Ethiopian and Kenyan rifts. Integration with new numerical models indicates that rift abandonment resulted from progressive focusing of the tectonic and magmatic activity into an oblique, throughgoing rift zone of near pure extension directly connecting the rift sectors.
Southern Patagonia is a prime example of ongoing oceanic ridge collision and slab-window formation sustained over several million years. The impact of these phenomena on the thermal structure and exhumation of the crust have been mainly assessed with low-temperature thermochronology of bedrock samples. Here, we infer thermal histories from new and existing thermochronological data from the region of most recent ridge collision. In particular, we evaluate the potential far-reaching thermal effects of the evolving slab window, which have previously been considered responsible for patterns of late Miocene reheating associated with back-arc alkaline volcanism. Our model results define protracted cooling since similar to 15 Ma and stepwise exhumation since the late Miocene. The pattern of stepwise exhumation closely matches the onset of Patagonian glaciation at 7 Ma and the successive pulse of glacial incision coeval with neotectonic activity since 3-4 Ma that are also documented by independent geological and geomorphological evidence in the region. Importantly, our findings challenge the recently suggested lack of glacial erosion and incision since 5 Ma in this region. Furthermore, in contrast to previous modelling studies, we find that the available data do not evidence a previously proposed northward-propagating heating event associated with alkaline volcanism. We hypothesize that the anomalous alkaline volcanism in the Patagonian back-arc might be related to trench-orthogonal tears aligned with transform faults in the subducting plate. The substantial differences from the previous modelling procedure on some of the same samples is demonstrated to result from an important lack of convergence in model runs. (C) 2019 Elsevier B.V. All rights reserved.