Fluid-rock interaction between syenites and marbles at Stephen Cross Quarry,Québec,Canada: petrological and stable isotope data

Grenville dolomitic marbles and calc-silicates at Stephen Cross Quarry, Québec, underwent contact metamorphism and metasomatism associated with the intrusion of the Wakefield syenite at ambient pressures of 0.4GPa at 1090–1070Ma. Fluid infiltration produced exoskarns, calcite+periclase+forsterite±diopside±orthoclase assemblages in the marbles, and quartz±calcite±wollastonite±diopside±anorthite assemblages in the calc-silicates. Phase-equilibria in the CaO–MgO–Al₂O₃–SiO₂–H₂O–CO₂ system suggest that fluid infiltration occurred close to the thermal peak of contact metamorphism (715–815°C) and that the fluids hadXCO₂0.15. In the metasediments, ¹⁸O values of calcite (Cc) are as low as 8.6, suggesting that the fluids were in isotopic equilibrium with the syenites (¹⁸O =8.8–10.2). Marble ¹³C(Cc) values are-0.1 to-3.2; the lack of correlation between ¹³C(Cc) and ¹⁸O(Cc) is consistent with the infiltration of water-rich fluids. The resetting of stable isotopes and the mineralogical changes can be explained by time-integrated fluid fluxes of up to 110 m³/m² (4×10⁶ mol/m²), corresponding to actual fluxes of 3×10^-11 to 3×10^-12 m³/m²-s and intrinsic permeabilities of 10^-18 to 10^-20 m² for fluid flow lasting 0.1-1Ma. Marble ¹⁸O(Cc) values do not correlate well with distance from the syenite, and fluids were probably channelled across lithological layering. The correlation between the degree of resetting of marble ¹⁸O(Cc) values with the abundance of submillimetre-wide veins, suggests that fluid focussing may have resulted from variations in fracture density. Late, lower temperature (<500°C), fluid flow formed serpentine (Serp) and brucite (Br) from periclase and forsterite. ¹⁸O(Br) and ¹⁸O(Scrp) values correlate with ¹⁸O(Cc), suggesting that retrogression involved only limited volumes of fluid. The observation that ¹⁸O(Cc-Br) and ¹⁸O(Cc-Serp) values are higher in marbles that have lower ¹⁸O(Cc) values is interpreted as indicating that fluid flow persisted to lower temperatures in those rocks due to higher intrinsic permeabilities. Calcite in the syenite was also formed by the influx of fluids during cooling. Syenite ¹⁸O(Cc) values are approximately in isotopic equilibrium with the high-temperature silicate minerals, suggesting that again only minor volumes of fluid were involved. In detail fluid flow was prolonged and complex, creating problems for the application of quantitative fluid flow models.This paper is a contribution to IGCP 304, Lower Crustal Processes     

A powder neutron diffraction study of the crystal structure of the fluoroperovskite NaMgF<Subscript>3</Subscript> (neighborite) from 300 to 3.6 K

The cell dimensions and crystal structures of the fluoroperovskite NaMgF₃ (neighborite), synthesized by solid state methods, have been determined by powder neutron diffraction and Rietveld refinement over the temperature range 300–3.6 K using Pt metal as an internal standard for calibration of the neutron wavelength. These data show that Pbnm NaMgF₃ does not undergo any phase transitions to structures of lower symmetry with decreasing temperature. The cell dimensions and atomic coordinates together with polyhedron volumes and distortion indices are given for Pbnm NaMgF₃ at 25 K intervals from 300 to 3.6 K. Decreases in the a and c cell dimensions reach a saturation point at 50 K, whereas the b dimension becomes saturated at 150 K. The distortion of the structure of Pbnm NaMgF₃ from the aristotype cubic structure is described in terms of the tilting of the MgF₆ octahedra according to the tilt scheme a ⁻ a ⁻ c ⁺ . With decreasing temperature the antiphase tilt (a ⁻) increases from 14.24° to 15.39°, whereas the in-phase tilt (c ⁺ ) remains effectively constant at ∼10.7°. Changes in the tilt angles are insufficient to cause changes in the coordination sphere of Na that might induce a low temperature phase transition. The structure of Pbnm NaMgF₃ is also described in terms of normal mode analysis and displacements of the condensed normal modes are compared with those of Pbnm KCaF₃.     

Geology and geochronology of the Betara region in south-southeastern Brazil: Evidence for possible Statherian (1.80–1.75 Ga) and Calymmian (1.50–1.45 Ga) extension events

The results of geological mapping, chemical analysis and radiometric dating of metabasic rocks of Betara Formation, and mapping and dating of those present in the Betara basement nucleus together with mylonitic granodiorite and syenogranite are reported here. U–Pb analysis of bulk zircon fractions from the metabasic rocks of the basement nucleus yielded a Statherian age of 1790 ± 22 Ma, while the metabasic rocks from the upper part of the Betara Formation yielded a Calymmian age between 1500 and 1450 Ma. This age is a minimum for the deposition of the Betara Formation. The older metabasic rocks are associated with post-tectonic, possibly anorogenic syenogranite, while the younger ones are gabbro or very porphyritic ankaramite whose REE patterns are consistent with crystallization from an N-MORB parent magma. The observations and data point to the probable events associated with extensional processes of the end of Paleoproterozoic and early Mesoproterozoic. Similar registers of Statherian (1.80–1.75 Ga) and Calymmian (1.50–1.45 Ga) extensional events are recorded in other parts of the South American and African continents. The Neoproterozoic witnessed the formation and junction of the tectonic slices which formed the Apiaí domain during the assemblage of western Gondwana.     

Paleogene of the Huanghekou Sag in the Bohai Bay Basin,NE China: deposition–erosion response to a slope break system of rift lacustrine basins

The sequence architecture and depositional systems of the Paleogene lacustrine rift succession in the Huanghekou Sag, Bohai Bay Basin, NE China were investigated based on seismic profiles, combined with well log and core data. Four second‐order or composite sequences and seven third‐order sequences were identified. The depositional systems identified in the basin include: fan delta, braid delta, meander fluvial delta, lacustrine and sublacustrine fan. Identification of the slope break was conducted combining the interpretation of faults of each sequence and the identification of syndepositional faults, based on the subdivision of sequence stratigraphy and analysis of depositional systems. Multiple geomorphologic units were recognized in the Paleogene of the Huanghekou Sag including faults, flexures, depositional slope break belts, ditch‐valleys and sub‐uplifts in the central sag. Using genetic division principles and taking into consideration tectonic features of the Paleogene of the Huanghekou Sag, the study area was divided into the Northern Steep Slope/Fault Slope Break System, the Southern Gentle Slope Break System and T10 Tectonic Slope Break System/T10 Tectonic Belt. Responses of slope break systems to deposition–erosion are shown as: (1) basin marginal slope break is the boundary of the eroded area and provenance area; (2) ditch‐valley formed by different kinds of slope break belts is a good transport bypass for source materials; (3) shape of the slope break belt of the slope break system controls sediments types; (4) the ditch‐valley and sub‐sag of a slope break system is an unloading area for sediments; and (5) due to their different origins, association characteristics and developing patterns, the Paleogene slope break belt systems in the Huanghekou Sag show different controls on depositional systems. The Northern Fault Slope Break system controls the deposition of a fan delta‐lacustrine‐subaqueous fan, the Southern Gentle Slope Break system controls the deposition of a fluvial–deltaic–shallow lacustrine and sublacustrine fan, and the T10 Tectonic Slope Break System controls the deposition of shallow lacustrine beach bar sandbodies. The existence of a slope break system is a necessary but not a sufficient condition for studying sandbody development. The formation of effective sandbodies along the slope break depends on the reasonable coupling of effective provenance, necessary association patterns of slope break belt, adequate unloading space and creation of definite accommodation space. Copyright © 2013 John Wiley & Sons, Ltd.     

Atmospheric deposition of Be and Be in New Zealand rain (1996-98)

Beryllium isotope concentrations were determined in monthly rainfall collections at three sites across New Zealand (36 to 45° S), from October 1996 to November 1998. At the northern sites of Leigh (near Auckland) and Gracefield (near Wellington), ⁷Be and ¹⁰Be concentrations are relatively constant at 1.2 to 1.4 × 10⁷ atoms kg⁻¹ rain and 2.1 to 2.6 × 10⁷ atoms kg⁻¹ rain, respectively. These concentrations correspond to annual flux rates of ∼15 and ∼27 × 10⁹ atoms m⁻² y⁻¹, respectively. At the southern site of Dunedin, concentrations are similar to those at the northern sites, but flux rates are significantly lower at ∼ 9 and ∼19 × 10⁹ atoms m⁻² y⁻¹, respectively, because of lower average rainfall east of the Southern Alps mountain range. These results are broadly similar to those reported by Brown et al. (1989) and Knies et al. (1994) for rain from midlatitude sites in the USA sampled from 1986 to 1994. Mean ⁷Be/¹⁰Be ratios for New Zealand (0.47 to 0.61) are, however, significantly lower than for the USA (0.69 to 0.78), due in part to the addition of ¹⁰Be from re-suspended dust. Subtraction of the dust component increases the New Zealand ⁷Be/¹⁰Be ratios to 0.70 (Leigh), 0.65 (Gracefield) and 0.50 (Dunedin). The adjusted results provide evidence for transfer of older stratospheric air to the troposphere in late-spring-summer, an effect which is strongest in the north. The overall reduction of ⁷Be/¹⁰Be from north to south implies an increase in residence time from ∼ 80 to ∼100 d for Be isotopes in the atmosphere above New Zealand.     

The stability of a remediated bed in Hamilton Harbour,Lake Ontario,Canada

ABSTRACT In situ measurements of lakebed sediment erodibility were made on three sites in Hamilton Harbour, Lake Ontario, using the benthic flume Sea Carousel. Three methods of estimating the surface erosion threshold (τ_c(0)) from a Carousel time series were evaluated: the first method fits measures of bed strength to eroded depth (the failure envelope) and evaluates threshold as the surface intercept; the second method regresses mean erosion rate (E_m) with bed shear stress and solves for the floc erosion rate (E_f) to derive the threshold for E_m = E_f = 1 × 10^?5 kg m^?2 s^?1; the third method extrapolates a regression of suspended sediment concentration (S) and fluid transmitted bed shear stress (τ₀) to ambient concentrations. The first field site was undisturbed (C) and acted as a control; the second (W) was disturbed through ploughing and water injection as part of lakebed treatment, whereas the third site (OIP) was disturbed and injected with an oxidant used for remediation of contaminated sediment. The main objectives of this study were: (1) to evaluate the three different methods of deriving erosion threshold; (2) to compare the physical behaviour of lacustrine sediments with their marine estuarine counterparts; and (3) to examine the effects of ploughing and chemical treatment of contaminated sediment on bed stability. Five deployments of Sea Carousel were carried out at the control site. Mean erosion thresholds for the three methods were: τ_c(0) = 0·5 (±0·06), 0·27 (±0·01) and 0·34 (±0·03) Pa respectively. Method 1 overpredicted bed strength as it was insensitive to effects in the surface 1–2 mm, and the fit of the failure envelope was also highly subjective. Method 2 exhibited a wide scatter in the data (low correlation coefficients), and definition of the baseline erosion rate (E_f) is largely arbitrary in the literature. Method 3 yielded stable (high correlation coefficients), reproducible and objective results and is thus recommended for evaluation of the erosion threshold. The results of this method correlated well with sediment bulk density and followed the same trend as marine counterparts from widely varying sites. Mass settling rates, expressed as a decay constant, k, of S(t), were strongly related to the maximum turbidity at the onset of settling (S_max) and were also in continuity with marine counterparts. Thus, it appears that differences in salinity had little effect on mass settling rates in the examples presented, and that biological activity dominated any effects normally attributable to changes in salinity. Bedload transport of eroded aggregates (2–4 mm in diameter) took place by rolling below a mean tangential flow velocity (U_y) of 0·32 ms^?1 and by saltation at higher velocities. Mass transport as bedload was a maximum at U_y = 0·4 ms^?1, although bedload never exceeded 1% of the suspended load. The proportion of material moving as bedload was greatest at the onset of erosion but decreased as flow competence increased. Given the low bulk density and strength of the lakebed sediment, the presence of a bedload component is notable. Bedload transport over eroding cohesive substrates should be greater in estuaries, where both sediment density and strength are usually higher. Significant differences between the ploughed and control sites were apparent in both the erosion rate and the friction coefficient (φ), and suggest that bed recovery after disruption is rapid (< 24 h). τ_c(0) increased linearly with time after ploughing and recovered to the control mean value within 3 days. The friction coefficient was reduced to zero by ploughing (diagnostic of fluidization), but increased linearly with time, regaining control values within 6 days. No long‐term reduction in bed strength due to remediation was apparent.     

Submarine channel levee shape and sediment waves from physical experiments

Submarine channel levees aggrade through repeated overspill events from the channel axis. The shape of the levees may therefore reflect some characteristic(s) of the overspilling flow. It has been noted that basin floor levees typically have a relatively low-relief and taper exponentially to their termination; in contrast slope channel levees may be much steeper close to the channel. A simple physical experiment was performed where a surge-like sediment-laden current flowed through a curved channel. Significant overspill occurred and generated a deposit flanking the channel on either side. The experiment was repeated 25 times to build up low-relief channel-levees. It was found that in proximal areas, levees were steep and characterised by power-law decays, a transitional zone of logarithmically thinning levee was found a little further down-channel, followed by exponential decays in medial to distal areas. The style of levee decay is a function of spatial variation in overbank sedimentation rates. Where flows rapidly lose momentum and deposit across the grain-size spectrum, i.e., in proximal areas, levees tend to be steep; farther down the channel, the steep levee slope gives way to a more gradually tapering deposit. In more distal parts of the channel, deposition is directly related to sediment settling velocity (rather than the suspended load exceeding flow transport capacity as is the case in proximal areas), the deposit reflects this with relatively simple exponential thickness decays. Additionally, small-scale sediment waves developed under lee wave conditions on the inner-bend overbank. The waves initially migrated slightly towards the channel, but as the style of overspill evolved due to intra-channel deposition, flows moved out of the lee wave window and sedimentation became out of phase with the wavelength of the features and the topography was healed.     

Mass-independent fractionation of oxygen isotopes in the mesostasis of a chondrule from the Semarkona LL3.0 ordinary chondrite

A large chondrule from Semarkona, the most primitive ordinary chondrite known, has been discovered to contain a record of mass transport during its formation. In most respects, it is a normal Type I, group A1, low-FeO chondrule that was produced by reduction and mass-loss during the unidentified flash-heating event that produced the chondrules, the most abundant structural component in primitive meteorites. We have previously measured elemental abundances and abundance profiles in this chondrule. We here report oxygen isotope ratio abundances and ratio abundance profiles. We have found that the mesostasis is zoned in oxygen isotope ratio, with the center of the chondrule containing isotopically heavier oxygen than the outer regions, the outer regions being volatile rich from the diffusion of volatiles into the chondrule during cooling. The δ¹⁷O values range from −2.0‰ to 9.9‰, while δ¹⁸O range from −1.9‰ to 9.6‰. More importantly, a plot of δ¹⁷O against δ¹⁸O has a slope of 1.1 ± 0.2 (1σ) and 0.88 ± 0.10 (1σ) when measured by two independent methods. Co-variation of δ¹⁷O with δ¹⁸O that does not follow mass-dependent fractionation has often been seen in primitive solar system materials and is usually ascribed to the mixing of different oxygen reservoirs. We argue that petrographic and compositional data indicate that this chondrule was completely melted at the time of its formation so that relic grains could not have survived. Furthermore, there is petrographic and compositional evidence that there was no aqueous alteration of this chondrule subsequent to its formation. Although it is possible to formulate a series of exchanges between the chondrule and external ¹⁶O-rich and ¹⁶O-poor reservoirs that may explain the detailed oxygen isotope systematics of this chondrule, such a sequence of events looks very contrived. We therefore hypothesize that reduction, devolatilization, and crystallization of the chondrule melt may have produced ¹⁶O-rich olivines and ¹⁶O-poor mesostasis plotting on a slope-one line as part of the chondrule-forming process in an analogous fashion to known chemical mass-independent isotopic fractionation mechanisms. During cooling, volatiles and oxygen near the terrestrial line in oxygen isotope composition produced the outer zone of volatile rich and ¹⁶O-rich mesostasis. The chondrule therefore not only retains a record of considerable mass transport accompanying formation, but also may indicate that the isotopes of oxygen underwent mass-independent fractionation during the process.     

Dolomitization of the Waulsortian Limestone (Lower Carboniferous) in the Irish Midlands

The Waulsortian Limestone (Lower Carboniferous) of the southern Irish Midlands is dolomitized pervasively over a much larger region than previous studies have documented. This study indicates a complex, multistage, multiple fluid history for regional dolomitization. Partially and completely dolomitized sections of Waulsortian Limestones are characterized by finely crystalline (0·01–0·3 mm) planar dolomite. Planar replacive dolomite is commonly followed by coarse (≥0·5 mm) nonplanar replacive dolomite, and pervasive void‐filling saddle dolomite cement is frequently associated with Zn–Pb mineralization. Planar dolomite has average δ¹⁸O and δ¹³C values (‰ PDB) of –4·8 and 3·9 respectively. These are lower oxygen and slightly higher carbon isotope values than averages for marine limestones in the Waulsortian (δ¹⁸O=–2·2, δ¹³C=3·7). Mean C and O isotope values of planar replacive dolomite are also distinct from those of nonplanar and saddle dolomite cement (–7·0 and 3·3; –7·4 and 2·4 respectively). Fluid inclusions indicate a complex history involving at least three chemically and thermally distinct fluids during dolomite cementation. The petrography and geochemistry of planar dolomites are consistent with an early diagenetic origin, possibly in equilibrium with modified Carboniferous sea water. Where the Waulsortian was exposed to hydrothermal fluids (70–280 °C), planar dolomite underwent a neomorphic recrystallization to a coarser crystalline, planar and nonplanar dolomite characterized by lower δ¹⁸O values. Void‐filling dolomite cement is isotopically similar to nonplanar, replacive dolomite and reflects a similar origin from hydrothermal fluids. This history of multiple stages of dolomitization is significantly more complex than earlier models proposed for the Irish Midlands and provides a framework upon which to test competing models of regional vs. localized fluid flow.     

Evidence for centennial-scale Lateglacial and early Holocene climatic complexity from Quoyloo Meadow,Orkney, Scotland

The influence of the North Atlantic on the margins of Europe means the region is particularly sensitive to changes in the ocean–atmospheric system. During the Last Glacial–Interglacial Transition (16–8 cal ka bp ) this system was repeatedly disrupted, leading to a series of abrupt and short-lived shifts in climate. Despite much research, the number and magnitude of these ‘centennial-scale’ events is not well understood. To address this, we expand upon investigations at Quoyloo Meadow, Orkney, Scotland, one of the best chronologically constrained palaeoclimate records in northern Britain. By coupling stable isotope and chironomid fossil analyses with existing data, this study identifies multiple phases of centennial-scale disturbance at: c. 14.0, 11.1, 10.8, 10.5, 10.45 and 10.3 cal ka bp , with the events at 14.0 and 10.3 exhibiting a particularly pronounced cold-climate signature. During the Holocene, the strongest response to climate forcing was at c. 10.3–10.0 cal ka bp , expressed as a two-stage drop in mean July temperatures, a shift in pollen spectra indicative of ‘less-stable’ climatic regimes, and a depletion in δ¹⁸O values. We interpret this as the first reliably dated incidence of the ‘10.3-ka event’ in the British Isles and consider the wider impact of this climatic reversal in other Holocene records.