The Stress Field in the Upper Crust as Determined from In Situ Measurements

This paper provides a synopsis of the state of stress in the upper parts of the earth's crust based uponin situ rock stress determinations. Despite the large scatter of the data, two dominant trends can be detected in the variations of average horizontal stress with depth in various geological environments. Basement rocks in ancient shields and deformed rocks in fold belts usually show horizontal stresses larger than the theoretical overburden pressure. Sedimentary cover rocks and fissured massive rocks show horizontal stresses smaller than the overburden pressure. The ratio of the maximum to the minimum horizontal stress exhibits a clear stress anisotropy in most cases. Directions of maximum horizontal compression are fairly consistent in areas where sufficient measurements are available (North America and Fennoscandia), although in many instances they do not conform to any simple predicted stress pattern. Many factors complicate the interpretation ofin situ stress determinations. There is no simple relationship between the stress trajectories and the free surface. Topographic features and erosional processes may cause horizontal stress concentrations. Remanent stresses of great age can be superimposed on current tectonic stresses, while sometimes current stresses no longer coincide with the stress systems that caused observable faulting and folding. Observation of remanent stresses in ancient rocks shows that rocks in the upper crust have finite strength even under geological time intervals. From the viewpoint of global tectonics,in situ stress determinations ought to be used with great caution, and in conjunction with focal mechanism solutions of earthquakes. Although the state of stress is everywhere compressive, the fact that undeformed sedimentary cover rocks often show no excess horizontal stress would seem to indicate that no active global horizontal compression is required. Horizontal stresses larger than the overburden pressure in regions of intense palaeodeformation may be due to remanent stress effects and to the influence of the local structure. It is premature to advance any general statement on tectogenesis on the basis ofin situ stress determinations. More measurements, and a quantitative evaluation of the factors affecting them, are required before further progress can be made.     

Evolution of Eolian arc volcanism (Southern Tyrrhenian Sea)

The Eolian arc is located at the boundary between the converging African and European plates. Its volcanism is characterized by a marked evolution in a restricted time (less than 1 my). A progressive transition is observed from typical calc-alkaline series toward shoshonitic rocks produced by past and present activity (Vulcano and Stromboli). By comparison with circum-Pacific island arcs, the evolution of Eolian volcanism can be related to a rapid deepening of the Benioff zone. The occurrence of shoshonitic rocks and the continental nature of the crust on both sides of the plate boundary suggest that the Eolian arc is in a senile stage of evolution. Petrology and Sr isotopic data indicate a probable mantle source for Eolian volcanism.     

Multiple intrusions in the Cima d'Asta pluton (Northern Italy)

The systematic modal analysis of four sectors of the Cima d'Asta pluton (North Italy) has shown the existence of multiple adjacent intrusions.The main data supporting this interpretation are reported. Some problems concerning the mechanism of adjacent intrusions and the plutonic magmatic evolution are also discussed in a preliminary way.This paper is published with the help of a grant from the Consiglio Nazionale delle Ricerche, Italia.     

Relative Sea Level Variations Caused by Subduction

By means of a stratified Earth model with viscoelastic rheology, we have studied the long-term global fluctuations of Relative Sea Level (RSL) induced by subducting slabs. We have computed RSL variations for both a single subduction and a realistic distribution of slabs by a numerical simulation based on a simplified model of the subduction process. RSL is determined by the offset between the geoid and the dynamic topography; our analysis demonstrates that the latter provides the prevailing contribution. We have studied, in addition, the effects of rheological stratification upon the amplitude and time-evolution of these two quantities and, consequently, of RSL fluctuations. According to our results, an upper bound for the rate of RSL associated with subduction is of the order of 0.1 mm/yr, in agreement with previous studies. This rate of sea level variation is comparable with that attributed to changes in the tectonic regime on a large scale. This preliminary result corroborates the suggestion by other authors to include subduction in the list of geophysical mechanisms which contribute to long-term RSL fluctuations.     

Postglacial rebound and lateral viscosity variations: a semi-analytical approach based on a spherical model with Maxwell rheology

A spectral method is employed to study the response to surface loads of a Maxwell earth including lateral viscosity variations. In particular, we focus on the effects of lithospheric cratons on the long-wavelength time-dependent displacement field for simple earth models. The viscosity contrast of the craton with respect to the surrounding mantle is kept fixed, whereas its thickness is allowed to vary. We show that the long-wavelength vertical displacement is not greatly affected by the presence of a lithospheric craton, while the tangential displacement is severely modified for the case of a homogeneous mantle. With increasing harmonic degree and thickness of the craton, the load-deformation coefficients deviate from those pertaining to a homogeneous mantle with a viscosity of 10²¹ Pa s. These deviations are particularly enhanced on timescales larger than a few hundred years. These findings indicate that the interpretation of the viscosity structure of the mantle inferred from postglacial rebound signatures based on radially stratified models is affected by the presence of lateral viscosity variations.     

PCBs and PAHs in surficial sediments from aquatic environments of Mexico City and the coastal states of Sonora,Sinaloa, Oaxaca and Veracruz (Mexico)

Concentrations of polychlorinated biphenyls (PCBs) and polycyclic aromatic hydrocarbons (PAHs) were measured in 13 surficial sediment samples collected at three lacustrine locations in the surroundings of Mexico City and four coastal areas of the States of Sinaloa, Sonora, Oaxaca and Veracruz. Total PCB concentrations span the interval 3.18–621 ng g⁻¹. The highest values (63.7–621 ng g⁻¹) were found in Mexico City, which is a highly anthropogenically impacted area, whereas low concentrations (3.18–12.9 ng g⁻¹) were characteristic of seven places, some of them subject to intense hydrodynamics. In these latter cases, values increase by 18–73 times if normalised against the fine fraction (silt plus clay) content in sediment. Two samples from Mexico City exceed the ERM (Effect Range Median) guidelines and are likely to cause adverse effects. Samples contain only lower chlorinated PCBs (mainly 2-, 3- and 4-CB), thus suggesting that the most used PCB commercial mixture was Aroclor 1242. The homologue composition of the sample taken close to the nuclear power plant of Laguna Verde is identical to this commercial mixture. PAHs in the same samples have relatively low concentrations (14.9–287 ng g⁻¹), well below ERL (Effect Range Low) guidelines. The composition of PAH mixtures accounts for the influence of both petrogenic and pyrolitic sources, with these latter prevailing at some places in Mexico City.     

Turbulence and Air Exchange in a Two-Dimensional Urban Street Canyon Between Gable Roof Buildings

We experimentally investigate the effect of a typical building covering: the gable roof, on the flow and air exchange in urban canyons. In general, the morphology of the urban canopy is very varied and complex, depending on a large number of factors, such as building arrangement, or the morphology of the terrain. Therefore we focus on a simple, prototypal shape, the two-dimensional canyon, with the aim of elucidating some fundamental phenomena driving the street-canyon ventilation. Experiments are performed in a water channel, over an array of identical prismatic obstacles representing an idealized urban canopy. The aspect ratio, i.e. canyon-width to building-height ratio, ranges from 1 to 6. Gable roof buildings with 1:1 pitch are compared with flat roofed buildings. Velocity is measured using a particle-image-velocimetry technique with flow dynamics discussed in terms of mean flow and second- and third-order statistical moments of the velocity. The ventilation is interpreted by means of a simple well-mixed box model and the outflow rate and mean residence time are computed. Results show that gable roofs tend to delay the transition from the skimming-flow to the wake-interference regime and promote the development of a deeper and more turbulent roughness layer. The presence of a gable roof significantly increases the momentum flux, especially for high packing density. The air exchange is improved compared to the flat roof buildings, and the beneficial effect is more significant for narrow canyons. Accordingly, for unit aspect ratio gable roofs reduce the mean residence time by a factor of 0.37 compared to flat roofs, whereas the decrease is only by a factor of 0.9 at the largest aspect ratio. Data analysis indicates that, for flat roof buildings, the mean residence time increases by 30% when the aspect ratio is decreased from 6 to 2, whereas this parameter is only weakly dependent on aspect ratio in the case of gable roofs.     

Water-Channel Estimation of Eulerian and Lagrangian Time Scales of the Turbulence in Idealized Two-Dimensional Urban Canopies

Lagrangian and Eulerian statistics are obtained from a water-channel experiment of an idealized two-dimensional urban canopy flow in neutral conditions. The objective is to quantify the Eulerian \((T^{\mathrm{E}})\) and Lagrangian \((T^{\mathrm{L}})\) time scales of the turbulence above the canopy layer as well as to investigate their dependence on the aspect ratio of the canopy, AR, as the latter is the ratio of the width (W) to the height (H) of the canyon. Experiments are also conducted for the case of flat terrain, which can be thought of as equivalent to a classical one-directional shear flow. The values found for the Eulerian time scales on flat terrain are in agreement with previous numerical results found in the literature. It is found that both the streamwise and vertical components of the Lagrangian time scale, \(T_\mathrm{u}^\mathrm{L} \) and \(T_\mathrm{w}^\mathrm{L} \), follow Raupach’s linear law within the constant-flux layer. The same holds true for \(T_\mathrm{w}^\mathrm{L} \) in both the canopies analyzed \((AR= 1\) and \(AR= 2\)) and also for \(T_\mathrm{u}^\mathrm{L} \) when \(AR = 1\). In contrast, for \(AR = 2\), \(T_\mathrm{u}^\mathrm{L} \) follows Raupach’s law only above \(z=2H\). Below that level, \(T_\mathrm{u}^\mathrm{L} \) is nearly constant with height, showing at \(z=H\) a value approximately one order of magnitude greater than that found for \(AR = 1\). It is shown that the assumption usually adopted for flat terrain, that \(\beta =T^{\mathrm{L}}/T^{\mathrm{E}}\) is proportional to the inverse of the turbulence intensity, also holds true even for the canopy flow in the constant-flux layer. In particular, \(\gamma /i_\mathrm{u} \) fits well \(\beta _\mathrm{u} =T_\mathrm{u}^\mathrm{L} /T_\mathrm{u}^\mathrm{E} \) in both the configurations by choosing \(\gamma \) to be 0.35 (here, \(i_\mathrm{u} =\sigma _\mathrm{u} / \bar{u} \), where \(\bar{u} \) and \(\sigma _\mathrm{u} \) are the mean and the root-mean-square of the streamwise velocity component, respectively). On the other hand, \(\beta _\mathrm{w} =T_\mathrm{w}^\mathrm{L} /T_\mathrm{w}^\mathrm{E} \) follows approximately \(\gamma /i_\mathrm{w} =0.65/\left( {\sigma _\mathrm{w} /\bar{u} } \right) \) for \(z > 2H\), irrespective of the AR value. The second main objective is to estimate other parameters of interest in dispersion studies, such as the eddy diffusivity of momentum \((K_\mathrm{{T}})\) and the Kolmogorov constant \((C_0)\). It is found that \(C_0\) depends appreciably on the velocity component both for the flat terrain and canopy flow, even though for the latter case it is insensitive to AR values. In all the three experimental configurations analyzed here, \(K_\mathrm{{T}}\) shows an overall linear growth with height in agreement with the linear trend predicted by Prandtl’s theory.     

Suspended sediment transport analysis in two Italian instrumented catchments

Suspended sediment transport in streams is an effective indicator of soil erosion at the watershed scale. This process was studied using the data obtained from two continuous monitoring stations installed in Italian watersheds, the Rio Cordon and Torrente Carapelle. The catchments have substantially different climatic, morphological and land use characteristics. The Rio Cordon, a small Alpine watershed, has been monitored since 1986, while in the Torrente Carapelle, a medium-size Mediterranean watershed, the monitoring station has provided data since 2007. Several flood events with suspended sediment transport were isolated in the two catchments, excluding those determined by snowmelt in the Rio Cordon as this process does not affect the Carapelle watershed. Analysis of the events showed different behavior in terms of hysteresis loop trends between water discharge (m³ s^-1) and suspended sediment concentration (g l^-1) values, as the Rio Cordon confirms clockwise relationships most frequently, while counter-clockwise trends represent the majority of cases for the Carapelle Torrent. The different behavior of the two watersheds was further analyzed by evaluating relationships between the main hydrological parameters related to suspended sediment transport. Runoff controls the total sediment budget in both catchments. In contrast, it was noted that the runoff process does not interfere in the magnitude of the suspended sediment transport in the Rio Cordon catchment, while variations due to the larger size of the catchment area characterize the events in the Carapelle watershed. Lastly, a flow peak threshold that determines the advance or delay of the hydrograph peak with respect to the sedigraph peak was established for the Carapelle Torrent, while it was not reported in the Rio Cordon.     

A new method for the interpretation of the constant-head well permeameter

A novel semi-analytical solution for the interpretation of the constant-head permeameter test is introduced, which accounts for the correct mixed-type boundary condition at the wellbore, unlike all published analytical solutions. Capillarity can also be accounted for. The simplifications are that flow from the bottom of the borehole is neglected (therefore the solution is applicable to slender boreholes, where the ponding depth is at least 10 times the radius) and capillarity can be modeled with a quasi-linear approach. The Green's function approach leads to an integral equation, the solution of which does not show significant ill-posedness. Two sub-cases are presented: the first neglects capillary effects (the all-saturated approximation) and the second (general solution) takes them into account. The all-saturated solution is successfully tested against finite element simulations. The corresponding values of the borehole shape factor C are slightly larger than the ones obtained with approximate analytical solutions from the literature. When capillarity is accounted for, C changes of a factor of 10 when the dimensionless sorptive number A goes from typical values for fine soils to typical values for coarse soils (about two orders of magnitude of variation for A). This range shifts to lower values of A as the dimensionless borehole depth increases. Consequently, the all-saturated solution is a good approximation of the soil behavior for boreholes with large ponding depth, and coarse soils. The proposed semi-analytical solution is fast to compute and thus it is possible to use it in an automated optimization technique to fit field data and estimate the field-saturated hydraulic conductivity and the sorptive number; this would not be feasible using a numerical solution.