Paleofloods and ancient fishing weirs in NW Iberian rivers

A 15-m-thick, fluvial sedimentary record of the NW Iberian lower Miño River was studied. Grain-size analyses were performed and twelve samples were dated using optically stimulated luminescence dating techniques, documenting a 1300-yr-old reconstructed fluvial record that does not match with known climate fluctuations in the area, but is linked instead to the construction of a series of ancient fishing weirs (pesqueiras). The sedimentation phases are in agreement with known episodes of increased population density, which suggests active use of the pesqueiras. A number of sedimentation hiatuses in the fluvial record point towards damage to the pesqueiras during large-scale flooding in the Miño River basin, and a sudden drop in population probably due to the arrival of the plague in the 13th century AD. The oldest sedimentation phases started just after 700 AD, and we infer that the first pesqueiras were constructed around this time. This timing coincides with the transition of the NW Iberian landscape towards a more intensively used agricultural landscape, as evidenced from other geo-archeological investigations. The results demonstrate that the pesqueiras are several hundreds of years older than known from historical records, but not so old as to date back to the Roman occupation.     

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.     

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.     

A detailed study of the 5-Hz quasi-periodic oscillations in the bright X-ray transient and black hole candidate GRS 1739–278

We present a detailed study of the 5-Hz quasi-periodic oscillation (QPO) recently discovered in the bright X-ray transient and black hole candidate (BHC) GRS     (Borozdin & Trudolyubov) during a Rossi X-ray Timing Explorer observation taken on 1996 March 31. In total 6.6 ksec of on-source data were obtained, divided in two data sets of 3.4 and 3.2 ksec which were separated by ∼2.6 ksec. The 5-Hz QPO was only present during the second data set. The QPO increased in strength from below 2 per cent rms amplitude for photon energies below 4 keV to ∼5 per cent rms amplitude for energies above 10 keV. The soft QPO photons (below 5 keV) lagged the hard ones (above 10 keV) by almost 1.5 rad. Besides the QPO fundamental, its first overtone was detected. The strength of the overtone increased with photon energy (from < 2 per cent rms below 5 keV to ∼8 per cent rms above 10 keV). Although limited statistics did not allow for an accurate determination of the lags of the first overtone, indications are that also for this QPO the soft photons lagged the hard ones. When the 5-Hz QPO was not detected (i.e., during the first part of the observation), a broad noise component was found for photon energies below 10 keV but it became almost a true QPO (with a Q value of ∼1.9) above that energy, with a frequency of ∼3 Hz. Its hard photons preceded the soft ones in a way reminiscent of the 5-Hz QPO, strongly suggesting that both features are physically related. We discuss our finding in the framework of low-frequency QPOs and their properties in BHCs.     

Bridging the Transition from Mesoscale to Microscale Turbulence in Numerical Weather Prediction Models

With a focus towards developing multiscale capabilities in numerical weather prediction models, the specific problem of the transition from the mesoscale to the microscale is investigated. For that purpose, idealized one-way nested mesoscale to large-eddy simulation (LES) experiments were carried out using the Weather Research and Forecasting model framework. It is demonstrated that switching from one-dimensional turbulent diffusion in the mesoscale model to three-dimensional LES mixing does not necessarily result in an instantaneous development of turbulence in the LES domain. On the contrary, very large fetches are needed for the natural transition to turbulence to occur. The computational burden imposed by these long fetches necessitates the development of methods to accelerate the generation of turbulence on a nested LES domain forced by a smooth mesoscale inflow. To that end, four new methods based upon finite amplitude perturbations of the potential temperature field along the LES inflow boundaries are developed, and investigated under convective conditions. Each method accelerated the development of turbulence within the LES domain, with two of the methods resulting in a rapid generation of production and inertial range energy content associated to microscales that is consistent with non-nested simulations using periodic boundary conditions. The cell perturbation approach, the simplest and most efficient of the best performing methods, was investigated further under neutral and stable conditions. Successful results were obtained in all the regimes, where satisfactory agreement of mean velocity, variances and turbulent fluxes, as well as velocity and temperature spectra, was achieved with reference non-nested simulations. In contrast, the non-perturbed LES solution exhibited important energy deficits associated to a delayed establishment of fully-developed turbulence. The cell perturbation method has negligible computational cost, significantly accelerates the generation of realistic turbulence, and requires minimal parameter tuning, with the necessary information relatable to mean inflow conditions provided by the mesoscale solution.     

Nesting Turbulence in an Offshore Convective Boundary Layer Using Large-Eddy Simulations

The applicability of the one-way nesting technique for numerical simulations of the heterogeneous atmospheric boundary layer using the large-eddy simulation (LES) framework of the Weather Research and Forecasting model is investigated. The focus of this study is on LES of offshore convective boundary layers. Simulations were carried out using two subgrid-scale models (linear and non-linear) with two different closures [diagnostic and prognostic subgrid-scale turbulent kinetic energy (TKE) equations]. We found that the non-linear backscatter and anisotropy model with a prognostic subgrid-scale TKE equation is capable of providing similar results when performing one-way nested LES to a stand-alone domain having the same grid resolution but using periodic lateral boundary conditions. A good agreement is obtained in terms of velocity shear and turbulent fluxes, while velocity variances are overestimated. A streamwise fetch of 14 km is needed following each domain transition in order for the solution to reach quasi-stationary results and for the velocity spectra to generate proper energy content at high wavelengths, however, a pile-up of energy is observed at the low-wavelength portion of the spectrum on the first nested domain. The inclusion of a second nest with higher resolution allows the solution to reach effective grid spacing well within the Kolmogorov inertial subrange of turbulence and develop an appropriate energy cascade that eliminates most of the pile-up of energy at low wavelengths. Consequently, the overestimation of velocity variances is substantially reduced and a considerably better agreement with respect to the stand-alone domain results is achieved.     

Quantifying the influence of sill intrusion on the thermal evolution of organic‐rich sedimentary rocks in nonvolcanic passive margins: an example from ODP 210‐1276, offshore Newfoundland,Canada

Intrusive magmatism is an integral and understudied component in both volcanic and nonvolcanic passive margins. Here, we investigate the thermal effects of widespread (ca. 20 000 km²) intrusive magmatism on the thermal evolution of organic‐rich sedimentary rocks on the nonvolcanic Newfoundland passive margin. ODP 210‐1276 (45.41°N, 44.79°W) intersects two sills: an older, upper sill and a younger, lower sill that are believed to correspond to the high amplitude ‘U‐reflector’ observed across the Newfoundland Basin. A compilation of previous work collectively provides; (1) emplacement depth constraints, (2) vitrinite reflectance data and (3) ⁴⁰Ar/³⁹Ar dates. Collectively, these data sets provide a unique opportunity to model the conductive cooling of the sills and how they affect thermal maturity of the sedimentary sequence. A finite differences method was used to model the cooling of the sills, with the model outputs then being entered into the EASY%R_o vitrinite reflectance model. The modelled maturation profile for ODP 210‐1276 shows a significant but localized effect on sediment maturity as a result of the intrusions. Our results suggest that even on nonvolcanic margins, intrusive magmatism can significantly influence the thermal evolution in the vicinity of igneous intrusions. In addition, the presence of widespread sills on nonvolcanic passive margins such as offshore Newfoundland may be indicative of regional‐scale thermal perturbations that should be considered in source rock maturation studies.