Global Rates of Free Hydrogen (H2) Production by Serpentinization and other Abiogenic Processes within Young Ocean Crust

The main conclusion of this dissertation is that global H2 production within young ocean crust (<10 Mya) is higher than currently recognized, in part because current estimates of H2 production accompanying the serpentinization of peridotite may be too low (Chapter 2) and in part because a number of abiogenic H2-producing processes have heretofore gone unquantified (Chapter 3). The importance of free H2 to a range of geochemical processes makes the quantitative understanding of H2 production advanced in this dissertation pertinent to an array of open research questions across the geosciences (e.g. the origin and evolution of life and the oxidation of the Earth’s atmosphere and oceans).

The first component of this dissertation (Chapter 2) examines H2 produced within young ocean crust [e.g. near the mid-ocean ridge (MOR)] by serpentinization. In the presence of water, olivine-rich rocks (peridotites) undergo serpentinization (hydration) at temperatures of up to ~500°C but only produce H2 at temperatures up to ~350°C. A simple analytical model is presented that mechanistically ties the process to seafloor spreading and explicitly accounts for the importance of temperature in H2 formation. The model suggests that H2 production increases with the rate of seafloor spreading and the net thickness of serpentinized peridotite (S-P) in a column of lithosphere. The model is applied globally to the MOR using conservative estimates for the net thickness of lithospheric S-P, our least certain model input. Despite the large uncertainties surrounding the amount of serpentinized peridotite within oceanic crust, conservative model parameters suggest a magnitude of H2 production (~1012 moles H2/y) that is larger than the most widely cited previous estimates (~1011 although previous estimates range from 1010-1012 moles H2/y). Certain model relationships are also consistent with what has been established through field studies, for example that the highest H2 fluxes (moles H2/km2 seafloor) are produced near slower-spreading ridges (<20 mm/y). Other modeled relationships are new and represent testable predictions. Principal among these is that about half of the H2 produced globally is produced off-axis beneath faster-spreading seafloor (>20 mm/y), a region where only one measurement of H2 has been made thus far and is ripe for future investigation.

In the second part of this dissertation (Chapter 3), I construct the first budget for free H2 in young ocean crust that quantifies and compares all currently recognized H2 sources and H2 sinks. First global estimates of budget components are proposed in instances where previous estimate(s) could not be located provided that the literature on that specific budget component was not too sparse to do so. Results suggest that the nine known H2 sources, listed in order of quantitative importance, are: Crystallization (6x1012 moles H2/y or 61% of total H2 production), serpentinization (2x1012 moles H2/y or 21%), magmatic degassing (7x1011 moles H2/y or 7%), lava-seawater interaction (5x1011 moles H2/y or 5%), low-temperature alteration of basalt (5x1011 moles H2/y or 5%), high-temperature alteration of basalt (3x1010 moles H2/y or <1%), catalysis (3x108 moles H2/y or <<1%), radiolysis (2x108 moles H2/y or <<1%), and pyrite formation (3x106 moles H2/y or <<1%). Next we consider two well-known H2 sinks, H2 lost to the ocean and H2 occluded within rock minerals, and our analysis suggests that both are of similar size (both are 6x1011 moles H2/y). Budgeting results suggest a large difference between H2 sources (total production = 1x1013 moles H2/y) and H2 sinks (total losses = 1x1011 moles H2/y). Assuming this large difference represents H2 consumed by microbes (total consumption = 9x1011 moles H2/y), we explore rates of primary production by the chemosynthetic, sub-seafloor biosphere. Although the numbers presented require further examination and future modifications, the analysis suggests that the sub-seafloor H2 budget is similar to the sub-seafloor CH4 budget in the sense that globally significant quantities of both of these reduced gases are produced beneath the seafloor but never escape the seafloor due to microbial consumption.

The third and final component of this dissertation (Chapter 4) explores the self-organization of barchan sand dune fields. In nature, barchan dunes typically exist as members of larger dune fields that display striking, enigmatic structures that cannot be readily explained by examining the dynamics at the scale of single dunes, or by appealing to patterns in external forcing. To explore the possibility that observed structures emerge spontaneously as a collective result of many dunes interacting with each other, we built a numerical model that treats barchans as discrete entities that interact with one another according to simplified rules derived from theoretical and numerical work, and from field observations: Dunes exchange sand through the fluxes that leak from the downwind side of each dune and are captured on their upstream sides; when dunes become sufficiently large, small dunes are born on their downwind sides (“calving”); and when dunes collide directly enough, they merge. Results show that these relatively simple interactions provide potential explanations for a range of field-scale phenomena including isolated patches of dunes and heterogeneous arrangements of similarly sized dunes in denser fields. The results also suggest that (1) dune field characteristics depend on the sand flux fed into the upwind boundary, although (2) moving downwind, the system approaches a common attracting state in which the memory of the upwind conditions vanishes. This work supports the hypothesis that calving exerts a first order control on field-scale phenomena; it prevents individual dunes from growing without bound, as single-dune analyses suggest, and allows the formation of roughly realistic, persistent dune field patterns.

Microplastics Alter the Properties and Sinking Rates of Zooplankton Faecal Pellets

Plastic debris is a widespread contaminant, prevalent in aquatic ecosystems across the globe. Zooplankton readily ingest microscopic plastic (microplastic, < 1 mm), which are later egested within their faecal pellets. These pellets are a source of food for marine organisms, and contribute to the oceanic vertical flux of particulate organic matter as part of the biological pump. The effects of microplastics on faecal pellet properties are currently unknown. Here we test the hypotheses that (1) faecal pellets are a vector for transport of microplastics, (2) polystyrene microplastics can alter the properties and sinking rates of zooplankton egests and, (3) faecal pellets can facilitate the transfer of plastics to coprophagous biota. Following exposure to 20.6 μm polystyrene microplastics (1000 microplastics mL–1) and natural prey (∼1650 algae mL–1) the copepod Calanus helgolandicus egested faecal pellets with significantly (P < 0.001) reduced densities, a 2.25-fold reduction in sinking rates, and a higher propensity for fragmentation. We further show that microplastics, encapsulated within egests of the copepod Centropages typicus, could be transferred to C. helgolandicus via coprophagy. Our results support the proposal that sinking faecal matter represents a mechanism by which floating plastics can be vertically transported away from surface waters.

Vitamin B-6 deficiency in rats reduces hepatic serine hydroxymethyltransferase and cystathionine beta-synthase activities and rates of in vivo protein turnover, homocysteine remethylation and transsulfuration

Vitamin B-6 deficiency causes mild elevation in plasma homocysteine, but the mechanism has not been clearly established. Serine is a substrate in one-carbon metabolism and in the transsulfuration pathway of homocysteine catabolism, and pyridoxal phosphate (PLP) plays a key role as coenzyme for serine hydroxymethyltransferase (SHMT) and enzymes of transsulfuration. In this study we used [H-2(3)]serine as a primary tracer to examine the remethylation pathway in adequately nourished and vitamin B-6-deficient rats pi and 0.1 mg pyridoxine (PN)/kg diet]. [H-2(3)]Leucine and [1-C-13]methionine were also used to examine turnover of protein and methionine pools, respectively, All tracers were injected intraperitoneally as a bolus dose, and then rats were killed (n = 4/time point) after 30, 60 and 120 min. Rats fed the low-PN diet had significantly lower growth and plasma and liver PLP concentrations, reduced liver SHMT activity, greater plasma and liver total homocysteine concentration, and reduced liver S-adenosylmethionine concentration. Hepatic and whole body protein turnover were reduced in vitamin B-6-deficient rats as evidenced by greater isotopic enrichment of [H-2(3)]leucine. Hepatic [H-2(2)]methionine production from [H-2(3)]serine via cytosolic SHMT and the remethylation pathway was reduced by 80.6% in vitamin B-6 deficiency. The deficiency did not significantly reduce hepatic cystathionine-beta-synthase activity, and in vivo hepatic transsulfuration flux shown by production of [H-2(3)]cysteine from the [H-2(3)]serine increased over twofold. In contrast, plasma appearance of [H-2(3)]cysteine was decreased by 89% in vitamin B-6 deficiency. The rate of hepatic homocysteine production shown by the ratio of [1-C-13]homocysteine/[1-C-13]methionine areas under enrichment vs. time curves was not affected by vitamin B-6 deficiency. Overall, these results indicate that vitamin B-6 deficiency substantially affects one-carbon metabolism by impairing both methyl group production for homocysteine remethylation and flux through whole-body transsulfuration.

Carbon flow in an upland grassland: effect of liming on the flux of recently photosynthesized carbon to rhizosphere soil

The effect of liming on the flow of recently photosynthesized carbon to rhizosphere soil was studied using (CO2)-C-13 pulse labelling, in an upland grassland ecosystem in Scotland. The use of C-13 enabled detection, in the field, of the effect of a 4-year liming period of selected soil plots on C allocation from plant biomass to soil, in comparison with unlimed plots. Photosynthetic rates and carbon turnover were higher in plants grown in limed soils than in those from unlimed plots. Higher delta(13)C% values were detected in shoots from limed plants than in those from unlimed plants in samples clipped within 15 days of the end of pulse labelling. Analysis of the aboveground plant production corresponding to the 4-year period of liming indicated that the standing biomass was higher in plots that received lime. Lower delta(13)C% values in limed roots compared with unlimed roots were found, whereas no significant difference was detected between soil samples. Extrapolation of our results indicated that more C has been lost through the soil than has been gained via photosynthetic assimilation because of pasture liming in Scotland during the period 1990-1998. However, the uncertainty associated with such extrapolation based on this single study is high and these estimates are provided only to set our findings in the broader context of national soil carbon emissions.

Concentrations et flux d'azote dans les sédiments hypoxiques de l'Estuaire Maritime du St-Laurent.

Les sédiments sont des sites importants d’élimination d’azote (N) puisqu’ils possèdent des gradients d’oxydoréduction leur conférant les conditions idéales pour les réactions microbiennes de transformation de N. L’eutrophisation des régions côtières peut altérer ces gradients, par des changements dans la concentration d’oxygène (O2) de l’eau interstitielle, et modifier l’importance relative des processus transformant le N. Afin de mieux comprendre comment l’O2 pourrait influencer les transformations de N, nous avons mesuré les flux diffusifs de diazote (N2), nitrate (NO3-), oxygène et ammonium (NH4+) dans les sédiments de l’Estuaire Maritime du St-Laurent (EMSL), et nous avons estimé les taux de dénitrification. L’importance du couple nitrification-dénitrification en fonction d’un gradient de concentrations d’O2 dans la zone d’hypoxie de l’EMSL fut aussi évaluée. La concentration des gaz dissous fut mesurée en utilisant une nouvelle approche développée dans cette étude. Les flux diffusifs de N2, O2, NO3- et NH4+ variaient de 5.5 à 8.8, de -37.1 à -84.8, de -4.0 à -5.8 et de 0.6 à 0.8 μmol N m-2 h-1 respectivement. Les concentrations de N2 et NO3- dans l’eau porale et les flux de NO3- et de N2 des sédiments, suggèrent que la diffusion de NO3- provenant de l’eau à la surface des sédiments ne peut pas expliquer par elle-même la production de N2 observée. En utilisant une approche stoichiométrique, les taux de nitrification potentielle estimés comptent pour 0.01 à 52% du flux total de NO3 nécessaire pour aboutir aux flux de N2 observés et diminuent avec l’augmentation de l’hypoxie.

Conservation of open solar magnetic flux and the floor in the heliospheric magnetic field

The near-Earth heliospheric magnetic field intensity, |B|, exhibits a strong solar cycle variation, but returns to the same ``floor'' value each solar minimum. The current minimum, however, has seen |B| drop below previous minima, bringing in to question the existence of a floor, or at the very least requiring a re-assessment of its value. In this study we assume heliospheric flux consists of a constant open flux component and a time-varying contribution from CMEs. In this scenario, the true floor is |B| with zero CME contribution. Using observed CME rates over the solar cycle, we estimate the ``no-CME'' |B| floor at ~4.0 +/- 0.3 nT, lower than previous floor estimates and below |B| observed this solar minimum. We speculate that the drop in |B| observed this minimum may be due to a persistently lower CME rate than the previous minimum, though there are large uncertainties in the supporting observational data.

Coronal mass ejections and magnetic flux buildup in the heliosphere

To test for magnetic flux buildup in the heliosphere from coronal mass ejections (CMEs), we simulate heliospheric flux as a constant background open flux with a time-varying interplanetary CME (ICME) contribution. As flux carried by ejecta can only contribute to the heliospheric flux budget while it remains closed, the ICME flux opening rate is an important factor. Two separate forms for the ICME flux opening rate are considered: (1) constant and (2) exponentially decaying with time. Coronagraph observations are used to determine the CME occurrence rates, while in situ observations are used to estimate the magnetic flux content of a typical ICME. Both static equilibrium and dynamic simulations, using the constant and exponential ICME flux opening models, require flux opening timescales of ∼50 days in order to match the observed doubling in the magnetic field intensity at 1 AU over the solar cycle. Such timescales are equivalent to a change in the ICME closed flux of only ∼7–12% between 1 and 5 AU, consistent with CSE signatures; no flux buildup results. The dynamic simulation yields a solar cycle flux variation with high variability that matches the overall variability of the observed magnetic field intensity remarkably well, including the double peak forming the Gnevyshev gap.

Growth of the lava dome and extrusion rates at Soufrière Hills Volcano, Montserrat, West Indies: 2005–2008

The third episode of lava dome growth at Soufrière Hills Volcano began 1 August 2005 and ended 20 April 2007. Volumes of the dome and talus produced were measured using a photo-based method with a calibrated camera for increased accuracy. The total dense rock equivalent (DRE) volume of extruded andesite magma (306 ± 51 Mm3) was similar within error to that produced in the earlier episodes but the average extrusion rate was 5.6 ± 0.9 m3s−1 (DRE), higher than the previous episodes. Extrusion rates varied in a pulsatory manner from <0.5 m3s−1 to ∼20 m3s−1. On 18 May 2006, the lava dome had reached a volume of 85 Mm3 DRE and it was removed in its entirety during a massive dome collapse on 20 May 2006. Extrusion began again almost immediately and built a dome of 170 Mm3 DRE with a summit height 1047 m above sea level by 4 April 2007. There were few moderate-sized dome collapses (1–10 Mm3) during this extrusive episode in contrast to the first episode of dome growth in 1995–8 when they were numerous. The first and third episodes of dome growth showed a similar pattern of low (<0.5 m3s−1) but increasing magma flux during the early stages, with steady high flux after extrusion of ∼25 Mm3

The evolution of the Sun's open magnetic flux: I. a single bipole

In this paper the origin and evolution of the Sun’s open magnetic flux are considered for single magnetic bipoles as they are transported across the Sun. The effects of magnetic flux transport on the radial field at the surface of the Sun are modeled numerically by developing earlier work by Wang, Sheeley, and Lean (2000). The paper considers how the initial tilt of the bipole axis (α) and its latitude of emergence affect the variation and magnitude of the surface and open magnetic flux. The amount of open magnetic flux is estimated by constructing potential coronal fields. It is found that the open flux may evolve independently from the surface field for certain ranges of the tilt angle. For a given tilt angle, the lower the latitude of emergence, the higher the magnitude of the surface and open flux at the end of the simulation. In addition, three types of behavior are found for the open flux depending on the initial tilt angle of the bipole axis. When the tilt is such that α ≥ 2◦ the open flux is independent of the surface flux and initially increases before decaying away. In contrast, for tilt angles in the range −16◦ < α < 2◦ the open flux follows the surface flux and continually decays. Finally, for α ≤ −16◦ the open flux first decays and then increases in magnitude towards a second maximum before decaying away. This behavior of the open flux can be explained in terms of two competing effects produced by differential rotation. Firstly, differential rotation may increase or decrease the open flux by rotating the centers of each polarity of the bipole at different rates when the axis has tilt. Secondly, it decreases the open flux by increasing the length of the polarity inversion line where flux cancellation occurs. The results suggest that, in order to reproduce a realistic model of the Sun’s open magnetic flux over a solar cycle, it is important to have accurate input data on the latitude of emergence of bipoles along with the variation of their tilt angles as the cycle progresses.

Reconfiguration and closure of lobe flux by reconnection during northward IMF: possible evidence for signatures in cusp/cleft auroral emissions

Observations are presented of the response of the dayside cusp/cleft aurora to changes in both the clock and elevation angles of the interplanetary magnetic field (IMF) vector, as monitored by the WIND spacecraft. The auroral observations are made in 630 nm light at the winter solstice near magnetic noon, using an all-sky camera and a meridian-scanning photometer on the island of Spitsbergen. The dominant change was the response to a northward turning of the IMF which caused a poleward retreat of the dayside aurora. A second, higher-latitude band of aurora was seen to form following the northward turning, which is interpreted as the effect of lobe reconnection which reconfigures open flux. We suggest that this was made possible in the winter hemisphere, despite the effect of the Earth's dipole tilt, by a relatively large negative X component of the IMF. A series of five events then formed in the poleward band and these propagated in a southwestward direction and faded at the equatorward edge of the lower-latitude band as it migrated poleward. It is shown that the auroral observations are consistent with overdraped lobe flux being generated by lobe reconnection in the winter hemisphere and subsequently being re-closed by lobe reconnection in the summer hemisphere. We propose that the balance between the reconnection rates at these two sites is modulated by the IMF elevation angle, such that when the IMF points more directly northward, the summer lobe reconnection site dominates, re-closing all overdraped lobe flux and eventually becoming disconnected from the Northern Hemisphere.

DIFFUSION OF MAGNETIC FIELD AND REMOVAL OF MAGNETIC FLUX FROM CLOUDS VIA TURBULENT RECONNECTION

The diffusion of astrophysical magnetic fields in conducting fluids in the presence of turbulence depends on whether magnetic fields can change their topology via reconnection in highly conducting media. Recent progress in understanding fast magnetic reconnection in the presence of turbulence reassures that the magnetic field behavior in computer simulations and turbulent astrophysical environments is similar, as far as magnetic reconnection is concerned. This makes it meaningful to perform MHD simulations of turbulent flows in order to understand the diffusion of magnetic field in astrophysical environments. Our studies of magnetic field diffusion in turbulent medium reveal interesting new phenomena. First of all, our three-dimensional MHD simulations initiated with anti-correlating magnetic field and gaseous density exhibit at later times a de-correlation of the magnetic field and density, which corresponds well to the observations of the interstellar media. While earlier studies stressed the role of either ambipolar diffusion or time-dependent turbulent fluctuations for de-correlating magnetic field and density, we get the effect of permanent de-correlation with one fluid code, i.e., without invoking ambipolar diffusion. In addition, in the presence of gravity and turbulence, our three-dimensional simulations show the decrease of the magnetic flux-to-mass ratio as the gaseous density at the center of the gravitational potential increases. We observe this effect both in the situations when we start with equilibrium distributions of gas and magnetic field and when we follow the evolution of collapsing dynamically unstable configurations. Thus, the process of turbulent magnetic field removal should be applicable both to quasi-static subcritical molecular clouds and cores and violently collapsing supercritical entities. The increase of the gravitational potential as well as the magnetization of the gas increases the segregation of the mass and magnetic flux in the saturated final state of the simulations, supporting the notion that the reconnection-enabled diffusivity relaxes the magnetic field + gas system in the gravitational field to its minimal energy state. This effect is expected to play an important role in star formation, from its initial stages of concentrating interstellar gas to the final stages of the accretion to the forming protostar. In addition, we benchmark our codes by studying the heat transfer in magnetized compressible fluids and confirm the high rates of turbulent advection of heat obtained in an earlier study.

Patterns of water and heat flux across a biome gradient from tropical forest to savanna in Brazil

We investigated the seasonal patterns of water vapor and sensible heat flux along a tropical biome gradient from forest to savanna. We analyzed data from a network of flux towers in Brazil that were operated within the Large-Scale Biosphere-Atmosphere Experiment in Amazonia (LBA). These tower sites included tropical humid and semideciduous forest, transitional forest, floodplain (with physiognomies of cerrado), and cerrado sensu stricto. The mean annual sensible heat flux at all sites ranged from 20 to 38 Wm(-2), and was generally reduced in the wet season and increased in the late dry season, coincident with seasonal variations of net radiation and soil moisture. The sites were easily divisible into two functional groups based on the seasonality of evaporation: tropical forest and savanna. At sites with an annual precipitation above 1900 mm and a dry season length less than 4 months (Manaus, Santarem and Rondonia), evaporation rates increased in the dry season, coincident with increased radiation. Evaporation rates were as high as 4.0 mm d(-1) in these evergreen or semidecidous forests. In contrast, ecosystems with precipitation less than 1700 mm and a longer dry season (Mato Grosso, Tocantins and Sao Paulo) showed clear evidence of reduced evaporation in the dry season. Evaporation rates were as low as 2.5 mm d(-1) in the transitional forests and 1 mm d(-1) in the cerrado. The controls on evapotranspiration seasonality changed along the biome gradient, with evaporative demand (especially net radiation) playing a more important role in the wetter forests, and soil moisture playing a more important role in the drier savannah sites.

Hydrochemistry and weathering rates on Corumbatai River basin, São Paulo State, Brazil

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Weathering rates and anthropogenic influences in a sedimentary basin, São Paulo State, Brazil

The weathering rate of rocks and chemical dynamics of the Corumbatai River basin, São Paulo State, Brazil, were evaluated using major elements as natural tracers. This basin has serious environmental problems in terms of quality of surface and rainwater, which affect the determination of weathering rate. The Corumbatai River, downstream from Rio Claro City, receives several elements/compounds through anthropogenic activities, with only K, SO42- and alkalinity yielding positive flux values. The negative flux of some anions/cations can be attributed to atmospheric loading mainly related to anthropogenic inputs, providing K a value of 16.7 ton/km(-2)a(-1) for the material removed by weathering in the Corumbatai River basin. This is equivalent to 26 x 10(6) kg of rock being removed each year by the Corumbatai River. The instantaneous flux was found to be a function of discharge, with the majority of dry residue (dissolved load) being transported during the summer (wet) months. The removed material in Corumbatai River basin derives mainly from two sub-basins (Cabegas River and Passa Cinco River), where the sandstones weather more easily than siltstones and claystones in the basin. (C) 2003 Elsevier Ltd. All rights reserved.

Spin accumulation encoded in electronic noise for mesoscopic billiards with finite tunneling rates

We study the effects of spin accumulation (inside reservoirs) on electronic transport with tunneling and reflections at the gates of a quantum dot. Within the stub model, the calculations focus on the current-current correlation function for the flux of electrons injected into the quantum dot. The linear response theory used allows us to obtain the noise power in the regime of thermal crossover as a function of parameters that reveal the spin polarization at the reservoirs. The calculation is performed employing diagrammatic integration within the universal groups (ensembles of Dyson) for a nonideal, nonequilibrium chaotic quantum dot. We show that changes in the spin distribution determine significant alterations in noise behavior at values of the tunneling rates close to zero, in the regime of strong reflection at the gates.