An approach for the determination of precipitation input for worst-case flood modelling

ABSTRACT

There is a lack of suitable methods for creating precipitation scenarios that can be used to realistically estimate peak discharges with very low probabilities. On the one hand, existing methods are methodically questionable when it comes to physical system boundaries. On the other hand, the spatio-temporal representativeness of precipitation patterns as system input is limited. In response, this paper proposes a method of deriving spatio-temporal precipitation patterns and presents a step towards making methodically correct estimations of infrequent floods by using a worst-case approach. A Monte Carlo approach allows for the generation of a wide range of different spatio-temporal distributions of an extreme precipitation event that can be tested with a rainfall–runoff model that generates a hydrograph for each of these distributions. Out of these numerous hydrographs and their corresponding peak discharges, the physically plausible spatio-temporal distributions that lead to the highest peak discharges are identified and can eventually be used for further investigations.

Editor A. Castellarin; Associate editor E. Volpi     

Geochemistry and quality of the groundwater from the karstic and coastal aquifer of Geropotamos River Basin at north-central Crete, Greece

In Geropotamos River Basin, located on the north-central part of Crete, Greece, two main factors were believed to be affecting the geochemistry of the groundwater with high salt contents: seawater intrusion and/or Miocene evaporates. To identify the origin of the high salinity in groundwater, a hydrogeochemical and isotopic study was performed. Water samples from 22 wells and 2 springs were analyzed for physico-chemical parameters, major ions analysis, as well as stable isotopes (??¹⁸O, ??D). From the present survey, in which detailed hydrogeochemical investigation was conducted, the uncertainty of the contamination sources was decreased in the northern part of Geropotamos Basin. The results complement the scenario in which seawater and the widespread human activities are the principal sources of groundwater contamination. Moreover, the results of the stable isotopes analyses (??¹⁸O and ??D) support the same hypothesis and make seawater intrusion the most probable cause for the highest salinity waters. It is indicated that saline intrusion is likely to occur along fractures in a fault zone through otherwise low-permeability phyllite?Cquartzite bedrock, which demonstrates the critical role of fracture pathways in salination problems of coastal aquifers.     

Evolution of anthropogenic and biomass burning emissions of air pollutants at global and regional scales during the 1980�C2010 period

Several different inventories of global and regional anthropogenic and biomass burning emissions are assessed for the 1980?C2010 period. The species considered in this study are carbon monoxide, nitrogen oxides, sulfur dioxide and black carbon. The inventories considered include the ACCMIP historical emissions developed in support of the simulations for the IPCC AR5 assessment. Emissions for 2005 and 2010 from the Representative Concentration Pathways (RCPs) are also included. Large discrepancies between the global and regional emissions are identified, which shows that there is still no consensus on the best estimates for surface emissions of atmospheric compounds. At the global scale, anthropogenic emissions of CO, NO_x and SO₂ show the best agreement for most years, although agreement does not necessarily mean that uncertainty is low. The agreement is low for BC emissions, particularly in the period prior to 2000. The best consensus is for NO_x emissions for all periods and all regions, except for China, where emissions in 1980 and 1990 need to be better defined. Emissions of CO need better quantification in the USA and India for all periods; in Central Europe, the evolution of emissions during the past two decades needs to be better determined. The agreement between the different SO₂ emissions datasets is rather good for the USA, but better quantification is needed elsewhere, particularly for Central Europe, India and China. The comparisons performed in this study show that the use of RCP8.5 for the extension of the ACCMIP inventory beyond 2000 is reasonable, until more global or regional estimates become available. Concerning biomass burning emissions, most inventories agree within 50?C80%, depending on the year and season. The large differences between biomass burning inventories are due to differences in the estimates of burned areas from the different available products, as well as in the amount of biomass burned.     

High-temperature emplacement of the Cerro Galán and Toconquis Group ignimbrites (Puna plateau,NW Argentina) determined by TRM analyses

Estimates of pyroclastic flow emplacement temperatures in the Cerro Galán ignimbrite and Toconquis Group ignimbrites were determined using thermal remanent magnetization of lithic clasts embedded within the deposits. These ignimbrites belong to the Cerro Galán volcanic system, one of the largest calderas in the world, in the Puna plateau, NW Argentina. Temperature estimates for the 2.08-Ma Cerro Galán ignimbrite are retrieved from 40 sites in 14 localities (176 measured clasts), distributed at different distances from the caldera and different stratigraphic heights. Additionally, temperature estimates were obtained from 27 sample sites (125 measured clasts) from seven ignimbrite units forming the older Toconquis Group (5.60–4.51 Ma), mainly outcropping along a type section at Rio Las Pitas, Vega Real Grande. The paleomagnetic data obtained by progressive thermal demagnetization show that the clasts of the Cerro Galán ignimbrite have one single magnetic component, oriented close to the expected geomagnetic field at the time of emplacement. Results show therefore that most of the clasts acquired a new magnetization oriented parallel to the magnetic field at the moment of the ignimbrite deposition, suggesting that the clasts were heated up to or above the highest blocking temperature (T _b) of the magnetic minerals (T _b = 580°C for magnetite; T _b = 600–630°C for hematite). We obtained similar emplacement temperature estimations for six out of the seven volcanic units belonging to the Toconquis Group, with the exception of one unit (Lower Merihuaca), where we found two distinct magnetic components. The estimation of emplacement temperatures in this latter case is constrained at 580–610°C, which are lower than the other ignimbrites. These estimations are also in agreement with the lowest pre-eruptive magma temperatures calculated for the same unit (i.e., 790°C; hornblende–plagioclase thermometer; Folkes et al. 2011b). We conclude that the Cerro Galán ignimbrite and Toconquis Group ignimbrites were emplaced at temperatures equal to or higher than 620°C, except for Lower Merihuaca unit emplaced at lower temperatures. The homogeneity of high temperatures from proximal to distal facies in the Cerro Galán ignimbrite provides constraints for the emplacement model, marked by a relatively low eruption column, low levels of turbulence, air entrainment, surface–water interaction, and a high level of topographic confinement, all ensuring minimal heat loss.     

Global mass wasting during the Middle Ordovician: Meteoritic trigger or plate-tectonic environment?

Mass wasting at continental margins on a global scale during the Middle Ordovician has recently been related to high meteorite influx. Although a high meteorite influx during the Ordovician should not be neglected, we challenge the idea that mass wasting was mainly produced by meteorite impacts over a period of almost 10 Ma. Having strong arguments against the impact-related hypothesis, we propose an alternative explanation, which is based on a re-evaluation of the mass wasting sites, considering their plate-tectonic distribution and the global sea level curve. A striking and important feature is the distribution of most of the mass wasting sites along continental margins characterised by periods of magmatism, terrane accretion and continental or back-arc rifting, respectively, related to subduction of oceanic lithosphere. Such processes are commonly connected with seismic activity causing earthquakes, which can cause downslope movement of sediment and rock. Considering all that, it seems more likely that most of this mass wasting was triggered by earthquakes related to plate-tectonic processes, which caused destabilisation of continental margins resulting in megabreccias and debris flows. Moreover, the period of mass wasting coincides with sea level drops during global sea level lowstand. In some cases, sea level drops can release pore-water overpressure reducing sediment strength and hence promoting instability of sediment at continental margins. Reduced pore-water overpressure can also destabilise gas hydrate-bearing sediment, causing slope failure, and thus resulting in submarine mass wasting. Overall, the global mass wasting during the Middle Ordovician does not need meteoritic trigger.     

Paleoecology and radiocarbon dating of the Pleistocene megafauna of the Brazilian Intertropical Region

During the Pleistocene a fauna composed of large (biomass > 44 kg) and giant mammals (biomass > 1000 kg) that are usually associated with open environments lived in the Brazilian Intertropical Region. We present here new information concerning the paleoecology and chronology of some species of this megafauna. Carbon isotope analyses were performed for a better understanding of the paleoecology of the species Eremotherium laurillardi (Lund, 1842), Notiomastodon platensis (Ameghino, 1888) and Toxodon platensis (Owen, 1849). The δ¹³C data allow attributing a generalist diet to these species, which varied according to the kind of habitat in which they lived. In more open habitats all species were grazers; in mixed habitats E. laurillardi and T. platensis were mixed feeders, and N. platensis was grazer; and in more closed habitats all species were mixed feeders.     

Impact of single inclined faults on the fluid flow and heat transport: results from 3-D finite element simulations

The impact of inclined faults on the hydrothermal field is assessed by adding simplified structural settings to synthetic models. This study is innovative in carrying out numerical simulations because it integrates the real 3-D nature of flow influenced by a fault in a porous medium, thereby providing a useful tool for complex geothermal modelling. The 3-D simulations for the coupled fluid flow and heat transport processes are based on the finite element method. In the model, one geological layer is dissected by a dipping fault. Sensitivity analyses are conducted to quantify the effects of the fault’s transmissivity on the fluid flow and thermal field. Different fault models are compared with a model where no fault is present to evaluate the effect of varying fault transmissivity. The results show that faults have a significant impact on the hydrothermal field. Varying either the fault zone width or the fault permeability will result in relevant differences in the pressure, velocity and temperature field. A linear relationship between fault zone width and fluid velocity is found, indicating that velocities increase with decreasing widths. The faults act as preferential pathways for advective heat transport in case of highly transmissive faults, whereas almost no fluid may be transported through poorly transmissive faults.     

Impact of Poroelastic Response of Sandstones on Geothermal Power Production

During geothermal power production using a borehole doublet consisting of a production and injection well, the reservoir conditions such as permeability k, porosity φ and Skempton coefficient B at the geothermal research site Gross Schoenebeck/Germany will change. Besides a temperature decrease at the injection well and a change of the chemical equilibrium, also the pore pressure p _p will vary in a range of approximately 44 MPa ± 10 MPa in our reservoir at ?3850 to ?4258 m depth. This leads to a poroelastic response of the reservoir rocks depending on effective pressure p _eff (difference between mean stress and pore pressure), resulting in a change in permeability k, porosity φ and the poroelastic parameter Skempton coefficient B. Hence, we investigated the effective pressure dependency of Flechtinger sandstone, an outcropping equivalent of the reservoir rock via laboratory experiments. The permeability decreased by 21% at an effective pressure range from 3 to 30 MPa, the porosity decreased by 11% (p _eff = 6 to 65 MPa) and the Skempton coefficient decreased by 24% (p _eff = 4 to 25 MPa). We will show which mechanisms lead to the change of the mentioned hydraulic and poroelastic parameters and the influence of these changes on the productivity of the reservoir. The most significant changes occur at low effective pressures until 15 to 20 MPa. For our in situ reservoir conditions p _eff = 43 MPa a change of 10 MPa effective pressure will result in a change in matrix permeability of less than 4% and in matrix porosity of less than 2%. Besides natural fracture systems, fault zones and induced hydraulic fractures, the rock matrix its only one part of geothermal systems. All components can be influenced by pressure, temperature and chemical reactions. Therefore, the determined small poroelastic response of rock matrix does not significantly influence the sustainability of the geothermal reservoir.     

The Holocene Secche di Lazzaro phreatomagmatic succession (Stromboli,Italy): evidence of pyroclastic density current origin deduced by facies analysis and AMS flow directions

The edifice of Stromboli volcano gravitationally collapsed several times during its volcanic history (>100 ka–present). The largest Holocene event occurred during the final stage of the Neostromboli activity (∼13–5 ka), and was accompanied by the emplacement of phreatomagmatic and lahar deposits, known as the Secche di Lazzaro succession. A stratigraphic and paleomagnetic study of the Secche di Lazzaro deposits allows the interpretation of the emplacement and the eruptive processes. We identify three main units within the succession that correspond to changing eruption conditions. The lower unit (UA) consists of accretionary lapilli-rich, thinly bedded, parallel- to cross-stratified ash deposits, interpreted to indicate the early stages of the eruption and emplacement of dilute pyroclastic density currents. Upward, the second unit (UB) of the deposit is more massive and the beds thicker, indicating an increase in the sedimentation rate from pyroclastic density currents. The upper unit (UC) caps the succession with thick, immediately post-eruptive lahars, which reworked ash deposited on the volcano’s slope. Flow directions obtained by Anisotropy of Magnetic Susceptibility (AMS) analysis of the basal bed of UA at the type locality suggest a provenance of pyroclastic currents from the sea. This is interpreted to be related to the initial base-surges associated with water–magma interaction that occurred immediately after the lateral collapse, which wrapped around the shoulder of the sector collapse scar. Upward in the stratigraphy (upper beds of UA and UB) paleoflow directions change and show a provenance from the summit vent, probably related to the multiple collapses of a vertical, pulsatory eruptive column.