Petrogenesis of the mafic igneous rocks of the Betic Cordilleras: A field, petrological and geochemical study

Mafic igneous rocks are widespread in the Nevado-Filábride Complex, the lowermost metamorphic unit of the internal zones of the Betic Cordilleras. They form intrusive, small, discontinuous bodies, predominantly dikes with subordinate small lava flows. The entire complex underwent alpine compressional metamorphism during the Paleogene continental collision, resulting in eclogites and blueschists in the mafic bodies and high-pressure assemblages in the intruded metasediments. Locally, weakly metamorphosed or completely unmetamorphosed igneous rocks with the same textural features occur as patches surrounded by eclogitized igneous rocks. The bulk rock chemistry of unmetamorphosed and completely metamorphosed mafic igneous rocks is consistent with an alkaline to transitional tholeiitic magmatism with typical within-plate geochemical characteristics. All but a few samples are nepheline normative and display REE and trace element characteristics typical of continental, rift-related magmatism. This conclusion is strongly supported by the mineral chemistry of the major constituents, in particular the calcic Ti-rich character of clinopyroxene, the lack of orthopyroxene, and the occurrence of kaersutitic amphibole. Incompatible trace element abundances and Sr and Nd isotopes support the provenance of these magmas from a variably metasomatized previously depleted sub-continental lithospheric mantle source. Received: 5 July 1999 / Accepted: 28 February 2000     

First study of fault trench stratigraphy at Mt. Etna volcano, Southern Italy: understanding Holocene surface faulting along the Moscarello fault

Paleoseismology, the study of past earthquakes based on their geological record in the stratigraphy and landscape, is a successful newly developing field of research. The application of fault trench studies in volcanic environments is one of the youngest branches of paleoseismology. In this paper, we present the results of the first exploratory trenches excavated at Mt. Etna in Sicily, the largest European volcano. Modern surface faulting at Etna is a very well known feature, which poses significant hazard to the local community, both in terms of ground displacement of essential lifelines and ground shacking from frequent damaging earthquakes. However, while the geomorphology and the seismicity of the active fault in the Etna region consistently show very high rates of tectonic activity, the Holocene cumulative throw and slip-rates, along with the nature (coseismic vs. creeping fault slip), dimension and timing of the displacement events, are still poorly constrained. For this purpose, we selected as a sample area the Moscarello fault, one of the most outstanding segment of the Timpe system of active normal faults in the volcano’s lower eastern flank. Displaced landforms and volcanic units at the Fondo Macchia basin, in the central sector of this fault, indicate some hundreds of meters of vertical offset in the last ca. 80 kyr, with a long-term slip-rate substantially higher than 1.5–2.0 mm/yr. According to the historical sources and instrumental observations, the Moscarello fault ruptured four times in the last 150 years during shallow (H < 5 km) and moderate magnitude (M < 4.8) earthquakes. These events were associated with severe damage in a narrow epicentral area (macroseismic intensities up to the IX–X grade of the MSK scale) and extensive surface faulting (end-to-end rupture length up to 6 km, vertical offsets up to 90 cm). This clearly indicates very high modern rates of deformation along this fault. We conducted trench investigations at the Fondo Macchia site, in a point where eyewitnesses observed ca. 20 cm of coseismic vertical displacement after the April 21, 1971, M_s=3.7, earthquake. The excavated sections provided direct stratigraphic evidence for a vertical slip-rate of 1.4–2.7 mm/yr in the last ca. 6 kyr. This should be regarded as a minimum slip-rate for the central section of the fault. We explored a single scarp at a single site, while we know from recent historical observations that several parallel scarps may rupture coseismically at Fondo Macchia. Thus, the relevant deformation rate documented for the modern period might be likely extended back in the past to a time-span of some thousands of years at least. As expected, for such a volcanic environment, the activity rates of the Moscarello fault are also significantly higher than for the Apennines normal faults, typically showing slip-rates lower than 1 mm/yr. The agriculturally reworked trench hangingwall stratigraphy did not allow to recognize individual displacement events. Nevertheless, the sedimentary structures observed in the trench footwall strongly suggest that, as for the last 150–200 years of detailed historical record, fault behavior at Fondo Macchia is governed by coseismic surface displacement rather than fault creep. This research confirms that paleoseismology techniques can be effectively applied also in active volcanic environments, typically characterized by rheology and, consequently, seismicity and fault dynamics very different from those of other tectonic environments in which paleoseismology has been firstly developed and is today extensively applied.     

Geomorphological and sedimentological features in Quaternary fluvial systems affected by solution-induced subsidence (Ebro Basin, NE-Spain)

The Quaternary evolution and the morpho-sedimentary features of some of the most important rivers in Spain (Ebro and Tagus rivers among others) have been controlled by subsidence due to alluvial karstification of the evaporitic bedrock. The subsidence mechanism may range from catastrophic collapse to slow sagging of the alluvium by passive bending. In the Ebro Basin, the mechanisms and processes involved in karstic subsidence were studied through the analysis of present-day closed depressions as well as through old subsidence depressions (palaeocollapses and solution-induced basins) and associated deformations recorded in the Quaternary alluvial sediments. The Gállego–Ebro river system is presented as a case study of channel adjustments and geomorphic and sedimentary evolution of fluvial systems in dissolution-induced subsidence areas. In this fluvial system, evaporite dissolution during particular Quaternary time intervals (namely early and middle Pleistocene) have lead to the development of a solution-induced basin, approximately 30 km-long by 8 km-wide, filled by Quaternary deposits with a total thickness in excess of 190 m. The main river response to balance the subsidence in the alluvial plain was aggradation in the central reach of the subsiding area, and degradation both in the upstream reach and in the valley sides where alluvial fans and covered pediments may prograde over the fluvial sediments. The main sinking areas are recognized in the sedimentary record by anomalous thickenings in the alluvial deposits and fine-grained sediments deposited in backswamp and ponded areas.     

Microseismicity and faulting geometry in the Gulf of Corinth (Greece)

During the summer of 1993, a network of seismological stations was installed over a period of 7 weeks around the eastern Gulf of Corinth where a sequence of strong earthquakes occurred during 1981. Seismicity lies between the Alepohori fault dipping north and the Kaparelli fault dipping south and is related to both of these antithetic faults. Focal mechanisms show normal faulting with the active fault plane dipping at about 45° for both faults. The aftershocks of the 1981 earthquake sequence recorded by King et al . (1985 ) were processed again and show similar results. In contrast, the observations collected near the western end of the Gulf of Corinth during an experiment conducted in 1991 ( Rigo et al . 1996 ), and during the aftershock studies of the 1992 Galaxidi and the 1995 Aigion earthquakes ( Hatzfeld et al . 1996 ; Bernard et al . 1997 ) show seismicity dipping at a very low angle (about 15°) northwards and normal faulting mechanisms with the active fault plane dipping northwards at about 30°. We suggest that the 8–12 km deep seismicity in the west is probably related to the seismic–aseismic transition and not to a possible almost horizontal active fault dipping north as previously proposed. The difference in the seismicity and focal mechanisms between east and west of the Gulf could be related to the difference in the recent extension rate between the western Gulf of Corinth and the eastern Gulf of Corinth, which rotated the faults dipping originally at 45° (as in the east of the Gulf) to 30° (as in the west of the Gulf).     

An Attempted Explanation for the Discrepancy between Internal and External Errors in Statistical Adjustments

It is shown that underparameterization of models used in statistical adjustments is at least in part responsible for the discrepancies between so-called internal and external rms errors (standard deviations), frequently encountered in derived sets of statistical parameters.     

Platinum-group element distribution in chromite ores from ophiolite complexes: implications for their exploration

Compilation of some new data on ophiolites for Greece and Yugoslavia, and published data from previous studies, indicate that platinum-group element (PGE) and gold concentrations in chromite ores are generally low, ranging from less than 100 ppb to a few hundred ppb. However, samples from several ophiolite complexes exhibit an enrichment (of a few ppm) (a) only in Os, Ir and Ru,(b) only in Pt and/or Pd or (c) in all PGE. This enrichment (up to 10s ppm) is mainly related with chromitites hosted in supra-Moho dunites and dunites of the uppermost stratigraphic levels of the mantle sequence and it seems to be local, independent of the chromitite major element composition and the chromite potential of the ophiolite complexes. The contents of PGE combined with less chalcophile elements (Ni, Co, Cu), the ratios of incompatible/compatible elements, and PGE-patterns provide evidence for discrimination between chromitites derived from primitive magmas and those derived from partially fractionated magmas, although they have a similar major element composition. Thus, they can be used for a stratigraphic orientation in the mantle sequence, and therefore for exploration targets. Moreover, PGE data offer valuable information for the evaluation of the chromite potential in ophiolite complexes. The most promising ophiolites seem to be those which apart from the petrological and geochemical characteristics indicating extensive degree of partial melting in the mantle source contain only one chromite ore type (the other type being only in small proportion) of limited compositional variation, in both major elements and PGE, low ratios of , while PGE-enriched chromitites in the mantle sequence are only occasionally present. In contrast, ophiolites which contain both high-Cr and -Al chromitites, and where their chalcophile element data implies relatively extensive fractionation trend are not good exploration targets for chromite ores, although they are related with a SSZ environment.     

The Role of Planetary Boundary-Layer Parameterizations in the Air Quality of an Urban Area with Complex Topography

The effect of different planetary boundary-layer (PBL) parameterization schemes on the spatial distribution of atmospheric pollution over the complex topography of the greater Athens area is investigated. Four PBL schemes originally implemented in a numerical meteorological model and a fifth one simulating the urban effect are examined. Two different atmospheric conditions are analyzed; a typical summer and a typical winter pollution episode. The relative importance of chemical and physical processes of the pollution predictions is discussed using process analysis. It is revealed that, for primary pollutants, a local scheme seems more adequate to represent the maximum observed concentrations while, completely different in structure, a non-local scheme reproduces the mean observed values in the basin. Concerning secondary pollutants, peak concentration differences, due to the different PBL schemes, are smoothed out. Nevertheless, the PBL scheme selection shapes the horizontal and the vertical extension of maximum values. The non-local and semi non-local schemes are superior to the others, favouring strong vertical mixing and transport towards the surface. The stronger turbulence accommodated effectively by the semi non-local urban scheme enhances ozone production along the sea-breeze axis and preserves the high ozone concentrations during the nighttime hours in the urban core.     

Analysis of Sediments and Soils by X-Ray Fluorescence Spectrometry Using Matrix Corrections Based on Fundamental Parameters

Wavelength dispersive X-ray fluorescence spectrometry (WD-XRF) is widely used for the analysis of soils and sediments using well characterised procedures. However, difficulties can occur with samples such as unknowns containing small amounts of ore materials and samples collected from contaminated sites where trace elemental concentrations can exceed the concentration range for routine analysis. We studied the performance of a commercially available method, based on fundamental parameters (FP) to correct matrix effects. The spectrometer was originally calibrated with elemental or simple compound calibrants. Samples were analysed as pressed powder pellets. Eighteen sediment and soil reference materials, three of them with certified values for some of their constituents, were used to evaluate accuracy, by comparing results with recommended values and their standard deviations (RV ± 2s) or certified values and their confidence intervals (CV ± Cl). When results fell systematically outside these intervals, calibrations were refined with geochemical reference materials. The best agreement of results with recommended and certified values was obtained when the contents of H₂O and C in each sample were included as matrix constituents during calculations. The detection limits of trace elements tended to be relatively high, because the measuring conditions employed were not maximised for sensitivity. The main advantage of the method tested was that it enabled the analysis of samples with high concentrations of trace elements and the determination of elements such as F, Bi, Sb and W, which are not commonly included in quantitative XRF analysis of geological samples.     

An Intercomparison of Large-Eddy Simulations of the Stable Boundary Layer

Results are presented from the first intercomparison of large-eddy simulation (LES) models for the stable boundary layer (SBL), as part of the Global Energy and Water Cycle Experiment Atmospheric Boundary Layer Study initiative. A moderately stable case is used, based on Arctic observations. All models produce successful simulations, in as much as they generate resolved turbulence and reflect many of the results from local scaling theory and observations. Simulations performed at 1-m and 2-m resolution show only small changes in the mean profiles compared to coarser resolutions. Also, sensitivity to subgrid models for individual models highlights their importance in SBL simulation at moderate resolution (6.25 m). Stability functions are derived from the LES using typical mixing lengths used in numerical weather prediction (NWP) and climate models. The functions have smaller values than those used in NWP. There is also support for the use of K-profile similarity in parametrizations. Thus, the results provide improved understanding and motivate future developments of the parametrization of the SBL.     

Future regional projections of extreme temperatures in Europe: a nonstationary seasonal approach

This paper analyzes changes of maximum temperatures in Europe, which are evaluated using two state-of-the-art regional climate models from the EU ENSEMBLES project. Extremes are expressed in terms of return values using a time-dependent generalized extreme value (GEV) model fitted to monthly maxima. Unlike the standard GEV method, this approach allows analyzing return periods at different time scales (monthly, seasonal, annual, etc). The study focuses on the end of the 20th century (1961?C2000), used as a calibration/validation period, and assesses the changes projected for the period 2061?C2100 considering the A1B emission scenario. The performance of the regional models is evaluated for each season of the calibration period against the high-resolution gridded E-OBS dataset, showing a similar South-North gradient with larger values over the Mediterranean basin. The inter-RCM changes in the bias pattern with respect to the E-OBS are larger than the bias resulting from a change in the boundary conditions from ERA-40 to ECHAM5 20c3m. The maximum temperature response to increased green house gases, as projected by the A1B scenario, is consistent for both RCMs. Under that scenario, results indicate that the increments for extremes (e.g. 40-year return values) will be two or three times higher than those for the mean seasonal temperatures, particularly during Spring and Summer in Southern Europe.