Inorganic nitrogen control in wastewater treatment ponds from a fish farm (Orbetello, Italy): denitrification versus Ulva uptake

The aim of this study was to quantify the N removal efficiency of an Ulva-based phytotreatment system receiving wastewaters from a land-based fish farm (Orbetello, Italy), to identify the main biogeochemical pathways involved and to provide basic guidelines for treatment implementation and management. Fluxes of O2 and nutrients in bare and in Ulva colonised sediments were assessed by light/dark core incubations; denitrification by the isotope pairing technique and Ulva growth by in situ incubation of macroalgal disks in cages. O2 and nutrient budgets were estimated as sum of individual processes and further verified by 24-h investigations of overall inlet and outlet loads. Ulva uptake (up to 7.8 mmol Nm(-2) h(-1)) represented a net sink for water column and regenerated NH4+ whilst N removal via denitrification (10-170 micromol Nm(-2) h(-1)) accounted for a small percentage of inorganic nitrogen load (<5%). Laboratory experiments demonstrated a high potential for denitrification (over 800 microM Nm(-2) h(-1)) indicating that N loss could be enhanced. The control of Ulva standing stocks by optimised harvesting of surplus biomass may represent an effective strategy to maximise DIN removal and could result in the assimilation of approximately 50% of produced inorganic nitrogen.     

Structural setting and tectonic evolution of the Apennine Units of northern Calabria

A new structural–stratigraphic synthesis of the Apennine units of northern Calabria is presented. The Meso-Cenozoic successions are grouped into two tectonic units, named Pollino–Ciagola Unit (PCU) and Lungro–Verbicaro Unit (LVU), comprising terrains formerly attributed to five different tectonic units. FeMg carpholite and blue amphibole record HP–LT metamorphism in the LVU, followed by progressive decompression leading to final greenschist facies re-equilibration during dominantly extensional deformation. Final tectonic emplacement of the LVU over the PCU post-dated the metamorphism of the LVU and was accompanied by intense ductile deformation along zones of strain localisation in footwall rocks. All of the units were later affected by folding and minor thrusting during subsequent Apennine tectonics. To cite this article: A. Iannace et al., C. R. Geoscience 337 (2005).     

Facing the challenge of shrinking cities in East Germany: The case of Leipzig

In the early 20^th century, the East German city of Leipzig seemed well on its way to become a metropolis of international importance. The city was expected to grow towards over one million inhabitants in 2000. Seventy years later, Leipzigs population has shrunk to less than 500,000 inhabitants instead. The German partition after World War II took away most of its national administrative and economic functions and much of its hinterland. The socialist GDR regime worsened the long-term development perspectives and living circumstances of the city. The German reunification brought new development chances, but like most East German cities, Leipzigs hopes soon became disappointed. The local politicians faced a difficult redevelopment task: apart from the question how to revive the local and regional economy, they also had to deal with a housing vacancy rate of 20%, a huge need for renovation in the older neighbourhoods as well as in the socialist high-rise areas, the negative effects of urban sprawl on the city core, and various environmental pollution problems. After briefly describing the development path of Leipzig until the 1990s, the paper will discuss the current attempts of the city government to give Leipzig a more positive post-industrial future. On the one hand, Leipzig is developing a strategy to downsize the citys built environment and infrastructure to adapt to a probably lastingly smaller population. On the other hand, many growth instruments well known from the international scientific and political debate are tried to put Leipzig back on the (inter)national map. The paper will discuss these development strategies in the light of the international debate on the question how to fight the shrinking city, with specific attention for post-socialist cities.     

Seismic activity along the central america volcanic arc: is it related to subduction of the Cocos plate?

We determine seismic strain rate of tectonic earthquakes along the Central America Volcanic Arc. We then compare this result to those obtained from earthquakes related to the convergence of the Cocos and Caribbean plates and to earthquakes in the back-arc region of northern Central America.

The seismic strain-rate tensor for shallow-focus earthquakes along the Central America volcanic arc since 1700, has a compressive eigenvector with a magnitude of 0.7 × 10⁻⁸ year⁻¹, and oriented in a 357° azimuth. The extensive eigenvector is oriented in a 86° azimuth, with a magnitude of 0.82 × 10⁻⁸ year⁻¹. When only Centroid Moment-tensor solutions (CMT) are considered, the respective eigenvectors are 1.2 × 10⁻⁸ year⁻¹ and 1.0 × 10⁻⁸ year⁻¹.

The compressive eigenvector from the seismic strain-rate tensor for earthquakes along the Cocos-Caribbean convergent margin is 2.0 × 10⁻⁸ year⁻¹, plunging at 25°, and oriented in a 29° azimuth. Its magnitude and direction are similar to those of the compressive eigenvector for earthquakes along the volcanic arc. The extensive eigenvector along the convergent margin, on the other hand, has a large vertical component. The compressive and extensive eigevenvectors are 4.9 × 10⁻⁸ year⁻¹ and 4.6 × 10⁻⁸ year⁻¹, using only CMTs as the database.

Earthquakes along the grabens of northern Central America yield a seismic strain-rate tensor whose extensive eigenvector has a magnitude of 2.4 × 10⁻⁸ year⁻¹, oriented in a 109° azimuth. Magnitude and direction are similar to those of the extensive eigenvector for earthquakes along the volcanic arc. The compressive eigenvector along the grabens is practically vertical.

Similarities in magnitudes and directions for compressive and extensive eigenvectors suggest to us that the strain field along the Central America volcanic arc is the result of compression along the convergent Cocos-Caribbean margin, and extension in the back-arc region, along the grabens of northern Central America. This field is resolved as strike-slip faulting along the arc.     

Regional Rainfall Depth–Duration–Frequency Equations for an Alpine Region

In this paper the scaling hypotheses are applied to annual maximum series of rainfall depth for different rainfall duration to derive the depth–duration–frequency (DDF) curve. It is shown that, based on the empirically observed scaling properties of rainfall and some general assumptions about the cumulative distribution function for the annual maximum of the rainfall depth, it is possible to derive a simple DDF relationship. This general framework provides a basis for the generation of maps that can be used to infer DDF curves at any point of a particular area. Data from a dense raingauge network in a mountainous region in north-eastern Italy (the Trentino Province) are used to clarify the methodology for the construction and regionalization of the DDF relationship. The geographical variation of short-duration (i.e., less than 60 xmin) rainfall extremes is also evaluated by using the same framework. It is found that depth–duration ratios, defined as the ratios of the t-min to the 60-min rainfall depths of the same return period, may be considered independent of return period and geographical location for any storm duration less than 60 min.     

High-resolution stratigraphy reveals repeated earthquake faulting in the Masada Fault Zone, Dead Sea Transform

A detailed study of the syndepositional Masada Fault Zone (MFZ) provides an example for fundamental characteristics of earthquakes, such as long term temporal clustering, repeated faulting on the same planes for a limited time of the order of a few thousands of years, and the formation of subaqueous breccia layers interpreted as seismites. The MFZ was studied in outcrops of 70–15 ka Lake Lisan sediments. Detailed columnar sections on both sides of well-exposed faults show that each individual fault exhibits a cluster, up to 4 ky long, with 3–5 slip events on the same plane. Each slip event is associated with the formation of widespread layers exhibiting soft sediment deformation, which are interpreted to be seismite layers. The uppermost part of the Lisan section, about 5 m, is not faulted, hence the last cluster of slip events ended about 25 ky ago. The clusters of activity of individual faults coalesce to form larger clusters. These are evident in the distribution of seismite layers throughout the entire Lisan section which shows earthquake clustering during periods of 10 ky. The clusters are separated by relatively quiescent periods of comparable duration.     

A physical analogue model to analyse interactions between tensile stresses and dissolution in carbonate slopes

A physical analogue model was developed to analyse the relationship between tensional states in rock masses, seepage and karst processes. Use was made of an experimental apparatus consisting of two hydraulic circuits realized by drilling two holes into each of two blocks of sampled limestone from the central Apennines (Italy). A static load (208.85 kg) was applied to one of the blocks in order to elicit tensile stresses within it. Physical and chemical monitoring data showed that the main process involved was temperature-dependent CaCO₃ dissolution. This process was more marked in the loaded block circuit, as Ca⁺⁺ concentration in circulated water reached 54 mg/L, whereas only 28 mg/L was reached in the unloaded one. The interaction between load and dissolution caused the observed opening of microcracks, as confirmed by further increase of water loss and by dilution in the loaded block circuit, resulting in a decrease of Ca⁺⁺ concentration. These findings were confirmed by recording additional water losses after increasing the load up to 445.05 kg. A finite difference numerical model showed that tensile stresses (max 20 kPa) within the loaded block were clustered at the intersection of the main joints with the flowpaths, thus representing points of preferential and accelerated dissolution.This study is part of a research project on analogical and numerical modelling for geological risk mitigation being conducted at CERI, Research Centre for Geological Risks of the University of Rome “La Sapienza” (Valmontone, Roma, Italy; http://www.ceri@uniroma1.it).     

Effective dispersion in a chemically heterogeneous medium under temporally fluctuating flow conditions

We investigate effective solute transport in a chemically heterogeneous medium subject to temporal fluctuations of the flow conditions. Focusing on spatial variations in the equilibrium adsorption properties, the corresponding fluctuating retardation factor is modeled as a stationary random space function. The temporal variability of the flow is represented by a stationary temporal random process. Solute spreading is quantified by effective dispersion coefficients, which are derived from the ensemble average of the second centered moments of the normalized solute distribution in a single disorder realization. Using first-order expansions in the variances of the respective random fields, we derive explicit compact expressions for the time behavior of the disorder induced contributions to the effective dispersion coefficients. Focusing on the contributions due to chemical heterogeneity and temporal fluctuations, we find enhanced transverse spreading characterized by a transverse effective dispersion coefficient that, in contrast to transport in steady flow fields, evolves to a disorder-induced macroscopic value (i.e., independent of local dispersion). At the same time, the asymptotic longitudinal dispersion coefficient can decrease. Under certain conditions the contribution to the longitudinal effective dispersion coefficient shows superdiffusive behavior, similar to that observed for transport in s stratified porous medium, before it decreases to its asymptotic value. The presented compact and easy to use expressions for the longitudinal and transverse effective dispersion coefficients can be used for the quantification of effective spreading and mixing in the context of the groundwater remediation based on hydraulic manipulation and for the effective modeling of reactive transport in heterogeneous media in general.     

Assessing the impact of large-scale dewatering on fault-controlled aquifer systems: a case study in the Acque Albule basin (Tivoli, central Italy)

The development of large-scale bedrock quarry operations often requires high-volume and long-term groundwater extraction to maintain a sustainable working environment. These dewatering activities often influence groundwater levels and flow patterns regionally. In the present study, the influence of the dewatering of the travertine quarry operations near the city of Tivoli, Italy, are quantitatively investigated through an integrated analysis of field data and numerical modeling. Lowering of regional groundwater levels in the vicinity of the quarry has led to destructive land subsidence and alterations to the flow system sustaining a hot-spring area. The study employs a finite element numerical model (FEFLOW) to evaluate and quantify the impact of the extensive dewatering on fault-controlled regional groundwater flow in the Acque Albule basin. By incorporating the physical field data and historical hydrologic information, the numerical model was calibrated against three groundwater scenarios, reproducing the effects of different exploitation activities, coupled with natural changes over the course of the quarry operation. The results indicate that groundwater withdrawals by the mining industry and by “Terme di Roma” spa resulted in the cessation of flow from the primary thermal spring and a drop in the phreatic level in the area consequently affected by land subsidence.     

The regularizing Levenberg–Marquardt scheme for history matching of petroleum reservoirs

In this paper, we study on a history matching approach that consists of finding stable approximations to the problem of minimizing the weighted least-squares functional that penalizes the misfit between the reservoir model predictions G(u) and noisy observations y ^η . In other words, we are interested in computing an approximation to the minimizer of $\frac {1}{2}\vert \vert \Gamma ^{-1/2}(y^{\eta }-G(u))\vert \vert _{Y}^{2} $ where Γ is the measurements error covariance, Y is the observation space, and X is a set of admissible parameters. This is an ill-posed nonlinear inverse problem that we address by means of the regularizing Levenberg–Marquardt scheme developed by Hanke (Inverse Probl. 13:79–95, 1997; J. Integr. Equ. Appl. 22(2):259–283, 2010). Under certain conditions on G, the theory of Hanke (Inverse Probl. 13:79–95, 1997; J. Integr. Equ. Appl. 22(2):259–283, 2010) ensures the convergence of the scheme to stable approximations to the inverse problem. We propose an implementation of the regularizing Levenberg–Marquardt scheme that enforces prior knowledge on the geologic properties. In particular, the prior mean $\overline {u}$ is incorporated in the initial guess of the algorithm, and the prior error covariance C is enforced through the definition of the parameter space X. Our main goal is to numerically show that the proposed implementation of the regularizing Levenberg–Marquardt scheme of Hanke is a robust method capable of providing accurate estimates of the geologic properties for small noise measurements. In addition, we provide numerical evidence of the convergence and regularizing results predicted by the theory of Hanke (Inverse Probl. 13:79–95, 1997; J. Integr. Equ. Appl. 22(2):259–283, 2010) for a prototypical oil–water reservoir model. The performance for recovering the true permeability with the regularizing Levenberg–Marquardt scheme is compared to the typical approach of computing the maximum a posteriori (MAP) estimator. In particular, we compare the proposed application of the regularizing Levenberg–Marquardt (LM) scheme against the standard LM approach of Li et al. (SPE J. 8(4):328–340, 2003) and Reynolds et al. (2008) for computing the MAP. Our numerical experiments suggest that the history matching approach based on iterative regularization is robust and could potentially be used to improve further on various methodologies already proposed as effective tools for history matching in petroleum reservoirs.