Quantifying turbulent mixing and oxygen fluxes in a Mediterranean-type,microtidal estuary

Experiments were carried out on an intermittent estuary during its closed (summer) and open (winter) states to identify the physical processes responsible for vertical mixing across the halocline, and to quantify vertical fluxes of oxygen and salt between water layers. During the blocked phase a two-layer structure was observed, with a brackish surface layer overlying old seawater. Within a deep basin the wind-driven turbulent mixing was consistent with the measured surface-layer turbulent dissipation, but the dissipation in the bottom layer appeared to be driven by internal seiching. In the shallow regions of the estuary vertical fluxes of dissolved oxygen were indicative of oxygen demand by respiration and remineralization of organic material in bottom water and sediments. During the estuary's open phase a three-layer structure was observed, having a fresh, river-derived surface layer, a middle layer of new seawater, and a bottom layer of old seawater. In the shallower regions surface-layer turbulent diffusion was consistent with the strong, gusty winds experienced at the time. The dissolved oxygen of the incoming seawater decreased to very low values by the time it reached the upstream deep basin as a result of the low cross-pycnocline oxygen flux being unable to compensate for the oxygen utilization. At least 50 % of the cross-pycnocline salt fluxes in the shallow reaches of the open estuary are suggested to be driven by Holmboe instabilities. Responsible Editor: Hans Burchard     

A test of the influence of tidal stream polarity on the structure of turbulent dissipation

Vertical mixing by the tides plays a key role in controlling water column structure over the seasonal cycle in shelf seas. The influence of tidal stirring is generally well represented as a competition between surface buoyancy input and the production of turbulent kinetic energy (TKE) by frictional stresses, a competition which is encapsulated in the Qh/u³ criterion. An alternative control mechanism arises from the limitation of the thickness of the bottom boundary layer due to the effects of rotation and the oscillation of the flow. Model studies indicate that, for conditions typical of the European shelf seas, the energy constraint exerts the dominant control but that for tidal streams with large positive polarisation (i.e. anti-clockwise rotation of velocity vector), some influence of rotation in limiting mixing should be detectable. We report here measurements of flow structure (with ADCPs) and turbulent dissipation (FLY Profiler) made at two similar locations in the Celtic Sea which differ principally in that the tidal currents rotate in opposite senses with approximately equal magnitude (polarity P=±0.6). A clear contrast was observed between the two sites in the vertical structure of the currents, the density profile and the rate of dissipation of TKE. At the positive polarity (PP) site (P≈+0.6), the bottom boundary layer in the tidal flow was limited to ∼20 mab (metre above the bed) and significant dissipation from bottom boundary friction was constrained within this layer. At the negative polarity (NP) site (P≈−0.6), the dominant clockwise rotary current component exhibited a velocity defect (i.e. reduction relative to the free stream) extending into the upper half of the water column while significant dissipation was observed to penetrate much further up the water column with dissipation levels ∼10^−4.5 W m⁻³ reaching to the base of the pycnocline at 70–80 mab. These contrasting features of the vertical distribution of dissipation are well reproduced by a 1-D model when run with windstress and tidal forcing and using the observed density profile. Model runs with reversed polarity at the two sites, support the conclusion that the observed contrast in the structure of tidal velocity, dissipation and stratification is due to the influence of tidal stream polarity. Increased positive polarity reduces the upward penetration of mixing which allows the development of stronger seasonal stratification, which, in turn, further inhibits vertical mixing.     

Mixing processes in lakes: mechanisms and ecological relevance

In Lake Baldegg, Switzerland (surface area 5.3 km², maximum depth 66 m) the analysis of data from moored instrument systems (atmospheric boundary layer, lake temperature distribution, bottom currents) was correlated to the long-term development of vertical mixing as seen from profiles of natural isotopes (radon-222, tritium and helium-3) and chemical species. The investigation shows: 1. Vertical mixing coefficients below 25 m are small. Consequently the vertical concentration distribution of sediment emanating species in the deep hypolimnion is controlled by the bottom topography. 2. Renewal of deep hypolimnic water is significant even during stratification. 3. Weakly damped internal waves characterize the internal dynamics during stratification. 4. Horizontal bottom currents play an important role in the hypolimnion mixing and can be correlated to internal waves during stratification.     

Physical conditions driving the spatial and temporal variability in aquatic metabolism of a subtropical coastal lake

We investigated the importance of meteorological and lake physical conditions for temporal, horizontal and vertical differences in the concentration of dissolved oxygen (DO) and water temperature, and the derived daily estimates of gross primary production (GPP), ecosystem respiration (R) and net ecosystem production (NEP). Our study was conducted in a subtropical and polymictic lake in Southern Brazil, during a spring–summer transition. Metabolic rates were determined from two sites using the open water oxygen technique. At the central deep site, oxygen sondes were deployed at three depths to assess patterns in vertical variability. During 10 days, an additional DO and temperature sonde was placed near the shoreline allowing us to compare metabolic differences in the surface layers between the central pelagic and littoral site. While GPP was similar, R was significantly higher at the shallower littoral site, causing NEP to be lower, although NEP was still positive. The littoral site had less diel changes in DO and higher daily variability in all metabolic rates. Variability in GPP and R at the littoral site was related to temperature, wind speed and rainfall suggesting that short-term variability in metabolic rates in shallow areas are sensitive to resuspension of sediments caused by a less stable water column. A clear vertical gradient was furthermore found for the metabolic rates at the central deep part of the lake, related to the light extinction, with highest GPP around 0.3 m and decreasing with depth, while respiration showed the inverse pattern. Below subsurface, respiration prevailed at 5.0 m depth and was uncoupled to primary production. Under conditions with high light and temperature, and low wind speeds, the mixing depth became shallower, in turn increasing the water column stability at the deep pelagic site, which resulted in higher mean light available and higher GPP in the water column. Our results confirm that deployment of sensors in different sites and depths allows for spatially, as well as temporally more representative estimates of lake metabolism.     

Darcian Flow Characteristics Upgradient of a Kettle Pond Determined by Direct Ground Water Flow Measurement

In glacial outwash deposits, the movement of ground water Is determined by small scale irregularities in the pattern of hydraulic conductivity. Permeability determinations on split spoon samples obtained from coring the site are not sufficient to predict the patchiness of flow since it cannot define continuity of the strata. The lattice work pattern can be determined by vertical profiling with direct ground water flow measurement. The rate and direction of flow is combined with head gradient changes to compute hydraulic conductivity changes across the site.
The results of the tests can be plotted on triangular graphs depicting the fundamental Darcy equation. The local conditions reflect a mathematical "patchiness" of hydraulic conductivity unique to outwash deposits.
The procedure was employed to determine flow characteristics and define the zone of contribution to porous bottom kettle lakes. The zone of contribution was defined by projecting backward from the vertical profiling and shallow measurements and taking into account the daily rain water recharge rate across the site.
For the unconfined aquifer north of the pond, shallow ground water flow measurements were necessary to define the recharge portion of the shoreline. Vertical profiling was required to define the recharge volume since the rate of flow was not even with depth. A simple differential equation for determining the recharge area is presented along with the calculations.     

Formation and Evolution of the Anaerobic Zone in the Black Sea from the Most Recent Salinization to the Present State: Retrospective Estimation Based on Mathematical Modeling

It is established that the formation and evolution of the anaerobic zone of the Black Sea are associated with the beginning and development of the most recent salinization of the sea and the formation of pycnocline at intermediate depths; a deterioration of deep-water aeration and the formation of anaerobic conditions in deeper layers; the rate of sulfate-reduction first in the near-bed layer and later, as oxygen is depleted, in the water mass. Formalization of these processes based on refined present-day data on water balance enabled the reproduction of profiles of water salinity and vertical-exchange coefficient for different formation stages of salinity regime in the Black Sea. The vertical distribution of oxygen and hydrogen sulfide is reconstructed, and the rates of oxygen consumption in water column and sulfate reduction in the near-bed layer and at the upper boundary of the anaerobic zone are evaluated in numerical experiments. The obtained data show the transformation of the vertical distribution of oxygen and hydrogen sulfide in the Black Sea from the beginning of its most recent salinization to the present-day state. It is shown that the anaerobic zone rises from 2000 to 200 m within 500–600 years (in the period 3.9–4.5 Ka from the beginning of water exchange through the Bosphorus), and next the upper boundary of the anaerobic zone slowly ascends up to its present-day position (130–180 m). Mathematical modeling was used to evaluate the fluxes of oxygen and hydrogen sulfide at different formation stages of the anaerobic zone in the sea.__________Translated from Vodnye Resursy, Vol. 32, No. 3, 2005, pp. 307–321.Original Russian Text Copyright © 2005 by Leonov, Shaporenko.     

An instrument system for the investigation of particle fluxes

We present the rationale, design, and use of an instrument system to measure the variability of vertical and horizontal particle fluxes. The system features a new sequentially sampling sediment trap which collects and seals 10 separate samples during a single deployment. Horizontal particle fluxes are simultaneously monitored with a beam transmissometer interfaced to a standard Aanderaa current meter. Results from a 10-week deployment of instruments at several depths in a deep fjord estuary indicate that the trapping rate increases from 0.5g m⁻² day⁻¹ at 20 m to 150g m⁻² day⁻¹ at 200 m (5 m above bottom) because of frequent erosion in the deep waters. Periodic flushing of the deep water by intrusions of marine water over the seaward sill markedly enhances erosion and causes an up-estuary particle transport comparable to the vertical particle flux originating at the surface.     

MELIMEX,an experimental heavy metal pollution study: Vertical mixing in a large limno-corral

In a limno-corral (diameter 12 m, depth to sediments 10 m), located in Baldeggersee (Switzerland), vertical mixing has been measured during more than one year and compared to the conditions in the open lake (maximum depth 65 m, surface area 5.3 km²). The temperature method by McEwen and Hutchinson yields K_z values between 5×10⁻² cm²s⁻¹ at the upper boundary of the thermocline and 2×10⁻³ cm²s⁻¹ at the bottom, a value near the molecular diffusion of heat at 4°C (1.36×10⁻³ cm²s⁻¹). K_z calculated from profiles of excess radon-222 generally agree with those from the temperature data. Compared to the open lake, the corral has a more shallow epilimnion. However, during calm meteorological conditions, vertical mixing in the upper 10 m is similar outside and inside the corral. Supported by the Swiss National Science Foundation within the framework for its National Research Program on ‘Lake Currents.’     

ANALYSIS OF AMPLIFICATION CHARACTERISTICS OF GROUND MOTIONS IN THE HEAVILY DAMAGED BELT ZONE DURING THE 1995 HYOGO-KEN NANBU EARTHQUAKE

To estimate the amplification characteristics of ground motions in the heavily damaged belt zone in Kobe City during the 1995 Hyogo-ken Nanbu earthquake, 3D wave propagation analyses of a 2D deep irregular underground structure model with a vertical discontinuity were performed at an early stage as a preliminary and qualitative study. The hyperelement method was applied to the analyses for incident plane waves expected from the wavefields due to the source mechanism. The observation records at Kobe University of the rock site were used as control motions. The ground motions on the engineering bedrock (assumed to be on the free surface of the Osaka group layers having a shear velocity of 500 m/s) and at ground surface were calculated. The effects of the deep irregular underground structure and shallow surface layers on the ground motion amplification are discussed. Although there are qualifications due to the uncertain characteristics of the input rock motion and shear wave velocities of the underground structure, the analytical results show that the ground motion in the heavily damaged belt zone were amplified due to the focusing effect of the deep irregular underground structure as well as the shallow surface layers, and that the calculated peak ground acceleration (PGA) distribution coincided closely with the distributions of structural damage. © 1997 by John Wiley & Sons, Ltd.     

The short term response to eutrophication abatement

We compare results of a new model for predicting the short term inter annual changes in chlorophyll-a (chl-a) in lakes after reductions in total phosphorus (TP) to predictions made by least squares regression models. In the new method, slopes of chl-a/TP graphs (both axes in mg · m^–3) are depicted in frequency diagrams and used to extract information on the expected, short term chl-a/TP response. The short term response for nine shallow (< 10 m deep) and nutrient rich lakes to changes in TP was found to be: Chl-a = 0.49 · TP + 17.3, and for nine deep, P-limited lakes: Chl-a = 0.08 · TP + 3.5. If the TP-reduction is known to be greater than 10 mg · m^–3, the expected slope increases to 0.58 for shallow lakes and to 0.26 for deep lakes. The slope, 0.58, is 8% lower than the slope for the long term response calculated by regression for the shallow lakes. For deep lakes the slope, 0.26, is 2 to 3 times higher than that calculated by regression, indicating that reductions in TP for deep lakes give greater effects than least squares regression equations suggest. We have also calculated the reduction in TP which will give about 80% probability that a reduction in chl-a will be observed next year. For shallow, P-limited lakes this reduction is about 30 mg · m^–3 (5% of average initial in-lake TP concentration), and for deep lakes about 14 mg · m^–3 (35% of average initial in-lake TP concentration).     

Inclined internal tide waves at a narrow Mexican Pacific shelf

The dynamics of a semidiurnal internal tidal wave at a narrow Mexican Pacific shelf is discussed using the data of temperature obtained by an anchored instrument and data of field surveys. The internal tide on the shelf is dominated by an inclined wave, which propagates upward and onshore along a continental slope. Despite its reflection from the bottom and from the surface of the ocean, they remain inclined and totally destroyed over the course of one wavelength. Due to wave reflection from the inclined bottom, the horizontal and vertical wave number increase threefold when the wave goes into shallow waters. The wave undergoes nonlinear transformation and overturns forming several homogeneous temperature layers up to 20 m thick. The most intense disturbances of water layers are observed near the bottom, where the slope angle approaches its critical value. Because of nonlinear effects, the wave carries cool deep water out to the shallow depth and causes coastal upwelling. Intense solar warming together with vertical mixing results in a rapid rise of temperature in the 130-m water column that was observed.     

Sedimentary inputs and oxygen uptake by the sediment in Lake Geneva (Switzerland)

Sediment cores and sediment traps were collected twice a month in two 35 m deep stations of Lake Geneva (Switzerland). The organic input sedimenting to the bottom is equal to 157 g C m⁻²y⁻¹ in station 1, to 214 g C in station 2. In spite of this difference, the oxygen uptake by the sediment (OUS) is similar in both locations (46–47 g C m⁻²y⁻¹). The oxygen uptake by the matter sedimenting to the bottom (OUSM) is respectively 45 g C m⁻²y⁻¹ and 41 g C in stations 1 and 2. The equivalence between OUS and OUSM implies that most of the sedimented matter arriving to the bottom is directly oxidized at the sediment surface. In station 1, OUS is positively correlated to OUSM, and OUSM is positively correlated to chlorophyll-a concentrations in the water column (0–20 m) one week before sediment sampling. In location 2, OUS is positively correlated to the percentage of organic carbon and nitrogen in the sedimented matter, negatively to its C:N ratio. Increasing allochthonous inputs have a negative influence on benthic respiration. At both sites, OUS is not directly related to macrobenthic biomass or to temperature of bottom water.     

Pb in the western Indian Ocean: distribution, disequilibrium, and partitioning between dissolved and particulate phases

²²⁶Ra profiles have been measured in the western Indian Ocean as part of the 1977–1978 Indian Ocean GEOSECS program. These profiles show a general increase in deep and bottom water Ra concentration from the Circumpolar region to the Arabian Sea. A deep Ra maximum which originates in the Arabian Sea and in the Somali basin at about 3000 m depth spreads southward into the Mascarene basin and remains discernible in the Madagascar and Crozet basins. In the western Indian Ocean, the cold Antarctic Bottom Water spreads northward under the possibly southward-flowing deep water, forming a clear benthic front along the Crozet basin across the Southwest Indian Ridge into the Madagascar and Mascarene basins. The Antarctic Bottom Water continues to spread farther north to the Somali basin through the Amirante Passage at 10°S as a western boundary current. The benthic front and other characteristic features in the western Indian Ocean are quite similar to those observed in the western Pacific where the benthic front as a distinctive feature was first described by Craig et al. [15]. Across the Mid-Indian Ridge toward the Ceylon abyssal plain near the triple junction, Ra profiles display a layered structure, reflecting the topographic effect of the mid-ocean ridge system on the mixing and circulation of the deep and bottom waters. Both Ra and Si show a deep maximum north of the Madagascar basin. Linear relationships between these two elements are observed in the deep and bottom water with slopes increasing northward. This suggests a preferential input of Ra over Si from the bottom sediments of the Arabian Sea and also from the flank sediments of the Somali basin.     

Horizontal ground deformation patterns and magma storage during the Puu Oo eruption of Kilauea volcano,Hawaii: episodes 22–42

Horizontal ground deformation measurements were made repeatedly with an electronic distance meter near the Puu Oo eruption site approximately perpendicular to Kilauea's east rift zone (ERZ) before and after eruptive episodes 22–42. Line lengths gradually extended during repose periods and rapidly contracted about the same amount following eruptions. The repeated extension and contraction of the measured lines are best explained by the elastic response of the country rock to the addition and subsequent eruption of magma from a local reservoir. The deformation patterns are modeled to constrain the geometry and location of the local reservoir near Puu Oo. The observed deformation is consistent with deformation patterns that would be produced by the expansion of a shallow, steeply dipping dike just uprift of Puu Oo striking parallel to the trend of the ERZ. The modeled dike is centered about 800 m uprift of Puu Oo. Its top is at a depth of 0.4 km, its bottom at about 2.9 km, and the length is about 1.6 km; the dike strikes N65° E and dips at about 87°SE. The model indicates that the dike expanded by 11 cm during repose periods, for an average volumetric expansion of nearly 500 000 m³. The volume of magma added to the dike during repose periods was variable but correlates positively with the volume of erupted lava of the subsequent eruption and represents about 8% of the new lava extruded. Dike geometry and expansion values are used to estimate the pressure increase near the eruption site due to the accumulation of magma during repose periods. On average, vent pressures increased by about 0.38 MPa during the repose periods, one-third of the pressure increase at the summit. The model indicates that the dikelike body below Puu Oo grew in volume from 3 million cubic meters (Mm³) to about 10–12 Mm³ during the series of eruptions. The width of this body was probably about 2.5–3.0 m. No net long-term deformation was detected along the measured deformation lines.     

Turbulence measurements in a submarine canyon

Profiles of velocity turbulence in Monterey Canyon, made with a recently developed expendable probe, show the existence of a very turbulent bottom boundary layer. The turbulent flow is up to 170 m thick and has peak microscale shears of 1 m s⁻¹ per meter. The rate of dissipation of kinetic energy, based on the observed shear variance, averaged over the depth of the turbulent boundary layer ranged from 70 to 500 × 10⁻⁶W m⁻³. Temperature measurements indicate that the flow was up canyon at a time of low tide. The upper bound for the vertical eddy viscosity is estimated to be17 × 10⁻⁴m²s⁻¹ and for the vertical eddy diffusivity is estimated to be 15 × 10⁻⁴m²s⁻¹. The large vertical scale and the intensity of the observed boundary layer suggest that the flow in Monterey Canyon may be important for the renewal and circulation of water over the continental shelf in the bay area.     

Swelling of a lava plug associated with a Vulcanian eruption at Sakurajima Volcano,Japan, as revealed by infrasound record: case study of the eruption on January 2, 2007

In order to clarify the time relation of the expansion of a gas pocket and failure of its overlying plug of lava during Vulcanian eruptions, infrasound records and video images of the Vulcanian eruption that occurred at Sakurajima volcano on January 2, 2007 were analyzed with respect to their origin times. Weak (≤3 Pa) and slowly increasing air pressure preceded the impulsive compression phase by 0.25–0.32 s, and a longer-period rarefaction phase of infrasound waves was recognized at all microphone stations. The velocity of the compression phase was assumed to be supersonic (ca. 400 m/s) up to 850 m above the crater bottom from other recent explosions. On the other hand, the propagation velocity of the preceding weak signal was regarded to be similar to the air sound velocity because the lack of impulsiveness is unlikely to be related to the main compression phase. Therefore, the estimated origin time of the main compression phase was delayed by 0.5–0.7 s from the preceding phase. The origin time of the preceding phase coincided with the onset of the isotropic expansion process of the pressurized gas pocket, which was obtained by the waveform inversion of the explosion earthquake. In contrast, the origin time of the main impulsive phase coincided with the time when the expansion rate reached its peak. This observation suggests that the volumetric increase of the gas pocket caused swelling of the surface of the crater bottom and its subsequent failure. When the expansion velocity exceeded a threshold level, the main impulsive compression phase radiated with a high velocity by the sudden releases of the pressurized gases. The volumetric change at the source was estimated to be 280–560 m³ from the preceding phase of the infrasound. This volume change indicates that the vertical displacement of the swelling ground was on the order of 1.0 m, assuming the radius of the lava plug was ca. 10 m.     

Test Of High-Resolution Seismic Reflection And Other Geophysical Techniques On The Boup Landslide In The Swiss Alps

Complementary geophysical surveys on large landslides help revealgeologic structures and processes, and thus can help devise mitigation strategies. The combined interpretation of these methods enhance the result of each data set interpretation and makes it possible to derive a geological model of the landslide.We chose a test site on the Boup landslide (Wallis, western SwissAlps) to test high-resolution seismic reflection surveyscombined with ground penetratingradar (GPR), electromagnetic (EM) and electrokinetic spontaneous potential (SP) measurements.The results of the high-resolution seismic surveys suggest thatthe sliding is within a gypsum layer at 50 m depth and not as previously believed along a deeper (70 m) gypsum-shale boundary, also mapped seismically. Inversion of electromagnetic profiles (EM-34) with constraints from seismic data provided a model cross-section of conductivity values of the landslide (20–25 mS/m) and of the surrounding stable ground (10–15 mS/m), and it helped outline their boundary at depth.The accurate surface location of the landslide limit could be detected withelectromagnetic measurements of shallower depths of investigation (EM-31). Positive PS anomalies revealed an upward flow of mineralised water interpreted to follow the lateral boundary of the Boup landslide on its east side. Limited success was obtained withGPR profiling. This method can be hampered by conductive shallow layers, and itssuccessful application on landslides is expected to be strongly site dependant.     

Influence of bedrock groundwater on streamflow characteristics in a volcanic catchment

Groundwater movements in volcanic mountains and their effects on streamflow discharge and representative elementary area (REA) have remained largely unclear. We surveyed the discharge and chemical composition of spring and stream water in two catchments: the Hontani river (NR) catchment (6.6 km²) and the Hosotani river (SR) catchment (4.0 km²) at the southern part of Daisen volcano, Japan. Daisen volcano is a young volcano (17 × 10³ years) at an early stage of erosion. Our study indicated that deep groundwater that moved through thick lava and pyroclastic flows and that could not be explained by shallow movements controlled by surface topography contributed dominantly to streamflow at larger catchment areas. At the NR catchment, the deep groundwater contribution clearly increased at a catchment boundary defined by an area of 3.0 km² and an elevation of 800 m. At the SR catchment, the contribution deep groundwater to the stream also increased suddenly at a boundary threshold of 2.0 and 700 m. Beyond these thresholds, the contributions of deep bedrock groundwater remained constant, indicating that the REA is between 2 and 3 km² at the observed area. These results indicate that the hydrological conditions of base flow were controlled mainly by the deep bedrock groundwater that moved through thick lava and pyroclastic flows in the undissected volcanic body of the upper part of the catchment. Our study demonstrates that deep and long groundwater movements via a deep bedrock layer including thick deposits of volcanic materials at the two catchments on Daisen volcano strongly determined streamflow discharge instead of the mixing of small‐scale hydrological conditions. Copyright © 2015 John Wiley & Sons, Ltd.     

Methane Occurrences in Aquifers Overlying the Barnett Shale Play with a Focus on Parker County,Texas

Clusters of elevated methane concentrations in aquifers overlying the Barnett Shale play have been the focus of recent national attention as they relate to impacts of hydraulic fracturing. The objective of this study was to assess the spatial extent of high dissolved methane previously observed on the western edge of the play (Parker County) and to evaluate its most likely source. A total of 509 well water samples from 12 counties (14,500 km²) were analyzed for methane, major ions, and carbon isotopes. Most samples were collected from the regional Trinity Aquifer and show only low levels of dissolved methane (85% of 457 unique locations <0.1 mg/L). Methane, when present is primarily thermogenic (δ¹³C 10th and 90th percentiles of ?57.54 and ?39.00‰ and C1/C2+C3 ratio 10th, 50th, and 90th percentiles of 5, 15, and 42). High methane concentrations (>20 mg/L) are limited to a few spatial clusters. The Parker County cluster area includes historical vertical oil and gas wells producing from relatively shallow formations and recent horizontal wells producing from the Barnett Shale (depth of ～1500 m). Lack of correlation with distance to Barnett Shale horizontal wells, with distance to conventional wells, and with well density suggests a natural origin of the dissolved methane. Known commercial very shallow gas accumulations (<200 m in places) and historical instances of water wells reaching gas pockets point to the underlying Strawn Group of Paleozoic age as the main natural source of the dissolved gas.     

Population dynamics of autotrophic picoplankton in relation to environmental factors in a productive lake

Spatial and seasonal fluctuations in autotrophic picoplankton (APP) abundance in a eutrophic, dimictic lake (Lake Aydat, France) were measured concurrently with a variety of environmental variables. Cell number ranged from 0.03 to 2.36×10⁶ cells·ml^–1 (highest concentrations were >5-fold higher than in oligotrophic lakes) and averaged 24 ± 7% of total picoplankton abundance (APP + heterotrophic bacteria). APP abundance (1) peaked in spring simultaneously with heterotrophic flagellate and ciliate densities, (2) decreased during the nitrogen-limited and summer stratification period, and (3) increased with fall turnover. In summer-autumn, the contribution of single-cell eukaryotic (up to 66%) and colonial prokaryotic (18%) forms to total abundance peaked in the bottom waters. Multivariate regression analyses suggest that >40% variance in APP number changes may be explained by ciliate abundance (at 0–4 m depth-range), heterotrophic flagellate number and oxygen concentration (5–9 m), and ciliate carbon biomass (10–14 m). The model accounting for changes in heterotrophic bacterial abundance (5–9 m) indicates chlorophylla concentration (r ²=58%) and ciliate abundance (r ²=34%) as dominant covariates. The data presented here suggest that micrograzers control APP abundance in Lake Aydat.