Determination of subcatchment and watershed boundaries in a complex and highly urbanized landscape

Urban development significantly alters the landscape by introducing widespread impervious surfaces, which quickly convey surface run‐off to streams via stormwater sewer networks, resulting in “flashy” hydrological responses. Here, we present the inadequacies of using raster‐based digital elevation models and flow‐direction algorithms to delineate large and highly urbanized watersheds and propose an alternative approach that accounts for the influence of anthropogenically modified land cover. We use a semi‐automated approach that incorporates conventional drainage networks into overland flow paths and define the maximal run‐off contributing area. In this approach, stormwater pipes are clustered according to their slope attributes, which define flow direction. Land areas drained by each cluster and contributing (or exporting) flow to a topographically delineated catchment were determined. These land masses were subsequently added or removed from the catchment, modifying both the shape and the size. Our results in a highly urbanized Toronto, Canada, area watershed indicate a moderate net increase in the directly connected watershed area by 3% relative to a topographically forced method; however, differences across three smaller scale subcatchments are greater. Compared to topographic delineation, the directly connected watershed areas of both the upper and middle subcatchments decrease by 5% and 8%, respectively, whereas the lower subcatchment area increases by 15%. This is directly related to subsurface storm sewer pipes that cross topographic boundaries. When directly connected subcatchment area is plotted against total streamflow and flashiness indices using this method, the coefficients of variation are greater (0.93 to 0.97) compared to the use of digital elevation model‐derived subcatchment areas (0.78 to 0.85). The accurate identification of watershed and subcatchment boundaries should incorporate ancillary data such as stormwater sewer networks and retention basin drainage areas to reduce water budget errors in urban systems.     

Bayesian estimation of seismic hazard for two sites in Switzerland

Seismic hazard analysis is based on data and models, which both are imprecise and uncertain. Especially the interpretation of historical information into earthquake parameters, e.g. earthquake size and location, yields ambiguous and imprecise data. Models based on probability distributions have been developed in order to quantify and represent these uncertainties. Nevertheless, the majority of the procedures applied in seismic hazard assessment do not take into account these uncertainties, nor do they show the variance of the results. Therefore, a procedure based on Bayesian statistics was developed to estimate return periods for different ground motion intensities (MSK scale).Bayesian techniques provide a mathematical model to estimate the distribution of random variables in presence of uncertainties. The developed method estimates the probability distribution of the number of occurrences in a Poisson process described by the parameter . The input data are the historical occurrences of intensities for a particular site, represented by a discrete probability distribution for each earthquake. The calculation of these historical occurrences requires a careful preparation of all input parameters, i.e. a modelling of their uncertainties. The obtained results show that the variance of the recurrence rate is smaller in regions with higher seismic activity than in less active regions. It can also be demonstrated that long return periods cannot be estimated with confidence, because the time period of observation is too short. This indicates that the long return periods obtained by seismic source methods only reflects the delineated seismic sources and the chosen earthquake size distribution law.     

Red dust rain within the Spanish Mediterranean area

Dust rain belongs to the climatological conditions of the Iberian Peninsula's Mediterranean seaboard. Traditionally known as muddy' or bloody rains, red dust rain has become a topical issue as a result of more frequent rainfall of this kind in recent years. In spite of the difficulties that studying this phenomenon involves, owing to the lack of systematically kept records at meteorological observatories, details are provided of the high frequency of such phenomena within the Spanish Mediterranean area, using records from other phenological sources. The study analyzes the chemical composition of red dust rain and the atmospheric factors giving rise to such phenomena. Whether there is a relationship betwen the higher frequency of dust rain and possible changes in global atmospheric circulation is as yet unclear.     

About the specification of erosion flux for soft stratified cohesive sediments

The present paper deals with the specification of bed erosion flux that accounts for the effects of sediment-induced stratification in the water column. Owing to difficulties in measuring the bed shear stress _b and the erosive shear strength _s, we suggest a series of methods that combine laboratory and numerical experiments. A simplified turbulent transport model that includes these effects helps to quantify _b and _s. Focusing on soft stratified beds, the present study considers erosion rate formulas of the form =_f exp {[T_b-T_s]} where is a model constant (=1 for Gularte's (1978) formula and =1/2 for Parchure's (1984) formula). First, the bed erosive strength profile _s(Z) is adjusted by forcing the turbulent transport model with measured erosion rates. Second, three procedures are suggested to determine the erosion rate formula coefficients _f and : a global procedure and two different layer-by-layer procedures. Each procedure is applied to an erosion experiment conducted in a rotating annular flume by Villaret and Paulic (1986). The use of the layer-by-layer procedure based on a least squares fitting technique provides a closer fit than the global procedure. The present study points out the complementarity of experimental and numerical approaches and also suggests possible improvements in laboratory test procedures.     

Molecular indicators for palaeoenvironmental change in a Messinian evaporitic sequence (Vena del Gesso, Italy). II: High-resolution variations in abundances and 13C contents of free and sulphur-bound carbon skeletons in a single marl bed

The extractable organic matter of 10 immature samples from a marl bed of one evaporitic cycle of the Vena del Gesso sediments (Gessoso-solfifera Fm., Messinian, Italy) was analyzed quantitatively for free hydrocarbons and organic sulphur compounds. Nickel boride was used as a desulphurizing agent to recover sulphur-bound lipids from the polar and asphaltene fractions. Carbon isotopic compositions (delta vs PDB) of free hydrocarbons and of S-bound hydrocarbons were also measured. Relationships between these carbon skeletons, precursor biolipids, and the organisms producing them could then be examined. Concentrations of S-bound lipids and free hydrocarbons and their delta values were plotted vs depth in the marl bed and the profiles were interpreted in terms of variations in source organisms, 13 C contents of the carbon source, and environmentally induced changes in isotopic fractionation. The overall range of delta values measured was 24.7%, from -11.6% for a component derived from green sulphur bacteria (Chlorobiaceae) to -36.3% for a lipid derived from purple sulphur bacteria (Chromatiaceae). Deconvolution of mixtures of components deriving from multiple sources (green and purple sulphur bacteria, coccolithophorids, microalgae and higher plants) was sometimes possible because both quantitative and isotopic data were available and because either the free or S-bound pool sometimes appeared to contain material from a single source. Several free n-alkanes and S-bound lipids appeared to be specific products of upper-water-column primary producers (i.e. algae and cyanobacteria). Others derived from anaerobic photoautotrophs and from heterotrophic protozoa (ciliates), which apparently fed partly on Chlorobiaceae. Four groups of n-alkanes produced by algae or cyanobacteria were also recognized based on systematic variations of abundance and isotopic composition with depth. For hydrocarbons probably derived from microalgae, isotopic variations are well correlated with those of total organic carbon. A resistant aliphatic biomacromolecule produced by microalgae is, therefore, probably an important component of the kerogen. These variations reflect changes in the depositional environment and early diagenetic transformations. Changes in the concentrations of S-bound lipids induced by variations in conditions favourable for sulphurization were discriminated from those related to variations in primary producer assemblages. The water column of the lagoonal basin was stratified and photic zone anoxia occurred during the early and middle stages of marl deposition. During the last stage of the marl deposition the stratification collapsed due to a significant shallowing of the water column. Contributions from anaerobic photoautotrophs were apparently associated with variations in depth of the chemocline.     

Dynamical evolution of NEAs: Close encounters, secular perturbations and resonances

We discuss the main mechanisms affecting the dynamical evolution of Near-Earth Asteroids (NEAs) by analyzing the results of three numerical integrations over 1 Myr of the NEA (4179) Toutatis. In the first integration the only perturbing planet is the Earth. So the evolution is dominated by close encounters and looks like a random walk in semimajor axis and a correlated random walk in eccentricity, keeping almost constant the perihelion distance and the Tisserand invariant. In the second integration Jupiter and Saturn are present instead of the Earth, and the 3/1 (mean motion) and v ₆ (secular) resonances substantially change the eccentricity but not the semimajor axis. The third, most realistic, integration including all the three planets together shows a complex interplay of effects, with close encounters switching the orbit between different resonant states and no approximate conservation of the Tisserand invariant. This shows that simplified 3-body or 4-body models cannot be used to predict the typical evolution patterns and time scales of NEAs, and in particular that resonances provide some fast-track dynamical routes from low-eccentricity to very eccentric, planet-crossing orbits.On leave from the Department of Mathematics, University of Pisa, Via Buonarroti 2, 56127 Pisa, Italy, thanks to the G. Colombo fellowships of the European Space Agency.     

Jupiter: New retrieved thermal profiles and ammonia distributions

New thermal profiles of Jupiter are retrieved from recent far infrared spectral measurements and for H₂ mixing ratios varying from 0.8 to 0.94. The effective temperature corresponding to the inferred thermal profile is 123.15 ± 0.35°K. Far-infrared brightness temperature spectra computed from these profiles are compared to experimental data including measurements made at high spectral resolution in the NH₃ν₂ band at 10 μm and in NH₃ pure rotational bands between 40 and 110 μm. It is found that a strong depletion of NH₃ does occur in the Jovian stratosphere and that ammonia seems to be undersaturated in the upper troposphere.     

Problems in relating the optical and X-ray emissions from a solar flare

We evaluate the possibility that the short-lived Balmer line emission at H9 3835 Å of the 1972, August 2 (1839 UT) solar flare is due to heating of the chromosphere by bombarding electrons. We point out some of the problems of comparing the time behavior and spatial distribution of simultaneous hard and soft X-ray emissions. It is concluded that the present data do not justify the attribution of the short-lived optical emission to the presumed hard X-ray producing electrons.     

Study of geothermal water intrusion due to groundwater exploitation in the Puebla Valley aquifer system, Mexico

Significant intrusion of geothermal water into fresh groundwater takes place in the Puebla Valley aquifer system, Mexico. The decline in the potentiometric surface due to the overexploitation of the groundwater induces this intrusion. This hydrological system comprises three aquifers located in Plio-Quaternary volcanic sediments and Mesozoic calcareous rocks. The hydraulic balance of the aquifer shows that the annual output exceeds the natural inputs by 12 million m³. Between 1973 and 2002, a drop in the potentiometric surface, with an 80 m cone of depression, was identified in a 5-km-wide area located southwest of the city of Puebla. Chemical analyses performed on water samples since 1990 have shown an increase in total dissolved solids (TDS) of more than 500 mg/L, coinciding with the region showing a cone of depression in the potentiometric surface. A three-dimensional flow and transport model, based on the hydrogeological and geophysical studies, was computed by using the MODFLOW and MT3D software. This model reproduces the evolution of the aquifer system during the last 30 years and predicts for 2010 an additional drawdown in the potentiometric surface of 15 m, and an increase in the geothermal water intrusion.     

3D numerical modelling of fluid flow in the Val-d’Or orogenic gold district: major crustal shear zones drain fluids from overpressured vein fields

Fluid flow patterns have been determined using oxygen isotope isopleths in the Val-d’Or orogenic gold district. 3D numerical modelling of fluid flow and oxygen isotope exchange in the vein field shows that the fluid flow patterns can be reproduced if the lower boundary of the model is permeable, which represents middle or lower crustal rocks that are infiltrated by a metamorphic fluid generated at deeper levels. This boundary condition implies that the major crustal faults so conspicuous in vein fields do not act as the only major channel for upward fluid flow. The upper model boundary is impermeable except along the trace of major crustal faults where fluids are allowed to drain out of the vein field. This upper impermeable boundary condition represents a low-permeability layer in the crust that separates the overpressured fluid from the overlying hydrostatic fluid pressure regime. We propose that the role of major crustal faults in overpressured vein fields, independent of tectonic setting, is to drain hydrothermal fluids out of the vein field along a breach across an impermeable layer higher in the crust and above the vein field. This breach is crucial to allow flow out of the vein field and accumulation of metals in the fractures, and this breach has major implications for exploration for mineral resources. We propose that tectonic events that cause episodic metamorphic dehydration create a short-lived pulse of metamorphic fluid to rise along zones of transient permeability. This results in a fluid wave that propagates upward carrying metals to the mineralized area. Earthquakes along crustal shear zones cause dilation near jogs that draw fluids and deposit metals in an interconnected network of subsidiary shear zones. Fluid flow is arrested by an impermeable barrier separating the hydrostatic and lithostatic fluid pressure regimes. Fluids flow through the evolving and interconnected network of shear zones and by advection through the rock matrix. Episodic breaches in the impermeable barrier along the crustal shear zones allow fluid flow out of the vein field.