Water transport from the North-east Irish Sea to western Scottish coastal waters: Further observations from time-trend matching of Sellafield radiocaesium

Radiocaesium isotopes, discharged into the North-east Irish Sea from the Sellafield (formerly Windscale) nuclear fuel reprocessing plant in Cumbria, have been employed as flow monitors to update and extend the record of coastal water movement from the Irish Sea to the Clyde Sea area and, further north, to Loch Etive. The temporal trends in radiocaesium levels have been used to determine the extent of water mixing en route and to define mean advection rates. Flow conditions from the Irish Sea have changed considerably since the mid-1970s, the residence time of northern Irish Sea waters being ～12 months during 1978–1980 inclusive. Average transport times of four and six months are estimated for the Sellafield to Clyde and Sellafield to Etive transects respectively. Sellafield ¹³⁷Cs levels in seawater were diluted by factors of 27 and 50 respectively during current movement to the Clyde and Etive areas. The decrease in salinity-corrected ¹³⁷Cs concentrations between the Clyde and Etive suggests that dilution by Atlantic water occurs, the latter mainly entering the Firth of Lorne from the west. The majority (～94%) of the radiocaesium supply to Loch Etive enters the Firth of Lorne via the portion of the coastal current circulating west of Islay, only ～6% arriving via the Sound of Jura.     

Fast iterative solution of large undrained soil‐structure interaction problems

In view of rapid developments in iterative solvers, it is timely to re‐examine the merits of using mixed formulation for incompressible problems. This paper presents extensive numerical studies to compare the accuracy of undrained solutions resulting from the standard displacement formulation with a penalty term and the two‐field mixed formulation. The standard displacement and two‐field mixed formulations are solved using both direct and iterative approaches to assess if it is cost‐effective to achieve more accurate solutions. Numerical studies of a simple footing problem show that the mixed formulation is able to solve the incompressible problem ‘exactly’, does not create pressure and stress instabilities, and obviate the need for an ad hoc penalty number. In addition, for large‐scale problems where it is not possible to perform direct solutions entirely within available random access memory, it turns out that the larger system of equations from mixed formulation also can be solved much more efficiently than the smaller system of equations arising from standard formulation by using the symmetric quasi‐minimal residual (SQMR) method with the generalized Jacobi (GJ) preconditioner. Iterative solution by SQMR with GJ preconditioning also is more elegant, faster, and more accurate than the popular Uzawa method. Copyright © 2003 John Wiley & Sons, Ltd.     

Subtidal sea level and current variations in the northern reach of San Francisco Bay

Analyses of sea level and current-meter data using digital filters and a variety of statistical methods show a variety of phenomena related to non-local coastal forcing and local tidal forcing in the northern reach of San Francisco Bay, a partially mixed estuary. Low-frequency variations in sea level are dominated by non-local variations in coastal sea level and also show a smaller influence from tidally induced fortnightly sea level variations. Low-frequency currents demonstrate a gravitational circulation which is modified by changes in tidal-current speed over the spring-neap tidal cycle. Transients in gravitational circulation induce internal oscillations with periods of two to four days.     

Numerical study of enhanced energy dissipation near a seamount

Using a two-dimensional primitive equation model, we examine nonlinear responses of a semidiurnal tidal flow impinging on a seamount with a background Garrett-Munk-like (GM-like) internal wavefield. It is found that horizontally elongated pancake-like structures of high vertical wavenumber near-inertial current shear are created both in the near-field (the region over the slope of the seamount) and far-field (the region over the flat bottom of the ocean). An important distinction is that the high vertical wavenumber near-inertial current shear is amplified only at mid-latitudes in the far-field (owing to a parametric subharmonic instability (PSI)), whereas it is amplified both at mid-and high-latitudes (above the latitude where PSI can occur) in the near-field. In order to clarify the generating mechanism for the strong shear in the near-field, additional numerical experiments are carried out with the GM-like background internal waves removed. The experiments show that the strong shear is also created, indicating that it is not caused by the interaction between the background GM-like internal waves and the semidiurnal internal tides. One possible explanation is proposed for the amplification of high vertical wavenumber near-inertial current shear in the near-field where tide residual flow resulting from tide-topography interaction plays an important role in transferring energy from high-mode internal tides to near-inertial internal waves.     

The Record of an Extreme Flood in the Sediments of Montane Lillooet Lake, British Columbia: Implications for Paleoenvironmental Assessment

Severe rainfall in mid October, 2003 produced the largest floods in almost a century of record on rivers in the Cordillera of southwestern British Columbia. Sediment deposited in Lillooet Lake as a result of this event is clearly distinguished by stratigraphy, colour, texture, magnetic properties, and organic content. Each of these physical properties is related to the lacustrine processes, especially turbid underflow, that distributed the sediment through the lake. The flood, which lasted less than a week, delivered 8–12 times the amount of sediment that accumulates in most entire years in the deepest, central parts of the lake. Recognition of events of this type in the stratigraphic record offers a means of assessing the changing nature of extreme hydroclimatic events, and their relation to more ubiquitous, lower-energy processes.     

Quantitative analysis of variations in depositional sequence thickness from submarine channel levees

Abstract Thickness variations across‐levee and downchannel in acoustically defined depositional sequences from six submarine channel‐levee systems show consistent and quantifiable patterns. The thickness of depositional sequences perpendicular to the channel trend, i.e. across the levee, decreases exponentially, as characterized by a spatial decay constant, k. Similarly, the thickness of sediment at the levee crest decreases exponentially down the upper reaches of submarine channels and can be characterized by a second spatial decay constant, λ. The inverse of these decay constants has units of length and defines depositional length scales such that k^?1 is a measure of levee width and λ^?1 is a measure of levee length. Quantification of levee architecture in this way allowed investigation of relationships between levee architecture and channel dimensions. It was found that these measures of levee e‐folding width and levee e‐folding length are directly related to channel width and relief. The dimensions of channels and levees are thus intimately related, thereby limiting the range of potential channel‐levee morphologies, regardless of allocyclic forcing. A simple sediment budget model relates the product of the levee e‐folding width and e‐folding length to through‐channel volume discharge. A classification system based on the quantitative downchannel behaviour of levee architecture allows identification of a ‘mid‐channel’ reach, where sediment is passively transferred from the through‐channel flow to the levees as an overspilling flow. Downstream from this reach, the channel gradually looses its control on guiding turbidity currents, and the resulting flow can be considered as an unconfined or spreading flow.     

Reconciling pyroclastic flow and surge: the multiphase physics of pyroclastic density currents

Two end-member types of pyroclastic density current are commonly recognized: pyroclastic surges are dilute currents in which particles are carried in turbulent suspension and pyroclastic flows are highly concentrated flows. We provide scaling relations that unify these end-members and derive a segregation mechanism into basal concentrated flow and overriding dilute cloud based on the Stokes number (S_T), the stability factor (Σ_T) and the dense-dilute condition (D_D). We recognize five types of particle behaviors within a fluid eddy as a function of S_T and Σ_T: (1) particles sediment from the eddy, (2) particles are preferentially settled out during the downward motion of the eddy, but can be carried during its upward motion, (3) particles concentrate on the periphery of the eddy, (4) particles settling can be delayed or ‘fast-tracked’ as a function of the eddy spatial distribution, and (5) particles remain homogeneously distributed within the eddy. We extend these concepts to a fully turbulent flow by using a prototype of kinetic energy distribution within a full eddy spectrum and demonstrate that the presence of different particle sizes leads to the density stratification of the current. This stratification may favor particle interactions in the basal part of the flow and D_D determines whether the flow is dense or dilute. Using only intrinsic characteristics of the current, our model explains the discontinuous features between pyroclastic flows and surges while conserving the concept of a continuous spectrum of density currents.     

Time-dependent, coupled, Ekman boundary layer solutions incorporating Stokes drift

Recent technological advances in current measuring devices has resulted in a large observational database related to wind-driven motions in the upper ocean mixed layer. This has served to highlight the fact that transient motions make up a substantial contribution of the resulting Ekman currents. At the same time, certain discrepancies have emerged between the observed angular deflections of the steady-state currents from the surface wind stress, both at the surface and at sub-surface depths, which cannot be reconciled using the classical Ekman model. This paper seeks to tackle these two issues.First a general analytical method is presented for solving the time dependent horizontal momentum Ekman equations. Analysis of the unsteady terms that arise from simple special cases shows how the evolution proceeds through three stages. At early times, the Coriolis acceleration is insignificant, and the current is unidirectional and deepens through downward diffusion of momentum. Later Coriolis acceleration deflects the current vectors in the upper layers, whilst downward diffusion of momentum continues to deepen the layer. Finally, once diffusion has penetrated down to the depth of the steady-state current, then the transients decay on the inertial or diffusive timescale, depending upon the boundary conditions of the particular problem.In the second half of the paper, a new steady-state model is developed that includes the effects of wind-generated waves, through the action of their Stokes drift on the planetary vorticity. Comparisons between observations and the theoretical predictions, demonstrate that inclusion of the Stokes drift is the key to reconciling the discrepancies in the angular deflections of the steady-state currents. This leads to the conclusion that Ekman layer currents are significantly influenced by the surface waves.