A plume tectonics model for the Tharsis province, Mars

Morphological and structural data from the whole Tharsis province suggest that a number of shallow grabens radially oriented about the Tharsis bulge on Mars are underlain by dykes, which define giant radiating swarms similar to, e.g. the Mackenzie dyke swarm of the Canadian shield. Mechanisms for graben formation are proposed, and the depth, width, and height of the associated dykes are estimated. Structural mapping leads to define successive stages of dyke emplacement, and provide stress-trajectory maps that indicate a steady source of the regional stress during the whole history of the Tharsis province. A new tectonic model of Tharsis is presented, based on an analogy with dyke swarms on the Earth that form inside hot spots. This model successfully matches the following features: (1) the geometry of the South Tharsis Ridge Belt, which may have been a consequence of the compressional stress field at the boundary between the uplifted and non-uplifted areas in the upper part of the lithosphere at the onset of hot spot activity; (2) extensive lava flooding, interpreted as a consequence of the high thermal anomaly at the onset of plume (hot spot) activity; (3) wrinkle ridge geometry in the Tharsis hemisphere, the formation of which is interpreted as a consequence of buoyant subsidence of the brittle crust in response to the lava load; (4) Valles Marineris limited stretching by preliminary passive rifting, and uplift, viewed as a necessary consequence of adiabatic mantle decompression induced by stretching. The geometrical analysis of dyke swarms suggests the existence of a large, Tharsis-independent extensional state of stress during all the period of tectonic activity, in which the minimum compressive stress is roughly N---S oriented. Although magmatism must have loaded the lithosphere significantly after the plume activity ceased and be responsible for additional surface deformations, there is no requirement for further loading stress to explain surficial features. Comparison with succession of magmatic and tectonic events related to hot spots on the Earth suggests that the total time required to produce all the surface deformation observed in the Tharsis province over the last 3.8 Ga does probably not exceed 10 or 15 Ma.     

Sensitivity analysis of a deterministic water temperature model to forest canopy and soil temperature in Catamaran Brook (New Brunswick,Canada)

A coupled deterministic hydrological and water temperature model, CEQUEAU, was modified to include soil temperature and crown closure in its calculation of local advective terms in the heat budget. The modified model was than tested to verify its sensitivity to these modifications. An analysis of the heat budget of a small forested catchment in eastern Canada revealed that the advective term related to interflow plays a significant role in the daily water heat budget, providing on average 28% of the local advective budget (which also includes advective heat terms from surface runoff and groundwater) and nearly 14% of the total heat budget (which includes all radiative terms at the water surface, convection and evaporation, as well as the local advective terms). Relative sensitivity indices (RSIs) were used to verify the impact of the newly introduced parameters and variables. Among them, parameters related to the forest cover (crown closure and leaf area index) have a maximum RSI of ?0·6; i.e. a 100% increase in value produces a 60% decrease in the local advective term. Parameters with the greatest influence are the volume of water contributing to interflow and the amplitude of the net radiative flux at the soil surface, which, if doubled, would double the contribution of the local interflow advective term to the heat budget. Copyright © 2003 John Wiley & Sons, Ltd.     

Relationship between footwall composition,crustal contamination,and fluid–rock interaction in the Platreef,Bushveld Complex,South Africa

In the northern limb of the 2.06-Ga Bushveld Complex, the Platreef is a platinum group elements (PGE)-, Cu-, and Ni-mineralized zone of pyroxenite that developed at the intrusion margin. From north to south, the footwall rocks of the Platreef change from Archaean granite to dolomite, hornfels, and quartzite. Where the footwall is granite, the Sr-isotope system is more strongly perturbed than where the footwall is Sr-poor dolomite, in which samples show an approximate isochron relationship. The Nd-isotope system for samples of pyroxenite and hanging wall norite shows an approximate isochron relationship with an implied age of 2.17 ± 0.2 Ga and initial Nd-isotope ratio of 0.5095. Assuming an age of 2.06 Ga, the ɛNd values range from −6.2 to −9.6 (ave. −7.8, n = 17) and on average are slightly more negative than the Main Zone of the Bushveld. These data are consistent with local contamination of an already contaminated magma of Main Zone composition. The similarity in isotope composition between the Platreef pyroxenites and the hanging wall norites suggests a common origin. Where the country rock is dolomite, the Platreef has generally higher plagioclase and pyroxene δ ¹⁸O values, and this indicates assimilation of the immediate footwall. Throughout the Platreef, there is considerable petrographic evidence for sub-solidus interaction with fluids, and the Δ _{plagioclase–pyroxene} values range from −2 to +6, which indicates interaction at both high and low temperatures. Whole-rock and mineral δD values suggest that the Platreef interacted with both magmatic and meteoric water, and the lack of disturbance to the Sr-isotope system suggests that fluid–rock interaction took place soon after emplacement. Where the footwall is granite, less negative δD values suggest a greater involvement of meteoric water. Consistently higher values of Δ _{plagioclase–pyroxene} in the Platreef pyroxenites and hanging wall norites in contact with dolomite suggest prolonged interaction with CO₂-rich fluid derived from decarbonation of the footwall rocks. The overprint of post crystallization fluid–rock interaction is the probable cause of the previously documented lack of correlation between PGE and sulfide content on the small scale. The Platreef in contact with dolomite is the focus of the highest PGE grades, and this suggests that dolomite contamination played a role in PGE concentration and deposition, but the exact link remains obscure. It is a possibility that the CO₂ produced by decarbonation of assimilated dolomite enhanced the process of PGE scavenging by sulfide precipitation.     

Controls on porphyroblast size along a regional metamorphic field gradient

Garnet-bearing schists from the Waterville Formation of south-central Maine provide an opportunity to examine the factors governing porphyroblast size over a range of metamorphic grade. Three-dimensional sizes and locations for all garnet porphyroblasts were determined for three samples along the metamorphic field gradient spanning lowest garnet through sillimanite grade, using high-resolution X-ray computed tomography. Comparison of crystal size distributions to previous data sets obtained by stereological methods for the same samples reveals significant differences in mode, mean, and shape of the distributions. Quantitative textural analysis shows that the garnets in each rock crystallized in a diffusion-controlled nucleation and growth regime. In contrast to the typical observation of a correlation between porphyroblast size and position along a metamorphic field gradient, porphyroblast size of the lowest-grade specimen is intermediate between the high- and middle-grade specimens’ sizes. Mean porphyroblast size does not correlate with peak temperatures from garnet-biotite Fe-Mg exchange thermometry, nor is post-crystallization annealing (Ostwald Ripening) required to produce the observed textures, as was previously proposed for these rocks. Robust pseudosection calculations fail to reproduce the observed garnet core compositions for two specimens, suggesting that these calc-pelites experienced metasomatism. For each of these two specimens, Monte Carlo calculations suggest potential pre-metasomatism bulk compositions that replicate garnet core compositions. Pseudosection analyses allow the estimation of the critical temperatures for garnet growth: ∼481, ∼477, and ∼485°C for the lowest-garnet-zone, middle-garnet-zone, and sillimanite-zone specimens, respectively. Porphyroblast size appears to be determined in this case by a combination of the heating rate during garnet crystallization, the critical temperature for the garnet-forming reaction and the kinetics of nucleation. Numerical simulations of thermally accelerated, diffusion-controlled nucleation, and growth for the three samples closely match measured crystal size distributions. These observations and simulations suggest that previous hypotheses linking the garnet size primarily to the temperature at the onset of porphyroblast nucleation can only partially explain the observed textures. Also important in determining porphyroblast size are the heating rate and the distribution of favorable nucleation sites.     

Evidence for crustal degassing of CF4 and SF6 in Mojave Desert groundwaters

Dissolved tetrafluoromethane (CF₄) and sulfur hexafluoride (SF₆) concentrations were measured in groundwater samples from the Eastern Morongo Basin (EMB) and Mojave River Basin (MRB) located in the southern Mojave Desert, California. Both CF₄ and SF₆ are supersaturated with respect to equilibrium with the preindustrial atmosphere at the recharge temperatures and elevations of the Mojave Desert. These observations provide the first in situ evidence for a flux of CF₄ from the lithosphere. A gradual basin-wide enhancement in dissolved CF₄ and SF₆ concentrations with groundwater age is consistent with release of these gases during weathering of the surrounding granitic alluvium. Dissolved CF₄ and SF₆ concentrations in these groundwaters also contain a deeper crustal component associated with a lithospheric flux entering the EMB and MRB through the underlying basement. The crustal flux of CF₄, but not of SF₆, is enhanced in the vicinity of local active fault systems due to release of crustal fluids during episodic fracture events driven by local tectonic activity. When fluxes of CF₄ and SF₆ into Mojave Desert groundwaters are extrapolated to the global scale they are consistent, within large uncertainties, with the fluxes required to sustain the preindustrial atmospheric abundances of CF₄ and SF₆.     

Pressures of Crystallization of Icelandic Magmas

Iceland lies astride the Mid-Atlantic Ridge and was createdby seafloor spreading that began about 55 Ma. The crust is anomalouslythick (20–40 km), indicating higher melt productivityin the underlying mantle compared with normal ridge segmentsas a result of the presence of a mantle plume or upwelling centeredbeneath the northwestern edge of the Vatnajökull ice sheet.Seismic and volcanic activity is concentrated in 50 km wideneovolcanic or rift zones, which mark the subaerial Mid-AtlanticRidge, and in three flank zones. Geodetic and geophysical studiesprovide evidence for magma chambers located over a range ofdepths (1·5–21 km) in the crust, with shallow magmachambers beneath some volcanic centers (Katla, Grimsvötn,Eyjafjallajökull), and both shallow and deep chambers beneathothers (e.g. Krafla and Askja). We have compiled analyses ofbasalt glass with geochemical characteristics indicating crystallizationof ol–plag–cpx from 28 volcanic centers in the Western,Northern and Eastern rift zones as well as from the SouthernFlank Zone. Pressures of crystallization were calculated forthese glasses, and confirm that Icelandic magmas crystallizeover a wide range of pressures (0·001 to 1 GPa), equivalentto depths of 0–35 km. This range partly reflects crystallizationof melts en route to the surface, probably in dikes and conduits,after they leave intracrustal chambers. We find no evidencefor a shallow chamber beneath Katla, which probably indicatesthat the shallow chamber identified in other studies containssilica-rich magma rather than basalt. There is reasonably goodcorrelation between the depths of deep chambers (> 17 km)and geophysical estimates of Moho depth, indicating that magmaponds at the crust–mantle boundary. Shallow chambers (<7·1 km) are located in the upper crust, and probablyform at a level of neutral buoyancy. There are also discretechambers at intermediate depths (11 km beneath the rift zones),and there is strong evidence for cooling and crystallizing magmabodies or pockets throughout the middle and lower crust thatmight resemble a crystal mush. The results suggest that themiddle and lower crust is relatively hot and porous. It is suggestedthat crustal accretion occurs over a range of depths similarto those in recent models for accretionary processes at mid-oceanridges. The presence of multiple stacked chambers and hot, porouscrust suggests that magma evolution is complex and involvespolybaric crystallization, magma mixing, and assimilation. KEY WORDS: Iceland rift zones; cotectic crystallization; pressure; depth; magma chamber; volcanic glass     

Structural investigation of the synthetic CaAn(PO4)2 (An = Th and Np) cheralite-like phosphates

The crystallographic structures of the synthetic cheralite, CaTh(PO₄)₂, and its homolog CaNp(PO₄)₂ have been investigated by X-ray diffraction at room temperature. Rietveld analyses showed that both compounds crystallize in the monoclinic system and are isostructural to monazite LnPO₄ (Ln = La to Gd). The space group is P2₁/n (I.T. = 14) with Z = 2. The refined lattice parameters of CaTh(PO₄)₂ are a = 6.7085(8) Å, b = 6.9160(6) Å, c = 6.4152(6) Å, and β = 103.71(1)° with best fit parameters R _wp = 4.87%, R _p = 3.69% and R _B = 3.99%. For CaNp(PO₄)₂, we obtained a = 6.6509(5) Å, b = 6.8390(3) Å, c = 6.3537(8) Å, and β = 104.12(6)° and R _wp = 6.74%, R _p = 5.23%, and R _B = 6.05%. The results indicate significant distortions of bond length and angles of the PO₄ tetrahedra in CaTh(PO₄)₂ and to a lesser extent in CaNp(PO₄)₂. The structural distortions were confirmed by Raman spectroscopy of CaTh(PO₄)₂. A comparison with the isostructural compounds LnPO₄ (Ln = Ce and Sm) confirmed that the substitution of the large rare earth trivalent cations with Ca²⁺ and Th⁴⁺ introduces a distortion of the PO₄ tetrahedra.     

In situ strength measurements on natural upper-mantle minerals

Using in situ strength measurements at pressures up to 10 GPa and at room temperature, 400, 600, and 700°C, we examined rheological properties of olivine, orthopyroxene, and chromian-spinel contained in a mantle-derived xenolith. Mineral strengths were estimated using widths of X-ray diffraction peaks as a function of pressure, temperature, and time. Differential stresses of all minerals increase with increasing pressure, but they decrease with increasing temperature because of elastic strain on compression and stress relaxation during heating. During compression at room temperature, all minerals deform plastically at differential stress of 4–6 GPa. During subsequent heating, thermally induced yielding is observed in olivine at 600°C. Neither orthopyroxene nor spinel shows complete stress relaxation, but both retain some stress even at 700°C. The strength of the minerals decreases in the order of chromian-spinel ≈ orthopyroxene > olivine for these conditions. This order of strength is consistent with the residual pressure of fluid inclusions in mantle xenoliths.     

Constraining the origin and evolution of confined turbidite systems: southern Cretan margin,Eastern Mediterranean Sea (34°30–36°N)

Bathymetric, 9.5-kHz long-range sidescan sonar (OKEAN), seismic reflection and sediment-core data are used in the analysis of two tectonic troughs south of Crete, Eastern Mediterranean Sea. Here, up to 1.2 s two-way travel time (TWTT) of strata have accumulated since the Middle Miocene in association with extension in the South Aegean region. The study area comprises >100-km- long by >25-km-wide basins filled by sediments subdivided into two seismic units: (1) an upper Unit 1 deposited in sub-basins which follow the present-day configuration of the southern Cretan margin; (2) a basal Unit 2, more than 500 ms (TWTT) thick, accumulated in deeper half-graben/grabens distinct from the present-day depocentres. Both units overlap a locally stratified Unit 3 comprising the pre-Neogene core complex of Crete and Gavdos. In this work, the interpreted seismic units are correlated with the onshore stratigraphy, demonstrating that denudation processes occurring on Crete and Gavdos in response to major tectonic events have been responsible for high sedimentation rates along the proximal southern Cretan margin. Consequently, topographically confined sedimentary units have been deposited south of Crete in the last 12 Ma, including turbidites and other mass-flow deposits fed by evolving transverse and axial channel systems. Surface processes controlling facies distribution include the direct inflow of sediment from alluvial-fan systems and incising mountain rivers onto the Cretan slope, where significant sediment instability processes occur at present. In this setting, seismic profiles reveal eight different types of stratigraphic contacts on basin-margin highs, and basinal areas show evidence of halokinesis and/or fluid escape. The acquired data also show that significant changes to the margin’s configuration occurred in association with the post-Alpine tectonic and eustatic episodes affecting the Eastern Mediterranean.