An integrated procedure to evaluate hydrological models

The growing concern for health‐related problems deriving from pollutants leaching is driving national and international administrations to support the development of tools for evaluating the effects of alternate management scenarios and identifying vulnerable areas. Cropping systems models are powerful tools for evaluating leachates under different environmental, social, and management conditions. As percolating water is the transport vehicle for pollutants transport in soil, a reliable evaluation of water balance models is a fundamental prerequisite for investigating pesticides and nitrate fate. As specific approaches for the evaluation of multi‐layer evolution of state variables are missing, we propose a fuzzy‐based, integrated indicator (I_SWC: 0, best; 1, worst) for a comprehensive evaluation of soil water content (SWC) simulations. We aggregated error metrics with others quantifying the homogeneity of errors across different soil layers, the capability of models to reproduce complex dynamics function of both time and soil depth, and model complexity. We tested I_SWC on a sample dataset where the models CropSyst and CERES‐Wheat were used to simulate SWC for winter wheat systems. I_SWC revealed that, in the explored conditions, the global assessment of the two models' performances allowed identification of CropSyst as the best (average I_SWC = 0·441, with a value of 0·537 obtained by CERES‐Wheat), although each model prevailed for some of the metrics. CropSyst presented the highest accuracy (average agreement module = 0·400), whereas CERES‐Wheat's accuracy was slightly worse, although achieved with a simplified modelling approach (average Akaike Information Criterion = − 230·44), thereby favouring large‐area applicability. The non‐univocal scores achieved by the models for the different metrics support the use of multi‐metric evaluation approaches for quantifying the different aspects of water balance model performances. Copyright © 2010 John Wiley & Sons, Ltd.     

Crystal-chemical examination of natural diopsides: Some geometrical indications of Si-Ti tetrahedral substitution

Summary The crystal structures of 20 natural diopsides have been refined. For some samples, microprobe analyses show a Si deficiency in the tetrahedron with no A1 present; this suggests that some Ti⁴⁺ may enter the tetrahedral site. Geometrical variations in these samples with respect to the others are discussed in the light of this Si Ti⁴⁺ substitution. The volumes of their M1 sites are larger because of the absence of R³⁺ cations, the octahedron being very regular; in the same way, the volumes of M2 sites of these samples are greater not because of different site occupancies but because of significant charge modifications around the O3 and O1 oxygens, caused by the tetrahedral site. The low distortion of M1 site and the lengthening of (M2-O1) and M2-O3 bond lengths are a clue to the presence of small amounts of Ti at the tetrahedral site.

---

Kristallchemische untersuchung natürlicher diopside: Geometrische hinweise für tetraedrische Si-Ti substitution

Zusammenfassung Die Kristallstrukturen von zwanzig natürlichen Diopsiden sind verfeinert worden. In einigen Proben zeigen Mikrosondenanalysen zu geringe Si-Werte im Tetraeder ohne daß A1 anwesend wäre. Dies weist darauf hin, daß Ti⁴⁺ Tetraederplätze einnehmen könnte. Geometrische Variationen in diesen Proben werden im Lichte der genannten Substitution diskutiert. Die Volumen von M1-Plätzen sind größer wegen der Abwesenheit von R³⁺ Kationen, wobei das Oktaeder sehr gleichmäßig ist. Die Volumina von M2-Plätzen sind noch größer nicht wegen verschiedener Besetzungen dieser Plätze sondern wegen signifikanter Ladungsmodifikationen im Bereich von O3 und O1 Sauerstoffen, die durch die tetraedrische Position bewirkt wird. Die geringe Distortion der M1-Plätze und die Verlängerung der (M2-O1) und (M2-O3) Abstände weisen auf die Anwesenheit von geringen Titangehalten auf Tetraederplätzen hin.

     

Petrogenesis of calc-alkaline,shoshonitic and associated ultrapotassic Oligocene volcanic rocks from the Northwestern Alps,Italy

Along the Western Alps there is geological evidence of late-Alpine (Oligocene) magmatic activity which clearly postdates the Lepontine (Eocene-early Oligocene) metamorphism and related deformation of the Alpine nappe pile. This magmatic activity was notably delayed in relation to the most important convergent processes and may be related to buoyancy of lithosphere, tensional tectonics and thermal updoming subsequent to the collision between the Eurasian and African plates. The geochemical features of the rocks and the geophysical characteristics of the Alpine chain, suggest that: (a) shoshonitic and calcalkaline melts may have been generated by partial melting of metasomatized peridotitic material and subsequent fractional crystallization and crustal contamination; silicic andesites and latites, however, could have been also derived from metasomatized eclogite or deep continental crust material; (b) the ultrapotassic lamprophyres with high K, P, LREE, Th, Zr, U and high ⁸⁷Sr/⁸⁶Sr ratios were generated by partial melting of strongly metasomatized mantle; the varied Sr-isotopic ratios may partially also reflect additional radiogenic component from the continental crust following magma segregation from the source.     

Potential merits for substorm research from imaging of charge-exchange neutral atoms

The in situ observations of the Earth magnetosphere performed over the past decades of space research have provided a rather good understanding of many partial localized processes of the magnetospheric substorm. The continuing lack of global observations inhibits the construction of a coherent picture of the substorm as a whole, which is actually determined by the coupling of the partial processes. In this context the importance of global observations for the advancement of magnetospheric substorm studies is critical. This paper presents briefly a promising technique of global observations, namely the imaging of charge exchange neutral atoms, or neutral atom imaging (NAI) of the magnetosphere. Model and theoretical estimates of charge-exchange neutral atom fluxes, as well as appropriate spacecraft orbit and instrumentation requirements are presented and discussed for specific regions of interest and vantage points. The potential merits of NAI for substorm research are presented along with possible combinations with other types of observational methods. Substorm issues that would benefit from NAI should include among others the assessment of the ionospheric contribution to the hot magnetospheric plasma, the relative importance of various ionospheric ion source regions, the resolution of spatial and temporal characteristics of substorm ion injections. NAI observations can be precious complements to local observations and lead to the understanding of how local processes, many of which are resolved quite well today, combine to form the global process of the magnetospheric substorm.     

On estimating frequency-magnitude relations from heterogeneous catalogs

The paper studies the effect of magnitude errors on heterogeneous catalogs, by applying the apparent magnitude theory (seeTinti andMulargia, 1985a), which proves to be the most natural and rigorous approach to the problem. Heterogeneities in seismic catalogs are due to a number of various sources and affect both instrumental as well as noninstrumental earthquake compilations.The most frequent basis of heterogeneity is certainly that the recent instrumental records are to be combined with the historic and prehistoric event listings to secure a time coverage, considerably longer than the recurrence time of the major earthquakes. Therefore the case which attracts the greatest attention in the present analysis is that of a catalog consisting of a subset of higher quality data, generallyS ₁, spanning the interval T ₁ (the instrumental catalog), and of a second subset of more uncertain magnitude determination, generallyS ₂, covering a vastly longer interval T ₂ (the historic and/or the geologic catalog). The magnitude threshold of the subcatalogS ₁ is supposedly smaller than that ofS ₂, which, as we will see, is one of the major causes of discrepancy between the apparent magnitude and the true magnitude distributions. We will further suppose that true magnitude occurrences conform to theGutenberg-Richter (GR) law, because the assumption simplified the analysis without reducing the relevancy of our findings.The main results are: 1) the apparent occurrence rate exceeds the true occurrence rate from a certain magnitude onward, saym _GR; 2) the apparent occurrence rate shows two distinct GR regimes separated by an intermediate transition region. The offset between the two regimes is the essential outcome ofS ₁ being heterogeneous with respect toS ₂. The most important consequences of this study are that: 1) it provides a basis to infer the parameters of the true magnitude distribution, by correcting the bias deriving from heterogeneous magnitude errors; 2) it demonstrates that the double GR decay, that several authors have taken as the incontestable proof of the failure of the GR law and of the experimental evidence of the characteristic earthquake theory, is instead perfectly consistent with a GR-type seismicity.     

Coupling of oceanic and continental crust during Eocene eclogite-facies metamorphism: evidence from the Monte Rosa nappe,western Alps

High precision U–Pb geochronology of rutile from quartz–carbonate–white mica–rutile veins that are hosted within eclogite and schist of the Monte Rosa nappe, western Alps, Italy, indicate that the Monte Rosa nappe was at eclogite-facies metamorphic conditions at 42.6 ± 0.6 Ma. The sample area [Indren glacier, Furgg zone; Dal Piaz (2001) Geology of the Monte Rosa massif: historical review and personal comments. SMPM] consists of eclogite boudins that are exposed inside a south-plunging overturned synform within micaceous schist. Associated with the eclogite and schist are quartz–carbonate–white mica–rutile veins that formed in tension cracks in the eclogite and along the contact between eclogite and surrounding schist. Intrusion of the veins at about 42.6 Ma occurred at eclogite-facies metamorphic conditions (480–570°C, >1.3–1.4 GPa) based on textural relations, oxygen isotope thermometry, and geothermobarometry. The timing of eclogite-facies metamorphism in the Monte Rosa nappe determined in this study is identical to that of the Gran Paradiso nappe [Meffan-Main et al. (2004) J Metamorphic Geol 22:261–281], confirming that these two units have shared the same Alpine metamorphic history. Furthermore, the Gran Paradiso and Monte Rosa nappes underwent eclogite-facies metamorphism within the same time interval as the structurally overlying Zermatt-Saas ophiolite [∼50–40 Ma; e.g., Amato et al. (1999) Earth Planet Sci Lett 171:425–438; Mayer et al. (1999) Eur Union Geosci 10:809 (abstract); Lapen et al. (2003) Earth Planet Sci Lett 215:57–72]. The nearly identical P–T–t histories of the Gran Paradiso, Monte Rosa, and Zermatt-Saas units suggest that these units shared a common Alpine tectonic and metamorphic history. The close spatial and temporal associations between high pressure (HP) ophiolite and continental crust during Alpine orogeny indicates that the HP internal basement nappes in the western Alps may have played a key role in exhumation and preservation of the ophiolitic rocks through buoyancy-driven uplift. Coupling of oceanic and continental crust may therefore be critical in preventing permanent loss of oceanic crust to the mantle.     

Modelling of mantle postglacial relaxation in axisymmetric geometry with a composite rheology and a glacial load interpolated by adjusted spherical harmonics analysis

Although studies on glacial isostatic adjustment usually assume a purely linear rheology, we have previously shown that mantle relaxation after the melting of Laurentide ice sheet is better described by a composite rheology including a non-linear term. This modelling is, however, based on axially symmetric geometry and glacial forcing derived from ICE-3G and suffers from a certain amount of arbitrariness in the definition of the ice load. In this work we apply adjusted spherical harmonics analysis to interpolate the ice thicknesses of ICE-3G and ICE-1 glaciological models. This filters out the non-axisymmetric components of the ice load by considering only the zonal terms in the spherical harmonics expansion. The resulting load function is used in finite-element simulation of postglacial rebound to compare composite versus purely linear rheology. Our results confirm that composite rheology can explain relative sea level (RSL) data in North America significantly better than a purely linear rheology. The performance of composite rheology suggests that in future investigations, it may be better to use this more physically realistic creep law for modelling mantle deformation induced by glacial forcing.     

Evaluation of the SEDD model for predicting sediment yield at the Sicilian experimental SPA2 basin

In this paper a spatially distributed model of the hillslope sediment delivery processes, named the sediment delivery distributed (SEDD) model, is initially reviewed; the model takes into account the sediment delivery processes due to both the hillslope sediment transport and the effects of slope curvature. Then the rainfall and sediment yield events measured at the experimental SPA2 basin, in Sicily, are used both to calibrate the SEDD model and to verify the predictive capability of the distributed sediment delivery approach at event scale. For the SPA2 basin discretized into morphological units and stream tubes, the SEDD model is calibrated at event scale using the measurements carried out at the outlet of the experimental basin in the period December 2000–January 2001. The model calibration is used to determine a relationship useful for estimating the unique coefficient β_e of the model by rainfall erosivity factor R_e at event scale. To test the predictive capability of the β_e = f(R_e) relationship, 20 events measured in the period September 2002–December 2005 are used; the comparison between measured sediment yield values and calculated ones for all monitored events shows that the sediment delivery distributed approach has a good predictive ability at event scale. The analysis also shows that estimating β_e by the relationship β_e = f(R_e) gives a better agreement between measured and calculated sediment yields than obtained with the median value β_e,m of all 27 calculated β_e values. Finally the analysis at annual scale, for the period 2000–2005, allows the estimation of the median value β_a,m representative of the annual behaviour. This analysis shows that the sediment delivery distributed approach also has a good predictive ability at annual scale. Copyright © 2006 John Wiley & Sons, Ltd.     

Effectiveness of simple approaches in mitigating residual deformations in buildings

Developments in performance‐based seismic design and assessment approaches have emphasized the importance of considering residual deformations. Recent investigations have also led to a proposed direct displacement‐based design (DDBD) approach which includes an explicit consideration of the expected residual deformations as an integral part of the design process. Having estimated the expected residual deformations in a structure, engineers are faced with the problem of reducing them to meet the targeted performance levels under pre‐defined seismic hazard levels. Previous studies have identified the post‐yield stiffness as a primary factor influencing the magnitude of residual deformations in single degree of freedom and multiple degree of freedom structures. In this paper, a series of simple approaches to increase the post‐yield stiffness of traditional framed and braced systems for the purpose of reducing residual deformations are investigated. These methods do not utilize recentring post‐tensioned technology. This contribution addresses the feasibility of altering the lateral post‐yield stiffness of structural systems by: (i) using different reinforcement materials with beneficial features in their stress–strain behaviour; (ii) re‐designing the section geometry and properties of primary seismic‐resisting elements; and (iii) introducing a secondary elastic frame to act in parallel with the primary system. The efficiency of each of these techniques is investigated through monotonic and cyclic moment‐curvature and non‐linear time‐history analyses. Of these approaches the design and introduction of an elastic secondary system was found to be most effective and consistent in reducing residual deformations. A simplified design approach for achieving the desired increase of a system's post‐yield stiffness is also presented. Copyright © 2007 John Wiley & Sons, Ltd.