Integrating process‐based flow and temperature models to assess riparian forests and temperature amelioration in salmon streams

The importance of riparian tree cover in reducing energy inputs to streams is increasingly recognized in schemes to mitigate climate change effects and protect freshwater ecosystems. Assessing different riparian management strategies requires catchment‐scale understanding of how different planting scenarios would affect the stream energy balance, coupled with a quantitative assessment of spatial patterns of streamflow generation. Here, we use the physically based MIKE SHE model to integrate simulations of catchment‐scale run‐off generation and in‐stream hydraulics with a heat transfer model. This was calibrated to model the spatio‐temporal distribution of hourly stream water temperature during warm low flow periods in a Scottish salmon stream. The model was explored as a “proof of concept” for a tool to investigate the effects of riparian management on high stream water temperatures that could affect juvenile Atlantic salmon. Uncertainty was incorporated into the assessment using the generalized likelihood uncertainty estimation approach. Results showed that by decreasing both the warming (daylight hours) and the cooling (night‐time hours) rates, forest cover leads to a reduction of the temperature range (with a delay of the time to peak by up to 2 hr) and can therefore be effectively used to moderate projected climate change effects. The modelling presented here facilitated the quantification of potential mitigating effects of alternative riparian management strategies and provided a valuable tool that has potential to be utilized as an evidence base for catchment management guidance.     

Axisymmetric Simulations of Pulsar Wind Nebulae with a New Central-Type Scheme for Relativistic MHD

A third-order shock-capturing numerical scheme for three-dimensional special relativistic magnetohydrodynamics (3-D RMHD) is presented. Due to the difficulties in developing exact or even approximate Riemann solvers in RMHD, a simple two-speed central-type solver that requires the knowledge of only the local fast magnetosonic velocities is employed. First results from axisymmetric simulations of Pulsar Wind Nebulae (PWNe) are finally presented. We show that when the wind magnetization is high enough, magnetic hoop stresses at work downstream of the termination shock are able to collimate a polar jet-like outflow, with velocities similar to those observed in the Crab and Vela Nebulae (v ≈ 0.5-0.7c).     

The Ferrara hard X-ray facility for testing/calibrating hard X-ray focusing telescopes

We will report on the current configuration of the X-ray facility of the University of Ferrara recently used to perform reflectivity tests of mosaic crystals and to calibrate the experiment JEM–X aboard Integral. The facility is now located in the technological campus of the University of Ferrara in a new building (named LARIX laboratory= ̳LARge ̳Italian ̳X-ray facility) that includes a tunnel 100 m long with, on the sides, two large experimental rooms. The facility is being improved for determining the optical axis of mosaic crystals in Laue configuration, for calibrating Laue lenses and hard X-ray mirror prototypes.     

Notes on Hori Method for Canonical Systems

In this paper a slightly different approach is proposed for the process of determining the functions S _m and H _m ^* of the algorithm of the canonical version of Hori method. This process will be referred to as integration theory of the mth order equation of the method. It will be shown that the ordinary differential equation with an auxiliary parameter t ^* as independent variable, introduced through Hori auxiliary system, can be replaced by a partial differential equation in the time t. In this way, the mth order equation of the algorithm assumes a form very similar to the one of other perturbation methods. In virtue of this new approach of the integration theory for Hori method, Lagrange's variational equations introduced by Sessin are revised. As an example, the Duffing equation is solved through this new approach.     

Local kinematics of OB stars

Analysis of observational data of OB stars show an, excellent agreement of the density distributions in space ?(x, y, z) as well as in velocity space \(\rho (\dot x,\dot y,\dot z)\) with the predictions of the density wave theory, the values for the density and velocity fluctuations are explained only by the non-linear theory. These theoretical calculations predict perturbations greater than ±10 km s^?1, consistent with the observations for the velocity field. Thus one should disregard analytical treatments of the linearized equations since they predict maximum perturbations of ±5km s^?1. Another consequence of this is the fact that the Gould's Belt is not a local anomaly, but a local feature of the density waves. The analysis of observational data show that the wave pattern is similar to that of the gas and dust.     

A seismic source zone model for the seismic hazard assessment of the Italian territory

We designed a new seismic source model for Italy to be used as an input for country-wide probabilistic seismic hazard assessment (PSHA) in the frame of the compilation of a new national reference map.

We started off by reviewing existing models available for Italy and for other European countries, then discussed the main open issues in the current practice of seismogenic zoning.

The new model, termed ZS9, is largely based on data collected in the past 10 years, including historical earthquakes and instrumental seismicity, active faults and their seismogenic potential, and seismotectonic evidence from recent earthquakes. This information allowed us to propose new interpretations for poorly understood areas where the new data are in conflict with assumptions made in designing the previous and widely used model ZS4.

ZS9 is made out of 36 zones where earthquakes with M_w > = 5 are expected. It also assumes that earthquakes with M_w up to 5 may occur anywhere outside the seismogenic zones, although the associated probability is rather low. Special care was taken to ensure that each zone sampled a large enough number of earthquakes so that we could compute reliable earthquake production rates.

Although it was drawn following criteria that are standard practice in PSHA, ZS9 is also innovative in that every zone is characterised also by its mean seismogenic depth (the depth of the crustal volume that will presumably release future earthquakes) and predominant focal mechanism (their most likely rupture mechanism). These properties were determined using instrumental data, and only in a limited number of cases we resorted to geologic constraints and expert judgment to cope with lack of data or conflicting indications. These attributes allow ZS9 to be used with more accurate regionalized depth-dependent attenuation relations, and are ultimately expected to increase significantly the reliability of seismic hazard estimates.     

Hydrogeochemical investigations and groundwater provinces of the Friuli Venezia Giulia Plain aquifers,northeastern Italy

Water resources are a key factor, particularly for the planning of the sustainable regional development of agriculture, as well as for socio-economic development in general. A hydrochemical investigation was conducted in the Friuli Venezia Giulia aquifer systems to identify groundwater evolution, recharge and extent of pollution. Temperature, pH, electric conductivity, total dissolved solids, alkalinity, total hardness, SAR, Ca²⁺, Na⁺, K⁺, Mg²⁺, Cl⁻, SO₄ ²⁻, NO₃ ⁻, HCO₃ ⁻ _, water quality and type, saturation indexes and the environmental stable isotope δ¹⁸O were determined in 149 sampling stations. The pattern of geochemical and oxygen stable isotope variations suggests that the sub-surface groundwater (from phreatic and shallow confined aquifers) is being recharged by modern precipitations and local river infiltrations. Four hydrogeological provinces have been recognised and mapped in the Friuli Venezia Giulia Plain having similar geochemical signatures. These provinces have different degrees of vulnerability to contamination. The deep confined groundwater samples are significantly less impacted by surface activities; and it appears that these important water resources have very low recharge rates and would, therefore, be severely impacted by overabstraction.     

Calcium isotope (δCa) fractionation along hydrothermal pathways, Logatchev field (Mid-Atlantic Ridge, 14°45′N)

We investigate the Logatchev Hydrothermal Field at the Mid-Atlantic Ridge, 14°45′N to constrain the calcium isotope hydrothermal flux into the ocean. During the transformation of seawater to a hydrothermal solution, the Ca concentration of pristine seawater ([Ca]_SW) increases from about 10 mM to about 32 mM in the hydrothermal fluid endmember ([Ca]_HydEnd) and thereby adopts a δ^44/40Ca_HydEnd of −0.95 ± 0.07‰ relative to seawater (SW) and a ⁸⁷Sr/⁸⁶Sr isotope ratio of 0.7034(4). We demonstrate that δ^44/40Ca_HydEnd is higher than that of the bedrock at the Logatchev field. From mass balance calculations, we deduce a δ^44/40Ca of −1.17 ± 0.04‰ (SW) for the host-rocks in the reaction zone and −1.45 ± 0.05‰ (SW) for the isotopic composition of the entire hydrothermal cell of the Logatchev field. The values are isotopically lighter than the currently assumed δ^44/40Ca for Bulk Earth of −0.92 ± 0.18‰ (SW) [Skulan J., DePaolo D. J. and Owens T. L. (1997) Biological control of calcium isotopic abundances in the global calcium cycle. Geochim. Cosmochim. Acta61,(12) 2505-2510] and challenge previous assumptions of no Ca isotope fractionation between hydrothermal fluid and the oceanic crust [Zhu P. and Macdougall J. D. (1998) Calcium isotopes in the marine environment and the oceanic calcium cycle. Geochim. Cosmochim. Acta62,(10) 1691-1698; Schmitt A. -D., Chabeaux F. and Stille P. (2003) The calcium riverine and hydrothermal isotopic fluxes and the oceanic calcium mass balance. Earth Planet. Sci. Lett. 6731, 1-16]. Here we propose that Ca isotope fractionation along the fluid flow pathway of the Logatchev field occurs during the precipitation of anhydrite. Two anhydrite samples from the Logatchev Hydrothermal Field show an average fractionation of about Δ^44/40Ca = −0.5‰ relative to their assumed parental solutions. Ca isotope ratios in aragonites from carbonate veins from ODP drill cores indicate aragonite precipitation directly from seawater at low temperatures with an average δ^44/40Ca of −1.54 ± 0.08‰ (SW). The relatively large fractionation between the aragonite precipitates and seawater in combination with their frequent abundance in weathered mafic and ultramafic rocks suggest a reconsideration of the marine Ca isotope budget, in particular with regard to ocean crust alteration.     

A model for the viscosity of rhyolite as a function of H2O-content and pressure: A calibration based on centrifuge piston cylinder experiments

The Newtonian viscosity of synthetic rhyolitic liquids with 0.15-5.24 wt% dissolved water was determined in the interval between 580 and 1640 °C and pressures of 1 atm and 5-25 kbar. Measurements were performed by combining static and accelerated (up to 1000g) falling sphere experiments on water-bearing samples, with high temperature concentric cylinder experiments on 0.15 wt% H₂O melts. These methods allowed viscosity determinations between 10² and 10⁷ Pa s, and cover the complete range of naturally occurring magmatic temperatures, pressures, and H₂O-contents for rhyolites.Our viscosity data, combined with those from previous studies, were modeled by an expression based on the empirical Vogel-Fulcher-Tammann equation, which describes viscosities and derivative properties (glass transition temperature T_g, fragility m, and activation volume of viscous flow Va) of silicic liquids as a function of P-T-X_(H2O). The fitted expressions do not account for composition-dependent parameters other than X_(H2O) and reproduce the entire viscosity database for silicic liquids to within 3.0% average relative error on log η (i.e. std. error of estimate of 0.26 log units).The results yield the expected strong decrease of viscosity with temperature and water content, but show variable pressure dependencies. Viscosity results to be strongly affected by pressure at low pressures; an effect amplified at low temperatures and water contents. Fragility, as a measure for the deviation from Arrhenian behavior, decreases with H₂O-content but is insensitive to pressure. Activation volumes are always largely negative (e.g., less than −10 cm³/mol) and increase strongly with H₂O-content. Variations in melt structure that may account for the observed property variations are discussed.