Innovating the IPCC review process—the potential of young talent

There is significant potential in young talent for enhancing the credibility of the scientific assessments such as the IPCC’s by contributing to quality assurance and quality control. In this essay, we reflect on an experiment that was done by the Dutch government as part of its government review of a contribution to the IPCC’s Fifth Assessment Report (AR5). In an effort to review the entire Working Group II contribution to the AR5 within the official review period for the Second Order Draft (SOD), the PBL Netherlands Environmental Assessment Agency turned to PhD students. This article shows that a systematic review focusing on transparency and errors of a large scientific assessment document using young talented scientists can be successful if certain conditions are met. The reviewers need to have intrinsic motivation to conduct the review. There needs to be a communication plan that fosters engagement and a clear methodology to guide the reviewers through their task. Based on this experiment in review, we reflect on the wider potential for openness and crowdsourcing in scientific assessment processes such as the IPCC’s.     

Weibull-k Revisited: “Tall” Profiles and Height Variation of Wind Statistics

The Weibull distribution is commonly used to describe climatological wind-speed distributions in the atmospheric boundary layer. While vertical profiles of mean wind speed in the atmospheric boundary layer have received significant attention, the variation of the shape of the wind distribution with height is less understood. Previously we derived a probabilistic model based on similarity theory for calculating the effects of stability and planetary boundary-layer depth upon long-term mean wind profiles. However, some applications (e.g. wind energy estimation) require the Weibull shape parameter (k), as well as mean wind speed. Towards the aim of improving predictions of the Weibull- \(k\) profile, we develop expressions for the profile of long-term variance of wind speed, including a method extending our probabilistic wind-profile theory; together these two profiles lead to a profile of Weibull-shape parameter. Further, an alternate model for the vertical profile of Weibull shape parameter is made, improving upon a basis set forth by Wieringa (Boundary-Layer Meteorol, 1989, Vol. 47, 85–110), and connecting with a newly-corrected corollary of the perturbed geostrophic-drag theory of Troen and Petersen (European Wind Atlas, 1989, Risø National Laboratory, Roskilde). Comparing the models for Weibull-k profiles, a new interpretation and explanation is given for the vertical variation of the shape of wind-speed distributions. Results of the modelling are shown for a number of sites, with a discussion of the models’ efficacy and applicability. The latter includes a comparative evaluation of Wieringa-type empirical models and perturbed-geostrophic forms with regard to surface-layer behaviour, as well as for heights where climatological wind-speed variability is not dominated by surface effects.     

Early Last Interglacial Greenland Ice Sheet melting and a sustained period of meridional overturning weakening: a model analysis of the uncertainties

Proxy-data suggest that the Last Interglacial (LIG; ~130–116 ka BP) climate was characterized by higher temperatures, a partially melted Greenland Ice Sheet (GIS) and a changed Atlantic meridional overturning circulation (AMOC). Notwithstanding the uncertainties in LIG palaeoclimatic reconstructions, this setting potentially provides an opportunity to evaluate the relation between GIS melt and the AMOC as simulated by climate models. However, first we need to assess the extent to which a causal relation between early LIG GIS melt and the weakened AMOC is plausible. With a series of transient LIG climate simulations with the LOVECLIM earth system model, we quantify the importance of the major known uncertainties involved in early LIG GIS melt scenarios. Based on this we construct a specific scenario that is within the parameter space of uncertainties and show that it is physically consistent that early LIG GIS melting kept the AMOC weakened. Notwithstanding, this scenario is at the extreme end of the parameter space. Assuming that proxy-based reconstructions of early LIG AMOC weakening offer a realistic representation of its past state, this indicates that either (1) the AMOC weakening was caused by other forcings than early LIG GIS melt or (2) the early LIG AMOC was less stable than indicated by our simulations and a small amount of GIS melt was sufficient to keep the AMOC in the weak state of a bi-stable regime. We argue that more intensive research is required because of the high potential of the early LIG to evaluate model performance in relation to the AMOC response to GIS melt.     

Comment on “Trends and low frequency variability of extra-tropical cyclone activity in the ensemble of twentieth century reanalysis” by Xiaolan L. Wang,Y. Feng,G. P. Compo,V. R. Swail,F. W. Zwiers,R. J. Allan,and P. D. Sardeshmukh,Climate Dynamics, 2012

The main subject of this article is to comment on the issue of storminess trends derived from the twentieth century reanalysis (20CR) and from observations in the North Atlantic region written about in Wang et al. (Clim Dyn 40(11–12):2775–2800, 2012). The statement that the 20CR estimates would be consistent with storminess derived from pressure-based proxies does not hold for the time prior to 1950.     

Shear Band Formation in Numerical Simulations Applying a Continuum Damage Rheology Model

In seismically active regions, faults nucleate, propagate, and form networks that evolve over time. To simulate crustal faulting processes, including the evolution of fault-zone properties, a rheological model must incorporate concepts such as damage rheology that describe the various stages of the seismic cycle (strain localization, subcritical crack growth and macroscopic failure) while accounting for material degradation and healing and off-fault deformation. Here we study the fundamental patterns of strain-localisation within the framework of a continuum damage rheology by performing a shear band analysis (linear instability analysis) and comparing predictions of shear band orientations with numerical results of the nonlinear problem. We find (analytically and numerically) that the angle between the shear band and the less compressive (transverse) direction is between 47° in compression tests with a strain ratio of 0.25 (highly confined compression test), and 60° for a strain ratio higher than 1.4 (axial compression and transverse extension). In addition we find that shear bands exhibit local dilation (I ₁ > 0) in a wide range of strain ratios excluding only simulations with highly confined compression (which yield compacting shear bands or non-localized deformation). Finally, we discuss the applicability of the damage model for simulating deformation in the seismogenic, brittle crust.     

Construction of Intensity–Duration–Frequency (IDF) curves for precipitation at Lubumbashi,Congo, under the hypothesis of inadequate data

Abstract

Intensity–Duration–Frequency (IDF) curves for precipitation constitute a probabilistic tool and have proven useful in water resources management. In particular, IDF curves for precipitation enable questions on the extreme character of precipitation to be answered. The construction of IDF curves for precipitation is difficult or impossible in tropical areas due to the lack of long-term extreme precipitation data. A technique is proposed to overcome this shortcoming by combining limited high-frequency information on rainfall extremes with long-term daily rainfall information. It may be regarded as an extension of Koutsoyiannis' approach. Using this technique, IDF curves for precipitation are produced for Lubumbashi in Congo.

Citation Van de Vyver, H. & Demarée, G. R. (2010) Construction of Intensity–Duration–Frequency (IDF) curves for precipitation at Lubumbashi, Congo, under the hypothesis of inadequate data. Hydrol. Sci. J. 55(4), 555–564.     

Ellipsoidal quadtrees for indexing of global geographical data

Systems for landscape visualization and geographical data handling require methods for efficient data access. Retrieval of data from large geographical databases, ten to thousands of Gbytes, is usually optimized with spatial indexing mechanisms. The simplest form of spatial indexing is achieved by dividing the database into congruent grid cells. The subsequent subdivision of the grid cells can be based on so-called quadtrees. Quadtrees for two-dimensional division and subdivision are appropriate for cartographical data. A geographical database, with objects stored in geocentric or geodetic (geographical) co-ordinates, requires indexing mechanisms that take into account the shape of the Earth. In this paper, we present a method for indexing of geographical data, named Ellipsoidal Quadtrees (EQT). In contrast to other global indexing methods, EQT is based on the Earth ellipsoid and not a spherical approximation. EQT division and subdivision make it possible to divide the Earth surface into a mesh of quadrangles with equal areas. We will demonstrate that EQT is flexible. It can be used for indexing databases of various sizes, including national and global databases. Tests on real data show that the performance of EQT is good.     

Comparing the adaptive Gaussian mixture filter with the ensemble Kalman filter on synthetic reservoir models

Over the last years, the ensemble Kalman filter (EnKF) has become a very popular tool for history matching petroleum reservoirs. EnKF is an alternative to more traditional history matching techniques as it is computationally fast and easy to implement. Instead of seeking one best model estimate, EnKF is a Monte Carlo method that represents the solution with an ensemble of state vectors. Lately, several ensemble-based methods have been proposed to improve upon the solution produced by EnKF. In this paper, we compare EnKF with one of the most recently proposed methods, the adaptive Gaussian mixture filter (AGM), on a 2D synthetic reservoir and the Punq-S3 test case. AGM was introduced to loosen up the requirement of a Gaussian prior distribution as implicitly formulated in EnKF. By combining ideas from particle filters with EnKF, AGM extends the low-rank kernel particle Kalman filter. The simulation study shows that while both methods match the historical data well, AGM is better at preserving the geostatistics of the prior distribution. Further, AGM also produces estimated fields that have a higher empirical correlation with the reference field than the corresponding fields obtained with EnKF.     

A parallel block preconditioner for large-scale poroelasticity with highly heterogeneous material parameters

Large-scale simulations of coupled flow in deformable porous media require iterative methods for solving the systems of linear algebraic equations. Construction of efficient iterative methods is particularly challenging in problems with large jumps in material properties, which is often the case in realistic geological applications, such as basin evolution at regional scales. The success of iterative methods for such problems depends strongly on finding effective preconditioners with good parallel scaling properties, which is the topic of the present paper. We present a parallel preconditioner for Biot’s equations of coupled elasticity and fluid flow in porous media. The preconditioner is based on an approximation of the exact inverse of the two-by-two block system arising from a finite element discretisation. The approximation relies on a highly scalable approximation of the global Schur complement of the coefficient matrix, combined with generally available state-of-the-art multilevel preconditioners for the individual blocks. This preconditioner is shown to be robust on problems with highly heterogeneous material parameters. We investigate the weak and strong parallel scaling of this preconditioner on up to 512 processors and demonstrate its ability on a realistic basin-scale problem in poroelasticity with over eight million tetrahedral elements.     

Eemian landscape response to climatic shifts and evidence for northerly Neanderthal occupation at a palaeolake margin in northern Germany

The prevailing view suggests that the Eemian interglacial on the European Plain was characterized by largely negligible geomorphic activity beyond the coastal areas. However, systematic geomorphological studies are sparse. Here we present a detailed reconstruction of Eemian to Early Weichselian landscape evolution in the vicinity of a small fingerlake on the northern margin of the Salzwedel Palaeolake in Lower Saxony (Germany). We apply a combination of seismics, sediment coring, pollen analysis and luminescence dating on a complex sequence of colluvial, paludal and lacustrine sediments. Results suggest two pronounced phases of geomorphic activity, directly before the onset and at the end of the Eemian period, with an intermediate period of pronounced landscape stability. The dynamic phases were largely driven by incomplete vegetation cover, but likely accentuated by fluvial incision in the neighbouring Elbe Valley. Furthermore, we discovered Neanderthal occupation at the lakeshore during Eemian pollen zone (PZ) E IV, which is chronologically in line with other known Eemian sites of central Europe. Our highly-resolved spatio-temporal data substantially contribute to the understanding of climate-induced geomorphic processes throughout and directly after the last interglacial period. It helps unraveling the landscape dynamics between the coastal areas to the north and the loess belt to the south.