Limnological investigation of biogenic silica sedimentation and silica biogeochemistry in Lake St. Moritz and Lake Zürich

Biogenic silica concentrations were determined from core samples of laminated sediments collected from Lake Zürich and Lake St. Moritz and used to calculate rates of biogenic silica sedimentation. In Lake Zürich biogenic silica sedimentation increased from 65 g SiO₂·^?2·yr^?1 in 1894 to nearly 900 g SiO₂·m^?2·yr^?1 in 1896 and in Lake St. Moritz biogenic silica sedimentation increased about 6-fold during the 1930s. Both periods of increased biogenic silica flux followed increased loadings of domestic sewage to the lakes. In Lake Zürich the period of increased flux only lasted for a few years whereas in Lake St. Moritz high fluxes were maintained to the top of the core. In Lake Zürich increased production and sedimentation of diatoms could be maintained only until the silca reservoir in the water mass was depleted whereas in Lake St. Moritz (0.10 year residence time and a mean depth of 25 m) diatom production and sedimentation was maintained at a high level by rapid replenishment of silica from tributary inputs. Although historical patterns of biogenic silica sedimentation differed for the two lakes, it is hypothesized that increased biogenic silica deposition in both lakes occurred because diatom production was stimulated by phosphorus enrichment of the water mass and that the silica biogeochemistry of both lakes was affected as a result.     

Holocene slip rate and paleoearthquake records of the Salt Lake segment of the Northern Zhongtiaoshan Fault,Shanxi Province

The Northern Zhongtiaoshan Fault is a major deep fault at the southern margin of the Yuncheng Basin. There have been few studies on the fault, and the historical earthquakes are few and weak. However, the intensity of activity on the fault should never be underestimated. Through interpretations of aerial images, topography measurements and excavation of trenches, this paper studied the fault distribution, the surface deformation and the activity of the normal fault south of Salt Lake near the city of Yuncheng. By tracing faults in the three trenches, it was found that there had been at least three large paleoseismic events, at 1–3.5, 3.6–4.4 and 7.4–8.8 ka BP. Employing 14 C dating, we determined the same gravel layers in the uplifted side and downthrown side. Making differential Global Positioning System measurements of the vertical difference and topographic profile, we obtained the mean slip rate of the Northern Zhongtiaoshan Fault since 24.7 ka BP(0.75±0.05 mm/a). Using the results of relevant studies, we calculated the possible vertical fault displacement of one earthquake(2.35 m) and obtained the recurrence interval of characteristic earthquakes as 2940–3360 a after dividing the displacement by the mean slip rate.     

A remote controlled freeze corer for sampling unconsolidated surface sediments

A new coring device is presented which allows the recovery of loose watery surface sediments and the water/sediment interface byin situ freezing, resulting in well preserved samples. The instrument consists of a tripod with adjustable legs, a hydraulic system, an insulated thermos (with two electrical pumps), and a wedge-shaped freeze box. Alcohol chilled with dry ice is the freezing agent. The corer can be deployed from a boat or a raft and permits recovery of undisturbed surface sediment sections up to 100 cm long in deep waters. The corer is particularly useful for sampling varved sediments.     

Reconnaissance surveys of near-event seismic activity in the volcanoes of the Cascade Range,Oregon

Surveys of near-event seismic activity were made at two principal locations in the Cascade Range in Oregon during the summers of 1969 and 1970. A tripartite array of ultrasensitive high frequency seismometers was deployed about 7.5 km north of the Mt. Hood summit with one of the 1-km legs oriented broadside to the dormant volcano. Seismometers were emplaced over olivine andesite flows associated with the Pinnacle, one of the parasitic cones formed on the flanks of the strato-volcano. During 16 days of operation on the north slope, 53 near events were recorded, most of which originated within the upper crust and were associated with the north-south trending zone of the Cascade Range. Event magnitudes for these near events range from ?1.7 to +1.8 and determination of b-values in the Gutenberg-Richter relationship was ?0.80, indicating a probable tectonic mechanism for the shocks. During the late summer of 1970, a four-station array was operated at Crater Lake Park about 13 km south-southwest of the caldera rim. In addition, an ultraportable outlier station was operated at two locations north of the caldera that resulted from the collapse of ancient Mt. Mazama some 6,600 years ago. Only a limited number of near events with S-P intervals of 4 sec or less were detected at Crater Lake; a larger number were recorded with S-P intervals longer than 4 sec. Event epicenters for the Crater Lake area are broadly distributed in azimuth, indicating the complex structure underlying the Cascade Range in southern Oregon. Crater Lake is located astride the broad upwarp of crystalline pre-Cenozoic rocks believed to extend northeast from the Klamath Mountains to the Ochoco Uplift of central Oregon. Major regional structural trends are also shown by the north-south trending belts of the Cascade volcanoes, probably related to deep fracture zones, and by the northeast-trending shear zones that exist in the Basin and Range province to the southeast of Crater Lake. Regional gravity and aeromagnetic surveys indicate that Crater Lake lies at the intersection of these zones that probably provided the conduits for the rise of magma that ultimately led to the collapse of Mt. Mazma and the formation of Crater Lake. Epicenters for near events recorded at this juncture do not reflect these linear trends and, indeed, a generally smaller incidence of near-event activity was recorded at Crater Lake than was recorded at Mt. Hood. Magnitudes for Crater Lake events with S-P intervals of 4 sec or more range from +0.25 to +2.19, and an examination of the relationship between cumulative frequency and magnitude for these events yields a b-value in the Gutenberg-Richter relationship of ?1.16, indicating the events at Crater Lake, like those detected at Mt. Hood, are associated with tectonic rather than volcanic sources. Events for which depth determinations were made show these sources to be within the crust, occurring in the upper 10 km of the earth’s crust. The relatively low incidence of small magnitude near events within the Oregon Cascade Range shows the aseismicity of the mountain chain which is consistent with the low incidence of earthquakes of a magnitude of 4.5 or greater detected for the volcanic range. The volcanoes of the Cascade Range in Oregon are dormant, and only small numbers of shocks are now being generated, probably from isostatic adjustments within the crust. The Cascade volcanic range, which once was a seismically active island are chain associated with subduction zones off the northwestern coast of America, has moved into a passive phase in which most seismic activity in western Oregon now occurs along the ridge and fracture zones offshore and within the Willamette Downwarp west of the dormant chain.     

Origin and size of hypolimnic mixing in Urnersee,the southern basin of Vierwaldstättersee (Lake Lucerne)

Urnersee and Gersauersee are two adjacent basins of Vierwaldstättersee (Lake Lucerne, Switzerland), seperated by a sill of 85 m depth, with similar topography (max. depth 195 and 213 m, respectively) but remarkably different exposure to “external forces”, such as wind and river input. Urnersee is exposed to diurnal winds and to occasional strong storms from the south (Föhn) whereas the wind over Gersauersee is moderate or weak. Two rivers, both having very large discharges during storms, replace the total water volume of Urnersee about once a year; in contrast, no large river flows directly into Gersauersee. Between March and October 1986, meteorological parameters, water temperatures and currents were measured quasi-continuously with the aim to quantify hypolimnic water exchange and mixing in Urnersee and to asses the relative importance of wind mixing versus river-induced water exchange for the renewal of the deep water layers. Three periods could be identified: (1) in April, weak stratification and strong episodic storms exchange about 50% of the deep hypolimnion (DH, defined as layer below 110 m depth) leading to a mean heat flux of 36 Wm². Because of the large wind mixing the water of the exposed Urnersee below about 20 m depth becomes lighter than in the sheltered Gersauersee. (2) In May and June, the horizontal density gradient causes about 65% renewal of the Urnersee DH by the heavier Gersauersee intermediate water but does not affect the heart content. (3) Simultaneously with these processes are the episodic river floods adding another 20% to the DH water exchange and causing a heat flux of about 6 Wm². During the rest of the summer, water exchange remains below 10% and is mainly due to episodic flood while wind mixing has little influence. Yet, during floods water input into the DH per unit time can still be very large and heat, fluxes reach 600 Wm² or more. The influence of lateral density currents between the two adjacent basins on hypolimnic mixing is of great ecological significance and explains the oxygen saturation found in the deep water of Urnersee compared to Gersauersee.     

Monitoring the water balance of Lake Victoria, East Africa, from space

Using satellite gravimetric and altimetric data, we examine trends in water storage and lake levels of multiple lakes in the Great Rift Valley region of East Africa for the years 2003–2008. GRACE total water storage estimates reveal that water storage declined in much of East Africa, by as much as , while altimetric data show that lake levels in some large lakes dropped by as much as 1–2 m. The largest declines occurred in Lake Victoria, the Earth’s second largest freshwater body. Because the discharge from the outlet of Lake Victoria is used to generate hydroelectric power, the role of human management in the lake’s decline has been questioned. By comparing catchment water storage trends to lake level trends, we confirm that climatic forcing explains only about 50decline. This analysis provides an independent means of assessing the relative impacts of climate and human management on the water balance of Lake Victoria that does not depend on observations of dam discharge, which may not be publically available. In the second part of the study, the individual components of the lake water balance are estimated. Satellite estimates of changes in lake level, precipitation, and evaporation are used with observed lake discharge to develop a parameterization for estimating subsurface inflows due to changes in groundwater storage estimated from satellite gravimetry. At seasonal timescales, this approach provides closure to Lake Victoria’s water balance to within . The third part of this study uses the water balance of a downstream water body, Lake Kyoga, to estimate the outflow from Lake Victoria remotely. Because Lake Kyoga is roughly 20 times smaller in area than Lake Victoria, its water balance is strongly influenced by inflow from Lake Victoria. Lake Kyoga has been shown to act as a linear reservoir, where its outflow is proportional to the height of the lake. This model can be used with satellite altimetric lake levels to estimate a time series of Lake Victoria discharge with an rms error of about .     

湖泊中溶解氧极大值之成因

本文对位于台湾南部两个次高海拔湖泊－大鬼湖及万山神池进行研究,并试图探讨此二无明显进,出水口之封闭型湖泊中,水体溶解氧垂直分布出现极大值之成因,此二次高海拔湖泊地处偏远,人烟罕至,因此较不受人为干扰,为研究湖泊水体中种种自然作用的良好对象,大鬼湖平均水深约14.8m,最深处约40m,降冬季外,水体均有分层现象,1988年夏季资料显示,水深16m以下水体趋于无氧状态,且于此深度以上的溶氧饱和值,均接近当地的大气饱和值（约78％）,经各种资料推断,此极大值的成因,除水团乃于春季时留下主溶氧值之外,应综合了季节增温效应下使表水向下温合的物理作用（尤其是山区明显的口,夜温差所引起的表水冷却向下混合作用,且其混合深度随季节增温而逐步变浅）,消耗溶解氧的生物作用及跃层存在等影响因素,而非单纯的物理或生物作用所造成,1991年4月万山神池观测资料显示,湖水平均深度约8m,最深可达14m ,其氧饱和程序分布在80％－104％之间,表水接近近当地饱和值（78％）,而于1.5m处往下增加,于2m处有溶解氧及和程度极大值,比当时大气饱合值高约20％,此极大值深度与叶绿素a极大值深度相吻合,主要由生物之光合作用造成,此外,于此深度pH值亦有明显增加现象,更证实了生物作用的存在,此二湖环境背景相似,且水泥中均出现溶氧极大值,但二者间极大值却有不同的成因,经曲两不同成因的比较,将可提供许多相关研究资料。     

Design and testing of a global climate prediction system based on a coupled climate model

A global climate prediction system(PCCSM4) was developed based on the Community Climate System Model, version 4.0, developed by the National Center for Atmospheric Research(NCAR), and an initialization scheme was designed by our group. Thirty-year(1981–2010) one-month-lead retrospective summer climate ensemble predictions were carried out and analyzed. The results showed that PCCSM4 can efficiently capture the main characteristics of JJA mean sea surface temperature(SST), sea level pressure(SLP), and precipitation. The prediction skill for SST is high, especially over the central and eastern Pacific where the influence of El Ni?o-Southern Oscillation(ENSO) is dominant. Temporal correlation coefficients between the predicted Ni?o3.4 index and observed Ni?o3.4 index over the 30 years reach 0.7, exceeding the 99% statistical significance level. The prediction of 500-hPa geopotential height, 850-hPa zonal wind and SLP shows greater skill than for precipitation. Overall, the predictability in PCCSM4 is much higher in the tropics than in global terms, or over East Asia. Furthermore, PCCSM4 can simulate the summer climate in typical ENSO years and the interannual variability of the Asian summer monsoon well. These preliminary results suggest that PCCSM4 can be applied to real-time prediction after further testing and improvement.     

A numerical-physical planetary boundary layer model

A numerical-physical model for the planetary boundary layer has been formulated for the purpose of predicting the winds, temperatures and humidities in the lowest 1600 m of the atmosphere. An application of the model to the synoptic situation of 30 August, 1972, demonstrates its ability to produce useful forecasts for a period of 24 h. Results are illustrated in terms of horizontal maps and time-height sections of winds and temperatures. The model is divided in the vertical direction into three layers that are governed, respectively, by different physical formulations. At the lowest level, which is the surface of the earth, forecasts of temperature and humidity are computed from empirical relations. In the first layer, the surface layer, application is made of the similarity theories of Monin-Obukhov, Monin-Kazanski and Businger’s form of the universal functions. The second layer, the Ekman layer, is 1550 m deep and is governed by diagnostic momentum and time-dependent thermodynamic and humidity equations. External input to the model are large-scale pressure gradients and middle-level cloudiness. Cressman’s objective analysis procedure is applied to conventional surface and upper air data over a horizontal region of about 2500 km by 2500 km, centered about Lake Ontario. With a grid distance of 127 km and a time interval of 30 min, the computer time required on Control Data Cyber 76 for a 24 h forecast for the case study is less than two minutes.     

Float‐transient electromagnetic method: in‐loop transient electromagnetic measurements on Lake Holzmaar,Germany

Lake sediments may serve as archives on paleoclimatic fluctuations, geomagnetic field variations and volcanic activities. Lake Holzmaar in Eifel/Germany is a maar lake and its lacustrine sediments provide paleoclimatic proxy data. Therefore, knowledge about the geometry and, especially, about the thickness of the sediments is very important for determining an optimum drilling location for paleoclimatic studies. We have developed a floating in‐loop transient electromagnetic method field set up (Float‐transient electromagnetic method) with a transmitter and receiver size of 18 × 18 m² and 6 × 6 m² respectively. This special set up enables in‐loop transient electromagnetic method measurements on the surface of freshwater lakes that define the geometry and the thickness of sediments beneath such lakes thus helping to determine optimum drilling locations. Due to the modular design of the new Float‐transient electromagnetic method field set up, this system can be handled by two operators and can easily be transported. Sixteen in‐loop soundings were carried out on the surface of Lake Holzmaar. The transient electromagnetic method data could not be interpreted by conventional 1D inversions because of the 3D distribution of subsurface conductivity caused by the lake's geometry. Three‐dimensional finite element modelling was applied to explain the observed transients and the 3D conductivity distribution beneath the lake was recovered by taking its geometry into account. The 3D interpretation revealed approximately 55 m thick sediments beneath 20 m deep water in the central part of the lake.     

内蒙古察干淖尔盐湖29-8ka时段的湖面波动

以内蒙古锡林郭勒盟苏尼特右旗的察干淖尔盐湖为研究对象,利用OSL(Optically Stimulated Luminescence)测年技术和DEM(Digital Elevation Model)数字高程模型,重建湖面波动历史,探讨湖泊形成与环境变化过程.通过对察干淖尔盐湖周边大量的野外考察,发现湖泊周围存在海拔高程为1020、978和973 m的三级古湖岸阶地,其OSL测年结果分别为29.2±1.3、18.4±0.8及8.2 8.0 ka.通过湖岸阶地高程恢复的上述3个时期的古湖面积分别为3600、500和400 km~2.与现今的干旱盐湖景观迥然不同.     

Geothermal energy

The potential of a geothermal area is primarily dependent on volume and temperature of the reservoir and adequacy of fluid supply. Inadequate fluid supply may be a more common limiting factor than inadequate heat supply, for heat stored in the upper 10,000 ft of many hot spring systems is 1,000 to 10,000 times their annual natural heat flow. Except in very porous reservoirs, most of this heat is stored in rocks rather than in pore fluids. Geothermal fields can be classified as hot spring systems or as deep insulated reservoirs with little surface expression; gradations also exist. Hot spring systems have high near-surface permeability, at least locally on faults and fractures, permitting fluids to escape at high rates. Owing to vigorous circulation and escaping fluids and heat, near-surface temperatures are high, but temperatures deep in the system are lower than would prevail with inhibited escape. Deep reservoirs with little surface expression require permeable reservoir rocks capped by insulating rocks of low permeability. Larderello, Italy, and Salton Sea, California, have slight leakage, but others may have no leakage. Liquid water, which can be at temperatures far above 100° C because of existing pressures, is generally the dominant fluid. Steam can form by boiling as hot water rises to levels of lower pressure. However, in several explored systems the heat supply is so high and rate of discharge of water so low that steam exists even deep in the system. Dry steam areas are probably rare. About 30 areas in the United States have been explored for geothermal energy, but dry steam has been proved only at « The Geysers ». Extensive utilisation of geothermal energy must therefore depend largely upon steam « flashed » from hot water with decrease in pressure. Problems that confront broad utilisation of geothermal energy include: 1) discovery of reservoirs with adequate supply of energy and natural fluids; 2) deposition of CaCO; or SiO₂; 3) chemical corrosion; 4) objectionable chemicals in some effluents; and 5) inapplicability of existing public laws. The optimum environment for a geothermal reservoir includes:

1. 1.
   Potent source of heat, such as a magma chamber. A depth of at least two miles provides enough pressure to insure water of high temperature; 5 miles may be too deep for effective transfer of heat to circulating water. Such heat sources are most likely to occur in regions of late Cenozoic volcanism.     
   
   

Interpretation of metal profiles in a sediment core from lake geneva: metal mobility or pollution

The transition elements Mn, Fe, Co, Zn, Cu and Cr together with Pb and P were determined in three sequential extractions of sediment and in the pore water of a 35-cm sediment core from Lake Geneva. A high percentage of the mobile Mn, Fe, P, Co, Pb and Zn was extracted in the first (inorganic) fraction, whereas the mobile Cu and Cr were mostly found in the second (organic and sulfide) fraction. Mn was subject to recycling with solubilization from the surface sediment to the oxic bottom water where it was precipitated. Highly enriched Co and P values were observed in the top few mm of the sediment resulting from post depositional solubilization with upward migration and precipitation with oxi-hydroxides at the sediment surface. The slower increase of Cu, Zn and Pb concentrations towards the sediment/water interface are the result of an increase of anthropogenic sources over the past decades. A further increase in Pb in the first cm of the core is believed to be related to the freshly precipitated Mn oxi-hydroxides.     

Generalization of the results of long-term strainmeter and GPS observations for intraplate regions

The measurements of contemporary deformations in Asia in the region of active seismic processes stretching from Northern Tien Shan through Altai and Sayan Mountains to Lake Baikal have been carried out for a few decades by strainmeters and tiltmeters. Recently, these works were expanded by the network measurements by the methods of space geodesy. Studying the contemporary displacements and deformations of the Earth’s crust is necessary for solving various problems on a wide range of spatial and time scales from exploring the contemporary deformations of the tectonic plates and present-day seismicity to estimating the man-made loads on the geological medium. The measurements of the deformations are conducted by different methods with short and long baselines, on the surface and at different depths. Based on the obtained data, it is possible to analyze the strain rates at different frequencies and to compare them with the level of seismic activity in the region. In recent decades, the seismic zone located in the south of the Siberian platforms experienced a few strong (M = 6.3–7.5) earthquakes, during which the deformation and displacement fields were captured. South of this region, the displacement fields during the earthquakes with the magnitude above 8 were recorded by the space geodetic methods in China. By excluding the coseismic effects, one can estimate the tectonic component of the displacements and deformations.     

A comparison of the recent eruptions of Mt. Pelée,Martinique and Soufrière,St. Vincent

The simultaneous eruption of Mt. Pelée, Martinique and Soufrière, St. Vincent are regarded as the first recognized examples of Pelean-type and St. Vincent-type pyroclastic eruptions. Both produced nuées ardentes, the former usually laterally directed because of the presence of a dome and the latter vertically directed from an open crater. Both volcanoes have subsequently erupted for a second time this century. The 1902–05 and 1929–32 eruptions of Mt. Pelée produced andesite lava of almost identical composition and mineralogy. Both contain two generations of plagioclase, orthopyroxene, Fe-Ti oxide, corroded brown amphibole and olivine rimmed by pyroxene. In contrast, the Soufrière material is more basic in composition varying from basaltic andesite to basalt in 1902–03 and basaltic andesite in 1971–72. The Soufrière material contains two generations of plagioclase (with those of 1971–72 having additional zones of labradorite), clinopyroxene, orthopyroxene, olivine and Fe-Ti oxide. The pyroclastic deposits are strikingly different, those from the Pelean-type eruption are termed «block and ash deposits» being characterised by poorly vesicular lava blocks up to 7 m in diameter, while the St. Vincent-type eruption produced «scoria and ash deposits» containing vesicular ropey blocks or bombs no larger than 1 m in diameter. The differences in styles of eruption are attributed to differences in viscosity and mechanism of eruption of the magmas. Stratigraphic studies of Mt. Pelée reveal that the volcano has produced basaltic andesite scoria and ash deposits from St. Vincent-type eruptions. It is concluded that the recent eruptions of Pelée tapped a deep level magma during both eruptions releasing magma of similar composition, while the 1971 Soufrière magma is thought to be a remnant of the 1903 basaltic magma which remained at a high level within the volcano where it underwent enrichment in plagioclase and loss of olivine and oxide.     

Michigan basin regional ground water flow discharge to three Great Lakes

Ground water discharge to the Great Lakes around the Lower Peninsula of Michigan is primarily from recharge in riparian basins and proximal upland areas that are especially important to the northern half of the Lake Michigan shoreline. A steady-state finite-difference model was developed to simulate ground water flow in four regional aquifers in Michigan's Lower Peninsula: the Glaciofluvial, Saginaw, Parma-Bayport, and Marshall aquifers interlayered with the Till/"red beds," Saginaw, and Michigan confining units, respectively. The model domain was laterally bound by a continuous specified-head boundary, formed from lakes Michigan, Huron, St. Clair, and Erie, with the St. Clair and Detroit River connecting channels. The model was developed to quantify regional ground water flow in the aquifer systems using independently determined recharge estimates. According to the flow model, local stream stages and discharges account for 95% of the overall model water budget; only 50% enters the lakes directly from the ground water system. Direct ground water discharge to the Great Lakes' shorelines was calculated at 36 m3/sec, accounting for 5% of the overall model water budget. Lowland areas contribute far less ground water discharge to the Great Lakes than upland areas. The model indicates that Saginaw Bay receives only approximately 1.13 m3/sec ground water; the southern half of the Lake Michigan shoreline receives only approximately 2.83 m3/sec. In contrast, the northern half of the Lake Michigan shoreline receives more than 17 m3/sec from upland areas.     

Major-element petrochemistry of some extrusive rocks from the volcanically active Mariana Islands

Volcanic rocks from six of the currently or recently active volcances of the Mariana Island are show little variation in major element abundances. SiO₂ content averages 51.5 wt.%. The flows are high in Al₂O (mean 17.7 wt.%) and Fe oxides (mean 10.1 wt.% calculated as FeO only), and moderate in MgO content (mean 4.7 wt.%), Na₂O (mean 2.7 wt.%), and K₂O (mean 0.7 wt.%). Only the rocks from Farallon de Pajaros, the northernmost of the Mariana Islands, deviate slightly from the average of the analyses. Three analyses from this island are slightly higher in SiO₂ (about 54 wt.%) and Al₂O₃, and are lower in total Fe oxides and MgO. According to preferred classification, the lavas of the Mariana Islands can be termed mela-andesites, high-alumina basalts, or calc-alkaline (orogenic) basalts. The K₂O values (mean 0.7 wt.%) obtained from lavas of the Mariana Islands are significantly higher than the K₂O values (about 0.33 wt.%) from volcanics of the Izu chain to the north. Inasmuch as the substantial scatter in location of earthquake foci beneath both arcs prevents accurate delineation of the upper boundary of the Benioff zone, it presently cannot be determined whether this discrepancy in K₂O values reflects a difference in depth from the volcanic are to the dipping seismic zone or relates to other phenomena. The older volcanic islands within the Mariana-Bonin island chain apparently defined an island arc system during Eocene to Miocene time. This indicates that the present plane of convergence between the Pacific plate and the Philippine Sea plate has defined the convergence between these plates since Eocene time.     

FEM × DEM modelling of cohesive granular materials: Numerical homogenisation and multi-scale simulations

The article presents a multi-scale modelling approach of cohesive granular materials, its numerical implementation and its results. At microscopic level, Discrete Element Method (DEM) is used to model dense grains packing. At the macroscopic level, the numerical solution is obtained by a Finite Element Method (FEM). In order to bridge the micro- and macro-scales, the concept of Representative Elementary Volume (REV) is applied, in which the average REV stress and the consistent tangent operators are obtained in each macroscopic integration point as the results of DEM’s simulation. In this way, the numerical constitutive law is determined through the detailed modelling of the microstructure, taking into account the nature of granular materials. We first elaborate the principle of the computation homogenisation (FEM × DEM), then demonstrate the features of our multiscale computation in terms of a biaxial compression test. Macroscopic strain location is observed and discussed.     

Petrochemistry of quaternary pumiceous pyroclastic products in southern Guadeloupe (F.W.I.)

Pumiceous pyroclastic products are present as flows and falls at several stages of the evolution of the southern Guadeloupe volcanic island. An understanding of this volcanism had to rely on detailed petrochemical data of these products to complement similar data for effusive rocks so as to yield complete stratigraphical coverage. On the other hand most pumiceous rocks are more or less conspicuousily banded suggesting that mixing phenomena occurred to different degrees in their genesis. Three major classes of pumiceous products are found: (1) the Axial Chain deposits (2.0–0.5 My) are characterized by An90, 75-55 + En55 + Wo42En37 + Usp35-37 ± Fo68 ± Hble, SiO₂ 60%, SiO₂/Th 35.6, and La/Th 3.9. Banded samples have components that differ in evolution indices by about 50%; (2) the Bouillante Chaine pyroclastics (0.3–0.1 My) consist of scattered deposits with variable mineralogical and geochemical compositions that seem to have erupted from a number of small eruptive centers. Qz-dacitic pumice is common with An90, 70-45 + En66-56 + Usp32-38 + Ilm94 + Hble ± Wo40-42En40-42, SiO₂ 62%, SiO₂/Th 22.9, and La/Th >4. Mixed pumice samples have highly contrasted evolution indices differing by up to 120%; (3) the Pintade pumice flows and falls correspond to the major pyroclastic event (approx. 10 km³) in the southern Basse Terre area. They are characterized by An85-70 + En66-56 + Usp32-37 ± Wo42En42, SiO₂/Th 18.7-22?6, and La/Th 3.0-4.0. Banded pumice lumps are scarce and show slight compositional contrasts; differences in evolution indices do not exceed 38%. Axial and Bouillante chain pyroclastics and Grande Découverte volcano pumice respectively, form two different families in trace element plots. Minor elements in pyroxenes also are distinctive. These trends are similar to those obtained for effusive rocks and define comagmatic series. Major and trace element data for the separated components of inhomogeneous pumice in each formation always plot in the corresponding series. These chemical discriminants can be used to attribute samples of unknown provenance to a given volcanic ensemble. An inverse relationship between differences in evolution indices in inhomegeneous pumice and the volume of any single eruptive sequence is noted. This is an indication that pumiceous pyroclastic rocks were erupted from a zoned magma chamber. We favor an interpretation where zonation is produced by influx of less envolved magma in superficial differentiated chambers which is a direct cause for eruption.     

五里湖富营养化过程中水生生物及生态环境的演变

五里湖是太湖西北部一个小型浅水湖湾,是无锡市的饮用水源和主要风景游览区。50年代时,该湖基本保持着原始状态,全湖以大型水生植物占优势。湖水清澈见底,水质为中营养水平,溶氧接近饱和,对外来的N、P污染冲击具有很强的缓冲能力。底泥的氧化程度较高,磷和有机物含量仅为0.023％和0.75％。浮游藻类受到了大型水生植物的强烈抑制,年均数量为26.7×10~4个/L,以硅藻和隐藻为主;从春季至秋季,随着大型水生植物的增长,浮游藻类数量大幅度减少。浮游动物多达190种,年均数量为5660ind./L。大型底栖动物较多,以日本沼虾和螺、蚌类为主。鱼类资源十分丰富,63种鱼中以凶猛性鱼类占优势,并有较多的底栖性鱼类。 自50年代以来,大约有1/2的湖面被围垦,沿岸带生态条件被破坏,失去了最适合于大型水生植物生长的浅水区。加之60年代后期在全湖放养草鱼,水生植被遭到彻底毁灭。外源污染加剧,引起了水质的严重富营养化。围垦和修建水闸隔断了五里湖与太湖间的通道,限制了两个水体间的水流交换,妨碍了污染物的稀释扩散,使得来自无锡市区的污水成了五里湖的主要补给水源,加速了富营养化的进程。五里湖水质已达重富营养水平,透明度小于0.5m,缺氧较为严重。营养物在底泥中大量积累,TP和TOC含量分别增高了4.17倍和1.87倍。在春末