Wetland development, permafrost history and nutrient cycling inferred from late Holocene peat and lake sediment records in subarctic Sweden

Permafrost in peatlands of subarctic Sweden is presently thawing at accelerated rates, which raises questions about the destiny of stored carbon and nutrients and impacts on adjacent freshwater ecosystems. In this study we use peat and lake sediment records from the Stordalen palsa mire in northern Sweden to address the late Holocene (5,000 cal BP-present) development of the mire as well as related changes in carbon and nutrient cycling. Formation, sediment accumulation and biogeochemistry of two studied lakes are suggested to be largely controlled by the development of the mire and its permafrost dynamics. Peat inception took place at ca. 4,700 cal BP as a result of terrestrialisation. Onset of organic sedimentation in the adjacent lakes occurred at ca. 3,400 and 2,650 cal BP in response to mire expansion and permafrost aggradation, respectively. Mire erosion, possibly due to permafrost decay, led to re-deposition of peat into one of the lakes after ca. 2,100 cal BP, and stimulated primary productivity in the other lake at ca. 1,900–1,800 cal BP. Carbonate precipitation appears to have been suppressed when acidic poor fen and bog (palsa) communities dominated the catchment mire, and permafrost-induced changes in hydrology may further have affected the inflow of alkaline water from the catchment. Elevated contents of biogenic silica and diatom pigments in lake sediments during periods of poor fen and bog expansion further indicate that terrestrial vegetation influenced the amount of nutrients entering the lake. Increased productivity in the lake likely caused bottom-water anoxia in the downstream lake and led to recycling of sediment phosphorous, bringing the lake into a state of self-sustained eutrophication during two centuries preceding the onset of twentieth century permafrost thaw. Our results give insight into nutrient and permafrost dynamics in a subarctic wetland and imply that continued permafrost decay and related vegetation changes towards minerotrophy may increase carbon and nutrient storage of mire deposits and reduce nutrient fluxes in runoff. Rapid permafrost degradation may on the other hand lead to widespread mire erosion and to relatively short periods of significantly increased nutrient loading in adjacent lakes.     

Toward Improved Solar Irradiance Forecasts: a Simulation of Deep Planetary Boundary Layer with Scattered Clouds Using the Weather Research and Forecasting Model

Accurate forecasts of solar irradiance are required for electric utilities to economically integrate substantial amounts of solar power into their power generation portfolios. A common failing of numerical weather models is the prediction of scattered clouds at the top of deep PBL which are generally difficult to be resolved due to complicated processes in the planetary boundary layer. We improved turbulence parameterization for better predicting solar irradiance during the scattered clouds’ events using the Weather Research and Forecasting model. Sensitivity tests show that increasing the exchange coefficient leads to enhanced vertical mixing and a deeper mixed layer. At the top of mixed layer, an adiabatically ascending air parcel achieved the water vapor saturation and finally scattered cloud is generated.     

Simulating Gravity Changes in Topologically Realistic Driven Earthquake Fault Systems: First Results

Currently, GPS and InSAR measurements are used to monitor deformation produced by slip on earthquake faults. It has been suggested that another method to accomplish many of the same objectives would be through satellite-based gravity measurements. The Gravity Recovery and Climate Experiment (GRACE) mission has shown that it is possible to make detailed gravity measurements from space for climate dynamics and other purposes. To build the groundwork for a more advanced satellite-based gravity survey, we must estimate the level of accuracy needed for precise estimation of fault slip in earthquakes. We turn to numerical simulations of earthquake fault systems and use these to estimate gravity changes. The current generation of Virtual California (VC) simulates faults of any orientation, dip, and rake. In this work, we discuss these computations and the implications they have for accuracies needed for a dedicated gravity monitoring mission. Preliminary results are in agreement with previous results calculated from an older and simpler version of VC. Computed gravity changes are in the range of tens of μGal over distances up to a few hundred kilometers, near the detection threshold for GRACE.     

Testing benthic foraminiferal distributions as a contemporary quantitative approach to biomonitoring estuarine heavy metal pollution

Biomonitoring of estuarine pollution is the subject of active research, and benthic foraminifera are an attractive group to use for these purposes due to their ubiquitous presence in saline water and wide diversity. Here, we describe a case study of biomonitoring using benthic foraminifera in the French Mediterranean lagoon, Bages-Sigean lagoon. In this case, the major pollutants of interest are heavy metals in the sediment, particularly contaminated by Cu and Cd derived from industrial and agricultural sources. The foraminiferal assemblages of the Bages-Sigean lagoon are typical of normal paralic environments, but unusually almost completely lack agglutinated forms. The density of benthic foraminifera was shown to be more influenced by the sediment characteristics rather than heavy metal pollution. However, the relative abundance of Quinqueloculina bicostata was shown to increase in the most polluted areas and we propose that this taxon may be used as an indicator of heavy metal pollution.     

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.     

Coseismic slip model of the 2007 August Pisco earthquake (Peru) as constrained by Wide Swath radar observations

The Pisco earthquake ( M _w 8.0; 2007 August 15) occurred offshore of Peru's southern coast at the subduction interface between the Nazca and South American plates. It ruptured a previously identified seismic gap along the Peruvian margin. We use Wide Swath InSAR observations acquired by the Envisat satellite in descending and ascending orbits to constrain coseismic slip distribution of this subduction earthquake. The data show movement of the coastal regions by as much as 85 cm in the line-of-sight of the satellite. Distributed-slip model indicates that the coseismic slip reaches values of about 5.5 m at a depth of ∼18–20 km. The slip is confined to less than 40 km depth, with most of the moment release located on the shallow parts of the interface above 30 km depth. The region with maximum coseismic slip in the InSAR model is located offshore, close to the seismic moment centroid location. The geodetic estimate of seismic moment is 1.23 × 10²¹ Nm ( M _w 8.06), consistent with seismic estimates. The slip model inferred from the InSAR observations suggests that the Pisco earthquake ruptured only a portion of the seismic gap zone in Peru between 13.5° S and 14.5° S, hence there is still a significant seismic gap to the south of the 2007 event that has not experienced a large earthquake since at least 1687.     

Estimation and application of near‐surface full waveform redatuming operators

Land seismic data quality can be severely affected by near‐surface anomalies. The imprint of a complex near‐surface can be removed by redatuming the data to a level below the surface, from where the subsurface structures are assumed to be relatively smooth. However, to derive a velocity‐depth model that explains the propagation effects of the near‐surface is a non‐trivial task. Therefore, an alternative approach has been proposed, where the redatuming operators are obtained in a data‐driven manner from the reflection event related to the datum. In the current implementation, the estimation of these redatuming operators is done in terms of traveltimes only, based on a high‐frequency approximation. The accompanying amplitudes are usually derived from a local homogeneous medium, which is obviously a simplification of reality. Such parametrization has produced encouraging results in the past but cannot completely remove the near‐surface complexities, leaving artefacts in the redatumed results. In this paper we propose a method that estimates the redatuming operators directly from the data, i.e., without using a velocity model, in a full waveform manner, such that detailed amplitude and phase variations are included. The method directly outputs the inverse propagation operators that are needed for true‐amplitude redatuming. Based on 2D synthetic data it is demonstrated that the resulting redatuming quality is improved and artefacts are reduced.     

Direct observation of cosmic strings via their strong gravitational lensing effect – I. Predictions for high-resolution imaging surveys

We use present theoretical estimates for the density of long cosmic strings to predict the number of strong gravitational lensing events in astronomical imaging surveys as a function of the angular resolution and survey area. We show that angular resolution is the single most important factor, and that interesting limits on the dimensionless string tension   G μ/ c ²  can be obtained by existing and planned surveys. At the resolution of the Hubble Space Telescope ( HST ) (0.14 arcsec), it is sufficient to survey of the order of a few square degrees – well within reach of the current HST archive – to probe the regime   G μ/ c ²∼ 10⁻⁷  . If lensing by cosmic strings is not detected, such a survey would improve the limit on the string tension by a factor of two over that available from the cosmic microwave background. Future high resolution imaging surveys, covering a few hundred square degrees or more, either from space in the optical or from large-format radio telescopes on the ground, would be able to further lower this limit to   G μ/ c ²∼ 10⁻⁸  . These limits will not be improved significantly by increasing the solid angle of the survey.