Advances in pan‐European flood hazard mapping

Flood hazard maps at trans‐national scale have potential for a large number of applications ranging from climate change studies, reinsurance products, aid to emergency operations for major flood crisis, among others. However, at continental scales, only few products are available, due to the difficulty of retrieving large consistent data sets. Moreover, these are produced at relatively coarse grid resolution, which limits their applications to qualitative assessments. At finer resolution, maps are often limited to country boundaries, due to limited data sharing at trans‐national level. The creation of a European flood hazard map would currently imply a collection of scattered regional maps, often lacking mutual consistency due to the variety of adopted approaches and quality of the underlying input data. In this work, we derive a pan‐European flood hazard map at 100 m resolution. The proposed approach is based on expanding a literature cascade model through a physically based approach. A combination of distributed hydrological and hydraulic models was set up for the European domain. Then, an observed meteorological data set is used to derive a long‐term streamflow simulation and subsequently coherent design flood hydrographs for a return period of 100 years along the pan‐European river network. Flood hydrographs are used to simulate areas at risk of flooding and output maps are merged into a pan‐European flood hazard map. The quality of this map is evaluated for selected areas in Germany and United Kingdom against national/regional hazard maps. Despite inherent limitations and model resolution issues, simulated maps are in good agreement with reference maps (hit rate between 59% and 78%, critical success index between 43% and 65%), suggesting strong potential for a number of applications at the European scale. Copyright © 2013 John Wiley & Sons, Ltd.     

Distributions-oriented wind forecast verification by a hidden Markov model for multivariate circular–linear data

Winds from the North–West quadrant and lack of precipitation are known to lead to an increase of PM10 concentrations over a residential neighborhood in the city of Taranto (Italy). In 2012 the local government prescribed a reduction of industrial emissions by 10% every time such meteorological conditions are forecasted 72 h in advance. Wind forecasting is addressed using the Weather Research and Forecasting (WRF) atmospheric simulation system by the Regional Environmental Protection Agency. In the context of distributions-oriented forecast verification, we propose a comprehensive model-based inferential approach to investigate the ability of the WRF system to forecast the local wind speed and direction allowing different performances for unknown weather regimes. Ground-observed and WRF-forecasted wind speed and direction at a relevant location are jointly modeled as a 4-dimensional time series with an unknown finite number of states characterized by homogeneous distributional behavior. The proposed model relies on a mixture of joint projected and skew normal distributions with time-dependent states, where the temporal evolution of the state membership follows a first order Markov process. Parameter estimates, including the number of states, are obtained by a Bayesian MCMC-based method. Results provide useful insights on the performance of WRF forecasts in relation to different combinations of wind speed and direction.     

Early Jurassic palaeoenvironments in the Surat Basin,Australia – marine incursion into eastern Gondwana

Interpretations of palaeodepositional environments are important for reconstructing Earth history. Only a few maps showing the Jurassic depositional environments in eastern Australia currently exist. Consequently, a detailed understanding of the setting of Australia in Gondwana is lacking. Core, wireline logs, two-dimensional and three-dimensional seismic from the Precipice Sandstone and Evergreen Formation in the Surat Basin have been used to construct maps showing the evolution of depositional environments through the Early Jurassic. The results indicate the succession consists of three third-order sequences (Sequence 1 to Sequence 3) that were controlled by eustatic sea level. The lowstand systems tract in Sequence 1 comprises braidplain deposits, confined to a fairway that parallels the basin centre. The strata were initially deposited in two sub-basins, with rivers flowing in different orientations in each sub-basin. The transgressive systems tract of Sequence 1 to lowstand systems tract of Sequence 3 is dominated by fluvio–deltaic systems infilling a single merged basin centre. Finally, the transgressive and highstand systems tracts of Sequence 3 show nearshore environments depositing sediment into a shallow marine basin. In the youngest part of this interval, ironstone shoals are the most conspicuous facies, the thickness and number of which increase towards the north and east. This study interprets a corridor to the open ocean through the Clarence–Moreton Basin, or the Carpentaria and Papuan basins, evidence of which has been eroded. These results challenge a commonly held view that eastern Australia was not influenced by eustasy, and propose a more dynamic palaeogeographic setting comprising a mixture of fluvial, deltaic and shallow marine sedimentary environments. This work can be used to unravel the stratigraphic relationships between Mesozoic eastern Australian basins, or in other basins globally as an analogue for understanding the complex interplay of paralic depositional systems in data poor areas.     

DISOLV: A Python Package for the Interpretation of Borehole Dilution Tests

Single borehole dilution tests (SBDTs) are an inexpensive but effective technique for hydrogeological characterization of hard-rock aquifers. We present a freely available, easy-to-use, open-source Python package, DISOLV, for plotting, analyzing, and modeling SBDT data. DISOLV can significantly reduce the time spent interpreting field data by helping to identify flowing fractures intersecting the borehole and estimate the corresponding flow rates. DISOLV is successfully benchmarked against two analytical solutions. We also present an example application to real data collected in a borehole in a crystalline basement aquifer in southern India.     

Sharp spatially constrained inversion with applications to transient electromagnetic data

Time‐domain electromagnetic data are conveniently inverted by using smoothly varying 1D models with fixed vertical discretization. The vertical smoothness of the obtained models stems from the application of Occam‐type regularization constraints, which are meant to address the ill‐posedness of the problem. An important side effect of such regularization, however, is that horizontal layer boundaries can no longer be accurately reproduced as the model is required to be smooth. This issue can be overcome by inverting for fewer layers with variable thicknesses; nevertheless, to decide on a particular and constant number of layers for the parameterization of a large survey inversion can be equally problematic. Here, we present a focusing regularization technique to obtain the best of both methodologies. The new focusing approach allows for accurate reconstruction of resistivity distributions using a fixed vertical discretization while preserving the capability to reproduce horizontal boundaries. The formulation is flexible and can be coupled with traditional lateral/spatial smoothness constraints in order to resolve interfaces in stratified soils with no additional hypothesis about the number of layers. The method relies on minimizing the number of layers of non‐vanishing resistivity gradient, instead of minimizing the norm of the model variation itself. This approach ensures that the results are consistent with the measured data while favouring, at the same time, the retrieval of horizontal abrupt changes. In addition, the focusing regularization can also be applied in the horizontal direction in order to promote the reconstruction of lateral boundaries such as faults. We present the theoretical framework of our regularization methodology and illustrate its capabilities by means of both synthetic and field data sets. We further demonstrate how the concept has been integrated in our existing spatially constrained inversion formalism and show its application to large‐scale time‐domain electromagnetic data inversions.     

The Seismotectonics of the Po Plain (Northern Italy): Tectonic Diversity in a Blind Faulting Domain

We present a systematic and updated overview of a seismotectonic model for the Po Plain (northern Italy). This flat and apparently quiet tectonic domain is, in fact, rather active as it comprises the shortened foreland and foredeep of both the Southern Alps and the Northern Apennines. Assessing its seismic hazard is crucial due to the concentration of population, industrial activities, and critical infrastructures, but it is also complicated because (a) the region is geologically very diverse, and (b) nearly all potential seismogenic faults are buried beneath a thick blanket of Pliocene–Pleistocene sediments, and thus can be investigated only indirectly. Identifying and parameterizing the potential seismogenic faults of the Po Plain requires proper consideration of their depth, geometry, kinematics, earthquake potential and location with respect to the two confronting orogens. To this end, we subdivided them into four main, homogeneous groups. Over the past 15 years we developed new strategies for coping with this diversity, resorting to different data and modeling approaches as required by each individual fault group. The most significant faults occur beneath the thrust fronts of the Ferrara-Romagna and Emilia arcs, which correspond to the most advanced and buried portions of the Northern Apennines and were the locus of the destructive May 2012 earthquake sequence. The largest known Po Plain earthquake, however, occurred on an elusive reactivated fault cutting the Alpine foreland south of Verona. Significant earthquakes are expected to be generated also by a set of transverse structures segmenting the thrust system, and by the deeper ramps of the Apennines thrusts. The new dataset is intended to be included in the next version of the Database of Individual Seismogenic Sources (DISS; http://diss.rm.ingv.it/diss/, version 3.2.0, developed and maintained by INGV) to improve completeness of potential sources for seismic hazard assessment.     

Towards a predictive model to assess the natural position of the Posidonia oceanica seagrass meadows upper limit

The upper portion of the meadows of the protected Mediterranean seagrass Posidonia oceanica occurs in the region of the seafloor mostly affected by surf-related effects. Evaluation of its status is part of monitoring programs, but proper conclusions are difficult to draw due to the lack of definite reference conditions. Comparing the position of the meadow upper limit with the beach morphodynamics (i.e. the distinctive type of beach produced by topography and wave climate) provided evidence that the natural landwards extension of meadows can be predicted. An innovative model was therefore developed in order to locate the region of the seafloor where the meadow upper limit should lie in natural conditions (i.e. those governed only by hydrodynamics, in absence of significant anthropogenic impact). This predictive model was validated in additional sites, which showed perfect agreement between predictions and observations. This makes the model a valuable tool for coastal management.     

Crustal thickness estimation beneath the southern central Andes at 30°S and 36°S from S wave receiver function analysis

We have used the S wave receiver function (SRF) technique to investigate the crustal thickness beneath two seismic profiles from the CHARGE project in the southern central Andes. A previous study employing the P wave receiver function method has observed the Moho interface beneath much of the profiles. They found, however, that the amplitude of the P to S conversion was diminished in the western part of the profiles and have attributed it to a reduction of the impedance contrast at the Moho due to lower crustal ecologitization. With SRF, we have successfully detected S to P converted waves from the Moho as well as possible conversions from other lithospheric boundaries. The continental South American crust reaches its maximum thickness of ∼70 km (along 30°S between 70°W and 68.5°W) beneath the Principal Cordillera and the Famatina system and becomes thinner towards the Sierras Pampeanas with a thickness of ∼40 km. Negative phases, possibly related to the base of the continental and oceanic lithosphere, can be recognized in the summation traces at different depths. By comparing our results with data obtained from previous investigations, we are able to further constrain the thickness of the crust and lithosphere beneath the central Andes.