Mercury biogeochemistry: Paradigm shifts,outstanding issues and research needs

Half a century of mercury research has provided scientists and policy makers with a detailed understanding of mercury toxicology, biogeochemical cycling and past and future impacts on human exposure. The complexity of the global biogeochemical mercury cycle has led to repeated and ongoing paradigm shifts in numerous mercury-related disciplines and outstanding questions remain. In this review, we highlight some of the paradigm shifts and questions on mercury toxicity, the risks and benefits of seafood consumption, the source of mercury in seafood, and the Arctic mercury cycle. We see a continued need for research on mercury toxicology and epidemiology, for marine mercury dynamics and ecology, and for a closer collaboration between observational mercury science and mercury modeling in general. As anthropogenic mercury emissions are closely tied to the energy cycle (in particular coal combustion), mercury exposure to humans and wildlife are likely to persist unless drastic emission reductions are put in place.     

VSX J074727.6 + 065050: A dwarf nova with a non-radially pulsating white dwarf primary

High speed photometric observations of the dwarf nova VSX J074727.6 + 065050 made in December 2009 during quiescence show an orbital modulation at P _orb =85.6 min. They also show that the star is a member of the relatively rare CV/ZZ group, i.e. the accreting white dwarf primary has non-radial pulsations. The two regions of oscillating power are at 684 s and 238 s. There is some evidence for hidden ∼1 μHz fine structure splitting, which has been seen in three other CV/ZZ stars.     

The effects of sulfur intercalation on the optical properties of artificial ‘hackmanite’, Na8[Al6Si6O24]Cl1.8S0.1; ‘sulfosodalite’, Na8[Al6Si6O24]S; and natural tugtupite, Na8[Be2Al2Si8O24](Cl,S)2?δ

The minerals ??hackmanite?? and tugtupite exhibit tenebrescence (reversible photochromism) and photoluminescence. These features are generally attributed to the presence of sulfide species within their structures. But how these optical properties might be affected by intercalating additional amounts of sulfur into their structures was until now unknown. Artificial ??hackmanite??, Na₈[Al₆Si₆O₂₄]Cl_1.8S_0.1, and ??sulfosodalite??, Na₈[Al₆Si₆O₂₄]S, were heated with sulfur in evacuated quartz-glass ampoules over the temperature range 450?C1,050°C. This work has shown that sulfur intercalation into Na₈[Al₆Si₆O₂₄]Cl_1.8S_0.1 destroys the tenebrescence and induces a permanently pale blue and, at higher temperature, a pale green coloration. The effect on Na₈[Al₆Si₆O₂₄]S induced similar colorations but of a deeper hue. Annealing tugtupite, Na₈[Be₂Al₂Si₈O₂₄](Cl,S)_2??? under a sulfur atmosphere over the range 600?C700°C, destroyed the tenebrescence and resulted in a colorless tugtupite; but did not effect the photoluminescence. This suggests that the chemical species responsible for the tenebrescence in tugtupite is unlikely to be the same as that for the luminescence.     

Probabilistic tsunami hazard assessment for Messina Strait Area (Sicily, Italy)

The general modular Bayesian procedure is applied to provide a probabilistic tsunami hazard assessment (PTHA) for the Messina Strait Area (MSA), Italy. This is the first study in an Italian area where the potential tsunamigenic events caused by both submarine seismic sources (SSSs) and submarine mass failures (SMFs) are examined in a probabilistic assessment. The SSSs are localized on active faults in MSA as indicated by the instrumental data of the catalogue of the Italian seismicity; the SMFs are spatially identified using their propensity to failure in the Ionian and Tyrrhenian Seas on the basis of mean slope and mean depth, and using marine geology background knowledge. In both cases the associated probability of occurrence is provided. The run-ups were calculated at key sites that are main cities and/or important sites along the Eastern Sicily and the Southern Calabria coasts where tsunami events were recorded in the past. The posterior probability distribution combines the prior probability and the likelihood calculated in the MSA. The prior probability is based on the physical model of the tsunami process, and the likelihood is based on the historical data collected by the historical catalogues, background knowledge, and marine geological information. The posterior SSSs and SMFs tsunami probabilities are comparable and are combined to produce a final probability for a full PTHA in MSA.     

Modelling sediment dynamics due to hillslope–river interactions: incorporating fluvial behaviour in landscape evolution model LAPSUS

Landscape evolution models (LEMs) simulate the three‐dimensional development of landscapes over time. Different LEMs have different foci, e.g. erosional behaviour, river dynamics, the fluvial domain, hillslopes or a combination. LEM LAPSUS is a relatively simple cellular model operating on timescales of centuries to millennia and using annual timesteps that has had a hillslope focus. Our objective was to incorporate fluvial behaviour in LAPSUS without changing the existing model equations. The model should be able to reproduce alternating aggradation and incision in the floodplains of catchments, depending on simulated conditions. Testing was done using an artificial digital elevation model (DEM) and a demonstration of the ability for fluvial simulation was performed for a real landscape (Torrealvilla catchment, southeast Spain). Model equations to calculate sediment dynamics and water routing were similar for both hillslope and fluvial conditions, but different parameter values were used for these domains, defined based on annual discharge. Parameters changing between the domains are convergence factor p, which is used in the multiple flow algorithm to route water, and discharge and gradient exponents m and n, used in transport capacity calculations. Erodibility and ‘sedimentability’ factors K and P were changed between cold (little vegetation, high erodibility) and warm conditions (more vegetation, lower erodibility). Results show that the adapted parameters reproduced alternating aggradation – due to divergent flow in the floodplain and sediment supply under cold conditions – and incision due to reduced sediment supply and resulting clean water erosion during simulated warm conditions. The simulated results are due to interactions between hillslopes and floodplains, as the former provide the sediments that are deposited in the latter. Similar behaviour was demonstrated when using the real DEM. Sensitivity and resolution analysis showed that the model is sensitive to changes in m, n and p and that model behaviour is influenced by DEM resolution. Copyright © 2012 John Wiley & Sons, Ltd.     

Estimates of the transition zone temperature in a mechanically mixed upper mantle

Partial melting of mantle peridotite generates a physically and chemically layered oceanic lithosphere that is cycled back into the mantle in subduction zones. Stirring times of the mantle are too long to allow for complete re-homogenization of subducted basalt and harzburgite, given the low chemical diffusivity of the solid mantle. This suggests that the Earth's mantle is a mechanical mixture of basaltic and harzburgitic components. Using a recently developed thermodynamic formulism we determine the phase equilibria and the seismic properties of a mantle comprised of a mechanical mixture of basalt and harzburgite (MM) and a homogeneous mantle (EA) with identical pyrolitic bulk chemistry. We use the theoretical shear velocity profiles as a new thermometer of the mantle below the magma-genetic zone by modeling the difference ΔT_410-660 between traveltimes of shear wave reflections off the 410-km and 660-km with the potential temperature T_P. ΔT_410-660 are measured from waveform stacks. They indicate that, over 1000+ km wave lengths, the temperature varies by about 200 K. Lowest and highest temperatures are resolved for the western Pacific subduction zones and the central Pacific, respectively. This variation is similar for the EA and MM and is in excellent agreement with estimates of transition zone thickness and shear velocity variations. The median value of T_P for the EA is 1720 K. It is about 1625 K for the MM, a value that is in better agreement with the Normal-MORB values of 1610 ± 40 K inferred from olivine-liquid equilibria given that our sampling region encompasses the Western Pacific subduction zones and the oldest parts of the Pacific Plate. We argue therefore that a mechanical mixed mantle, with generally higher velocities and steeper velocities gradients, represents a better physical reference model than a model based on a fully equilibrated assemblage.     

No evidence for thermogenic methane release in coal from the Karoo-Ferrar large igneous province

The Early Jurassic Toarcian oceanic anoxic event (T-OAE) and concurrent negative carbon-isotope (δ¹³C) excursion have recently been attributed to either the release of methane (CH₄) clathrates or thermogenic CH₄ gas associated with the Karoo-Ferrar large igneous province (LIP) into coals and organic-rich shales. ¹²C-enriched thermogenic CH₄ production associated with the Karoo-Ferrar would result in residual material being ¹²C-depleted nearer the intrusions. In this study, geochemical analyses (carbon isotopes, volatile matter (VM), vitrinite reflectance (R_o)) are reported for two coal transects associated with dykes intruding the No. 4L coal in the Highveld Coalfield, Karoo Basin, South Africa. VM decreases from over 35% to around 15% in one transect, and the second transect shows a less pronounced decrease (from > 25% to ~ 16%). Accompanying the decrease in VM content is an increase in R_o from background levels of around 0.7% to over 4% adjacent to the dyke; used as a palaeo-geothermometer, R_o values indicate background temperatures of ~ 100 °C increasing to > 300 °C close to the contact. Despite changes in VM and R_o, there are no significant changes in δ¹³C, certainly not of the magnitude that would be expected associated with large-scale thermogenic CH₄ generation. These and other Gondwanan coals have low vitrinite and liptinite contents (components more prone to CH₄ generation), in part explaining the modest decreases in VM adjacent to the dykes. This, combined with the relatively narrow metamorphic aureole surrounding the intrusions and the likelihood that at least some of the volatiles generated by the intrusion were trapped as coalbed CH₄ or condensed as pyrolytic carbon, suggests only limited CH₄ release. In addition, based on original estimates of moisture contents in these coals and the depth at time of intrusion (1,000–2,000 m) the dykes would have lost most of their energy heating and evaporating water, thus having very little remaining energy to generate thermogenic CH₄.     

Focal depths and mechanisms of Tohoku-Oki aftershocks from teleseismic P wave modeling

Aftershocks of the 2011 Tohoku-Oki great earthquake have a wide range of focal depths and fault plane mechanisms. We constrain the focal depths and focal mechanisms of 69 aftershocks with M _w > 5.4 by modeling the waveforms of teleseismic P and its trailing near-surface reflections pP and sP. We find that the “thrust events” are within 10 km from the plate interface. The dip angles of these thrust events increase with depth from ~5° to ~25°. The “non-thrust events” vary from 60 km above to 40 km below the plate interface. Normal and strike-slip events within the overriding plate point to redistribution of stress following the primary great earthquake; however, due to the spatially variable stress change in the Tohoku-Oki earthquake, an understanding of how the mainshock affected the stresses that led to the aftershocks requires accurate knowledge of the aftershock location.     

Storm-driven sediment transport in Massachusetts Bay

Massachusetts Bay is a semi-enclosed embayment in the western Gulf of Maine about 50 km wide and 100 km long. Bottom sediment resuspension is controlled predominately by storm-induced surface waves and transport by the tidal- and wind-driven circulation. Because the Bay is open to the northeast, winds from the northeast (‘Northeasters’) generate the largest surface waves and are thus the most effective in resuspending sediments. The three-dimensional oceanographic circulation model Regional Ocean Modeling System (ROMS) is used to explore the resuspension, transport, and deposition of sediment caused by Northeasters. The model transports multiple sediment classes and tracks the evolution of a multilevel sediment bed. The surficial sediment characteristics of the bed are coupled to one of several bottom-boundary layer modules that calculate enhanced bottom roughness due to wave–current interaction. The wave field is calculated from the model Simulating WAves Nearshore (SWAN). Two idealized simulations were carried out to explore the effects of Northeasters on the transport and fate of sediments. In one simulation, an initially spatially uniform bed of mixed sediments exposed to a series of Northeasters evolved to a pattern similar to the existing surficial sediment distribution. A second set of simulations explored sediment-transport pathways caused by storms with winds from the northeast quadrant by simulating release of sediment at selected locations. Storms with winds from the north cause transport southward along the western shore of Massachusetts Bay, while storms with winds from the east and southeast drive northerly nearshore flow. The simulations show that Northeasters can effectively transport sediments from Boston Harbor and the area offshore of the harbor to the southeast into Cape Cod Bay and offshore into Stellwagen Basin. This transport pattern is consistent with Boston Harbor as the source of silver found in the surficial sediments of Cape Cod Bay and Stellwagen Basin.