Sub‐arctic Holocene climatic and oceanographic variability in Stjernsund,northern Norway: evidence from benthic foraminifera and stable isotopes

A high‐resolution record, covering 9.3–0.2 ka BP, from the sub‐arctic Stjernsund (70°N) was studied for benthic foraminiferal faunas and stable isotopes, revealing three informally named main phases during the Holocene. The Early‐ to Mid‐Holocene (9.3–5.0 ka BP) was characterized by the strong influence of the North Atlantic Current (NAC), which prevented the reflection of the Holocene Climatic Optimum (HCO) in the bottom‐water temperature. During the Mid‐Holocene Transition (5.0–2.5 ka BP), a turnover of benthic foraminiferal faunas occurred, Atlantic Water species decreased while Arctic‐Polar species increased, and the oxygen isotope record showed larger fluctuations. Those variations correspond to a period of global climate change, to spatially more heterogeneous benthic foraminiferal faunas in the Nordic Seas region, and to regionally diverging terrestrial temperatures. The Cool Late Holocene (2.5–0.2 ka BP) was characterized by increased abundances of Arctic‐Polar species and a steady cooling trend reflected in the oxygen isotopes. In this period, our record differs considerably from those on the SW Barents Sea shelf and locations farther south. Therefore, we argue that regional atmospheric cooling triggered the late Holocene cooling trend. Several cold episodes centred at ～8.3, ～7.8, ～6.5, ～4.9, ～3.9 and ～3.3 ka BP were identified from the benthic foraminiferal faunas and the δ¹⁸O record, which correlated with marine and atmospherically driven proxy records. This suggests that short‐term cold events may result from reduced heat advection via the NAC or from colder air temperatures.     

Comparison of numerical methods for simulating strongly nonlinear and heterogeneous reactive transport problems—the MoMaS benchmark case

Although multicomponent reactive transport modeling is gaining wider application in various geoscience fields, it continues to present significant mathematical and computational challenges. There is a need to solve and compare the solutions to complex benchmark problems, using a variety of codes, because such intercomparisons can reveal promising numerical solution approaches and increase confidence in the application of reactive transport codes. In this contribution, the results and performance of five current reactive transport codes are compared for the 1D and 2D subproblems of the so-called easy test case of the MoMaS benchmark (Carrayrou et al., Comput Geosci, 2009, this issue). This benchmark presents a simple fictitious reactive transport problem that highlights the main numerical difficulties encountered in real reactive transport problems. As a group, the codes include iterative and noniterative operator splitting and global implicit solution approaches. The 1D easy advective and 1D easy diffusive scenarios were solved using all codes, and, in general, there was a good agreement, with solution discrepancies limited to regions with rapid concentration changes. Computational demands were typically consistent with what was expected for the various solution approaches. The differences between solutions given by the three codes solving the 2D problem are more important. The very high computing effort required by the 2D problem illustrates the importance of parallel computations. The most important outcome of the benchmark exercise is that all codes are able to generate comparable results for problems of significant complexity and computational difficulty.     

Atmospheric profiling with the UAS SUMO: a new perspective for the evaluation of fine-scale atmospheric models

For the first time, unmanned aerial system measurements collected by the small unmanned meteorological observer (SUMO) are used to evaluate atmospheric boundary layer (ABL) parameterization schemes embedded in the Advanced Weather Research and Forecasting model (AR-WRF). Observation sites were located in the vicinity of the almost idealized shaped mountain Hofsj?kull, Central Iceland. SUMO profiles provided temperature, relative humidity and wind data to maximum heights of 3?km above ground. Two cases are investigated, one with calm wind conditions and development of a convective ABL and one with moderate winds and gravity waves over Hofsj?kull. For the high-resolution simulation with AR-WRF, three two-way nested domains are chosen with a grid size of 9, 3 and 1?km. During its first meteorological test, SUMO has proved its great value for the investigation of the diurnal evolution of the ABL and the identification of mesoscale features residing above the ABL, such as subsidence.     

The ion mass loading rate at Io

The Io plasma torus, composed of mostly heavy ions of oxygen and sulfur, is sustained by an Iogenic mass loading rate of ∼10³⁰ amu s⁻¹ = 1.6 × 10²⁸ SO₂ s⁻¹ or approximately 10³ kg s⁻¹(A.L. Broadfoot et al., 1979, Science 204, 979-982). We argue on the basis of available power sources, reanalysis of F. Bagenal (1997, Geophys. Res. Lett. 24, 2111-2114), HST UV remote sensing, and detailed model calculations that at most 20% of this mass leaves Io in the form of ions, i.e., ≤3 × 10²⁷ × (n_e,0/3600 cm⁻³) ions s⁻¹, where n_e,0 is the average torus electron density. For the Galileo spacecraft Io pass in December 1995, the ion mass loading rate was ≤3 × 10²⁷ ions s⁻¹, whereas for the Voyager epoch with lower n_e,0 (=2000 cm⁻³), this rate would be ≤1.7 × 10²⁷ ions s⁻¹, consistent with the D.E. Shemansky (1980, Astrophys. J. 242, 1266-1277) mass loading limit of ≤1 × 10²⁷ ions s⁻¹. We investigate the processes that control Io’s large scale electrodynamic interaction and find that the elastic collision rate exceeds the ionization/pickup rate by at least a factor of 5 for all atmospheric column densities considered (10¹⁶-10²¹ m⁻²) and by a factor of ∼100 for the most realistic column density. Consequently, elastic collisions are mostly responsible for Io’s high conductances and thus generate Io’s large scale electrodynamic interaction such as the generation of Io’s electric current system and the slowing of the plasma flow. The electrodynamic part of Io’s interaction is thus best described as an ionosphere-like interaction rather than a comet-like interaction. An analytic expression for total electron impact rates is derived for Io’s atmosphere, which is independent of any particular model for the 3D interaction of torus electrons with its atmosphere.     

Tube wave modeling by the finite-difference method with varying grid spacing

A finite-difference approach of aP-SV modeling scheme is applied to compute seismic wave propagation in heterogeneous isotropic media, including fluid-filled boreholes. The discrete formulation of the equation of motion requires the definition of the material parameters at the grid points of the numerical mesh. The grid spacing is chosen as coarse as possible with respect to the accurate representation of the shortest wavelength. If we assume frequencies lower than 250 Hz then the grid spacing is usually chosen in the range of a few meters. One encounters difficulties because of the large-scale difference between the grid spacing and the size of the borehole, usually several centimeters.These difficulties can be overcome by a grid refinement technique. This technique provides the construction of grids with varying grid spacing. The grid spacing in the vicinity of the borehole is chosen such that the borehole is properly represented. An example demonstrates the accuracy of this technique by comparisons with other methods. Unlike many analytical methods, the FD method can handle complex subsurface geometries. Further numerical examples of walk-awayVSP configurations show tube wave propagation within fluid-filled boreholes of realistic diameters.     

Vacuum energy in cosmic dynamics

The analysis of the Th/U ratio in meteorites and the evolutionary ages of globular clusters favour values of the cosmic age of (19±5)×10⁹ yr. This evidence together with a Hubble parameterH ₀>70 km s^–1 Mpc^–1=(14×10⁹ yr)^–1 cannot be reconciled in a Friedmann model with =0. It requires a cosmological constant in the order of 10^–56 cm^–2, equivalent to a vacuum density _v =10^–29 g cm^–3 The Friedmann-Lemaître models (>0) with a hot big-bang have been calculated. They are based on a present value of the baryonic matter density of ₀=0.5×10^–30 g cm^–3 as derived from the primordial⁴He and²H abundances.For a Hubble parameter ofH ₀=75 km s^–1 Mpc^–1, our analysis favours a set of models which can be represented by a model with Euclidean metric (density parameter ₀=1.0, deceleration parameterq ₀=–0.93, aget ₀=19.7×10⁹ yr) and by a closed model with perpetual expansion (₀=1.072,q ₀=–1.0, aget ₀=21.4×10⁹ yr). A present density parameter close to one can indeed be expected if the conjecture of an exponential inflation of the very early universe is correct.The possible behaviour of the vacuum density is demonstrated with the help of Streeruwitz' formula in the context of the closed model with an inflationary phase at very early times.     

Monsoon-induced upwelling off the Vietnamese coast

During the southwest monsoon from July 8 to 28, 2003, an interdisciplinary cruise took place in the central area of Vietnamese upwelling with “MV Nghien Cuu Bien” in the South China Sea. Physical observations in the upwelling area are analyzed with respect to local/regional wind forcing and far field forcing. Nutrients and phytoplankton measurements are discussed with respect to exchange processes between different water masses. The wind-induced coastal upwelling by local wind forcing is much weaker than in the previous years due to weaker-than-normal winds. This can be attributed to the far field forcing of the 2002/2003 El Niño event which modulates the upwelling intensity. The atmospheric conditions reflect the typical situation after an El Niño event which weakens the wind-induced coastal upwelling, reduces the latent heat flux, and results in higher-than-normal sea-surface temperatures. The general circulation pattern during SW monsoon is driven by the spatial asymmetry in the monsoon forcing. The flow pattern is characterized by an upwelling-induced northward undercurrent and a recently detected southward countercurrent. The resulting stretching deformation of this flow pattern forms an offshore jet between ～12°N and 12.5°N and causes a local enhancement of the upwelling intensity. The upwelling due to stretching deformation is a peculiarity, which makes the Vietnamese upwelling area different to other upwelling areas. A budget of the upwelling components is presented: the strongest contribution in 2003 to the Vietnamese upwelling is the dynamical upwelling due to the clockwise rotation of the northward undercurrent. The internal radius of deformation separates the upwelling area from the offshore area as well as different water masses. Mekong River and the Gulf of Thailand waters which are offshore show nutrient depletion. Therefore, high chlorophyll maxima cannot be explained by nutrient supply from river runoff. The dynamical upwelling brings in nutrient-rich Maximum Salinity Water into the euphotic zone. This causes a subsurface chlorophyll maximum between 20 and 40 m water depth along the northward undercurrent. Deflection from the Redfield ratio in the C:N ratio and negative excess nitrogen identifies the region as nitrogen-limited which may favor cyanobacteria blooms. The consequence is a unique feature in new production: in the upwelling area, new production is based on upwelled nitrate, whereas offshore in the nutrient-depleted Mekong and Gulf of Thailand water, new production is based in addition on nitrogen fixation.