Pollution history and recovery of a boreal lake exposed to a heavy bleached pulping effluent load

We examined the effects of heavy pulp mill discharges on the Lake Lievestuoreenjärvi ecosystem and the later recovery of diatom and chironomid communities from age-dated short core samples. Beginning in 1927 the lake received a heavy effluent load from a sulphite pulp mill. Except for the recession during the Second World War and the temporary closure of the mill from 1967 to 1971, the industrial load, containing large quantities of nutrients, organic matter and toxic compounds, increased continuously. In the early 1980s, laboratory documents were falsified by the directors of the mill and the systematic illegal effluent overload led to a collapse of the whole lake ecosystem. In 1985, the outdated plant was finally closed down. Based on the assessment of chemical properties and biological remains of the sediment, we distinguished five developmental phases in the ecological state of the lake. In the pre-industrial phase, the pelagic and profundal benthic communities were dominated by species preferring ultraoligotrophic or oligotrophic lakes. Concomitant with the increasing discharge and deposition of chlorine compounds, resin acids, and mercury, as well as strong acidity and hypolimnetic and epilimnetic anoxia, the ecological status changed in a short period from excellent to bad. Finally, in the early 1960s, the majority of the lake was virtually dead and the aquatic life survived only in the uppermost littoral zone. Since 1985, a fast recovery in the water quality has led to a strong, but temporary eutrophy in pelagic communities. The main peak of eutrophication was caused by the invasion of a species new to the lake,Aulacoseira granulata var.angustissima. Later, the pelagic communities shifted towards oligotrophy, but the original, pre-industrial status has not been re-established. The profundal benthic communities have not achieved the pre-industrial structure, but at present indicate mesotrophy.     

利用SPOT和ETM+影像进行县级土地利用调查--以河北省万全县为例

本文以SPOT和ETM 影像作为主要信息源，经过对图像进行几何纠正、配准、增强和融合等一系列处理后，在ARCVIEW GIS基础上，成功完成了对河北省万全县土地资源利用情况的调查分析。同时还对应用遥感技术成功获取土地资源利用信息的方法和技术进行了探讨。     

Record of Late Pleistocene Glaciation and Deglaciation in the Southern Cascade Range. I. Petrological Evidence from Lacustrine Sediment in Upper Klamath Lake, Southern Oregon

Petrological and textural properties of lacustrine sediments from Upper Klamath Lake, Oregon, reflect changing input volumes of glacial flour and thus reveal a detailed glacial history for the southern Cascade Range between about 37 and 15 ka. Magnetic properties vary as a result of mixing different amounts of the highly magnetic, glacially generated detritus with less magnetic, more weathered detritus derived from unglaciated parts of the large catchment. Evidence that the magnetic properties record glacial flour input is based mainly on the strong correlation between bulk sediment particle size and parameters that measure the magnetite content and magnetic mineral freshness. High magnetization corresponds to relatively fine particle size and lower magnetization to coarser particle size. This relation is not found in the Buck Lake core in a nearby, unglaciated catchment. Angular silt-sized volcanic rock fragments containing unaltered magnetite dominate the magnetic fraction in the late Pleistocene sediments but are absent in younger, low magnetization sediments. The finer grained, highly magnetic sediments contain high proportions of planktic diatoms indicative of cold, oligotrophic limnic conditions. Sediment with lower magnetite content contains populations of diatoms indicative of warmer, eutrophic limnic conditions. During the latter part of oxygen isotope stage 3 (about 37–25 ka), the magnetic properties record millennial-scale variations in glacial-flour content. The input of glacial flour was uniformly high during the Last Glacial Maximum, between about 21 and 16 ka. At about 16 ka, magnetite input, both absolute and relative to hematite, decreased abruptly, reflecting a rapid decline in glacially derived detritus. The decrease in magnetite transport into the lake preceded declines in pollen from both grass and sagebrush. A more gradual decrease in heavy mineral content over this interval records sediment starvation with the growth of marshes at the margins of the lake and dilution of detrital material by biogenic silica and other organic matter.     

Long-term bed load transport rate based on aerial-photo and ground penetrating radar surveys of fan-delta growth, Coast Mountains, British Columbia

Sequential aerial photography, sonar bathymetry, ground-penetrating radar (GPR), and sediment sampling and analysis provide the basis for calculating the volumetric and mass rate of progradation of the delta of Fitzsimmons Creek, a steep, high-energy, debris-flow-dominated channel draining about 100 km² of the southern Coast Mountains of British Columbia. Fitzsimmons Creek is typical of small mountain rivers in the region. GPR imaging is used to define the pre-depositional morphology of the receiving basin, a technique that improves the accuracy of the volumetric survey. The 52-year record (1947–1999) of progradation yielded an average annual volumetric transport rate of 1.00±0.16×10⁴ m³ year⁻¹ for bed load, corresponding to a mass transport rate of 1.60±0.28×10⁴ Mg year⁻¹. Bed load yields are consistent with those obtained in hydrogeomorphically similar basins in the region and elsewhere. Decade-based annual rates, which vary from 0.64±0.11×10⁴ to 2.85±0.38×10⁴ Mg year⁻¹, provide poor estimates of the 52-year average. Indeed, the 52-year record may also not be long enough to fully integrate the significant fluctuations in the sediment efflux from Fitzsimmons Creek. The methodology proposed in this paper can be transferred to other comparable mountain environments worldwide.     

Sequence of faulting and deformation during the 2000 eruptions of the Usu Volcano, Hokkaido, Japan—interpretation and image analyses of aerial photographs

Significant faulting and deformation of the ground surface has been rarely known during volcanic eruptions. Usu Volcano, Hokkaido, Japan, is a unique example of deformation due to felsic magma intrusion. Usu Volcano has a history of such types of eruptions as phreatic, pumice eruption (Plinian type), pyroclastic flowing and lava doming since 1663. On March 31, 2000, phreatomagmatic to phreatic eruptions took place after 23 years of dormancy in the western piedmont, followed by explosions on the western flank of Usu Volcano. They were associated with significant deformation including faulting and uplift. The eruptions and deformation were continuing up to the end of May 2000.We identified the faulting using total nine sets of aerial photographs taken from before the eruption (March 31, 2000) to more than 1 year (April 27, 2001) after the end of the activity, and traced deformation processes through image processing using aerial photographs. We found that some of the new faults and the associated phreatic eruptions were related to old faults formed during the 1977–1981 eruptive episode.The image processing has revealed that the surface deformation is coincident with the area of faulting forming small grabens and the phreatic explosion vents. However, the faulting and main explosive eruptions did not take place in the highest uplift area, but along the margin. This suggests that the faulting and explosive activities were affected by small feeder channels diverging from the main magma body which caused the highest uplift.     

A parameterisation for the threshold shear velocity to initiate deflation of dry and wet sediment

A new parameterisation for the threshold shear velocity to initiate deflation of dry and wet particles is presented. It is based on the balance of moments acting on particles at the instant of particle motion. The model hence includes a term for the aerodynamic forces, including the drag force, the lift force and the aerodynamic-moment force, and a term for the interparticle forces. The effect of gravitation is incorporated in both terms. Rather than using an implicit function for the effect of the aerodynamic forces as reported earlier in literature, a constant aerodynamic coefficient was introduced. From consideration of the van der Waals force between two particles, it was further shown that the effect of the interparticle cohesion force between two dry particles on the deflation threshold should be inversely proportional to the particle diameter squared. The interparticle force was further extended to include wet bonding forces. The latter were considered as the sum of capillary forces and adhesive forces. A model that expresses the capillary force as a function of particle diameter squared and the inverse of capillary potential was deduced from consideration of the well-known model of Fisher and the Young–Laplace equation. The adhesive force was assumed to be equal to tensile strength, and a function which is proportional to particle diameter squared and the inverse of the potential due to adhesive forces was derived. By combining the capillary-force model and the adhesive force model, the interparticle force due to wet bonding was simplified and written as a function of particle diameter squared and the inverse of matric potential. The latter was loglinearly related to the gravimetric moisture content, a relationship that is valid in the low-moisture content range that is important in the light of deflation of sediment by wind. By introducing a correction to force the relationship to converge to zero moisture content at oven dryness, the matric potential–moisture content relationship contained only one unknown model parameter, viz. moisture content at −1.5 MPa. Working out the model led to a rather simple parameterisation containing only three coefficients. Two parameters were incorporated in the term that applies to dry sediment and were determined by using experimental data as reported by Iversen and White [Sedimentology 29 (1982) 111]. The third parameter for the wet-sediment part of the model was determined from wind-tunnel experiments on prewetted sand and sandy loam aggregates. The model was validated using data from wind-tunnel experiments on the same but dry sediment, and on data obtained from simulations with the model of Chepil [Soil Sci. Soc. Am. Proc. 20 (1956) 288]. The experiments showed that soil aggregates should be treated as individual particles with a density equal to their bulk density. Furthermore, it was shown that the surface had to dry to a moisture content of about 75% of the moisture content at −1.5 MPa before deflation became sustained. The threshold shear velocities simulated with our model were found to be in good agreement with own observations and with simulations using Chepil's model.     

Potential toxicity of resuspended particulate matter and sediments: Environmental samples from the Bay of Banyuls-sur-Mer and Thermaikos Gulf

This contribution attempts to determine whether, and to what degree, environmental samples of resuspended particulate matter and sediments exert a toxicological impact. Further, an attempt is made also to screen the toxic level of potentially hazardous sites, based upon established sediment toxicity criteria. Therefore, a rapid, cost-effective and highly sensitive biotest (bioluminescence assay, based upon marine bacteria) has been applied on: biological fluid extracts; bottom sediments; and sediment trap samples. Samples were taken either from the Bay of Banyuls-sur-Mer (northwestern Mediterranean, France) or Thermaikos Gulf (northeastern Mediterranean, Greece). Biological fluid extracts and sediment trap samples corresponded to periods of resuspension events, or preceding and following such events.The results have revealed that the sampling strategy and biotest implemented in this study might be a useful tool for screening the toxicity of resuspended matter and sediments. Resuspension events appear to be able to exert an influence on the chemical forms of the micro-pollutants; thus on their bioavailability and toxicity. Nevertheless, based upon chemical analysis combined with the bioassay, the toxic level of the sediment samples could be determined, with the level at potentially hazardous sites being characterised.     

Image processing of 2D resistivity data for imaging faults

A methodology to locate automatically limits or boundaries between different geological bodies in 2D electrical tomography is proposed, using a crest line extraction process in gradient images. This method is applied on several synthetic models and on field data set acquired on three experimental sites during the European project PALEOSIS where trenches were dug. The results presented in this work are valid for electrical tomographies data collected with a Wenner-alpha array and computed with an l₁ norm (blocky inversion) as optimization method. For the synthetic cases, three geometric contexts are modelled: a vertical and a dipping fault juxtaposing two different geological formations and a step-like structure. A superficial layer can cover each geological structure. In these three situations, the method locates the synthetic faults and layer boundaries, and determines fault displacement but with several limitations. The estimated fault positions correlate exactly with the synthetic ones if a conductive (or no superficial) layer overlies the studied structure. When a resistive layer with a thickness of 6 m covers the model, faults are positioned with a maximum error of 1 m. Moreover, when a resistive and/or a thick top layer is present, the resolution significantly decreases for the fault displacement estimation (error up to 150%). The tests with the synthetic models for surveys using the Wenner-alpha array indicate that the proposed methodology is best suited to vertical and horizontal contacts. Application of the methodology to real data sets shows that a lateral resistivity contrast of 1:5–1:10 leads to exact faults location. A fault contact with a resistivity contrast of 1:0.75 and overlaid by a resistive layer with a thickness of 1 m gives an error location ranging from 1 to 3 m. Moreover, no result is obtained for a contact with very low contrasts (1:0.85) overlaid by a resistive soil. The method shows poor results when vertical gradients are greater than horizontal ones. This kind of image processing technique should be systematically used for improving the objectiveness of tomography interpretation when looking for limits between geological objects.     

Distribution of trace metals in the Odra River system: Water–suspended matter–sediments

Extensive investigations of trace metals concentrations in water, suspended particulate matter (SPM) and bottom sediments of the whole Odra River system were carried out over the years 1997–2000. The vertical distribution of selected metals and their mobility were also studied in the sediment cores from upper and middle river sections. Significant levels of metal contamination were found. Median concentrations (Cd, Pb, Cu, Zn and As) in the SPM and sediments were (mg kg⁻¹) 7.1 and 8.9 Cd, 128 and 146 Pb, 81 and 119 Cu, 1198 and 1204 Zn, 48 and 54 As, respectively. The highest metal pollution of the Odra River solids was found with cadmium, zinc, lead and arsenic, showing high similarity in their frequency distributions in both SPM and sediments. Cd, Zn and As appear to be of particular concern because of the high levels, that appear to be bioavailable, and their high mobility. The exchangeable and carbonate chemical forms of Cd and Zn reached up to 50% of their total amount. Besides the determination of total metal concentration, the metal chemical forms in river solids were investigated. The results of very wide studies of the Odra River system through 4 years suggest that metal pollution decreased, especially for Zn, Pb and Cu. Among all metals studied in the Odra River sediments, substantial reductions of Cd contamination were observed neither in the period after ’97 flood, nor if compared with the earlier results obtained before ’97. No essential differences of the metal contents were observed among the samples for the same river compartment, from the same locality, taken within the five sampling campaigns. The pattern of spatial and vertical metal distributions in the river solids indicates that a variety of sources might be responsible for the contamination; very intensive, historical and current mining and smelting activities probably are the most important ones.     

Data assimilation (4D-VAR) to forecast flood in shallow-waters with sediment erosion

In this paper, the four-dimensional variational data assimilation technique (4D-VAR) is presented as a tool to forecast floods. Our study is limited to purely hydrological flows and supposes that the weather, here a big rain, has been already forecasted by meteorological services. The technique consists in minimizing, in the sense of Lagrange, the cost function: a measure of the difference between calculated data and available observations, here the water level. This is done under constraints that are the equations of the physical model. In our case, we modified the shallow-water equations to include a simplified sediment transport model. The steepest descent algorithm is then used to find the minimum. This is made possible because we can compute analytically the gradient of the cost function by using the adjoint equations of the model. As an application of the 4D-VAR technique, the overflowing of the Chicoutimi River at the Chute-Garneau dam, during the 1996 flood, is investigated. It is found that the 4D-VAR method reduces the error in the water height forecast even when the erosion model is not activated. In terms of Lyapunov exponents, we estimate the predictability horizon of such an event to be about half-an-hour after a big rain. However, this limit of predictability can be increased by using more observations or by using a finer computational grid.