Hässeldala – a key site for Last Termination climate events in northern Europe

The Last Termination (19 000–11 000 a BP) with its rapid and distinct climate shifts provides a perfect laboratory to study the nature and regional impact of climate variability. The sedimentary succession from the ancient lake at Hässeldala Port in southern Sweden with its distinct Lateglacial/early Holocene stratigraphy (>14.1–9.5 cal. ka BP) is one of the few chronologically well‐constrained, multi‐proxy sites in Europe that capture a variety of local and regional climatic and environmental signals. Here we present Hässeldala's multi‐proxy records (lithology, geochemistry, pollen, diatoms, chironomids, biomarkers, hydrogen isotopes) in a refined age model and place the observed changes in lake status, catchment vegetation, summer temperatures and hydroclimate in a wider regional context. Reconstructed mean July temperatures increased between c. 14.1 and c. 13.1 cal. ka BP and subsequently declined. This latter cooling coincided with drier hydroclimatic conditions that were probably associated with a freshening of the Nordic Seas and started a few hundred years before the onset of Greenland Stadial 1 (c. 12.9 cal. ka BP). Our proxies suggest a further shift towards colder and drier conditions as late as c. 12.7 cal. ka BP, which was followed by the establishment of a stadial climate regime (c. 12.5–11.8 cal. ka BP). The onset of warmer and wetter conditions preceded the Holocene warming over Greenland by c. 200 years. Hässeldala's proxies thus highlight the complexity of environmental and hydrological responses across abrupt climate transitions in northern Europe.     

Why the <Emphasis Type="Italic">Euler</Emphasis> scheme in particle tracking is not enough: the shallow-sea pycnocline test case

During the last decades, the Euler scheme was the common “workhorse” in particle tracking, although it is the lowest-order approximation of the underlying stochastic differential equation. To convince the modelling community of the need for better methods, we have constructed a new test case that will show the shortcomings of the Euler scheme. We use an idealised shallow-water diffusivity profile that mimics the presence of a sharp pycnocline and thus a quasi-impermeable barrier to vertical diffusion. In this context, we study the transport of passive particles with or without negative buoyancy. A semi-analytic solutions is used to assess the performance of various numerical particle-tracking schemes (first- and second-order accuracy), to treat the variations in the diffusivity profile properly. We show that the commonly used Euler scheme exhibits a poor performance and that widely used particle-tracking codes shall be updated to either the Milstein scheme or second-order schemes. It is further seen that the order of convergence is not the only relevant factor, the absolute value of the error also is.     

Mars Express and Venus Express multi-point observations of geoeffective solar flare events in December 2006

In December 2006, a single active region produced a series of proton solar flares, with X-ray class up to the X9.0 level, starting on 5 December 2006 at 10:35 UT. A feature of this X9.0 flare is that associated MeV particles were observed at Venus and Mars by Venus Express (VEX) and Mars Express (MEX), which were ∼80° and ∼125° east of the flare site, respectively, in addition to the Earth, which was ∼79° west of the flare site. On December 5, 2006, the plasma instruments ASPERA-3 and ASPERA-4 on board MEX and VEX detected a large enhancement in their respective background count levels. This is a typical signature of solar energetic particle (SEP) events, i.e., intensive MeV particle fluxes. The timings of these enhancements were consistent with the estimated field-aligned travel time of particles associated with the X9.0 flare that followed the Parker spiral to reach Venus and Mars. Coronal mass ejection (CME) signatures that might be related to the proton flare were twice identified at Venus within <43 and <67 h after the flare. Although these CMEs did not necessarily originate from the X9.0 flare on December 5, 2006, they most likely originated from the same active region because these characteristics are very similar to flare-associated CMEs observed on the Earth. These observations indicate that CME and flare activities on the invisible side of the Sun may affect terrestrial space weather as a result of traveling more than 90° in both azimuthal directions in the heliosphere. We would also like to emphasize that during the SEP activity, MEX data indicate an approximately one-order of magnitude enhancement in the heavy ion outflow flux from the Martian atmosphere. This is the first observation of the increase of escaping ion flux from Martian atmosphere during an intensive SEP event. This suggests that the solar EUV flux levels significantly affect the atmospheric loss from unmagnetized planets.     

Superfluid signatures in magnetar seismology

We investigate the role of neutron star superfluidity for magnetar oscillations. Using a plane-wave analysis, we estimate the effects of a neutron superfluid in the elastic crust region. We demonstrate that the superfluid imprint is likely to be more significant than the effects of the crustal magnetic field. We also consider the region immediately beneath the crust, where superfluid neutrons are thought to coexist with a type II proton superconductor. Since the magnetic field in the latter is carried by an array of fluxtubes, the dynamics of this region differ from standard magnetohydrodynamics. We show that the presence of the neutron superfluid (again) leaves a clear imprint on the oscillations of the system. Taken together, our estimates show that the superfluid components cannot be ignored in efforts to carry out 'magnetar seismology'. This increases the level of complexity of the modelling problem, but also points to the possibility of using observations to probe the superfluid nature of supranuclear matter.     

Influence of the boundary layer height on the global air–sea surface fluxes

Results from large-eddy simulations and field measurements have previously shown that the velocity field is influenced by the boundary layer height, z _i , during close to neutral, slightly unstable, atmospheric stratification. During such conditions the non-dimensional wind profile, φ _m , has been found to be a function of both z/L and z _i /L. At constant z/L, φ _m decreases with decreasing boundary layer height. Since φ _m is directly related to the parameterizations of the air–sea surface fluxes, these results will have an influence when calculating the surface fluxes in weather and climate models. The global impact of this was estimated using re-analysis data from 1979 to 2001 and bulk parameterizations. The results show that the sum of the global latent and sensible mean heat fluxes increase by 0.77 W m⁻² or about 1% and the mean surface stress increase by 1.4 mN m⁻² or 1.8% when including the effects of the boundary layer height in the parameterizations. However, some regions show a larger response. The greatest impact is found over the tropical oceans between 30°S and 30°N. In this region the boundary layer height influences the non-dimensional wind profile during extended periods of time. In the mid Indian Ocean this results in an increase of the mean annual heat fluxes by 2.0 W m⁻² and an increase of the mean annual surface stress by 2.6 mN m⁻².     

Potential climate change impacts on the probability of wind damage in a south Swedish forest

We estimated how the possible changes in wind climate and state of the forest due to climate change may affect the probability of exceeding critical wind speeds expected to cause wind damage within a forest management unit located in Southern Sweden. The topography of the management unit was relatively gentle and the forests were dominated by Norway spruce (Picea abies (L.) Karst.). We incorporated a model relating the site index (SI) to the site productivity into the forest projection model FTM. Using estimated changes in the net primary production (NPP) due to climate change and assuming a relative change in NPP equal to a relative change in the site productivity, we simulated possible future states of the forest under gradual adjustment of SI in response to climate change. We estimated changes in NPP by combining the boreal-adapted BIOMASS model with four regional climate change scenarios calculated using the RCAO model for the period 2071–2100 and two control period scenarios for the period 1961–1990. The modified WINDA model was used to calculate the probability of wind damage for individual forest stands in simulated future states of the forest. The climate change scenarios used represent non-extreme projections on a 100-year time scale in terms of global mean warming. A 15–40% increase in NPP was estimated to result from climate change until the period 2071–2100. Increasing sensitivity of the forest to wind was indicated when the management rules of today were applied. A greater proportion of the calculated change in probability of wind damage was due to changes in wind climate than to changes in the sensitivity of the forest to wind. While regional climate scenarios based on the HadAM3H general circulation model (GCM) indicated no change (SRES A2 emission scenario) or a slightly reduced (SRES B2 emission scenario) probability of wind damage, scenarios based on the ECHAM4/OPYC3 GCM indicated increased probability of wind damage. The assessment should, however, be reviewed as the simulation of forest growth under climate change as well as climate change scenarios are refined.     

Spectra of CO<Subscript>2</Subscript> and Water Vapour in the Marine Atmospheric Surface Layer

Spectra of CO₂ and water vapour fluctuations from measurements made in the marine atmospheric surface layer have been analyzed. A normalization of spectra based on Monin–Obukhov similarity theory, originally developed for wind speed and temperature, has been successfully extended also to CO₂ and humidity spectra. The normalized CO₂ spectra were observed to have somewhat larger contributions from low frequencies compared to humidity spectra during unstable stratification. However, overall, the CO₂ and humidity spectra showed good agreement as did the cospectra of vertical velocity with water vapour and CO₂ respectively. During stable stratification the spectra and cospectra displayed a well-defined spectral gap separating the mesoscale and small-scale turbulent fluctuations. Two-dimensional turbulence was suggested as a possible source for the mesoscale fluctuations, which in combination with wave activity in the vertical wind is likely to explain the increase in the cospectral energy for the corresponding frequency range. Prior to the analysis the turbulence time series of the density measurements were converted to time series of mixing ratios relative to dry air. Some differences were observed when the spectra based on the original density measurements were compared to the spectra based on the mixing ratio time series. It is thus recommended to always convert the density time series to mixing ratio before performing spectral analysis.     

Storm surges in the Western Baltic Sea: the present and a possible future

Globally-coupled climate models are generally capable of reproducing the observed trends in the globally averaged atmospheric temperature or mean sea level. However, the global models do not perform as well on regional/local scales. Here, we present results from four 100-year ocean model experiments for the Western Baltic Sea. In order to simulate storm surges in this region, we have used the General Estuarine Transport Model (GETM) as a high-resolution local model (spatial resolution ≈ 1?km), nested into a regional atmospheric and regional oceanic model in a fully baroclinic downscaling approach. The downscaling is based on the global model ECHAM5/MPI-OM. The projections are imbedded into two greenhouse-gas emission scenarios, A1B and B1, for the period 2000–2100, each with two realisations. Two control runs from 1960 to 2000 are used for validation. We use this modelling system to statistically reproduce the present distribution of surge extremes. The usage of the high-resolution local model leads to an improvement in surge heights of at least 10% compared to the driving model. To quantify uncertainties associated with climate projections, we investigate the impact of enhanced wind velocities and changes in mean sea levels. The analysis revealed a linear dependence of surge height and mean sea level, although the slope parameter is spatially varying. Furthermore, the modelling system is used to project possible changes within the next century. The results show that the sea level rise has greater potential to increase surge levels than does increased wind speed. The simulations further indicate that the changes in storm surge height in the scenarios can be consistently explained by the increase in mean sea level and variation in wind speed.     

The distribution of neodymium isotopes in Arctic Ocean basins

Nd concentration and isotope data have been obtained for the Canada, Amundsen, and Makarov Basins of the Arctic Ocean. A pattern of high Nd concentrations (up to 58 pM) at shallow depths is seen throughout the Arctic, and is distinct from that generally seen in other oceans where surface waters are relatively depleted. A range of isotopic variations across the Arctic and within individual depth profiles reflects the different sources of waters. The dominant source of water, and so Nd, is the Atlantic Ocean, with lesser contributions from the Pacific and Arctic Rivers. Radiogenic isotope Nd signatures (up to ε_Nd = −6.5) can be traced in Pacific water flowing into the Canada Basin. Waters from rivers draining older terrains provide very unradiogenic Nd (down to ε_Nd = −14.2) that can be traced in surface waters across much of the Eurasian Basin. A distinct feature of the Arctic is the general influence of the shelves on the Nd concentrations of waters flowing into the basins, either from the Pacific across the Chukchi Sea, or from across the extensive Siberian shelves. Water-shelf interaction results in an increase in Nd concentration without significant changes in salinity in essentially all waters in the Arctic, through processes that are not yet well understood. In estuarine regions other processes modify the Nd signal of freshwater components supplied into the Arctic Basin, and possibly also contribute to sedimentary Nd that may be subsequently involved in sediment-water interactions. Mixing relationships indicate that in estuaries, Nd is removed from major river waters to different degrees. Deep waters in the Arctic are higher in Nd than the inflowing Atlantic waters, apparently through enrichments of waters on the shelves that are involved in ventilating the deep basins. These enrichments generally have not resulted in major shifts in the isotopic compositions of the deep waters in the Makarov Basin (ε_Nd ∼ −10.5), but have created distinctive Nd isotope signatures that were found near the margin of the Canada Basin (with ε_Nd ∼ −9.0). The deep waters of the Amundsen Basin are also distinct from the Atlantic waters (with ε_Nd = −12.3), indicating that there has been limited inflow from the adjacent Makarov Basin through the Lomonosov Ridge.     

Precise electron microprobe age determination in altered uraninite: consequences on the intrusion age and the metallogenic significance of the Kirchberg granite (Erzgebirge,Germany)

An integral method of data evaluation for chemical dating of accessory uraninite by electron microprobe analysis and electron microscopy is proposed. The internal structure of uraninite grains, the contents of Si, Ca, and Fe, and tests for loss or gain of Pb and U are suitable criteria to decide whether or not a given crystal area was affected by secondary alteration. The method was applied to uraninite from the Kirchberg granite, Erzgebirge, Germany, and an age of 330LJ Ma (DŽ†) was constrained. This result supports earlier assumptions that there were at least two events of Hercynian granite magmatism separated by a time gap of about 40 Ma, and that the magmatic evolution of the Kirchberg granite is not related to the formation of the W mineralisation developed in the endocontact of this granite. Secondary alteration of uraninite is probably closely related to the formation of U deposits in the vicinity of the Kirchberg massif which was extensively mined in the past.