Some remarks on high-temperature—low-pressure metamorphism in convergent orogens

Many high-temperature–low-pressure (high- T –low- P ) metamorphic terranes show evidence for peak mineral growth during crustal thickening strain increments at pressures near the maximum attained during the heating–cooling cycle. Such terranes are not readily explained as the conductive response to crustal thickening since the resulting Moho temperatures would greatly exceed the crustal liquidus and because heating due to conductive equilibration on length scales appropriate to lithospheric-scale strains must greatly outlast the deformation. Consequently, high- T –low- P metamorphism may be generated during crustal thickening only when significant heat is advected within the crust, as for example may occur during the segregation of granitic melts. We show that without the addition of asthenospheric melts and at strain rates appropriate to continental deformation the conditions required for significant lower crustal melting during deformation are only likely to be attained if heat flow into the lower crust during crustal thickening is increased substantially, for example, by removing the mantle part of the lithosphere. A simple parameterization of lithospheric deformation involving the vertical strain on the scale of the crust, _c, and the lithosphere, ₁ respectively, allows the potential energy of the evolving orogen to be readily evaluated. Using this parameterization we show that an important isostatic consequence of the deformation geometries capable of generating such high- T –low- P metamorphism during crustal thickening (with _c1) is an imposed upper limit to crustal thicknesses which is much lower than for homogeneous deformations (f_c= f₁) for the same initial lithospheric configuration.     

Extrusion vs. accretion at the frictional–viscous décollement transition in experimental thrust wedges: the role of convergence velocity

In many cases, thrust wedges accreted at shallow crustal levels show an across‐strike rheological variability along the basal décollement, notably from brittle to ductile behaviour. In this paper, we illustrate the results of sandbox analogue modelling research devoted to studying the influence of convergence velocity on wedge architecture when laterally juxtaposed frictional and viscous materials occur along the basal décollement of accreting thrust wedges. Our results show that slow convergence favours a near symmetrical distribution of thrust vergence within wedge sectors accreted above viscous décollement material, whereas fast convergence favours vergence asymmetry. In particular, at fast convergence rates the hinterlandward extrusion of viscous décollement material at the toe of the frictional wedge is favoured and contributed to accommodate a significant amount of the total contraction. Terra Nova, 18, 241–247, 2006     

Boundary effects on stream meandering and river morphology

ABSTRACT An account of stream meandering and river morphology is given for five experiments, each having different conditions at the wetted perimeter of the channel.
On the basis of the present experiments, it is proposed that a necessary (but insufficient) condition for the formation of stream meandering is an appropriately small ratio of total bank shear force to total bed shear force, both per unit length downstream . In mathematical form, the criterion is expressed by τ_s and τ_s are an average bank stress and bed stress, respectively, and p _s and p _bb are an average bank wetted perimeter and bed wetted perimeter of a half-channel, respectively. The value of α is evaluated by Kinoshita's (1957) experiment and others.
It is also found that the flow pattern, channel morphology and sand waves on the bed are dependent upon the boundary conditions at the wetted perimeter.     

Magnification of undamped seismographs and the analysis of the 1906 San Francisco earthquake record obtained on the Vicentini seismograph in Zagreb (Croatia)

A method to estimate the dynamic magnification of nonlinear seismographs is presented. The method was tested by applying it to the record of the 1906 San Francisco earthquake seismogram obtained on the Vicentini seismograph in Zagreb, in order to reconstruct the ground motion amplitudes. The magnitudes inferred from the corrected Vicentini recording of P- and surface waves are internally consistent ( m _b= 7.74, M _s= 8.48) and agree with the values published in various catalogues, within the error limits of magnitude determination for individual events. The corrected P-waveform is compatible with a rupture slip profile exhibiting several distinct maxima along the 300 km-long section of the San Andreas fault, in accordance with recent studies and contemporary geodetic measurements.     

Local displacement rate cycles in the life of a fold-thrust belt

The evolution of minor structures during the growth of major folds and thrusts, in the Chartreuse district of the French Subalpine thrust belt, shows that each thrust evolved through a phase of distributed faulting, major thrust propagation and displacement, followed by distributed shear modification of the hanging-wall fold. Microstructural studies suggest that the distributed faulting phases, early and late in the history, were characterized by strain rates limited by diffusive mass transfer processes ( c . 10^-15-10^-16 s^-1). Faulting whose rate is limited by DMT is too slow on its own to accommodate the regional time-averaged shortening rates for the thrust belt as a whole, implying that the slow thrusts operated in tandem with those major, fast thrusts where deformation was primarily cataclastic. Consequently each thrust anticline experienced a displacement rate cycle and an array of thrust anticlines must amplify simultaneously. These interpretations raise important issues for the dynamics of fault populations, the evolution of thrust wedges and the history of fluid migration in thrust belts.     

Sensitivity of shear zones in orogenic wedges to surface processes and strain softening

We investigate the internal deformation of orogenic wedges growing by frontal accretion with a two-dimensional numerical model. Our models are limited to crustal deformation and assume a horizontal detachment as observed for various natural orogens (e.g. Alaska and Costa Rica). The model wedges develop as a result of convergence of a brittle sediment layer in front of a strong backstop. We find that our reference model develops in-sequence forward-thrusts which propagate upward from the basal detachment. For this reference model we investigate the sensitivity of shear zone activity to surface processes and strain softening. Model results show that diffusive or slope dependent erosion enhances material transport across the wedge and slows down forward propagation of the deformation front. Frictional strain softening focuses deformation into narrow shear zones and enhances displacement along them. This has also been postulated for natural thrusts such as the Glarus thrust in the Swiss Alps and the Moine thrust in the Scottish Caledonides. A second series of models investigates the effects of regularly spaced weak inclusions within the sediment layer which simulate remnants of previous deformation phases. These inclusions facilitate and focus internal deformation, influence the thrust dip and thrust vergence and enable thrust reactivation in the internal part of the wedge. Our results show that inactive thrusts in the internal part of the wedges may be reactivated in models with diffusive surface processes, strain softening or weak inclusions. Thrust reactivation occurs as models seek to maintain their critical taper angle. First order characteristics of our numerical models agree well with natural orogenic wedges and results from other numerical and analogue models.     

CO2 transfer between the upper mantle and the atmosphere: temporary storage in the lower continental crust

The CO₂ atmospheric content has shown large variations over geological times. High contents (up to one order of magnitude more than present-day values) ultimately correspond to discrete episodes of mantle degassing, either juvenile, or subduction-related (carbon recycling). A number of arguments (e.g. the continuous volume increase of carbonate-bearing sediments with time) suggest that, throughout the Earth's history, juvenile CO₂ has formed a major contribution to the global carbon budget of the Earth.
The absence of a direct relationship between major volcanic episodes and the average CO₂ atmospheric content suggests that volcanoes might not be the only way by which mantle CO₂ is transported to the surface. It is proposed that large quantities of juvenile CO₂ could temporarily be stored in the lower continental crust during major episodes of granulite formation. These are primarily caused by magmatic underplating and they result in a vertical accretion of the crust by accumulation of CO₂-bearing, mantle-derived magmas. Most of the CO₂ migrates through the crust during post-metamorphic evolution and isostatic restoration of the normal continental thickness. However, large quantities of CO₂ can still be present in some areas, notably as high-density fluids enclosed in minerals.     

Major Element Data (1966–1978) for The Six "Nimroc" Reference Samples

The results are given for all the determinations made by the co-operating laboratories of major elements in the six NIMROC rock samples (granite, syenite, lujavrite, norite, pyroxenite and dunite) prepared by the National Institute for Metallurgy in 1966. Relevant statistical data are given for the sets of results for each major constituent, and recommended values for all constituents except Al₂0₃, Na₂0, K₂0and CO₂ in dunite, Fe₂O₃, MgO and CO₂ in the granite, Fe₂₃ and CO₂ in the norite and CO₂ in the pyroxenite.
This report and that on the trace and minor elements issued in 1978 complete the revision of the recommended values. It is suggested that analysts should concentrate rather on those constituents for which the results have shown such a wide scatter that they can be of no usc for reference purposes, than on those for which the ualues are fairly well established.     

Petrographic Characteristics and Alteration Geochemistry of Granite-hosted Tungsten Mineralization at Degana, NW India

Abstract. Vein type tungsten mineralization at Degana is genetically and spatially associated with the Degana Granite. The deposit is characterized by pervasive wall rock alteration around the mineralized quartz veins. Laterally three different alteration zones, greisen, silicification and potassic zones, are marked based on the field features, mineral assemblages and geo-chemical characteristics. In the present paper, systematic mineralogical and chemical variation in these alteration zones is reported. Thick mono-mineralic (zinnwaldite) selvages around the veins characterize the deposit. Plagioclase and alkali feldspar are low in the greisen zones while K-feldspar shows more increase than plagioclase in the potassic zone. Quartz is uniformly high in all the alteration zones, but it shows an anomalous value in the silicification zone. Al₂O₃ concentration shows initial depletion in greisen zone with gradual increase away from the contact. MgO and FeO are higher in greisen zone than silicification and potassic zones. The potassic zone is characterized by the depletion of Na₂O and higher value of K₂O.
The common presence of topaz and fluorite as both primary and secondary minerals and fluorine-bearing micas suggest fluorine partitioning in substantial amount between granitic melt and coexisting aqueous fluid phase and higher HF activity during the evolution of hydrothermal fluid. The mutual relationship of the fluorine minerals (topaz and fluorite) in the different alteration zones suggests an increase in the Ca²⁺ activity and decrease of H+ activity during the fluid evolution from greisenization towards alkali-metasomatised granite and the fluid is assumed to change from low to high activity ratio of Ca²⁺/H+.     

Oxygen isotope exchange and closure temperatures in cooling rocks

Retrograde exchange of oxygen isotopes between minerals in igneous and metamorphic rocks by means of diffusion is explored using a finite difference computer model, which predicts both the zonation profile of δ¹⁸O within grains, and the bulk δ¹⁸O value of each mineral in the rock. Apparent oxygen isotope equilibrium temperatures that would be observed in these rocks are calculated from the δ¹⁸O values of each mineral pair within the rock. In systems which cool linearly from a sufficiently high temperature or at a low enough cooling rate, such that the final oxygen isotope values are not dependent upon the initial oxygen isotope values ('slow cooling'), the apparent oxygen isotope temperature derived for a rock composed of a single mineral pair can be shown to be simply related to the Dodson closure temperatures ( T _c) for the two phases and the mode of the rock. Adding a third phase into a system which undergoes 'slow' cooling will cause the apparent temperature derived for the two minerals already present to differ from the simple relationship for a two-phase system. In some systems oxygen isotope reversals can be developed. If cooling is not 'slow', then the mineral δ¹⁸O values resulting from cooling will be partly dependent upon the initial temperature of the system concerned. The model successfully simulates the mineral δ¹⁸O values that are often observed in granitic rocks. Application of the model will help in assessing the validity of oxygen isotope thermometry in different geological settings, and allows quantitative prediction of the oxygen isotope fractionations that are developed in cooling closed systems.     

Refinement of thrust faulting models for the Central New Madrid seismic zone

We present the results of the joint relocation of events recorded during 1989–1992 by the PANDA network in the central New Madrid seismic zone. The near-surface material in the study area is a gently-dipping layer of poorly consolidated sediments with low P-wave velocity and high V_p/V_s (estimated values: 1.8 km s⁻¹ and 3). The sediments are underlain by high-velocity Paleozoic rocks. Under the network the difference in sediment thickness is only 0.6 km, but because of the low velocities the location of the events using layered models is affected by errors. Application of the joint hypocentral determination (JHD) technique to a subset of 580 events shows that the single-event locations may be in error by as much as 1 km in depth, depending on where the events are located. Analysis of synthetic data generated for a realistic 3-D velocity model supports the JHD results. The analysis of synthetic data also suggests that a V_p/V_s≤ 2.3 is more appropriate for the post-Paleozoic Mississippi embayment sediments. Based on the JHD locations we present a new interpretation of the seismicity, with two en-echelon SW-dipping thrust faults connected by a west-dipping thrust fault. These faults appear associated with the Reelfoot scarp and its northern extension, the Kentucky bend scarp.     

Impact of surface processes on the growth of orogenic wedges: Insights from analog models and case studies

Interaction between surface processes and deep tectonic processes plays a key role in the structural evolution, kinematics and exhumation of rocks in orogenic wedges. The deformation patterns observed in analogue models applied to natural cases of present active or ancient mountain belts reflect several first order processes that result of these interactions. Internal strain partitioning due to mechanical behaviour of a thrust wedge has a strong impact on the vertical component of displacement of tectonic units that in return favour erosion in domains of important uplift. Such strain partitioning is first controlled by tectonic processes, but surface processes exert a strong feed back on wedge dynamics. Indeed, material transfer in thrust wedges not only depends on its internal dynamics, it is also influenced by climate controlled surface processes involving erosion and sedimentation. Effects of erosion are multiple: they allow long term localization of deformed domains, they favour important exhumation above areas of deep underplating and combined with sedimentation in the foreland they contribute to maintain the wedge in a critical state for long time periods. The simple models illustrate well how mountain belts structure, kinematics of tectonic units and exhumation are determined by these complex interactions.     

Additional Provisional Elemental Values for LKSD-1, LKSD-2, LKSD-3, LKSD-4, STSD-1, STSD-2, STSD-3 and STSD-4

Provisional values related to element concentration derived from partial extractions as well as some additional total concentration data are reported for eight CCRMP sediment reference materials. These values are intended to supplement previously published data for these samples. The partial extractions for which data are reported include dilute HNO₃ - dilute HCl, concentrated HNO₃ -concentrated HCl and concentrated HNO₃ -concentrated HClO₄.     

On the mechanics of basin formation in the Pannonian basin: Inferences from analogue and numerical modelling

We present the results of a thrust fault reactivation study that has been carried out using analogue (sandbox) and numerical modelling techniques. The basement of the Pannonian basin is built up of Cretaceous nappe piles. Reactivation of these compressional structures and connected weakness zones is one of the prime agents governing Miocene formation and Quaternary deformation of the basin system. However, reactivation on thrust fault planes (average dip of ca. 30°) in normal or transtensional stress regimes is a problematic process in terms of rock mechanics. The aim of the investigation was to analyse how the different stress regimes (extension or strike-slip), and the geometrical as well as the mechanical parameters (dip and strike of the faults, frictional coefficients) effect the reactivation potential of pre-existing faults.

Results of analogue modelling predict that thrust fault reactivation under pure extension is possible for fault dip angle larger than 45° with normal friction value (sand on sand) of the fault plane. By making the fault plane weaker, reactivation is possible down to 35° dip angle. These values are confirmed by the results of numerical modelling. Reactivation in transtensional manner can occur in a broad range of fault dip angle (from 35° to 20°) and strike angle (from 30° to 5° with respect to the direction of compression) when keeping the maximum horizontal stress magnitude approximately three times bigger than the vertical or the minimum horizontal stress values.

Our research focussed on two selected study areas in the Pannonian basin system: the Danube basin and the Derecske trough in its western and eastern part, respectively. Their Miocene tectonic evolution and their fault reactivation pattern show considerable differences. The dominance of pure extension in the Danube basin vs. strike-slip faulting (transtension) in the Derecske trough is interpreted as a consequence of their different geodynamic position in the evolving Pannonian basin system. In addition, orientation of the pre-existing thrust fault systems with respect to the Early to Middle Miocene paleostress fields had a major influence on reactivation kinematics.

As part of the collapsing east Alpine orogen, the area of the Danube basin was characterised by elevated topography and increased crustal thickness during the onset of rifting in the Pannonian basin. Consequently, an excess of gravitational potential energy resulted in extension (σ_v > σ_H) during Early Miocene basin formation. By the time topography and related crustal thickness variation relaxed (Middle Miocene), the stress field had rotated and the minimum horizontal stress axes (σ_h) became perpendicular to the main strike of the thrusts. The high topography and the rotation of σ_h could induce nearly pure extension (dip-slip faulting) along the pre-existing low-angle thrusts. On the contrary, the Derecske trough was situated near the Carpathian subduction belt, with lower crustal thickness and no pronounced topography. This resulted in much lower σ_v value than in the Danube basin. Moreover, the proximity of the retreating subduction slab provided low values of σ_h and the oblique orientation of the paleostress fields with respect to the master faults of the trough. This led to the dominance of strike-slip faulting in combination with extension and basin subsidence (transtension).     

Preliminary Investigation into the Use of a High Resolution Inductively Coupled Plasma-Mass Spectrometer with Laser Ablation for Bulk Analysis of Geological Materials Fused with Li2B4O7

The paper presents preliminary results of the use of a high resolution double-focussing, magnetic sector inductively coupled plasma-mass spectrometer (HR-ICP-MS) with ultraviolet laser ablation (LA) for the bulk analysis of geological materials fused with Li₂B₄O₇. Detection limits are based on data from precision measurements of a fused SiO₂ sample of high purity, and sensitivity data (cps/μg g^-1) obtained on the Reference Material (RM) Syenite SY-2. For many trace elements, the detection limits are better than 0.05 μg g^-1 using a sample to flux weight ratio of 1:7.
Calibration curves, which are based entirely on RMs, are established for Hf, Ta, Tb, Tm and Lu. They indicate that, even at this early stage in the development of the technique, data accurate to ˜ 25% can be collected. It is concluded that the method may prove to be a valuable supplement to XRF for low level element concentration measurements; it is also very practical, as the same sample discs can be used for both XRF and LA-ICP-MS analyses.     

Hydrogen, Oxygen and Sulfur Isotope Studies of Seafloor Hydrothermal System at the Desmos Caldera, Manus Back-arc Basin, Papua New Guinea: An Analogue of Terrestrial Acid Hot Crater-lake

Abstract. The Onsen site is an active submarine hydrothermal system hosted by the Desmos caldera in the Eastern Manus Basin, Papua New Guinea. The hydrothermal fluid is very acidic (pH=1.5) and abundant native sulfur is deposited around the vent. The δ³⁴S values of native sulfur range from -6.5 to -9.3 %o. δ³⁴S values of H₂S and SO₄ in the hydrothermal fluid are -4.3 to -9.9 %o and +18.6 to +20.0 %o, respectively. These δ³⁴S values are significantly lower than those of the other hydrothermal systems so far reported. These low δ³⁴S values and the acidic nature of the vent fluids suggest that volcanic SO₂ gas plays an important role on the sulfur isotope systematic of the Onsen hydrothermal system. Relationship among the δ³⁴S values of S-bearing species can be successively explained by the model based on the disproportionation reaction starting from the volcanic SO₂ gas. The predicted δ³⁴S values of SO₂ agree with the measured whole rock δ³⁴S values. δD and δ¹⁸O values of clay minerals separated from the altered rock samples also suggest the contribution of the magmatic fluid to the hydrothermal system. Present stable isotopic study strongly suggests that the Onsen hydrothermal site in the Desmos caldera is a magmatic submarine hydrothermal system.     

The Determination of Total Hydrogen, Carbon, Nitrogen and Sulfur in Silicates, Silicate Rocks, Soils and Sediments

Determinations of the volatile elements carbon, hydrogen, sulfur and nitrogen in many geological RM, performed with the LECO CHN and SC analysers, are presented. The method allowed the determination of S in concentrations from a few % m/m to 0.001% m/m or less, of C from % m/m to 0.01% m/m and of H from % m/m to 0.004% m/m. Accuracy was usually better than the XRF method (for S). All obtained values passed the Sutarno-Steger test, which establishes that |(mean_analysed - mean_certified)|/ S_certified < 2, for the cases with an appropriate number of determinations (n > 10 for each element). It was possible to perform routine determination of C, H and S with the instrumentation, coupled with the determination of major and minor elements in geological materials. Determination of nitrogen could also be performed on an exploratory basis, with improvements in the method dependent on the future availability of more reference materials with reliable composition of this element.     

Climate variability in the SE Alps of Italy over the past 17 000 years reconstructed from a stalagmite record

Stalagmite SV1 from Grotta Savi, located at the SE margin of the European Alps (Italy), is the first Alpine speleothem that continuously spans the past c . 17kyr. Extension rate and δ¹⁸O_c record for the Lateglacial probably reflect a combination of temperature and rainfall, with rainfall exerting the dominant effect. Low speleothem calcite δ¹⁸ O_c values were recorded from c . 14.5 and 12.35 kyr, during GI-1 (Bølling— Allerød) interstadial, which in our interpretation, was warm and wet. The GS-1 (Younger Dryas) was characterized by a shift to heavier δ¹⁸ O_c, coinciding with δ¹³C_c enrichment and extremely low extension rate (<8 μm/year). These characteristics indicate that GS-1 climate was cool and dry in the SE Alps. Calibration using historical data revealed that there is a positive δ¹⁸O_c/dT relationship. A 1°C rise in mean annual temperature should correspond to c . 2.85% increase of SV-1 δc¹⁸O_c. We reconstructed a slow and steady temperature rise of c . 0.5°C since 10 kyr BP, in broad agreement with reconstructions from pollen data for SE Europe. Stalagmite SV1 indicates that climate variability in the SE Alps has been influenced by the Mediterranean Sea for the past c . 17 kyr.     

A Simple Method for Precise Determination of 23 Trace Elements in Granitic Rocks by ICP-MS after Lithium Tetraborate Fusion

Abstract. An improved alkali fusion method followed by HF-HNO₃-HC1O₄ treatment is performed for simultaneous determination of 23 trace elements (Sr, Y, Zr, Nb, Ba, Hf, Ta, Th, U, and REE) by ICP-MS in rock reference materials: basaltic rocks (JB-2, JB-3) and granitic rocks (JG-la, JG-2, JG-3). Our improved method offers several advantages including: (1) suppression of whitish precipitates probably composed of insoluble fluorides by addition of HCIO4, (2) simple and reliable preparation procedure, (3) instrument calibration which enables straightforward simultaneous multi-elemental analysis, and (4) the very low background levels by using pure lithium tetraborate flux. We obtained the analytical results with a reproducibility of mostly <2 % (1) for the basaltic rocks and <7 % for the granitic rocks. The higher relative standard deviation (RSD) values for granitic rocks may be attributed to sample heterogeneity of coarse-grained granitic rocks. The analytical results of the granitic rocks demonstrate that Zr and Hf abundances are consistent with the compiled values and that REE concentrations agree well with recently published data, suggesting that the Li₂B₄O₇ fusion method applied in the present study is suitable for the analysis of the granitic rocks.     

Empirical garnet–muscovite geothermometry in low-grade metapelites, Selwyn Range (Canadian Rockies)

Abstract Partitioning of Fe and Mg between garnet and phengitic muscovite was calibrated as a geothermometer by Green & Hellman (1982) using experimental data at 25–30 kbar. When the thermometer is applied to pelites regionally metamorphosed at pressures of between 3 and 7 kbar it yields temperatures much higher than those from the garnet–biotite thermometer. A new empirical calibration is proposed for use with such rocks, with particular application where garnet occurs at lower grades than biotite. The new calibration is where K is given by: In K = In K _d and X _ii are mole fractions in the garnets.
The calibration was derived from comparison with the garnet–biotite thermometer of Ferry & Spear (1978), assuming no pressure-dependence for the partitioning between garnet and muscovite, no ferric iron partitioning, ideal mixing in muscovite, and the garnet mixing model of Ganguly & Saxena (1984) modified for a non-linear Ca effect. This latter garnet mixing model was selected because it gave the geologically most reasonable results. It has not proved possible to distinguish a pressure effect from a ferric-iron effect.
Despite the simplifying assumptions used to derive the calibration, it yields temperatures generally within 15°C of those given by the garnet–biotite thermometer, and has been used to supply thermometric data in a low-grade region of the Canadian Rockies.