Lower lithospheric structure beneath the Trans-European Suture Zone from POLONAISE''97 seismic profiles

The large-scale POLONAISE'97 seismic experiment investigated the velocity structure of the lithosphere in the Trans-European Suture Zone (TESZ) region between the Precambrian East European Craton (EEC) and Palaeozoic Platform (PP). In the area of the Polish Basin, the P-wave velocity is very low (Vp <6.1 km/s) down to depths of 15–20 km, and the consolidated basement (Vp5.7–5.8 km/s) is 5–12 km deep. The thickness of the crust is 30 km beneath the Palaeozoic Platform, 40–45 km beneath the TESZ, and 40–50 km beneath the EEC. The compressional wave velocity of the sub-Moho mantle is >8.25 km/s in the Palaeozoic Platform and 8.1 km/s in the Precambrian Platform. Good quality record sections were obtained to the longest offsets of about 600 km from the shot points, with clear first arrivals and later phases of waves reflected/refracted in the lower lithosphere. Two-dimensional interpretation of the reversed system of travel times constrains a series of reflectors in the depth range of 50–90 km. A seismic reflector appears as a general feature at around 10 km depth below Moho in the area, independent of the actual depth to the Moho and sub-Moho seismic velocity. “Ringing reflections” are explained by relatively small-scale heterogeneities beneath the depth interval from 90 to 110 km. Qualitative interpretation of the observed wave field shows a differentiation of the reflectivity in the lower lithosphere. The seismic reflectivity of the uppermost mantle is stronger beneath the Palaeozoic Platform and TESZ than the East European Platform. The deepest interpreted seismic reflector with zone of high reflectivity may mark a change in upper mantle structure from an upper zone characterised by seismic scatterers of small vertical dimension to a lower zone with vertically larger seismic scatterers, possible caused by inclusions of partial melt.     

Equilibria in the system Au–Ag–S–fluid: computed and experimental data

We present results of computations on the interaction of solid-phase electrum–argentite–pyrite (weight ratios 210⁻⁵/ 210⁻³/1 and 210⁻⁵/410⁻²/1) association with Cl-containing aqueous moderately acid solutions (0.5m NaCl, pH = 3.08) at 300 °C and 500 bars. These data are a physicochemical basis for predicting the geochemical behavior of Au and Ag during the hydrothermal-metasomatic transformation of Au-Ag-pyrite. We also propose a technique of study of this process based on the phase equilibria of the subsystem Au–Ag–S with the aqueous solution at different liquid/solid (l/s) ratios, with the use of new graphic diagrams. The relationship of the composition of the solid-phase association with l/s ratio in real boundary conditions (Au = 17 ppm, mAu/mAg = 10^–3.57–10^–2.28) is shown. The maximum l/s values for complete leaching of gold and silver (l/smax = 200–800) are estimated. It has been established that argentite is the first to dissolve when mAu/mAg(s) > mAu/mAg(sol), and electrum, when mAu/mAg(s) < mAu/mAg(sol).

The experimental results showed that at 300 °C, the conversion of electrum (N_Au = 300‰) nonequilibrated with pyrite into an Au-richer form (N_Au = 730‰) and argentite follows an intricate kinetic scheme. Using the Pilling-Bedwords kinetic equation for processing data yielded the process rate constant K = 2.8(±0.5)10⁻⁵ g²cm⁻⁴day⁻¹. With this equation, the time of the complete conversion of 200 μm thick flat gold grains is 604 days. These data evidence a significant role of kinetic factors in hydrothermal-metasomatic processes involving native gold, which requires combination of thermodynamic and kinetic approaches on the construction of geologo-genetic models for hydrothermal sulfide formation.     

Main features of crustal structure in western and southern Europe based on data of explosion seismology

During the past two decades deep seismic sounding measurements have been carried out in western and southern Europe, mainly using the refraction method. These investigations were performed partly on a national basis but as well within international cooperative programs under the sponsorship of the European Seismological Commission.

In France, a systematic study has been executed to determine the main feature of deep structures under the Central Massif and the Paris Basin. In the Forez and Margeride regions, the sub-crustal velocity is lower (7.2 km/sec) than the normal value (8.0 km/sec) observed in the adjacent areas.

The central and southern part of Western Germany is covered by an extensive network of refraction profiles. The crustal thickness varies, similarly to France, from 25 to 35 km. A great amount of deep reflection data was obtained by commercial and special reflection work. The crust beneath the Rhinegraben area shows the typical “rift system” structure with a low subcrustal velocity (7.4–7.7 km/sec).

Very intensive refraction work has been carried out in the Alpine area. The maximum crustal thickness found near the axis of the negative gravity anomaly is about 55–60 km. Furthermore, a clear lowvelocity layer at a depth between 10 and 30 km has been detected. A key position with regard to the geotectonic structure of the Alps is held by the zone of Ivrea characterized by a pronounced gravity high. From the refraction work it may be concluded that there material of the lower crust and the upper mantle (7.2–7.5 km/sec) is overlying a layer of extremely low velocity (5.0 km/sec) which is interpreted as sialic crust.

Three years ago, a systematic study of crustal structure of the Italian peninsula has been started. Reversed profiles were observed on Sicily, in Calabria, and in Puglia. On Sicily, the structure is very complicated; the crust of the western part looks like a transition between a continental and oceanic structure whereas the eastern side shows a continental-type crust. In Calabria and Puglia, the crustal thickness has been determined to be about 25–35 km.     

VRANCEA99—the crustal structure beneath the southeastern Carpathians and the Moesian Platform from a seismic refraction profile in Romania

The VRANCEA99 seismic refraction experiment is part of an international and multidisciplinary project to study the intermediate depth earthquakes of the Eastern Carpathians in Romania. As part of the seismic experiment, a 300-km-long refraction profile was recorded between the cities of Bacau and Bucharest, traversing the Vrancea epicentral region in NNE–SSW direction.

The results deduced using forward and inverse ray trace modelling indicate a multi-layered crust. The sedimentary succession comprises two to four seismic layers of variable thickness and with velocities ranging from 2.0 to 5.8 km/s. The seismic basement coincides with a velocity step up to 5.9 km/s. Velocities in the upper crystalline crust are 5.9–6.2 km/s. An intra-crustal discontinuity at 18–31 km divides the crust into an upper and a lower layer. Velocities within the lower crust are 6.7–7.0 km/s. Strong wide-angle PmP reflections indicate the existence of a first-order Moho at a depth of 30 km near the southern end of the line and 41 km near the centre. Constraints on upper mantle seismic velocities (7.9 km/s) are provided by Pn arrival times from two shot points only. Within the upper mantle a low velocity zone is interpreted. Travel times of a PLP reflection define the bottom of this low velocity layer at a depth of 55 km. The velocity beneath this interface must be at least 8.5 km/s.

Geologic interpretation of the seismic data suggests that the Neogene tectonic convergence of the Eastern Carpathians resulted in thin-skinned shortening of the sedimentary cover and in thick-skinned shortening in the crystalline crust. On the autochthonous cover of the Moesian platform several blocks can be recognised which are characterised by different lithological compositions. This could indicate a pre-structuring of the platform at Mesozoic and/or Palaeozoic times with a probable active involvement of the Intramoesian and the Capidava–Ovidiu faults. Especially the Intramoesian fault is clearly recognisable on the refraction line. No clear indications of the important Trotus fault in the north of the profile could be found. In the central part of the seismic line a thinned lower crust and the low velocity zone in the uppermost mantle point to the possibility of crustal delamination and partial melting in the upper mantle.     

Velocity structure of the continental upper mantle: evidence from southern Africa

The velocity model for southern Africa of Qiu et al. [Qiu, X., Priestley, K., McKenzie, D., 1996. Average lithospheric structure of southern Africa. Geophys. J. Int. 127, 563–587] is revised so as to satisfy both the regional seismic waveform data and the fundamental mode Rayleigh wave phase velocity data for the region. The revised S-wave model is similar to the original model of Qiu et al. except that the high velocity, upper mantle lid extends to 160 km depth in the revised model rather than to 120 km in the original model. Sensitivity tests of the regional seismic data show that the minimum velocity in the S-wave low velocity zone can be as high as 4.45 km s⁻¹ compared to 4.32 km s⁻¹ in the Qiu et al. model. The vertical S-wave travel time for the revised south African model is compared with the vertical S-wave travel times for the global tomographic models S12WM13 and S16B30, and they are found to be similar.     

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.     

Crustal structure across the TESZ along POLONAISE''97 seismic profile P2 in NW Poland

The POLONAISE'97 (POlish Lithospheric ONset—An International Seismic Experiment, 1997) seismic experiment in Poland targeted the deep structure of the Trans-European Suture Zone (TESZ) and the complex series of upper crustal features around the Polish Basin. One of the seismic profiles was the 300-km-long profile P2 in northwestern Poland across the TESZ. Results of 2D modelling show that the crustal thickness varies considerably along the profile: 29 km below the Palaeozoic Platform; 35–47 km at the crustal keel at the Teisseyre–Tornquist Zone (TTZ), slightly displaced to the northeast of the geologic inversion zone; and 42 km below the Precambrian Craton. In the Polish Basin and further to the south, the depth down to the consolidated basement is 6–14 km, as characterised by a velocity of 5.8–5.9 km/s. The low basement velocities, less than 6.0 km/s, extend to a depth of 16–22 km. In the middle crust, with a thickness of ca. 4–14 km, the velocity changes from 6.2 km/s in the southwestern to 6.8 km/s in the northeastern parts of the profile. The lower crust also differs between the southwestern and northeastern parts of the profile: from 8 km thickness, with a velocity of 6.8–7.0 km/s at a depth of 22 km, to ca.12 km thickness with a velocity of 7.0–7.2 km/s at a depth of 30 km. In the lowermost crust, a body with a velocity of 7.20–7.25 km/s was found above Moho at a depth of 33–45 km in the central part of the profile. Sub-Moho velocities are 8.2–8.3 km/s beneath the Palaeozoic Platform and TTZ, and about 8.1 km/s beneath the Precambrian Platform. Seismic reflectors in the upper mantle were interpreted at 45-km depth beneath the Palaeozoic Platform and 55-km depth beneath the TTZ.

The Polish Basin is an up to 14-km-thick asymmetric graben feature. The basement beneath the Palaeozoic Platform in the southwest is similar to other areas that were subject to Caledonian deformation (Avalonia) such that the Variscan basement has only been imaged at a shallow depth along the profile. At northeastern end of the profile, the velocity structure is comparable to the crustal structure found in other portions of the East European Craton (EEC). The crustal keel may be related to the geologic inversion processes or to magmatic underplating during the Carboniferous–Permian extension and volcanic activity.     

Structure of the crust and upper mantle beneath the central european rift system

The crustal structure of the central European rift system has been investigated by seismic methods with varying success. Only a few investigations deal with the upper-mantle structure. Beneath the Rhinegraben the Moho is elevated, with a minimum depth of 25 km. Below the flanks it is a first-order discontinuity, while within the graben it is replaced by a transition zone with the strongest velocity gradient at 20–22 km depth. An anomalously high velocity of up to 8.6 km/s seems to exist within the underlying upper mantle at 40–50 km depth. A similar structure is also found beneath the Limagnegraben and the young volcanic zones within the Massif Central of France, but the velocity within the upper mantle at 40–50 km depth seems to be slightly lower. Here, the total crustal thickness reaches only 25 km. The crystalline crust becomes extremely thin beneath the southern Rhônegraben, where the sediments reach a thickness of about 10 km while the Moho is found at 24 km depth. The pronounced crustal thinning does not continue along the entire graben system. North of the Rhinegraben in particular the typical graben structure is interrupted by the Rhenohercynian zone with a “normal” West-European crust of 30 km thickness evident beneath the north-trending Hessische Senke. A single-ended profile again indicates a graben-like crustal structure west of the Leinegraben north of the Rhenohercynian zone. No details are available for the North German Plain where the central European rift system disappears beneath a sedimentary sequence of more than 10 km thickness.     

Deep Crustal Electrical Signatures of Eastern Dharwar Craton, India

Wide band magnetotelluric (MT) investigations were carried out along a profile from Kavali in the east to Anantapur towards west across the Eastern Ghat Granulite Terrain (EGGT), Eastern Dhanvar Craton (EDC) and a Proterozoic Cuddapah Basin. This 300 km long profile was covered with 20 stations at an interval of 12–18 km. The MT data is subjected to robust processing, decomposition and static shift correction before deriving a 2-D model. The model shows a resistive crust (−10,000–30,000 ohm-m) to a depth of 8–10 km towards west of the Cuddapah basin. The mid crust is less resistive (about 500 ohm-m) and the lower crust with a slight increase in resistivity (about 1,500 ohm-m) in the depth range of 20–22 km. The resistivity picture to the east of the Cuddapah basin also showed a different deep crustal structure. The resistivity of upper crust is about 5,000 ohm-m and about 200 ohm-m for mid and lower crust. The sediment resistivity of Cuddapah basin is of the order of 15–20 ohm-m. MT model has shown good correlation with results from other geophysical studies like deep seismic sounding (DSS), gravity and magnetics. The results indicate that the lower crustal layers are of intermediate type showing hydrous composition in Eastern Dhanvar Craton.     

Crustal structure beneath the Faroe Islands and the Faroe–Iceland Ridge

We have mapped the transition from the continental Faroe block (the Faroe Islands and surrounding shelf) to the thickened oceanic crust of the Faroe–Iceland Ridge in the North Atlantic using the results of a detailed sea-to-land seismic profile with wide-angle to normal-incidence recordings of explosive and airgun shots fired at sea along the Faroe–Iceland Ridge. Interpretation of all available seismic and gravity data indicates that this aseismic ridge is composed of 30±3-km-thick oceanic crust, with a gradual transition to ancient continental crust from 100 to 40 km northwest of the Faroe Islands, close to the shelf edge. This confirms that the crust beneath the Faroe Islands, which may be up to 46 km thick, comprises continental material in agreement with previous seismic and geochemical results. Results suggest that the upper 5.2±0.7 km of the Faroe crust consists of Tertiary basalts generated during continental breakup, overlying the continental crust beneath. The lower crust, where seismic constraint is poor, may exhibit high seismic velocities (7.1–7.6 km s⁻¹) which we attribute to underplating or intrusion by mafic melts during continental breakup in the early Tertiary.     

Seismic, gravity and petrological evidence for partial melt beneath the thickened Central Andean crust (21–23°S)

On the basis of seismic refraction investigations and gravimetric data we have modelled the crustal structure of the southern Central Andes (21–23°S). A pronounced variation in crustal parameters is seen in N-S- and W-E-crossing seismic profiles over the entire Andean orogene, characterized by a crustal thickness of up to 70 km under the magmatic arc and backarc, strongly reduced seismic velocities and a Bouguer minimum of −450 mGal. Anomalously low velocities of 5.9–6.0 km/s in the deeper crust of the Western Cordillera and Altiplano regions lead to an over-compensation of the Bouguer minima resulting in values of crustal densities higher than estimates based purely on seismic velocity measurements. In an attempt to reconcile these differences, the behavior of crystalline rocks based on published laboratory data was studied under varying pressure and temperature conditions up to the range of partial melting. If the temperature is increased above the melting point, a rapid decrease in seismic velocity is accompanied by a slow decrease in density. For the Central Andes, a good fit of the observed and calculated Bouguer anomalies is obtained if the densities of the rocks from the low-velocity zone (LVZ) beneath the Western Cordillera and the Altiplano are varied. Model calculations lead to a velocity-density relation for partial molten rocks that allows the melt proportions of rocks to be estimated. Model calculations indicate that 15–20 vol.% of basaltic to andesitic melt at depth is necessary to explain the LVZ and Bouguer anomaly beneath the arc and parts of the backarc. High heat flow values (100 mW/m²) support the idea that large areas of the deeper Andean crust are strongly weakened by the presence of partially molten rocks, resulting in reduced seismic velocities, with the Western Cordillera, the active volcanic arc of the Andean mountain range, acting as a ductile buffer between the two more rigid crustal blocks of the forearc and backarc regions.     

Variation in element release rate from different mineral size fractions from the B horizon of a granitic podzol

Far-from-equilibrium batch dissolution experiments were carried out on the 2000–500, 500–250, 250–53 and 53–2 μm size fractions of the mineral component of the B horizon of a granitic iron humus podzol after removal of organic matter and secondary precipitates. The different size fractions were mineralogically and chemically similar, the main minerals present being quartz, alkali and plagioclase feldspar, biotite and chlorite. Specific surface area increased with decreasing grain size. The measured element release rates decreased in the order 53–2>>>2000–500>500–250>250–53 μm. Surface area normalised element release rates from the 2000–500, 500–250 and 250–53 μm size fractions (0.6–77×10⁻¹⁴ mol/m²/s) were intermediate between literature reported surface area normalised dissolution rates for monomineralic powders of feldspar (0.1–0.01×10⁻¹⁴ mol/m²/s) and sheet silicates (100×10⁻¹⁴ mol/m²/s) dissolving under similar conditions. Element release rates from the 53–2 μm fraction (400–3000×10⁻¹⁴ mol/m²/s) were a factor of 4–30 larger than literature reported values for sheet silicates. The large element release rate of the 53–2 μm fraction means that, despite the small mass fraction of 53–2 μm sized particles present in the soil, dissolution of this fraction is the most important for element release into the soil. A theoretical model predicted similar (within a factor of <2) bulk element release rates for all the mineral powders if observed thicknesses of sheet silicate grains were used as input parameters. Decreasing element release rates with decreasing grain size were only predicted if the thickness of sheet silicates in the powders was held constant. A significantly larger release rate for the 53–2 μm fraction relative to the other size fractions was only predicted if either surface roughness was set several orders of magnitude higher for sheet silicates and several orders of magnitude lower for quartz and feldspars in the 53–2 μm fraction compared to the other size fractions or if the sheet silicate thickness input in the 53–2 μm fraction was set unrealistically low. It is therefore hypothesised that the reason for the unpredicted large release rate from the 52–3 μm size fraction is due to one or more of the following reasons: (1) the greater reactivity of the smaller particles due to surface free energy effects, (2) the lack of proportionality between the BET surface area used to normalise the release rates and the actual reactive surface area of the grains and, (3) the presence of traces quantities of reactive minerals which were undetected in the 53–2 μm fraction but were entirely absent in the coarser fractions.     

The distribution of radiogenic heat production as a function of depth in the Sierra Nevada Batholith, California

Geochemical analyses and geobarometric determinations have been combined to create a depth vs. radiogenic heat production database for the Sierra Nevada batholith, California. This database shows that mean heat production values first increase, then decrease, with increasing depth. Heat production is 2 μW/m³ within the 3-km-thick volcanic pile at the top of the batholith, below which it increases to an average value of 3.5 μW/m³ at 5.5 km depth, then decreases to 0.5–1 μW/m³ at 15 km depth and remains at these values through the entire crust below 15 km. Below the crust, from depths of 40–125 km, the batholith's root and mantle wedge that coevolved beneath the batholith appears to have an average radiogenic heat production rate of 0.14 μW/m³. This is higher than the rates from most published xenolith studies, but reasonable given the presence of crustal components in the arc root assemblages. The pattern of radiogenic heat production interpreted from the depth vs. heat production database is not consistent with the downward-decreasing exponential distribution predicted from modeling of surface heat flow data. The interpreted distribution predicts a reasonable range of geothermal gradients and shows that essentially all of the present day surface heat flow from the Sierra Nevada could be generated within the 35 km thick crust. This requires a very low heat flux from the mantle, which is consistent with a model of cessation of Sierran magmatism during Laramide flat-slab subduction, followed by conductive cooling of the upper mantle for 70 m.y. The heat production variation with depth is principally due to large variations in uranium and thorium concentration; potassium is less variable in concentration within the Sierran crust, and produces relatively little of the heat in high heat production rocks. Because silica content is relatively constant through the upper 30 km of the Sierran batholith, while U, Th, and K concentrations are highly variable, radiogenic heat production does not vary directly with silica content.     

Cooling history of lunar Mg-suite gabbronorite 76255, troctolite 76535 and Stillwater pyroxenite SC-936: The record in exsolution and ordering in pyroxenes

We have determined cooling rates of orthopyroxene crystals from two Mg-suite lunar samples (gabbronorite 76255 and troctolite 76535) and one terrestrial sample (orthopyroxenite SC-936 from the Stillwater Complex), on the basis of their Fe–Mg ordering states. In addition, a cooling rate of 76255 was determined by modeling the formation of exsolution lamellae in pyroxenes. The M1–M2 site occupancies of the orthopyroxene crystals were determined by single crystal X-ray diffraction and the rate constant for the ordering reaction was used along with calibrations of the equilibrium intracrystalline fractionation of Fe and Mg as a function of temperature to calculate cooling rates. The closure temperatures (T_C) of cation ordering are 525 °C for 76255, 500 °C for 76535 and 350 °C for SC-936 corresponding to cooling rates of 4 × 10⁻² °C/year at the closure temperature for the lunar samples and 10⁻⁶ °C/year for the Stillwater sample. A cooling rate for 76255, determined by simulating the exsolution process, is 1.7 × 10⁻² °C/year at a closure temperature for exsolution of 700 °C. The Fe–Mg ordering cooling rate determined for 76535 reflects a complex thermal history superimposed on the initial plutonic provenance established for this sample [McCallum, I.S., Schwartz, J.M., 2001. Lunar Mg suite: thermobarometry and petrogenesis of parental magmas. J. Geophys. Res. 106, 27969–27983]. The preservation of a crystallization age of 4.51 Ga and a metamorphic age of 4.25 Ga for 76535 is consistent with a model in which excavation of this sample from the lower lunar crust took place while the sample was at a temperature above the closure temperatures for the Sm–Nd, U–Pb and Ar–Ar isotopic systems. Temperatures in excess of the isotopic closure temperatures (i.e., >600 °C) in the lower lunar crust were maintained by heat diffusing from concentrations of U- and Th-rich KREEP material at the base of the crust. On the other hand, 76255 formed at a much shallower depth in the lunar crust (2 km) and was well below its isotopic closure temperatures at the time of excavation, most likely during the Serenitatis basin-forming impact event. Both lunar samples were reheated during transport to the surface and deposition in hot ejecta blankets. The reheating was short lived but apparently sufficient to redistribute Fe and Mg in M sites in orthopyroxenes. For the lunar samples, the cooling rates based on Fe–Mg ordering represent final stage cooling within an ejecta blanket.     

A review and analysis of silicate mineral dissolution experiments in natural silicate melts

We present a database and a graphical analysis of published experimental results for dissolution rates of olivine, quartz plagioclase, clinopyroxene, orthopyroxene, spinel, and garnet in basaltic and andesitic melts covering a range of experimental temperatures (1100–1500°C) and pressures (10⁵ Pa-3.0 GPa). The published datasets of Donaldson (1985, 1990) and Brearly and Scarfe (1986) are the most complete. Experimental dissolution rates from all datasets are recalculated and normalized to a constant oxygen basis to allow for direct comparison of dissolution rates between different minerals. Dissolution rates (ν) range from 5·10⁻¹⁰ oxygen equivalent moles (o.e.m.) cm⁻² s⁻¹ for olivine in a basaltic melt to 1.3·10⁻⁵ o.e.m. cm⁻² s⁻¹ for garnet in a basaltic melt. Values of ln ν are Arthenian for the experiments examined and activation energies range from 118 to 1800 kJ/o.e.m. for quartz and clinopyroxene, respectively.

The relationship between calculated A/RT for the dissolution reactions, where A is the thermodynamic potential affinity, and values of ν is linear for olivine, plagioclase, and quartz. We interpret this as strong evidence in support of using calculated A as a predictor of ν for, at least, superliquidus melt conditions.     

Studying the migration behaviour of selenate in Boom Clay by electromigration

For the disposal of high-level waste (HLW) in a deep geological formation as Boom Clay, safety assessment studies have shown that long lived ⁷⁹Se is one of the more critical fission products. Therefore, the knowledge of its migration properties (diffusion, retention) through the geological barrier (Boom Clay) is of paramount importance. The migration behaviour of selenium strongly depends on its speciation. Under reducing conditions, selenide would be the dominant species and selenium migration would mainly be controlled by the low solubility of Se(−II)-bearing minerals. However Se species are often found in redox disequilibrium and more oxidized species might also coexist. Therefore, the study of selenate migration requires attention, as it might be the most mobile selenium species in the host rock. Electromigration experiments performed with a ⁷⁵Se-labeled selenate in Boom Clay indicate a high mobility for this species. The apparent diffusion coefficient (D_app) of selenate in Boom Clay is estimated from electromigration experiments performed under different electric fields. Using two independent approaches, the value of D_app for selenate is shown to fall in the range from 1.7×10⁻¹¹ to 6.2×10⁻¹¹ m² s⁻¹. Moreover, no reduction of selenate in Boom Clay was observed.     

Crustal structure in the southwestern part of the Iberian Peninsula

Seismic investigations to determine the crustal structure in the southwestern part of the Iberian Peninsula have been initiated in 1970. First experiments were carried out during July 1970, when a series of ten shots was fired off Cabo de Sines (Portugal) in shallow water and recorded up to distances of 185 km along a SE-profile towards Huelva (Spain). The profile was reversed in December 1970, when a series of twelve shots was fired off the south coast near Fuzeta (east of Faro) and recorded up to distances of about 260 km along a NW-profile towards Cabo da Roca west of Lisboa. A considerable increase in the seismic efficiency of the explosions could be achieved by generating standing waves in the water.

The structure deduced exhibits some peculiar features. Below the Palaeozoic sediments a fairly high velocity of 6.4 km/sec is found for the dome-shaped basement in that area. The lower crust, which is separated from the upper crust by a distinct velocity inversion (with a minimum velocity of about 5.3–5.6 km/sec), is characterized by a velocity of 7.1 km/sec. From the geological evidence and the sequence of seismic velocities it must be concluded that the upper crustal block in the southwestern part of the Iberian Peninsula has been uplifted by about 2–5 km since Permo-Triassic time, thus emphasizing the significance of vertical movement in tectonic activity.

The top of the upper mantle (8.15 km/sec) was detected at a depth of 30 km close to the Atlantic coast in the west, while near the Algarve coast in the south the depth to the M-discontinuity is about 34–35 km. This result in conjunction with studies of earthquake focal mechanisms confirms the suggestion that the Iberian block is being underthrust under the African plate.     

Determining mineral weathering rates based on solid and solute weathering gradients and velocities: application to biotite weathering in saprolites

Chemical weathering gradients are defined by the changes in the measured elemental concentrations in solids and pore waters with depth in soils and regoliths. An increase in the mineral weathering rate increases the change in these concentrations with depth while increases in the weathering velocity decrease the change. The solid-state weathering velocity is the rate at which the weathering front propagates through the regolith and the solute weathering velocity is equivalent to the rate of pore water infiltration. These relationships provide a unifying approach to calculating both solid and solute weathering rates from the respective ratios of the weathering velocities and gradients. Contemporary weathering rates based on solute residence times can be directly compared to long-term past weathering based on changes in regolith composition. Both rates incorporate identical parameters describing mineral abundance, stoichiometry, and surface area.

Weathering gradients were used to calculate biotite weathering rates in saprolitic regoliths in the Piedmont of Northern Georgia, USA and in Luquillo Mountains of Puerto Rico. Solid-state weathering gradients for Mg and K at Panola produced reaction rates of 3 to 6×10⁻¹⁷ mol m⁻² s⁻¹ for biotite. Faster weathering rates of 1.8 to 3.6×10⁻¹⁶ mol m⁻² s⁻¹ are calculated based on Mg and K pore water gradients in the Rio Icacos regolith. The relative rates are in agreement with a warmer and wetter tropical climate in Puerto Rico. Both natural rates are three to six orders of magnitude slower than reported experimental rates of biotite weathering.     

Crustal properties from S-wave and gravity data along a seismic refraction profile in Romania

VRANCEA'99 is a seismic refraction line that was carried out in 1999 to investigate the deep structure and physical properties of the upper lithosphere of the southeastern Carpathians and its foreland. It runs from the city of Bacau to the Danube River, traversing the Vrancea epicentral area of strong intermediate-depth seismicity and the city of Bucharest.

Interpretation of P-wave arrivals led to a velocity model that displays a multi-layered crust with velocities increasing with depth. The range of P wave velocities in the sedimentary cover increases from N to S and a structuring of the autochthonous basement of the Moesian Platform is observed. The crystalline crust displays thickness variations, but at the same time the lateral velocity structure along the seismic line remains almost constant. An intra-crustal boundary separates an upper crust from the lower crust. Within the upper mantle a low velocity zone is detected at a depth of about 55-km.

The interpretation of observable S-waves resulted in a velocity model that shows the same multi-layered crust, with S-velocities increasing similarly with depth as the P-waves. The corresponding Poisson's ratio is highly variable throughout the crust and ranges from 0.20–0.35 for the sedimentary cover to 0.22–0.25 for the crystalline crust. The interpretation of the V_p, V_s and Poisson's ratio in petrological terms suggests a large variety of rocks from sand and clay to sandstone, limestone and dolomite within the sedimentary cover. Within the crystalline crust the most probably rock types are granite, granodiorite, granite–gneiss and/or felsic amphibolite–gneiss in the upper part and gneiss and /or amphibolite in the lower part.

Based on the 2-D seismic velocity model, a density model is developed. Density values are assigned to each layer in agreement with the P-wave velocity model and with values accepted for the geological units in the area. After several iterations a good fit between the computed and observed Bouguer anomalies was obtained along the seismic line.     

Ray path interpretation of the crustal structure beneath Saudi Arabia

Interpretation of a long-range seismic refraction line in Saudi Arabia has shown that beneath the Arabian Shield velocity generally increases with depth, from about 6 km s⁻¹ at the surface to about 7 km s⁻¹ at the top of the crust-mantle transition zone. The base of this transition zone (Moho) occurs at 37–44 km in depth. Intracrustal discontinuities can also be recognized, the most important being in the 10–20 km-depth range and separating the upper from the lower crust. Laterally, the variations in the intracrustal discontinuities and the total crustal thickness can be correlated with previously defined tectonic regions. Beneath the Red Sea shelf and coastal plain the crust, including 4 km of sediments, is only 15–17.5 km thick. With the aid of both seismic and gravity data an abrupt, steeply dipping transition from the crust of the Red Sea shelf and coastal plain to that of the Arabian Shield has been derived. With a jump of more than 20 km in Moho depth, this appears to be the major discontinuity between the Red Sea depression and the Arabian continental shield.