Monitoring volcanic hazard using eddy covariance at Solfatara volcano, Naples, Italy

An eddy covariance (EC) station was deployed at Solfatara crater, Italy, June 8–25, 2001 to assess if EC could reliably monitor CO₂ fluxes continuously at this site. Deployment at six different locations within the crater allowed areas of focused gas venting to be variably included in the measured flux. Turbulent (EC) fluxes calculated in 30-min averages varied between 950 and 4460 g CO₂ m⁻² d⁻¹; the highest measurements were made downwind of degassing pools. Comparing turbulent fluxes with chamber measurements of surface fluxes using footprint models in diffuse degassing regions yielded an average difference of 0% (±4%), indicating that EC measurements are representative of surface fluxes at this volcanic site. Similar comparisons made downwind of degassing pools yielded emission rates from 12 to 27 t CO₂ d⁻¹ for these features. Reliable EC measurements (i.e. measurements with sufficient and stationary turbulence) were obtained primarily during daytime hours (08:00 and 20:00 local time) when the wind speed exceeded 2 m s⁻¹. Daily average EC fluxes varied by ±50% and variations were likely correlated to changes in atmospheric pressure. Variations in CO₂ emissions due to volcanic processes at depth would have to be on the same order of magnitude as the measured diurnal variability in order to be useful in predicting volcanic hazard. First-order models of magma emplacement suggest that emissions could exceed this rate for reasonable assumptions of magma movement. EC therefore provides a useful method of monitoring volcanic hazard at Solfatara. Further, EC can monitor significantly larger areas than can be monitored by previous methods.     

Modelling the effects of temperature and food on individual growth and reproduction of Daphnia and their consequences on the population level

Individual based simulations of population dynamics require the availability of growth models with adequate complexity. For this purpose a simple-to-use model (non-linear multiple regression approach) is presented describing somatic growth and reproduction of Daphnia as a function of time, temperature and food quantity. The model showed a good agreement with published observations of somatic growth (r² = 0.954, n = 88) and egg production (r² = 0.898, n = 35). Temperature is the main determinant of initial somatic growth and food concentration is the main determinant of maximal body length and clutch size. An individual based simulation was used to demonstrate the simultaneous effects of food and temperature on the population level. Evidently, both temperature and food supply affected the population growth rate but at food concentrations above approximately 0.4 mg Cl⁻¹ Scenedesmus acutus temperature appeared as the main determinant of population growth.

Four simulation examples are given to show the wide applicability of the model: (1) analysis of the correlation between population birth rate and somatic growth rate, (2) contribution of egg development time and delayed somatic growth to temperature-effects on population growth, (3) comparison of population birth rate in simulations with constant vs. decreasing size at maturity with declining food concentrations and (4) costs of diel vertical migration. Due to its plausible behaviour over a broad range of temperature (2–20 °C) and food conditions (0.1–4 mg Cl⁻¹) the model can be used as a module for more detailed simulations of Daphnia population dynamics under realistic environmental conditions.     

哈尔滨市磨盘山水利枢纽区三维渗流数值模拟模型研究

哈尔滨市磨盘山水库防渗型式拟采用垂直砼防渗墙及帷幕灌浆方案。为了分析左岸不同长度防渗对水库渗漏量及坝基渗透稳定性的影响。对磨盘山水利枢纽区渗流场进行了三维数值模拟研究。结果表明，随着左岸防渗长度的增大。渗漏量及下游渗流出口渗透比降基本呈线性趋势减小。当左岸防渗长度为200m时。渗漏量及各渗透比降均小于允许值，为推荐的渗控方案。     

地下水及溶质运移数值模拟系统

介绍了地下水及溶质运移数值模拟的研究进展、数值模型的分类和数值模拟的步骤。根据数值模拟中涉及到的地质信息资源的特点，提出了基于地理信息系统(GIS)技术的地下水及溶质运移数值模拟系统的设计思路。该系统能够对模拟区域的基础地质资料、水文地质勘探资料以及数值模型的输入、输出资料进行动态管理，而且随着勘探、科研和生产的不断进行，通过对数据库的及时更新以及对数值模型的不断校正，可为不同的科研和生产方案的预测与评价提供信息资源。     

The design of post‐mining landscapes using geomorphic principles

Nature can provide analogues for post‐mining landscapes in terms of landscape stability and also in terms of the rehabilitated structure ‘blending in’ with the surrounding undisturbed landscape. In soil‐mantled landscapes, hillslopes typically have a characteristic pro?le that has a convex upper hillslope pro?le with a concave pro?le lower down the slope. In this paper hillslope characteristic form is derived using the area–slope relationship from pre‐mining topography at two sites in Western Australia. Using this relationship, concave hillslope pro?les are constructed and compared to linear hillslopes in terms of sediment loss using the SIBERIA erosion model. It is found that concave hillslopes can reduce sediment loss by up to ?ve times that of linear slopes. Concave slopes can therefore provide an alternative method for the construction of post‐mining landscapes. An understanding of landscape geomorphological properties and the use of erosion models can greatly assist in the design of post‐mining landscapes. Copyright © 2003 John Wiley & Sons, Ltd.     

基于GIS和SD方法的社会经济发展与生态环境响应动态模拟预测研究

在着重于空间现象分析处理的GIS与动态模拟的SD方法结合的模式下，动态模拟预测了不同社会经济发展策略对生态环境的影响，揭示了生态经济系统社会经济发展与生态环境相互影响关系，为区域可持续发展提供决策支持服务。     

水下无源导航系统仿真匹配算法研究

针对水下无源导航(地形辅助导航)系统的特点，充分利用各种海洋地理特征传感器获取信息，进行数据预处理、算法分析、性能评估，然后进行智能数据融合，进而完成水下无源导航系统模块化的仿真设计，为水下无源导航系统的进一步研究提供可行的方法和技术。仿真结果表明，采用水下地形匹配辅助导航，能够有效地修正惯导误差达到设计要求，为真实海洋环境下的精确导航提供可靠的技术依据。     

Atomistic simulation of the structure and elastic properties of pyrite (FeS2) as a function of pressure

Interatomic potential parameters have been derived at simulated temperatures of 0 K and 300 K to model pyrite FeS₂. The predicted pyrite structures are within 1% of those determined experimentally, while the calculated bulk modulus is within 7%. The model is also able to simulate the properties of marcasite, even though no data for this phase were included in the fitting procedure. There is almost no difference in results obtained for pyrite using the two potential sets; however, when used to model FeS₂ marcasite, the potential fitted at 0 K performs better. The potentials have also been used to study the high-pressure behaviour of pyrite up to 44 GPa. The calculated equation of state gives good agreement with experiment and shows that the Fe–S bonds shorten more rapidly that the S–S dimer bonds. The behaviour of marcasite at high pressure is found to be similar to that of pyrite.     

Are catch-up reefs an artefact of coring?

Drill cores through modern coral reefs commonly show a time lag in reef initiation followed by a phase of rapid accretion to sea level from submerged foundations – the so-called ‘catch-up response’. But because of the difficulty of drilling in these environments, core distribution is usually restricted to accessible areas that may not fully represent reef history, especially if the reef initiated in patches or developed with a prograde or retrograde geometry. As a consequence, core data have the potential to give a misleading impression of reef development, particularly with respect to the timing of initiation and response to sea-level rise. Here, we use computer models to simulate keep-up reef development and, from them, quantify variations in the timing of reef initiation and accretion rate using mock cores taken through the completed simulations. The results demonstrate that cores consistently underestimate the timing of reef initiation and overestimate the reef accretion rate so that, statistically, a core through a keep-up reef will most likely produce a catch-up pattern – an initiation lag followed by a phase of rapid accretion to sea level. This implies that catch-up signatures may be an artefact of coring and that keep-up reefs are significantly more common than previous core studies claim.     

Pull-down basin in the central part of Japan due to subduction-induced mantle flow

The driving force for the basin subsiding against isostatic balance in and around Lake Biwa in the Kinki district, Japan is discussed. The lake region is characterized by strong negative Bouguer anomalies, especially by a steep horizontal gradient zone of gravity anomaly running along the western margin of the lake. The large negative anomaly (>50 mgal) cannot be explained by low-density sediments beneath it. A down-warping structure extending to the Moho depth should be taken into account. This conjecture has been strongly supported by a short-period receiver function imaging, which shows a clear offset of about 8 km for the Moho discontinuity under the steep gravity gradient zone.A question arises as to what is the driving force to create such a large down-warping structure. We consider that the subduction of the shallow-dipping slab under the region (Philippine Sea Slab) may cause crustal deformation by dragging the viscous mantle downward. In order to verify this model, we simulated the induced mantle flow due to the subduction of the Philippine Sea Slab and the pressure distribution on the crust–mantle boundary. This numerical experiment showed that the induced flow makes a strong negative pressure zone under the lake region if the slab has a vertical offset along the direction of subduction. This offset of the slab is consistent with plate models deduced from hypocentral distributions and Sp phases of the deep-focus earthquakes.