Investigating the mechanisms leading to the deglaciation of past continental northern hemisphere ice sheets with the CLIMBER–GREMLINS coupled model

A coupling procedure between a climate model of intermediate complexity (CLIMBER-2.3) and a 3-dimensional thermo-mechanical ice-sheet model (GREMLINS) has been elaborated. The resulting coupled model describes the evolution of atmosphere, ocean, biosphere, cryosphere and their mutual interactions. It is used to perform several simulations of the Last Deglaciation period to identify the physical mechanisms at the origin of the deglaciation process. Our baseline experiment, forced by insolation and atmospheric CO₂, produces almost complete deglaciation of past northern hemisphere continental ice sheets, although ice remains over the Cordilleran region at the end of the simulation and also in Alaska and Eastern Siberia. Results clearly demonstrate that, in this study, the melting of the North American ice sheet is critically dependent on the deglaciation of Fennoscandia through processes involving switches of the thermohaline circulation from a glacial mode to a modern one and associated warming of the northern hemisphere. A set of sensitivity experiments has been carried out to test the relative importance of both forcing factors and internal processes in the deglaciation mechanism. It appears that the deglaciation is primarily driven by insolation. However, the atmospheric CO₂ modulates the timing of the melting of the Fennoscandian ice sheet, and results relative to Laurentide illustrate the existence of threshold CO₂ values, that can be translated in terms of critical temperature, below which the deglaciation is impeded. Finally, we show that the beginning of the deglaciation process of the Laurentide ice sheet may be influenced by the time at which the shift of the thermohaline circulation from one mode to the other occurs.     

Static compression of α-MnS at 298 K to 21 GPa

To investigate the equation of state of -MnS at high pressure and the possibility of a phase transition, the compression curve was measured at 298 K from 0 to 21 GPa using powder x-ray diffraction with a diamond anvil cell. The compression data are fit to a thirdorder Birch-Murnaghan equation of state, with parameters K ₀ = 72(2) GPa and K ₀ = 4.2(13). To compare present results with previous work, the data sets from three previous investigations (Clendenen and Drickamer 1966; Wakabayashi et al. 1968; Kraft and Greuling 1988) are refit to a Birch-Murnaghan equation of state. In the low pressure region (P < 10=" gpa),=" the=" results=" of=" clendenen=" and=" drickamer=" (1966)=" agree=" with=" the=" present=" data;=" however=" the=" results=" of=" wakbayashi=" et=" al.=" (1968)=" differ=" by=" more=" than=" 10%.=" a=" greater=" discrepancy=" between=" the=" present=" and=" previous=" results=" occurs=" above=" 10=" gpa.=" kraft=" and=" greuling=" (1988)=" reported=" a=" structure=" transition=" at=" 7=" gpa,=" and=" clendenen=" and=" drickamer=" (1966)=" observed=" a=" structure=" distortion=" at=" approximately=" 10=" gpa;=" the=" present=" data=" show=" no=" evidence=" of=" either=" transition,=" and=" are=" well=" fit=" by=" a=" single=" equation=" of=" state=" from=" 0=" to=" 21=" gpa.=" nonhydrostatic=" stress=" is=" discussed=" as=" one=" possibility=" for=" the=">     

Particulate barium fluxes and export production in the northwestern Mediterranean

Biogenic barium, mostly in the barite (BaSO₄) form, has been proposed as a tracer for export production in the ocean. Here we report on biogenic barium (Ba_xs) and particulate organic carbon (POC) fluxes from sediment traps deployed at the DYFAMED site in the Northwestern Mediterranean Sea. Ba_xs fluxes display average values of 37 ± 45 and 50 ± 58 μg/m²/d at 200 and 1000 m respectively, and are linearly correlated to POC fluxes (mean values of 7.9 ± 9.3 and 6.8 ± 6.8 mg C/m²/d at 200 and 1000 m). Export production estimates, calculated using published Ba_xs- or POC-based algorithms, all fall below or close to the lower limit of potential export values proposed in the literature. This work clearly demonstrates the usefulness of Ba_xs as a tracer of oceanic export production in the Northwestern Mediterranean Sea. However, development of a quantitative export production proxy requires a clear understanding of the underlying cause(s) for the observed spatial variations in the relationship between Ba_xs and POC fluxes. The present study confirms that the processes leading to barite formation differ between margin and open-ocean sites and probably account for much of the regional variability in the POC/Ba_xs ratio.     

Authigenic carbonates from methane seeps of the Congo deep-sea fan

Submersible investigations with the ROV Victor 6000 of some pockmark structures on the seafloor of the Congo deep-sea fan have shown that they are active venting sites of methane-rich fluids, associated with abundant fauna and carbonate crusts. Moreover, methane hydrates have been observed both outcropping and deep in the sediments in the centre of the “Regab” giant pockmark. Authigenic carbonates, mostly calcite sometimes mixed with aragonite, are cementing the sedimentary matrix components and fauna; diatoms are abundant but only as moulds, indicating that biogenic silica dissolution occurred in situ synchronous with carbonate precipitation. The occurrence of diagenetic barite and pyrite in some carbonate crusts demonstrates that they can be formed either within the sulphate/methane transition zone or deeper in sulphate-depleted sediments. The oxygen isotopic compositions of the diagenetic carbonates (3.17–6.01‰ V-PDB) indicate that precipitation occurred with bottom seawater mixed with a variable contribution of water from gas hydrate decomposition. The very low carbon isotopic compositions of the diagenetic carbonates (−57.1 to −27.75‰ V-PDB) demonstrate that carbon derives mostly from the microbial oxidation of methane.     

Relative microbial activity and mineralization associated with water masses in the lower Cook Inlet,Alaska

A study of the microbial uptake and respiration (mineralization) of radioactive glutamic acid was made in 141 surface water samples in the Cook Inlet, Alaska during three sampling periods (October, 1976; April, 1977; and November, 1977). The rate at which the test substrate was incorporated into cell material plus that respired as CO₂ was used to calculate the relative microbial activity. A northern water mass with salinities less than 31 showed high rates of relative microbial activity but low mineralization rates as indicated by respiration of¹⁴CO₂. The reverse pattern was seen in the water mass (salinity higher than 31) to the south which was more typical of coastal water. These data suggest that the organic nutrients in these two water masses are quite different either qualitatively or quantitatively. This study illustrates how the observed measurement of relative microbial activity might be helpful in defining specific water masses.(The publication of this article was unintentionally delayed owing to a fault in secretarial routine.)Published as Technical Paper No. 6009, Oregon Agricultural Experiment Station.     

Lake Sediment Core Records of Sulphur Accumulation and Sulphur Isotopic Composition in Central Ontario, Canada Lakes

Stable isotopic compositions and concentrations of total sedimentary sulphur (S) were determined in cores from 6 lakes in the acid-sensitive Muskoka-Haliburton region of south-central Ontario. The isotopic composition of S in deep sediment (> ~ 20 cm) was approximately constant in all lakes, and indicated a pre-industrial δ ³⁴S value between +4.0 and +5.3‰, which is similar to current bulk deposition. Similarly, total S concentrations in deep sediment were relatively low (1.9–5 mg S g⁻¹ dwt) and approximately constant with depth within cores. All lakes exhibited up-core increases in total S and decreases in δ ³⁴S at a depth corresponding to the beginning of industrialization in the Great Lakes region ( ~ 1900), resulting in a generally reciprocal depth pattern between total S concentration and δ ³⁴S ratios. While initial shifts in total S and δ ³⁴S were likely due to enhanced SO₄ reduction of newly available anthropogenic SO₄, both the magnitude and pattern of up-core S enrichment and shifts in δ ³⁴S varied greatly among lakes, and did not match changes in S deposition post 1900. Differences between lakes in total S and δ ³⁴S were not related to any single hydrologic (e.g., residence time) or physical (e.g., catchment-area-to-lake area ratio) lake characteristic. This work indicates that sediment cores do not provide consistent records of changes in post-industrial S deposition in this region, likely due to redox-related mobility of S in upper sediment.     

3D numerical modelling of fluid flow in the Val-d’Or orogenic gold district: major crustal shear zones drain fluids from overpressured vein fields

Fluid flow patterns have been determined using oxygen isotope isopleths in the Val-d’Or orogenic gold district. 3D numerical modelling of fluid flow and oxygen isotope exchange in the vein field shows that the fluid flow patterns can be reproduced if the lower boundary of the model is permeable, which represents middle or lower crustal rocks that are infiltrated by a metamorphic fluid generated at deeper levels. This boundary condition implies that the major crustal faults so conspicuous in vein fields do not act as the only major channel for upward fluid flow. The upper model boundary is impermeable except along the trace of major crustal faults where fluids are allowed to drain out of the vein field. This upper impermeable boundary condition represents a low-permeability layer in the crust that separates the overpressured fluid from the overlying hydrostatic fluid pressure regime. We propose that the role of major crustal faults in overpressured vein fields, independent of tectonic setting, is to drain hydrothermal fluids out of the vein field along a breach across an impermeable layer higher in the crust and above the vein field. This breach is crucial to allow flow out of the vein field and accumulation of metals in the fractures, and this breach has major implications for exploration for mineral resources. We propose that tectonic events that cause episodic metamorphic dehydration create a short-lived pulse of metamorphic fluid to rise along zones of transient permeability. This results in a fluid wave that propagates upward carrying metals to the mineralized area. Earthquakes along crustal shear zones cause dilation near jogs that draw fluids and deposit metals in an interconnected network of subsidiary shear zones. Fluid flow is arrested by an impermeable barrier separating the hydrostatic and lithostatic fluid pressure regimes. Fluids flow through the evolving and interconnected network of shear zones and by advection through the rock matrix. Episodic breaches in the impermeable barrier along the crustal shear zones allow fluid flow out of the vein field.     

Validation of a fish farm waste resuspension model by use of a particulate tracer discharged from a point source in a coastal environment

To validate a resuspension model of particulate material (salmonid farm wastes), a UV fluorescent particle tracer was selected with similar settling characteristics. Tracer was introduced to the seabed (water depth ≈30 m) and sediment samples taken on days 0, 3, 10, 17 and 30 to measure the horizontal and vertical distribution of tracer in sediments. A concentric sampling grid was established at radii of 25, 50, 100, 150, 200, 400, 700 and 1, 000 m from the source on transects 30° apart. The bulk of the deployed tracer was initially concentrated in an area 25 m radius from the release point; tracer was observed to steadily decrease to zero over a period of 30 days. In a 200 m region measured from the release point in the direction of the residual current, the redeposition of tracer was low. A Lagrangian particle tracking model was validated using these observed data by varying resuspension model parameters within limits to obtain the best agreement between spatial and temporal distributions. The validated model generally gave good predictions of total mass budgets (±7% of total tracer released), particulary where tracer concentrations were high near the release point. Best fit model parameters (critical erosion shear stress=0.018 N m⁻², erodibility constan=60 g m⁻² d⁻¹) are at the low end of reported parameters for coastal resuspension models. Such a low critical erosion shear stress indicates that the frequency of resuspension and deposition events for freshly deposited material is high.     

Comparison of the pseudo-static and dynamic behaviour of gravity retaining walls

Summary Pseudo-static and dynamic non-linear finite element analyses have been performed to assess the dynamic behaviour of gravity retaining walls subjected to horizontal earthquake loading. In the pseudo-static analysis, the peak ground acceleration is converted into a pseudo-static inertia force and applied as a horizontal incremental gravity load. In the dynamic analysis, an actual measured earthquake acceleration time history has been scaled to provide peak ground acceleration values of 0.1 g and 0.3 g. Good agreement is obtained between the pseudo-static analysis and analytical methods for the calculation of the active coefficient of earth pressure. However, the results from the dynamic analysis require careful interpretation. In the pseudo-static analysis, the increase in the point of application of the resultant active force with the horizontal earthquake coefficient k _h from the one-third point to the mid-height of the wall is clearly observed. In the dynamic analysis, the variation in the point of application is shown to be a function of the type of wall deformation. Both finite element analyses indicate the importance of determining the magnitude of the predicted displacements when assessing the behaviour of the wall to seismic loading.