A simulation experiment for optimal design hyetograph selection

The aim of this work is to assess the accuracy of literature design hyetographs for the evaluation of peak discharges during flood events. Five design hyetographs are examined in a set of simulations, based upon the following steps: (i) an ideal river basin is defined, characterized by a Beta distribution shaped unit hydrograph (UH); (ii) 1000 years of synthetic rainfall are artificially generated; (iii) a discharge time‐series is obtained from the convolution of the rainfall time‐series and the UH, and the reference T‐years flood is computed from this series; (iv) for the same return period T, the parameters of the intensity–duration–frequency (IDF) curve are estimated from the 1000 years of synthetic rainfall; (v) five design hyetographs are determined from the IDF curves and are convolved with the discrete UH to find the corresponding design hydrographs; (vi) the hydrograph peaks are compared with the reference T‐years flood and the advantages and drawbacks of each of the five approaches are evaluated. The rainfall and UH parameters are varied, and the whole procedure is repeated to assess the sensitivity of results to the system configuration. We found that all design hyetographs produce flood peak estimates that are consistently biased in most of the climatic and hydrologic conditions considered. In particular, significant underestimation of the design flood results from the adoption of any rectangular hyetograph used in the context of the rational formula. In contrast, the Chicago hyetograph tends to overestimate peak flows. In two cases it is sufficient to multiply the result by a constant scaling factor to obtain robust and nearly unbiased estimates of the design floods. Copyright © 2007 John Wiley & Sons, Ltd.     

The Breccia Museo formation, Campi Flegrei, southern Italy: geochronology, chemostratigraphy and relationship with the Campanian Ignimbrite eruption

The Breccia Museo is one of the most debated volcanic formations of the Campi Flegrei volcanic district. The deposit, made up of six distinctive stratigraphic units, has been interpreted by some as the proximal facies of the major caldera-forming Campanian Ignimbrite eruption, and by others as the product of several, more recent, independent and localized events. New geochemical and chemostratigraphical data and Ar–Ar age determinations for several units of the Breccia Museo deposits (~39 ka), correlate well with the Campanian Ignimbrite-forming eruption. The chemical zoning of the Breccia Museo deposits is interpreted here to be a consequence of a three-stage event that tapped a vertically zoned trachytic magma chamber. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

The 2D Steady Hydraulic Head Field Surrounding a Pumping Well in a Finite Heterogeneous Confined Aquifer

We present a second-order analytic solution [in terms of a heterogeneous log-transmissivity Y(r) = ln T(r)] for the hydraulic head field in a finite 2D confined heterogeneous aquifer under steady radial flow conditions assuming fixed head boundary conditions at the well and at a circular exterior boundary. The solution may be used to obtain the gradient used in calculation of solute transport to a well in a heterogeneous transmissivity field. The solution, obtained using perturbation methods coupled with Green's function techniques, leads us to postulate a more general form of the head for arbitrarily large-variance fields and may be used to obtain moment relations between the log-transmissivity and head under convergent flow conditions when Y(r) is expressed as a random space function. We present expressions for the mean head field when the log-transmissivity is Gaussian and conditioned on the transmissivity value at the well for an arbitrary ln T covariance. Finally, we look at the effect of parameter variations on the mean head behavior and present numerical simulations verifying the second-order mean head expressions.     

Petrology of an ophiolitic cumulate sequence from Pindos,Greece

The ophiolitic sequence which crops out along the Aspropotamos Valley, Northern Pindos, Greece is composed from the bottom to the top of cumulates, dolerites, basaltic lavas, upper pillow lavas with basaltic/andesitic composition, and scarce basaltic dykes. The intrusive sequence, which is the subject of the present paper, exhibits magmatic layering more pronounced at the bottom than at the top where isotropic gabbros occur; they grade into the overlying dolerites. Troctolites with rare ultramafites prevail in the lower section and olivine gabbros in the upper section; at the top two-pyroxene gabbros appear. The rocks are mainly adcumulates and mesocumulates with subordinate heteradcumulates. The cumulus phases separated in the order: olivine and Cr-spinel, plagioclase, clinopyroxene, orthopyroxene. Olivine, plagioclase and pyroxenes frequently exhibit adeumulus overgrowth. Intercumulus phases may be plagioclase, clinopyroxene, orthopyroxene, pale brown amphibole and magnetite. Where pore material is present, it is composed of plagioclase, clinopyroxene, orthopyroxene, hornblende and ores. Cr-spinel occurs mainly at the bottom of the sequence (Cr₂O₃ between 30·5 and 39·8 per cent), while magnetite appears as a very rare phase in the upper section. Olivine, orthopyroxene, clinopyroxene exhibit slight cryptic variation (Mg × 100/(Mg + Fe) in the range 90–79, 90–70, 93–72 respectively). The investigated dolerites are non-cumulus rocks where clinopyroxene may be more magnesian than in the uppermost gabbros. The cumulate sequence and dolerites underwent variable but generally slight spilitization, in contrast to the overlying lavas. The sequence was generated through crystal accumulation probably from periodic pulses of tholeiitic magma; newly injected magma batches mixing with magma fractions already differentiated in the magma chamber. The high fluid pressure evidenced by the fluid inclusions in plagioclase and the whole chemical trend of the cumulate sequence are consistent with a genesis above a subduction zone, as already hypothesized for the overlying lavas.     

Linkages between the North Pacific Oscillation and central tropical Pacific SSTs at low frequencies

The North Pacific Oscillation (NPO) recently (re-)emerged in the literature as a key atmospheric mode in Northern Hemisphere climate variability, especially in the Pacific sector. Defined as a dipole of sea level pressure (SLP) between, roughly, Alaska and Hawaii, the NPO is connected with downstream weather conditions over North America, serves as the atmospheric forcing pattern of the North Pacific Gyre Oscillation (NPGO), and is a potential mechanism linking extratropical atmospheric variability to El Ni?o events in the tropical Pacific. This paper explores further the forcing dynamics of the NPO and, in particular, that of its individual poles. Using observational data and experiments with a simple atmospheric general circulation model (AGCM), we illustrate that the southern pole of the NPO (i.e., the one near Hawaii) contains significant power at low frequencies (7–10?years), while the northern pole (i.e., the one near Alaska) has no dominant frequencies. When examining the low-frequency content of the NPO and its poles separately, we discover that low-frequency variations (periods >7?years) of the NPO (particularly its subtropical node) are intimately tied to variability in central equatorial Pacific sea surface temperatures (SSTs) associated with the El Ni?o-Modoki/Central Pacific Warming (CPW) phenomenon. This result suggests that fluctuations in subtropical North Pacific SLP are important to monitor for Pacific low-frequency climate change. Using the simple AGCM, we also illustrate that variability in central tropical Pacific SSTs drives a significant fraction of variability of the southern node of the NPO. Taken together, the results highlight important links between secondary modes (i.e., CPW-NPO-NPGO) in Pacific decadal variability, akin to already established relationships between the primary modes of Pacific climate variability (i.e., canonical El Ni?o, the Aleutian Low, and the Pacific Decadal Oscillation).     

Field application of calcite Dispersed Alkaline Substrate (calcite-DAS) for passive treatment of acid mine drainage with high Al and metal concentrations

Passive treatment systems are widely used for remediation of acid mine drainage (AMD), but existing designs are prone to clogging or loss of reactivity due to Al- and Fe-precipitates when treating water with high Al and heavy metal concentrations. Dispersed alkaline substrate (DAS) mixed from a fine-grained alkaline reagent (e.g. calcite sand) and a coarse inert matrix (e.g. wood chips) had shown high reactivity and good hydraulic properties in previous laboratory column tests. In the present study, DAS was tested at pilot field scale in the Iberian Pyrite Belt (SW Spain) on metal mine drainage with pH near 3.3, net acidity 1400–1650 mg/L as CaCO₃, and mean concentrations of 317 mg/L Fe (95% Fe(II)), 311 mg/L Zn, 74 mg/L Al, 20 mg/L Mn, and 1.5–0.1 mg/L Cu, Co, Ni, Cd, As and Pb. The DAS-tank removed an average of 870 mg/L net acidity as CaCO₃ (56% of inflow), 25% Fe, 93% Al, 5% Zn, 95% Cu, 99% As, 98% Pb, and 14% Cd, but no Mn, Ni or Co. Average gross drain pipe alkalinity was 181 mg/L as CaCO₃, which increased total Fe removal to 153 mg/L (48%) in subsequent sedimentation ponds. Unfortunately, the tank suffered clogging problems due to the formation of a hardpan of Al-rich precipitates. DAS lifetime could probably be increased by lowering Al-loads.     

Salt Waters of the Northern Apennine Foredeep Basin (Italy): Origin and Evolution

The salt waters from the Emilia-Romagna sector of the Northern Apennine Foredeep have been investigated using major and trace element and stable isotope (δ²H, δ¹⁸O, δ³⁷Cl, δ⁸¹Br and ⁸⁷Sr/⁸⁶Sr ratio). Ca, Mg, Na, K, Sr, Li, B, I, Br and SO₄ vs. Cl diagrams suggest the subaerial evaporation of seawater beyond gypsum and before halite precipitation as primary process to explain the brine’s salinity, whereas saline to brackish waters were formed by mixing of evaporated seawater and water of meteoric origin. A diagenetic end-member may be a third component for mud volcanoes and some brackish waters. Salinization by dissolution of (Triassic) evaporites has been detected only in samples from the Tuscan side of the Apennines and/or interacting with the Tuscan Nappe. In comparison with the seawater evaporation path, Ca–Sr enrichment and Na–K–Mg depletion of the foredeep waters reveal the presence of secondary processes such as dolomitization–chloritization, zeolitization–albitization and illitization. Sulfate concentration, formerly buffered by gypsum-anhydrite deposition, is heavily lowered by bacterial and locally by thermochemical reduction during burial diagenesis. From an isotopic point of view, data of the water molecule confirm mixing between seawater and meteoric end-members. Local ¹⁸O-shift up to +11‰ at Salsomaggiore is related to water–rock interaction at high temperature (≈150°C) as confirmed by chemical (Mg, Li, Ca distribution) and isotopic (SO₄–H₂O) geothermometers. ³⁷Cl/³⁵Cl and ⁸¹Br/⁷⁹Br ratios corroborate the marine origin of the brines and evidence the diffusion of halogens from the deepest and most saline aquifers toward the surface. The ⁸⁷Sr/⁸⁶Sr ratio suggests a Miocene origin of Sr and rule out the hypothesis of a Triassic provenance of the dissolved components for the analyzed waters issuing from the Emilia-Romagna sector of the foredeep. Waters issuing from the Tuscan side of the Apennines and from the Marche sector of the foredeep show higher ⁸⁷Sr/⁸⁶Sr ratios because of the interaction with siliciclastic rocks.     

Numerical identification of soil-structure interaction pressures

A numerical approach is presented for the ‘identification’ (or back calculation) of the earth pressure acting on embedded or retaining structures. The procedure is applicable to structures of any shape and requires a set of in situ measurements that may include displacements of points on the structure, values of concentrated forces, values of distributed loads at some locations, etc. Possible limiting values of the unknown loads, non-linear structural behaviour, varying accuracies of the input data are accounted for in the problem formulation. Depending upon the type of problem, the solution is reached by means of the unconstrained or constrained minimization of a suitably defined error function. As an example, the proposed approach is applied to the identification of the earth pressure acting on some typical geotechnical engineering structures.