Climatic changes in central and eastern Canada inferred from deep borehole temperature data

We analyzed data from 23 boreholes at 19 sites in central and eastern Canada, for the purpose of estimating ground surface temperature (GST) histories. These boreholes were logged down to at least 550 m depth with thermistor probes. Thermal conductivity measurements had been previously made at small depth intervals for the entire depth ranges of most of the boreholes. The temperature profiles of these boreholes do not indicate water disturbance. We estimated terrain effects for each borehole using a time dependent solid-angle method. The thermal perturbations caused by lakes or deforestation near the borehole sites are insignificant in most cases. However, four of the holes were found to be severely influenced by terrain effects. GSTs estimated from the borehole data less influenced by the terraineffects form two groups. The first group, which are generally from data of better quality, show a cold period near the end of the last century before the recent warming trend; the second show it 80–100 years earlier. We consider the former typical of the climate of the Boreal climatic region of Canada. The difference between the two groups may reflect the spacial variability of the climate. Four GST estimates do not belong to either type, and the reasons are discussed.     

Rift-related Jurassic basaltic phreatomagmatic volcanism in the central Transantarctic Mountains: precursory stage to flood-basalt effusion

 The Middle Jurassic Kirkpatrick flood basalts and comagmatic Ferrar intrusions in the Transantarctic Mountains represent a major pulse of tholeiitic magmatism related to early stages in the breakup of Gondwana. A record of the volcano-tectonic events leading to formation of this continental flood-basalt province is provided by strata underlying and only slightly predating the Kirkpatrick lavas. In the central Transantarctic Mountains, the lavas rest on widespread (≥7500 km²) tholeiitic pyroclastic deposits of the Prebble Formation. The Prebble Formation is dominated by lahar deposits and is an unusual example of a regionally developed basaltic lahar field. Related, partly fault-controlled pyroclastic intrusions cut underlying strata, and vents are represented by the preserved flanks of two small tephra cones associated with a volcanic neck. Lahar and air-fall deposits typically contain 50–60% accidental lithic fragments and sand grains derived from underlying Triassic – Lower Jurassic strata in the upper part of the Beacon Supergroup. Juvenile basaltic ash and fine lapilli consist of nonvesicular to scoriaceous tachylite, sideromelane, and palagonite, and have characteristics indicating derivation from hydrovolcanic eruptions. The abundance of accidental debris from underlying Beacon strata points to explosive phreatomagmatic interaction of basaltic magma with wet sediment and groundwater, which appears to have occurred in particular where rising magma intersected upper Beacon sand aquifers. Composite clasts in the lahar deposits exhibit complex peperitic textures formed during fine-scale intermixing of basaltic magma with wet sand and record steps in subsurface fuel-coolant interactions leading to explosive eruption. The widespread, sustained phreatomagmatic activity is inferred to have occurred in a groundwater-rich topographic basin linked to an evolving Jurassic rift zone in the Transantarctic Mountains. Coeval basaltic phreatomagmatic deposits of the Mawson and Exposure Hill Formations, which underlie exposures of the Kirkpatrick Basalt up to 1500 km to the north along strike in Victoria Land, appear to represent other parts of a regional, extension-related Middle Jurassic phreatomagmatic province which developed immediately prior to rapid outpouring of the flood basalts. This is consistent with models which assign an important role to lithospheric stretching in the generation of flood-basalt provinces. Received: 28 August 1995 / Accepted: 18 April 1996     

Trace Element Partitioning and Accessory Phase Saturation during H2O-Saturated Melting of Basalt with Implications for Subduction Zone Chemical Fluxes

Experimental phase equilibrium and trace element partitioningdata are reported for H₂O-saturated mid-ocean ridge basalt at2·5 GPa, 750–900°C and oxygen fugacities atthe nickel–nickel oxide buffer. Garnet, omphacite andrutile are present at all temperatures. Amphibole and epidotedisappear as residual phases above 800°C; allanite appearsabove 750°C. The Na–Al-rich silicate glass presentin all run products is likely to have quenched from a supercriticalliquid. Trace element analyses of glasses demonstrate the importantcontrol exerted by residual minerals on liquid chemistry. Inaddition to garnet, which controls heavy rare earth elements(HREE) and Sc, and rutile, which controls Ti, Nb and Ta, allanitebuffers the light REE (LREE; La–Sm) contents of liquidsto relatively low levels and preferentially holds back Th relativeto U. In agreement with previous experimental and metamorphicstudies we propose that residual allanite plays a key role inselectively retaining trace elements in the slab during subduction.Experimental data and analyses of allanite-bearing volcanicrocks are used to derive a model for allanite solubility inliquids as a function of pressure, temperature, anhydrous liquidcomposition and LREE content. The large temperature dependenceof allanite solubility is very similar to that previously determinedfor monazite. Our model, fitted to 48 datapoints, retrievesLREE solubility (in ppm) to within a factor of 1· 40over a pressure range of 0–4 GPa, temperature range of700–1200°C and for liquids with anhydrous SiO₂ contentsof 50–84 wt %. This uncertainty in LREE content is equivalentto a temperature uncertainty of only ± 27°C at 1000K, indicating the potential of allanite as a geothermometer.Silicic liquids from either basaltic or sedimentary protolithswill be saturated in allanite except for Ca-poor protolithsor at very high temperatures. For conventional subduction geothermsthe low solubility of LREE (+ Th) in liquids raises questionsabout the mechanism of LREE + Th transport from slab to wedge.It is suggested either that, locally, temperatures experiencedby the slab are high enough to eliminate allanite in the residueor that substantial volumes of H₂O-rich fluids must pass throughthe mantle wedge prior to melting. The solubility of accessoryphases in fluids derived from subducted rocks can provide importantconstraints on subduction zone thermal structure. KEY WORDS: subduction; experimental petrology; allanite; solubility; supercritical liquid; eclogite     

Photoadaptation of carboxylating enzymes and photosynthesis during a spring bloom

Properties of the light saturation curve of photosynthesis and ribulose-1,5-bisphosphate carboxylase (RuBPC) activity are shown to change qualitatively in a natural population of marine phytoplankton during a spring bloom. Evidence is presented to show that these changes constitute photoadapative responses to increasing irradiance. As irradiance increased during the bloom, both the level of light-saturated photosynthesis (P_m) and the initial slope of the light saturation curve (α = photosynthetic efficiency) increased whether those parameters were normalized to chlorophyll a concentration (P_m^b, α^b) or to cell numbers (P_m^c, α^c). The magnitudes of these changes were such that I_k (= P_m/α, the photoadaptation parameter) did not change, but I_m, the light intensity at which photosynthesis becomes saturated, increased. RuBPC activity, both chlorophyll a (RuBPC^b) and cell number normalized (RuBPC^c), also increased during the bloom. We suggest that these adaptations were achieved by simultaneously increasing the number of photosynthetic units, proportionately decreasing the photosynthetic unit size, and increasing both the concentrations of the enzymes of the dark reactions and possibly also of photosynthetic electron transport components.We also observed diminished levels of photoinhibition in the high light adapted cells late in the bloom and have suggested that this was a consequence of the same suite of physiological changes.In situ carbon fixation per cell increased during the bloom whereas no change occurred in this parameter when normalized to chlorophyll a concentration. Although these photoadaptive responses thus permitted carbon to be fixed in situ more rapidly per cell, at a constant efficiency with respect to investment of energy in the photosynthetic apparatus, they did not result in a change in growth rate. Based on consideratios of the role of time scale in physiological adaptation, however, it is suggested that the observed alterations in photosynthesis with increasing irradiance might permit a cell to more rapidly fill an energy quota for division, possibly an advantage in a mixing environment in which energy is patchily distributed, both spatially and temporalyy.Phosphoenolpyruvate carboxylase activity when normalized to chlorophyll a (PEPC^b) did not change during the bloom while chlorophyll a normalized dark carbon fixation decreased sharply and was quantitatively small compared to PEPC^b. On this basis and considering that RuBPC^b increased during the bloom, it is suggested that, although PEPC may be involved in dark carbon fixation, its most important quantitative role is probably an indirect one in light dependent photosynthesis.We have also considered the relevance of laboratory results on photoadaptation to interpretations of field studies and have suggested that batch culture studies must be treated with caution but that turbidistat and semi-continuous methods provide reasonable simulations of natural conditions.     

Patterns of biomass-size spectra from oligotrophic waters of the Northwest Atlantic

The study of the distribution of biomass by size provides an ataxonomic approach for analyzing the structure of the pelagic ecosystem. However, empirical data regarding planktonic size-structure in offshore areas are scarce. Here, we report the results of a study of the planktonic biomass size-distribution at several stations located within two regions (Sargasso Sea and New England Seamounts Area) of the Northwest Atlantic. The biomass size-spectra covered a body-size range from bacteria to mesozooplankton and a depth range from the surface to 400 m. It is shown that the slope of the normalized biomass size-spectrum (NBS-spectrum) varies depending on whether volume or carbon units are used. The transformation from volume to carbon units makes the slope of the NBS-spectrum approximately 0.15 units more negative. The distribution of normalized-biomass by size was linear (plotted on a log–log scale) at all stations. The slopes of the NBS-spectra (volume scale) ranged from −0.96 to −1.01. There were no significant differences among the slopes of the NBS-spectra within either of the two areas studied. In addition, no significant differences were detected between the stations in the Sargasso Sea and those located in the New England Seamounts area. Apart from a tendency towards a decrease in the intercept of the normalized-biomass axis of the size-spectra in deeper waters, the NBS-spectra were also very similar through depth. The slopes of the NBS-spectra in biovolume units are not significantly different from −1.0 (p<0.01), and therefore are in agreement with Sheldon’s Linear Biomass Hypothesis. By contrast, the slopes of the NBS-spectra in carbon units are significantly different (p<0.01) from −1.0 (range=–1.09 to –1.17; all stations TOGETHER=–1.14) and their numerical values are in close agreement with Platt and Denman’s model and with the findings of [Rodríguez and Mullin, 1986a]) for the North Pacific Central Gyre. The results of this study support the hypothesis that the planktonic size-structure of offshore oligotrophic systems is a conservative property.     

Double-diffusive convection with a non-linear equation of state: Part I. The accurate conservation of properties in a two-layer system

If the equation of state is nonlinear, a given flux of heat across a double diffusive interface causes different buoyancy fluxes in the upper and lower layers. This results in different convective activity in the two layers and can lead to preferential entrainment across the interface in one direction (i.e. a migration of the interface). In this paper we derive the conservation equations for properties (e.g. heat and a solute) across a double diffusive interface between two well-mixed layers. A nondimensional measure of the entrainment across an interface and the most suitable choice for the buoyancy flux ratio are presented. Some surprising facts emerge. First, even for a linear equation of state and in the absence of direct entrainment across the interface, the flux of water across a finger interface is shown to be important. Second, for the heat-solute system, the heat balance equations for each well-mixed layer contain terms proportional to the heat of solution of the solute and the partial specific enthalpy of pure water in a seawater solution. Third, the rate of change of gravitational potential energy of the two-layer system is shown to have several extra terms in addition to the two commonly quoted major terms.     

An evaluation of systems for the integrated assessment of capture fisheries

Integrated fishery assessments are intended to take an integrated and wholistic view of fisheries management, to provide comprehensive systems-based evaluations of the sustainability of individual fisheries and their products. Using guidance derived from decision theory and management science, we develop 18 criteria to evaluate the likely effectiveness of nine integrated fishery assessment systems. The assessment systems do not consistently consider four key criteria: the socio-economic impacts of fisheries, the provision of food security, the local relevance of the assessment to the fishery under assessment, and the independent peer review of assessment outcomes. Other important areas of weakness include the lack of appropriate consultation processes relating to both the design of the assessment system and the outcomes of fishery assessments, and limitations in the quality of data and information used in assessments. Improvements will be needed in these areas across all the systems we evaluated if fishers, consumers, managers and other stakeholders are to be provided with consistently reliable assessment outcomes, and an assurance of the sustainability of capture fisheries.     

On matching the spectrum of Io: Variations in the photometric properties of sulfur-containing mixtures

The problem of comparing laboratory spectra of sulfur-containing binary mixtures with the spectrum of Io is discussed. For the satellite, the observable is the geometric albedo as a function of wavelength, whereas in the laboratory one often measures some other type of albedo. In a previous paper we demonstrated that for pure sulfur the multiplicative factor which converts the laboratory albedos to geometric albedos can be strongly wavelength dependent. The present paper demonstrates that this is also true for binary sulfur-containing mixtures. Furthermore, there is no universal conversion factor applicable to all binary mixtures, nor can the factor be interpolated for a particular mixture from the conversion factors of the two end members. The conversion factor is a function not only of the specific composition of a binary mixture, but of the relative particle size distributions of the two components, and must be measured specifically for each individual sample if a quantitative comparison between a laboratory sample and Io's surface is desired.     

Analysis of stellar occultation data: Effects of photon noise and initial conditions

A new inversion technique for obtaining temperature, pressure, and number density profiles of a planetary atmosphere from an occultation light curve is described. This technique employs an improved boundary condition to begin the numerical inversion and permits the computation of errors in the profiles caused by photon noise in the light curve. We present our assumptions about the atmosphere, optics, and noise and develop the equations for temperature, pressure, and number density and their associated errors. By inverting in equal increments of altitude, Δh, rather than in equal increments of time, Δt, the inversion need not be halted at the first negative point on the light curve as required by previous methods. The importance of the boundary condition is stressed, and a new initial condition is given. Numerical results are presented for the special case of inversion of a noisy isothermal light curve. From these results, simple relations are developed which can be used to predict the noise quality of an occultation. It is found that fractional errors in temperature profiles are comparable to those of pressure and number density profiles. An example of the inversion method is shown. Finally, we discuss the validity of our assumptions. In an appendix we demonstrate that minimum fractional errors in scale height determined from the inversion are comparable to those from an isothermal fit to a noisy isothermal light curve.