Matching hydrologic response to measured effective hydraulic conductivity

The objective of this study was to test the practicability of defining hydrologic response units as combinations of soil, land use and topography for modelling infiltration at the hillslope and catchment scales. In an experimental catchment in the East African Highlands (Kwalei, Tanzania), three methods of measuring infiltration were compared for their ability to capture the spatial variability of effective hydraulic conductivity: the constant head (CH) method; the tension infiltration (TI) method; and the mini‐rainfall simulation (RS) method. The three methods yielded different probability distributions of effective hydraulic conductivity and suggested different types of hydrologic response units. Independently from these measurements, the occurrence of infiltration‐excess overland flow was monitored over an area of 6 ha by means of overland flow detectors. The observed pattern of overland flow occurrence did not match any of the patterns suggested by the infiltration measurements. Instead, clusters of spots with overland flow were practically independent from field borders. Geostatistical analysis of the overland flow confirmed the absence of spatial correlation for distances over 40 m. The RS method yielded the pattern closest to the observations, probably because the method simulated better the processes that trigger infiltration‐excess overland flow, i.e. soil sealing and infiltration through macroporosity. The RS hydrologic response unit correlated significantly with observed overland flow frequency. However, the location of clusters and ‘hot spots’ of overland flow remained largely unexplained by land use, soil and topographic variables. It is concluded that using such landscape variables to define hydrologic units may create artificial boundaries that do no correspond to physical realities, especially if the stochastic component within hydrologic units is neglected. Copyright © 2005 John Wiley & Sons, Ltd.     

Mn-Cr systematics of carbonates in CM chondrites: Implications for the timing and duration of aqueous alteration

CM chondrites contain carbonates and other secondary minerals such as phyllosilicates, sulfides, sulfates, oxides and hydroxides that are believed to have formed by aqueous alteration reactions on their parent asteroid. We report in situ Mn-Cr isotope measurements in the highly aqueously altered CM2.1 chondrites QUE 93005 and ALH 83100 using secondary ion mass spectrometry (Cameca ims-1270 ion microprobe). The ⁵³Cr excesses are correlated with the ⁵³Mn/⁵⁵Mn ratio and result from the in situ decay of ⁵³Mn, a short-lived radioisotope with a half-life of 3.7 Ma. If we assume that carbonate grains in samples QUE 93005 and ALH 83100 are cogenetic, then the excesses define initial ⁵³Mn/⁵⁵Mn ratios ((⁵³Mn/⁵⁵Mn)₀) of (4.1 ± 1.2) × 10⁻⁶ and (5.1 ± 1.7) × 10⁻⁶, respectively. These values are comparable to those in carbonates from other CM chondrites as reported in the literature. Initial ⁵³Mn/⁵⁵Mn ratios for calculated model isochrones for individual carbonate grains range from (3.8 ± 1.4) × 10⁻⁶ to (4.8 ± 2.1) × 10⁻⁶ for QUE 93005 and from (3.1 ± 1.6) × 10⁻⁶ to (1.3 ± 0.5) × 10⁻⁵ for ALH 83100. A possible interpretation for the ranges in (⁵³Mn/⁵⁵Mn)₀ could be that alteration in individual CM chondrites was episodic and occurred over an extended period of time. However, isochrones based on the entire set of carbonate grains in each of the CM chondrites imply that the degree of aqueous alteration is roughly correlated with the age of carbonate formation in CM chondrites of different subtypes and that alteration on the CM parent asteroid started contemporaneously with or shortly after CAI formation and lasted at least 4 Ma.     

Rates of hydrothermal cooling of new oceanic upper crust derived from lithium-geospeedometry

Episodic emplacement and cooling of lavas and dikes at mid-ocean ridges leads to large fluctuations in hydrothermal fluxes and biological activity. However, the processes operating beneath the seafloor during these transient events such as permeability creation and dike cooling are poorly understood. We have developed a new approach to determine the cooling rate of the sheeted dike complex based on the extent of diffusion of lithium from plagioclase into clinopyroxene during cooling. We have calibrated this Li-geospeedometer using new high-temperature experiments to determine both the temperature dependence of the partitioning of Li between plagioclase and clinopyroxene and the diffusion coefficient for Li in clinopyroxene. Application of this method to lavas and dikes from ODP Hole 504B shows that cooling rates vary dramatically with depth in the upper oceanic crust. Extremely rapid cooling rates (> 450 °C hr^− 1) in the upper part of the sheeted dike complex are sufficient to power hydrothermal megaplume formation within the overlying water column.     

Chronology of Vesuvius’ activity from A.D. 79 to 1631 based on archeomagnetism of lavas and historical sources

The activity of Vesuvius between A.D. 79 and 1631 has been investigated by means of precise archaeomagnetic dating of primary volcanic deposits and taking into account the stratigraphy of lavas and tephra, historical written accounts, archaeological evidence related to the developing urbanisation, and radiocarbon ages. We found that the historical records are highly useful in constraining the timing of the main events, even if the data are often too scarce and imprecise for ascertaining the details of all phases of activity, especially their magnitude and emplacement of all the deposit types. In addition, some eruptions that took place in the 9^th and 10^th centuries appear to be unnoticed by historians. The archaeomagnetic study involved 26 sites of different lavas and 2 pyroclastic deposits. It shows that within the 15 centuries which elapsed between A.D. 79 and 1631, the effusive activity of Vesuvius clustered in the relatively short period of time between A.D. 787 and 1139 and was followed by a 5-century-long repose period. During this time Vesuvius prepared itself for the violent explosive eruption of 1631. The huge lavas shaping the morphology of the coast occurred largely through parasitic vents located outside the Mount Somma caldera. One of these parasitic vents is located at low elevation, very close to the densely inhabited town of Torre Annunziata. Among the various investigated lavas, a number of which were previously attributed to the 1631 eruption, none is actually younger than the 12^th century. Therefore it is definitively concluded that the destructive 1631 event was exclusively explosive.Editorial handling: J. McPhie     

Meiofaunal colonisation on artificial substrates: a tool for biomonitoring the environmental quality on coastal marine systems

The use of artificial substrates for monitoring environmental quality and contaminant impact has been repeatedly suggested, but little information is yet available on the applicability and the sensitivity of this approach. One of the crucial points to be addressed for using artificial substrates, as tools mirroring actual benthic community state, is assessing whether, and how, they reflect in situ conditions characterising natural sediments. In this study we investigated short-term dynamics of meiofaunal colonisation on artificial substrates (bottle brushes) deployed in impacted and "pristine" (control) sites. Meiofaunal colonisation reached a threshold level after 4-6 days in both impacted and control sites but with significantly higher densities in artificial substrates deployed in the control site. Meiofaunal parameters in artificial substrates reflected those of meiofaunal assemblages inhabiting natural sediments where artificial substrates were deployed. Colonised artificial substrates were then transplanted from impacted to control sites. This caused a significant increase of meiofaunal abundance, which after 5 days reached values indistinguishable from the non-impacted control. Given the rapid colonisation time, meiofaunal sensitivity to changing environmental conditions, it can be concluded that artificial substrates can represent an useful tool to be further developed for routine and low-cost monitoring studies aiming at integrating biological indicators of environmental quality.     

Benthic input rates predict seagrass (Posidonia oceanica) fish farm-induced decline

Fish farms represent a growing source of anthropogenic disturbance to benthic communities, and efficient predictors of such impacts are urgently needed. We explored the effects of fish farm benthic organic and nutrient inputs on the population dynamics of a key seagrass species (Posidonia oceanica) in four Mediterranean deep meadows adjacent to sea bream and sea bass farms. We performed two annual plant censuses on permanent plots at increasing distance from farms and measured benthic sedimentation rates around plots. High shoot mortality rates were recorded near the cages, up to 20 times greater than at control sites. Recruitment rates increased in variability but could not compensate mortality, leading to rapid seagrass decline within the first 100 m from cages. Seagrass mortality increased with total sedimentation rates (K=0.55, p<0.0002), and with organic matter (K=0.50, p=0.001), total nitrogen (K=0.46, p=0.002) and total phosphorus (K=0.56, p<3.10(-5)) inputs. P. oceanica decline accelerated above a phosphorus loading threshold of 50mg m(-2)day(-1). Phosphorus benthic sedimentation rate seems a powerful predictor of seagrass mortality from fish farming. Coupling direct measurements of benthic sedimentation rates with dynamics of key benthic species is proposed as an efficient strategy to predict fish farm impacts to benthic communities.     

A strong motion network in northern Italy: detection capabilities and first analysis

The necessity of a dense network in Northern Italy started from the lack of available data after the occurrence of the 24th November 2004, Ml 5.2, Salò earthquake. Since 2006, many efforts have been made by the INGV (Italian National Institute for Geophysic and Vulcanology), Department of Milano-Pavia (hereinafter INGV MI-PV), to improve the strong-motion monitoring of the Northern Italy regions. This activity led to the installation of a strong-motion network composed by 20 accelerometers, 4 coupled with 20-bits Lennartz Mars88 recorders, 12 coupled with 24-bits Reftek 130 recorders and 4 coupled with 24-bits Gaia2 recorders. The network allow us to reduce, in the area under study, the average inter-distances between strong-motion stations from about 40 km (at November 2004) to 15 km. At present the network includes nine 6-channels stations where velocity sensors work together the strong-motion ones. The data transmission is assured by modem-gsm, with the exception of four stations that send data in real time through a TCP/IP protocol. In order to evaluate different site responses, the stations have been installed both in free field and near (or inside) public buildings, located in the center of small villages. From June 2006 to December 2008 a dataset of 94 events with local magnitude range from 0.7 to 5.1 has been collected. An ad hoc data-processing system have been created in order to provide, after each recorded event, engineering parameters such as peak ground acceleration (PGA) and velocity (PGV), response spectra (SA and PSV), Arias and Housner intensities. Data dissemination is achieved through the web site , while the waveforms are distributed through the Italian strong motion database ().     

Calibration of rainfall-runoff models in ungauged basins: A regional maximum likelihood approach

Parameter estimation for rainfall-runoff models in ungauged basins is a challenging task that is receiving significant attention by the scientific community. In fact, many practical applications suffer from problems induced by data scarcity, given that hydrological observations are often sparse or unavailable. This study focuses on regional calibration for a generic rainfall-runoff model. The maximum likelihood function in the spectral domain proposed by Whittle [40] is approximated in the time domain by maximising the fit of selected statistics of the river flow process, with the aim to propose a calibration procedure that can be applied at regional scale. Accordingly, the statistics above are related to the dominant climate and catchment characteristics, through regional regression relationships. The proposed technique is applied to the case study of 4 catchments located in central Italy, which are treated as ungauged and are located in a region where detailed hydrological, as well as geomorphologic and climatic information, is available. The results obtained with the regional calibration are compared with those provided by a classical least squares calibration in the time domain. The outcomes of the analysis confirm the potential of the proposed methodology and show that regional information can be very effective for setting up hydrological models.     

Discrimination of induced seismicity by full moment tensor inversion and decomposition

Human activities, including operations related to mining and reservoir exploitation, may induce seismicity and pose a risk for population and infrastructures. While different observations are commonly used to assess the origin of earthquakes, there is a lack of rules and methods for the discrimination between natural and induced seismicity. The inversion and decomposition of the full moment tensor and the observation of relevant deviation from a pure double couple (DC) model may be an indicator for induced seismicity. We establish here a common procedure to analyse a set of natural and induced events of similar magnitude, which occurred in Germany and neighbouring regions. The procedure is based on an inversion method and on a consistent velocity model and recording network. Induced seismicity is recorded during different mining and/or reservoir exploitations. Moment tensors are inverted using a multi-step inversion approach. This method, which was successfully applied in previous studies at regional and teleseismic distances, is further developed here to account for full moment tensor analysis. We first find a best DC solution and then perform a full moment tensor inversion, fitting full waveforms amplitude spectra at regional distances. The moment tensor solution is decomposed into DC, compensated linear vector dipole and isotropic terms. The discrimination problem is then investigated through the evaluation of distributions of non-DC source components for natural and induced data sets. Results illustrate the potential of the inversion and discrimination approach. Additional detailed analyses are carried out for the two most significant induced earthquakes, and rupture models are compared with the full moment tensor solutions.