Robert H. Spigel

The classification of Mixed-Layer Dynamics of Lakes of Small to Medium Size

Abstract An analysis of the time scales of processes relevant to wind mixing in lakes indicates that the response of a lake to wind may be classified into four regimes with respect to thermocline deepening behavior, depending on the relative sizes of the parameters describing wind strength, basin size and stratification. The dependence is described in terms of a mixed layer Richardson number and the aspect ratio of the mixed layer thickness to length. The classification is used to explain the diversity of phenomena reported in the literature for wind events in a number of different lakes and laboratory tanks which are either short enough or narrow enough for rotational effects to be unimportant. The classification is derived with reference to a two-layer, rectangular basin in the absence of Coriolis forces and surface heating. The classification is exended in a simple way to more realistic stratifications and basin shapes to predict the overall mixing features of a wind event. Response to wind varies from...

Mixing processes relevant to phytoplankton dynamics in lakes

Abstract Active turbulence in lakes is confined to the surface mixed layer, to boundary layers on the lake sides and bottom, and to turbulent patches in the interior. The density stratification present in most lakes fundamentally alters the pathways connecting external mechanical energy inputs, for example by wind, with its ultimate fate as dissipation to heat; the density stratification supports internal waves and intrusions that distribute the input energy throughout the lake. Intrusions may be viewed as internal waves with zero horizontal wavenumber and are formed each time localised mixing occurs in a stratified fluid. Intrusions are also formed in the epilimnion by differential heating or cooling and by differential deepening. The fraction of lake volume below the diurnal mixed layer that is subject to active turbulence is very small, probably of the order of 1% or less in small to medium‐sized lakes. By contrast, in the surface mixed layer, turbulence is less intermittent and maintains phytoplankton...

Hydrodynamic implications of large offshore mussel farms

Around 3000 Ha of New Zealands coastal waters are presently allocated to mussel farming of the Greenshell mussel Perna canaliculus mainly in small inshore farms of 30 000 to 80 000 m/sup 2/ each. The farms are constructed from blocks of parallel moored long-lines, with loops of spat-impregnated line hanging vertically in the upper 10-15 m of the water column (provided sufficient water depth). When fully laden, the mussel loops can present a rough diameter of 0.20 m or more to the flow. Industry expansion is likely to be in the form of larger offshore farms. Some of the larger proposed farms will have more than 1000 km of mussel line. This presents a potentially significant drag to currents and waves. In this paper, observations are used from one of the first large offshore farms constructed in New Zealand, in conjunction with scaling estimates of energy loss, to investigate the effect of the farm on waves, currents, and stratification. The farm, consisting of 230 long-lines arranged into 20 blocks, measured 650 m /spl times/ 2450 m. The hanging mussel loops at this site reached an average depth of 8 m, and water depth at midtide varied across the site from 10-12 m. Transects were taken through the farm using both conductivity, temperature, depth (CTD), and Acoustic Doppler Profiler (ADP) instruments. A second ADP was moored inside the farm providing velocity time series. The average current was reduced within the farm by 36%-63%. The moored ADP detected an undercurrent beneath the farm with velocities nearly twice that within the farm. Sharp changes in stratification coincide with the start of fully stocked mussel long-lines. Furthermore, other than at the very upstream end of the farm, stratification appears to have prevented significant vertical mixing between the fluid within the farm and that flowing under the farm. Wave energy was low during the 7-d deployment, with significant wave heights peaking at 0.25 m. The attenuation of transitional waves through the farm is examined by comparing simultaneous pressure sensor measurements from offshore and inshore of the farm. The observations show that wave energy attenuation was frequency dependent. Observed wave energy attenuation across the 650 m wide farm was approximately 5%, 10%, and 17% at wave frequencies of 0.1, 0.2, and 0.25 Hz, respectively.

Evolution of temperature and salt structure of Lake Bonney, a chemically stratified Antarctic lake

A resurgence of interest in the ecology of perennially ice-covered lakes in the McMurdo dry valleys has necessitated a review of our knowledge of the physical and chemical properties of these unusual lakes. Salinities in the ice-covered lakes cover a range from freshwater to hypersaline brines. Recent measurements of salt composition and concentrations in Lake Bonney reveal little change below the chemocline since extensive measurements made in 1960–1961, although lake level has risen by approximately 5 m since that time. The rise in lake level has resulted in a thickening of the freshwater layer above the chemocline. Temperature structure has adjusted to the effects of increased lake level on heat transfer processes such as transmission and absorption of solar radiation in the water column.Questions about how water-column stability affects biology in Lake Bonney have motivated the formulation of a method to compute density from in situ measurements of temperature, conductivity and pressure. Owing to high salt concentration and unique ion ratios, we modified the UNESCO Equation of State for seawater to predict density at salinities greater than 42. The modifications merge smoothly with the UNESCO equations at a salinity of 42. At salinities below 42 the UNESCO equations give excellent predictions of density.

Chaotic characteristics of the Southern Oscillation Index time series

Abstract The monthly time series of the Southern Oscillation Index (SOI) is analysed to examine its chaotic characteristics. Three schemes, moving average, low-pass filter and nonlinear smoothing, were used to reduce noise and enhance chaotic properties. Autocorrelation and spectral characteristics, as well as three chaos-oriented properties — phase space trajectory, the largest Lyapunov exponent and correlation dimension — were then examined. No significant signs of chaotic behaviour were found for either the noise-reduced SOI time series or the raw one. Although it contains long-term periodicity, the SOI time series is considered to be stochastic rather than chaotic from the viewpoint of dynamical systems theory.

Eutrophication processes regulated by a plunging river inflow

Lake Rotoiti (North Island, New Zealand) is a deep mesotrophic lake that has declined in water quality over the last 30 years. The main river entering the lake was identified as the primary enrichment source, but its interaction with the surrounding lakewater varied with season and time of day. During winter the river was colder than the lake and penetrated 6–8 km into the main basin as an underflow. In summer the river often entered as a plunging inflow during the early morning, but it warmed during the day, and in the afternoon entered the lake as a buoyant jet that flowed directly to a nearby outlet river. From continuous temperature measurements in the inflow and lake surface it was estimated that the river plunged and penetrated the lake as an interflow or underflow for 60.2% of the year. This translated into 31% of the N and 64% of the P loading on the lake. The river also injected phytoplankton into the main basin of Lake Rotoiti, including populations of bloom-forming cyanobacteria. The underflow was the dominant term in calculating the hydraulic flushing time of the main basin, and also made a large contribution of dissolved oxygen to the subsurface waters. These complex interactions between Lake Rotoiti and its inflowing river were controlled by small temperature differences ( < 3 °C, and had wide-ranging implications in the eutrophication process.

A microbiologically clean strategy for access to the Whillans Ice Stream subglacial environment

Abstract The Whillans Ice Stream Subglacial Access Research Drilling (WISSARD) project will test the overarching hypothesis that an active hydrological system exists beneath a West Antarctic ice stream that exerts a major control on ice dynamics, and the metabolic and phylogenetic diversity of the microbial community in subglacial water and sediment. WISSARD will explore Subglacial Lake Whillans (SLW, unofficial name) and its outflow toward the grounding line where it is thought to enter the Ross Ice Shelf seawater cavity. Introducing microbial contamination to the subglacial environment during drilling operations could compromise environmental stewardship and the science objectives of the project, consequently we developed a set of tools and procedures to directly address these issues. WISSARD hot water drilling efforts will include a custom water treatment system designed to remove micron and sub-micron sized particles (biotic and abiotic), irradiate the drilling water with germicidal ultraviolet (UV) radiation, and pasteurize the water to reduce the viability of persisting microbial contamination. Our clean access protocols also include methods to reduce microbial contamination on the surfaces of cables/hoses and down-borehole equipment using germicidal UV exposure and chemical disinfection. This paper presents experimental data showing that our protocols will meet expectations established by international agreement between participating Antarctic nations.

Model application to assess effects of urbanisation and distributed flow controls on erosion potential and baseflow hydraulic habitat

The catchment model SWMM was modified to include on-site flow-control devices, and then linked with models of hydraulic habitat suitability and erosive potential for specific reaches in a study catchment. Urbanisation decreased the modelled baseflow by up to 33% and reduced the area of suitable hydraulic habitat by 0.5–13.3%, depending on the reach and species. Hypothetical infiltration devices increased the baseflow to pre-development levels, with an associated increase in habitat, but other measures were not effective. Urbanisation increased the erosion potential index by a factor of 1.58–9.32, depending on the reach. Erosion-control ponds decreased the erosion potential to pre-development level in some reaches. Detention tanks reduced the erosion potential significantly, but not to pre-development levels. The poor predictive ability of the models for baseflow compromised the utility of the hydraulic habitat assessment. Predictions of the effects of urbanisation on baseflow habitat should therefore be treated cautiously.

Field calibration of a formula for entrance mixing of river inflows to lakes: Lake Taupo, North Island, New Zealand

Abstract Field measurements were used to validate predictions for the initial dilution of negatively buoyant, cold‐water inflows to Lake Taupo, as part of a study to quantify mixing processes associated with the two largest inflows to the lake. The predictions were made using a formulation originally derived for positively buoyant, warm‐water inflows to cooling ponds. The formulation predicts the total dilution of an inflow during its inertia‐dominated phase between its entrance to the lake and the point where buoyancy forces are great enough to cause the inflow to plunge and form a submerged density current. In one of the measured inflows, the inflowing jet was free to entrain from both sides; in the other, entrainment was restricted on one side by attachment of the inflowing jet to the shoreline of a bay just upstream of the plunge point. In the former example, the unmodified coefficients from the cooling pond formulation provided an excellent prediction of initial dilution. In the latter example, entrainment was reduced and different coefficients were derived. In both examples the inflows remained attached to the lake bed throughout their course until their liftoff at depths of 45–55 m to form interflows. The difference between coefficients for the two inflows indicates that the coefficient values should be considered site‐specific. The formulation is not valid for inflows that separate from the bottom of the inflow channel before plunging. The entrance mixing formulation was incorporated in a more general model of lake stratification, DYRESM, which already includes a well‐documented routine for routing underflows down submerged channels on the bed of a lake. Application of the model to the inflows measured in Lake Taupo gave good results for two model outputs that were not involved in the calibration of the entrance mixing formulation, but that are affected by the result of the initial dilution calculation—the temperatures in the river plume after it has plunged, and the insertion depth.

Coupling of Internal Wave Motion with Entrainment at the Density Interface of a Two-Layer Lake

Abstract A solution is presented to the problem of a two-layer rectangular basin subject to a suddenly applied, uniform wind stress; Coriolis effects are ignored. The solution is obtained for the case in which the time scales of internal wave motion, wave decay and entrainment are widely separated; in this range the entrainment across the interface is a perturbation on the mean motions. The solution includes an oscillatory initial response, followed by wave decay and a steady-state interface setup with baroclinic circulation, all superimposed on a slow deepening of the top layer by entrainment. The entrainment in turn affects the frequency of interfacial waves as the mixed layer deepens—the deepening alters top and bottom layer thickness and the density jump between them. For the range of mixed-layer Richardson number considered, entrainment is energized by wind stirring at the water surface, which does work at a rate proportional to u*3.