Nutrient transports in a Swedish estuary

Himmerfjärd is a Swedish estuary bordering on the Baltic. The estuary lacks astronomical tides and its circulation is driven by winds and freshwater runoff. Because of a tertiary sewage treatment plant located at its inner end, the estuary is becoming increasingly eutrophic. A field study was carried out for a 78-day period in late summer and early fall of 1977 to determine rates of nutrient transport and to construct nutrient budgets. Since physical parameters (current velocity, temperature, salinity, winds and water level changes) were measured more frequently than nutrients (phosphate, nitrate, nitrite and ammonium) it was necessary to develop a suitable method to calculate nutrient flux time series and net nutrient fluxes. Over the study period, Himmerfjärd imported phosphorus and exported nitrogen. Direction of nutrient fluxes and changes in flux direction were consistent with the structure of the baroclinic currents.     

Mezen Bay—a key area for understanding Weichselian glaciations in northern Russia

Sediment successions in coastal cliffs around Mezen Bay, southeastern White Sea, record an unusually detailed history of former glaciations, interstadial marine and fluvial events from the Weichselian. A regional glaciation model for the Weichselian is based on new data from the Mezen Bay area and previously published data from adjacent areas. Following the Mikulinian (Eemian) interglacial a shelf‐centred glaciation in the Kara Sea is reflected in proglacial conditions at 100–90 ka. A local ice‐cap over the Timan ridge existed between 75 and 65 ka. Renewed glaciation in the Kara Sea spread southwestwards around 60 ka only, interrupted by a marine inundation, before it advanced to its maximum position at about 55–50 ka. After a prolonged ice‐free period, the Scandinavian ice‐sheet invaded the area from the west and terminated east of Mezen Bay about 17 ka. The previously published evidence of a large ice‐dammed lake in the central Arkhangelsk region, Lake Komi, finds no support in this study. Copyright © 2003 John Wiley & Sons, Ltd.     

On the structure of oscillatory boundary layers

An empirical analysis is performed on the most detailed, recent measurements of turbulent oscillatory boundary layer flow. The measurements show that throughout elevations where the flow can be considered horizontally uniform, one deficit model is sufficient for describing the fundamental mode. Some general properties of the non dimensional velocity deficit D₁(z) appear with striking consistency. First of all the identity $\text{ln}\text{¦D}{}_{\mathrm{<mtext>1</mtext>}}\text{¦ ≡}\text{Arg}\text{D}{}_{\mathrm{<mtext>1</mtext>}}$, which is a theoretical result for smooth laminar flow, seems to hold with great accuracy for a large range of turbulent flow conditions as well. This is of principal theoretical interest because all previous analytical eddy viscosity models as well as numerical mixing length models predict a consistent and fairly large difference between Arg D₁ and $\text{ln}\text{¦D}{}_{\mathrm{<mtext>1</mtext>}}\text{¦}$. If the identity between $\text{ln}\text{¦D}{}_{\mathrm{<mtext>1</mtext>}}\text{¦}$ and Arg D₁ extends all the way to the bed, it means that the bed shear stress leads the free stream velocity by 45 degrees. It is also found that the structure of smooth turbulent oscillatory flows as measured by Kalkanis (1964) corresponds to a sharp maximum in the normalized energy dissipation rate.     

Weak-contrast approximation of the elastic scattering matrix in anisotropic media

The plane-wave reflection and transmission coefficients at a plane interface between two anisotropic media constitute the elements of the elastic scattering matrix. For a 1-D anisotropic medium the eigenvector decomposition of the system matrix of the transformed elasto-dynamic equations is used to derive a general expression for the scattering matrix. Depending on the normalization of the eigenvectors, the expressions give scattering coefficients for amplitudes or for vertical energy flux.Computing the vertical slownesses and the corresponding polarizations, the eigenvector matrix and its inverse can be found. We give a simple formula for the inverse, regardless of the normalization of the eigenvectors. When the eigenvectors are normalized with respect to amplitudes of displacement (or velocity), the calculation of the scattering matrix for amplitudes is simplified.When the relative changes in all parameters are small, a weak-contrast approximation of the scattering matrix, based on the exactly determined polarization vectors in an average medium, is obtained. The same approximation is also derived directly from the transformed elasto-dynamic equations for a smooth vertically inhomogeneous medium, proving the consistency of the approximation.For monoclinic media, with the mirror symmetry plane parallel to the interface, the approximative scattering matrix is given in terms of analytic expressions for the non-normalized eigenvectors and vertical slownesses. For transversely isotropic media with a vertical axis of symmetry (VTI) and isotropic media, explicit solutions for the weak-contrast approximations of the scattering matrices have been obtained. The scattering matrix for amplitudes for isotropic media is well known. The scattering matrix for vertical energy flux may have applications in AVO analysis and inversion due to the reciprocity of the reflection coefficients for converted waves.Numerical examples for monoclinic and VTI media provide good agreement between the approximative and the exact reflection matrices. It is, however, expected that the approximations cannot be used when the symmetry properties of the two media are very different. This is because the approximation relies on a small relative contrast between the eigenvectors in the two media.Presented at the Workshop Meeting on Seismic Waves in Laterally Inhomogeneous Media, Castle of Trest, Czech Republic, May 22–27, 1995.     

Post-mid-Cretaceous eastern North Sea evolution inferred from 3D thermo-mechanical modelling

The Late Cretaceous–Cenozoic evolution of the eastern North Sea region is investigated by 3D thermo-mechanical modelling. The model quantifies the integrated effects on basin evolution of large-scale lithospheric processes, rheology, strength heterogeneities, tectonics, eustasy, sedimentation and erosion.

The evolution of the area is influenced by a number of factors: (1) thermal subsidence centred in the central North Sea providing accommodation space for thick sediment deposits; (2) 250-m eustatic fall from the Late Cretaceous to present, which causes exhumation of the North Sea Basin margins; (3) varying sediment supply; (4) isostatic adjustments following erosion and sedimentation; (5) Late Cretaceous–early Cenozoic Alpine compressional phases causing tectonic inversion of the Sorgenfrei–Tornquist Zone (STZ) and other weak zones.

The stress field and the lateral variations in lithospheric strength control lithospheric deformation under compression. The lithosphere is relatively weak in areas where Moho is deep and the upper mantle warm and weak. In these areas the lithosphere is thickened during compression producing surface uplift and erosion (e.g., at the Ringkøbing–Fyn High and in the southern part of Sweden). Observed late Cretaceous–early Cenozoic shallow water depths at the Ringkøbing–Fyn High as well as Cenozoic surface uplift in southern Sweden (the South Swedish Dome (SSD)) are explained by this mechanism.

The STZ is a prominent crustal structural weakness zone. Under compression, this zone is inverted and its surface uplifted and eroded. Contemporaneously, marginal depositional troughs develop. Post-compressional relaxation causes a regional uplift of this zone.

The model predicts sediment distributions and paleo-water depths in accordance with observations. Sediment truncation and exhumation at the North Sea Basin margins are explained by fall in global sea level, isostatic adjustments to exhumation, and uplift of the inverted STZ. This underlines the importance of the mechanisms dealt with in this paper for the evolution of intra-cratonic sedimentary basins.     

Phytoplankton, nutrients, and transparency in Danish coastal waters

We present a comparative analysis of 1400 data series of water chemistry (particularly nitrogen and phosphorus concentrations), phytoplankton biomass as chlorophylla (chla) concentrations, concentrations of suspended matter and Secchi depth transparency collected from the mid-1980s to the mid-1990s from 162 stations in 27 Danish fjords and coastal waters. The results demonstrate that Danish coastal waters were heavily eutrophied and had high particle concentrations and turbid waters. Median values were 5.1 μg chla 1⁻¹, 10.0 mg DW 1⁻¹ of suspended particles, and Secchi depth of 3.6 m. Chlorophyll concentration was strongly linked to the total-nitrogen concentration. The strength of this relationship increased from spring to summer as the concentration of total nitrogen declined. During summer, total nitrogen concentrations accounted for about 60% of the variability in chlorophyll concentrations among the different coastal systems. The relationship between chlorophyll and total phosphorus was more consistant over the year and correlations were much weaker than encountered for total nitrogen. Secchi depth could be predicted with good precision from measurements of chlorophyll and suspended matter. In a multiple stepwise regression model with In-transformed values the two variables accounted for most of the variability in water transparency for the different seasons and the period March–October as a whole (c. 80%). We were able to demonstrate a significant relationship between total nitrogen and Secchi depth, with important implications for management purposes.     

Neotectonics, sea-level changes and biological evolution in the Fennoscandian Border Zone of the southern Kattegat Sea

Shallow seismic data and vibrocore information, sequence stratigraphic and faunal evidence have been used for documentation of Late Weichselian reactivation of faulting in the south central Kattegat, southern Scandinavia. The study area is situated on the Fennoscandian Border Zone, where tectonic activity has been recurrent since Early Palaeozoic time and still occurs, as shown by present earthquake activity. New data from the area south of the island of Anholt show that after deglaciation fast isostatic rebound resulted in reactivation of a NW-SE striking normal fault system. This tectonic episode is dated to a period starting shortly before 15.0 cal. ka BP and ending around 13.5 cal. ka BP, after regression had already reached a level of about 30 m b.s.l. The vertical displacement associated with the faulting was in the order of 20 m. More generally, the results support the previously reported late Weichselian sea-level highstand, which was followed by forced regression until the eustatic sea-level rise surpassed the rate of glacio-isostatic rebound in early Preboreal. Our findings further imply that drainage of the Baltic Ice Lake through the Øresund at c. 15 cal. ka BP (Bergsten & Nordberg 1992) may have been triggered by tectonic activity in this region.     

Erratum: Comparison of Seismic Dispersion and Attenuation Models: Stud. Geophys. Geod., 46(2002), 293-320

Studia Geophysica et Geodaetica -     

Reduction of ground vibration by means of barriers or soil improvement along a railway track

Trains running in built-up areas are a source to ground-borne noise. A careful design of the track may be one way of minimizing the vibrations in the surroundings. For example, open or infilled trenches may be constructed along the track, or the soil underneath the track may be improved. In this work, the influence of the track design and properties on the level of ground vibration due to a vehicle moving with subsonic speed is examined. A coupled finite element-boundary element model of the track and subsoil is employed, adopting a formulation in the moving frame of reference following the vehicle. The computations are carried out in the frequency domain for various combinations of the vehicle speed and the excitation frequency. The analyses indicate that open trenches are more efficient than infilled trenches or soil stiffening–even at low frequencies. However, the direction of the load is of paramount importance. For example, the response outside a shallow open trench may change dramatically when horizontal load is applied instead of vertical load.