龙滩水库诱发地震的孔隙压力扩散特征

利用孔隙压力扩散机制,分析了广西省龙滩水库与蓄水有关的小震精定位数据,获得了5个分区的流体孔隙压力扩散特征.龙滩水库与蓄水有关的小震活动表现出明显分区现象,使用多源孔隙压力扩散过程触发机理,计算出5个分区流体孔隙压扩散系数值为0.04～0.65 m2/s.根据5个分区流体孔隙压力扩散程度,将这5区划分为I、Ⅱ、Ⅲ三种类...     

REE diffusion in calcite

Chemical diffusion of four rare-earth elements (La, Nd, Dy and Yb) has been measured in natural calcite under anhydrous conditions, using rare-earth carbonate powders as the source of diffusants. Experiments were run in sealed silica capsules along with finely ground calcite to ensure stability of the single-crystal samples during diffusion anneals. Rutherford backscattering spectroscopy (RBS) was used to measure diffusion profiles. The following Arrhenius relations were obtained over the temperature range 600–850°C: D_La =2.6×10⁻¹⁴ exp(−147±14 kJ mol⁻¹/RT) m² s⁻¹, D_Nd =2.4×10⁻¹⁴ exp(−150±13 kJ mol⁻¹/RT) m² s⁻¹, D_Dy =2.9×10⁻¹⁴ exp(−145±25 kJ mol⁻¹/RT) m² s⁻¹, D_Yb =3.9×10⁻¹² exp(−186±23 kJ mol⁻¹/RT) m² s⁻¹. In contrast to previous findings for refractory silicates (e.g. zircon), differences in transport rates among the REE are not pronounced over the range of temperature conditions investigated in this study. Diffusion of the REE is significantly slower than diffusion of the divalent cations Sr and Pb and slower than transport of Ca and C at temperatures above 650°C. Fine-scale zoning and isotopic and REE chemical signatures may be retained in calcites under many conditions if diffusion is the dominant process affecting alteration.     

Self-diffusion in liquid Fe at high pressure

The self-diffusivity of liquid iron has been measured between 2 and 20 GPa and 1883 and 2393 K. There is a large pressure dependence across the pressure range studied and simple Arrhenius extrapolation to outer core conditions yields unrealistically low diffusivities. The large pressure and temperature range covered, however, allows different a priori diffusion models to be tested using the experimental results. The free volume model [J. Chem. Phys. 31 (1959) 1164] is found to be in excellent agreement with the experimental data and provides a basis for extrapolation of the results to outer core conditions. An extrapolation to core–mantle boundary (CMB) pressure at 4300 K yields a self-diffusivity for liquid iron of 2.5×10⁻⁵ cm² s⁻¹, which corresponds to a viscosity of 11 mPa s, in reasonable agreement with previous, theoretical, estimates.     

Strain localisation in bimineralic rocks: Experimental deformation of synthetic calcite–anhydrite aggregates

Deformation of synthetic calcite–anhydrite aggregates to large shear strains (up to γ = 12.4 at 600 °C, 300 MPa confining pressure and a constant angular displacement rate corresponding to a shear strain rate of 10^− 3 s^− 1) resulted in the first experimental observation of strain localisation from initially homogeneous rocks. In contrast to experiments on pure calcite and anhydrite, which deformed homogeneously to large strains (γ ≥ 5), all experiments on calcite–anhydrite mixtures resulted in heterogeneous deformation at γ > 1 and the formation of narrow localised bands in the microstructures at γ > 4. In these bands, the amount of strain is at least twice as large as in the rest of the sample and individual grains of the same phase cluster and align, thereby forming microstructural layering similar to planar fabrics in natural mylonites. A switch in deformation mechanism in anhydrite from dislocation creep to diffusion creep and/or grain boundary sliding occurs simultaneously with strain localisation. It is concluded that deformation-induced heterogeneous phase distributions cause local strength differences initiating strain localisation in the calcite–anhydrite mixtures. The study suggests that the presence of two phases in combination with a change in deformation mechanism may be responsible for strain localisation in natural poly-mineralic mylonites.     

An empirical test of helium diffusion in apatite: borehole data from the Otway basin, Australia

We have analyzed helium ages of apatites from several boreholes in the Otway basin, Australia, to evaluate whether laboratory helium diffusivity can be accurately extrapolated to conditions relevant in nature. Downhole apatite helium ages define a broad swath of values from 78–71 Ma at the surface (15°C) to nearly zero at depths corresponding to ambient temperatures of 80°C. The width of the swath results from uncertainties in corrected borehole temperatures, differences in the thermal history experienced by the various boreholes, and possibly from slightly different helium diffusivities among the detrital apatite grains studied. In the eastern Otway basin, the shape and position of the helium age profile is in good agreement with predictions based on the extrapolation of laboratory diffusivity data for Durango apatite coupled with published thermal histories for this part of the basin. In contrast, helium ages are much younger than predicted in the western Otway basin. Based on measured ages from Otway sediments, which have been essentially isothermal over the last few million years, an empirical diffusivity (D/a²) of 2×10⁻¹⁵ s⁻¹ can be assigned to apatites residing at downhole temperatures of 67–97°C. This empirical diffusivity is consistent with laboratory diffusion measurements, demonstrating that such measurements are reasonably accurate and can be applied with confidence to natural geologic settings. Given this confirmation of the laboratory diffusivity data, the discrepancy between the observed and modeled helium age profiles in the western Otway basin suggests that these sediments recently experienced higher temperatures than presently supposed.     

基于镭同位素分布的黄海和东海垂直混合速率计算

在黄海和东海采样测定了水体中的镭同位素分布,用平流扩散模型描述镭同位素分布,最小二乘方法计算了垂直涡动扩散系数和上升流或下降流流速.结果给出北黄海中部、南黄海中部、浙江沿岸和台湾北部海域存在上升流,流速分别为0.46×10^-3cm·s^-1、0.17×10^-3~1.39×10^-3cm·s^-1、2.02×10^-3~3.04×10^-3 cm·s^-1和1.06×10^-3~2.51×10^-3 cm·s^-1.北黄海中部和东海东北部存在下降流.流速分别为-2.30×10^-3 cm·s^-1和-0.61×10^-3~-2.10×10^-3 cm·s^-1.计算同时给出的垂直涡动扩散系数为5.84~48.2 cm²·s^-1,平均值为22.3 cm²·s^-1.北黄海和浙江沿岸上升流流速与文献的结果一致;北黄海中部存在下降流与文献的结论一致.本研究结果与文献结果一致是对所建立的方法的肯定,也是对文献研究结果的支持.     

2013年岷漳地震和2014年景谷地震的能矩比对比研究

2013年岷漳地震和2014年景谷地震两个地震的矩震级均为M_W6.1,但相同等级烈度区域的面积差异显著,前者均明显大于后者.地面震动的强烈程度与震源破裂过程中辐射的地震波高频成分密切相关,作者系统地测定与研究了这两个地震的震源谱和反映其震源特征的能矩比等震源参量.结果表明,虽然岷漳地震和景谷地震释放的地震矩基本相同,分别为1.6×10¹⁸Nm和1.8×10¹⁸Nm,而且前者还比后者略小些,但其释放的地震波辐射能却差异显著,分别为1.3×10¹⁴J和0.6×10¹⁴J;因而前者的能矩比（8.1×10^-5）是后者能矩比（3.3×10^-5）的大约2.5倍.相应地,岷漳地震的视应力和应力降也都明显大于景谷地震的视应力和应力降.此外,这两个地震震源谱结果均显示出其高频趋势均按角频率ω的-2幂次衰减,符合"ω²模式";采用理论震源谱公式计算结果分别约占实际结果的54%和50%,两个地震的地震波辐射能解析解之间的差异依然约2.3倍.震源破裂特征的显著差异导致岷漳地震单位地震矩释放的地震波辐射能更多,即岷漳地震能矩比更大,是导致这两个矩震级均为M_W6.1地震的地震烈度差异显著的一个重要原因.     

Al depth profile in Apollo 15 drill core

Cosmic-ray-produced²⁶Al (t_1/2 = 7.05 × 10⁵ years) has been measured in the Apollo 15 long core (surface to 390 g/cm²—218 cm) for study of galactic cosmic ray production profiles, using accelerator mass spectrometry. The results are in general accord with non-destructive counting data obtained earlier, but systematically lower, and significantly higher precision. From this experiment the half-attenuation length for²⁶Al production can be calculated to be 122 g/cm² (150–400 g/cm² region) after normalizing the data to average chemical composition. The⁵³Mn (t_1/2 = 3.7 × 10⁶ years) production profile in deep cores was also compiled to date. The half-attenuation length for⁵³Mn production was calculated to be 123 g/cm² (150–400 g/cm² region).     

Diffusive Oxygen Emitters for Enhancement of Aerobic In Situ Treatment

Aerobic biodegradation can be enhanced within contaminant plumes by elevating typically low dissolved oxygen (D.O.) levels using materials or devices that passively release oxygen. We have developed passive devices that provide a uniform, steady, long-term source of oxygen by diffusion from pressurized polymeric tubing and report test results under lab and field conditions. Lab flow-through reactor tests were conducted to determine the diffusion coefficient (D) of oxygen through four readily available tubing materials. Oxygen diffusion was greatest through Tygon® 3350 platinum-cured silicone (D = 6.67 ± 10^-7 cm²/sec), followed by 2075 Ultra Chemical Resistant Tygon (1.59 ± 10^-7 cm²/sec), 2275 High Purity Tygon (5.11 ± 10^-8 cm²/sec), and low-density polyethylene (LDPE; 1.73 ± 10^-8 cm²/sec). Variable-pressure release tests with LDPE resulted in very close estimates of D, which confirmed that mass transfer is controlled by diffusion and that the concentration gradient is a valid approximation of the chemical potential driving diffusion. LDPE emitter devices were designed and installed in seven 8-inch-diameter well screens across a portion of a gasoline plume at a former service station. With the devices pressurized to 620.5 kPag (kilopascals gauge) late in the test, steady-state D.O. concentrations reached as high as 25 mg/L, comparing favorably to the value predicted using the mass-transfer coefficient estimated from the lab test (26.3 mg/L). The method can also be used to release other gases for other reasons: gaseous tracers (i.e., sulphur hexafluoride, helium, and argon), hydrogen (for reductive dechlorination), or light alkanes (for cometabolic biodegradation of methyl tertiary butyl ether [MTBE] or chlorinated solvents).     

Uranium and thorium diffusion in diopside

This paper presents new experimental data on the tracer diffusion rates of U and Th in diopside at 1 atm and 1150–1300°C. Diffusion couples were prepared by depositing a thin layer of U–Th oxide onto the polished surface of a natural diopside single crystal, and diffusion profiles were measured by ion microprobe depth profiling. For diffusion parallel to [001] the following Arrhenius relations were obtained: log₁₀D_U=(−5.75±0.98)−(418±28 kJ/mol)/2.303RT log₁₀D_Th=(−7.77±0.92)−(356±26 kJ/mol)/2.303RT. The diffusion data are used to assess the extent to which equilibrium is obtained during near fractional melting of a high-Ca pyroxene bearing mantle peridotite. We find that the diffusion rates for both elements are slow and that disequilibrium between solid and melt will occur under certain melting conditions. For near-fractional adiabatic decompression melting at ascent rates >3 cm/yr, high-Ca pyroxene will exhibit disequilibrium effects. High-Ca pyroxene will become zoned in U and Th and the melts extracted will be depleted in these incompatible elements relative to melts produced by equilibrium fractional melting. U and Th diffusivities in high-Ca pyroxene are similar, and diffusive fractionation of these elements will be limited. Numerical solutions to a dynamic melting model with diffusion-controlled chemical equilibration indicate that the activity ratio [²³⁰Th/²³⁸U] in a partial melt of spinel peridotite will be slightly less than 1 for a broad range of melting parameters. This result reinforces the already widely accepted conclusion that melting of spinel peridotite cannot account for ²³⁰Th excesses in mid-ocean ridge and ocean island basalts, and that garnet must therefore be present over part of the melting column.     

Ecological–economic modeling as a tool for watershed management: A case study of Lake Qionghai watershed, China

This paper presents an ecological–economic model for a lake and its watershed systems. We describe the linkage between the watershed system and the lake aquatic ecosystem and the modeling process. The lake–watershed system was divided into six subsystems: social system, economic system, terrestrial ecosystem, lake water system, pollutant system, and lake aquatic ecosystem. The model equations were constructed based on five main assumptions. The Lake Qionghai watershed in southwestern China, which is undergoing rapid eutrophication, was used as a case study. The targeted goals for total phosphorus (TP) and chlorophyll a (Chl a) concentrations in the lake in 2015 are 0.025 and 10.0 mg m⁻³, respectively. We present two scenarios from 2004 to 2015 based on the ecological–economic model. In both scenarios, the TP and Chl a concentrations in the lake are predicted to increase under the effects of watershed pressures and the targeted goals cannot be met. The application of techniques to reduce pollutants loading and the corresponding pollutants reductions are reflected again in the constructed model. The model predicts that TP and Chl a concentrations will decrease to 0.024 and 7.71 mg m⁻³, respectively, which meet the targeted thresholds. The model results provide directions for local government management of watersheds and lake aquatic ecosystem restoration.     

He diffusion in olivine

Helium diffusion in olivine (dunite xenolith) has been measured in the temperature range 1180–1460°C; a linear Arrhenius function was obtained with an activation energy of 120₋₂₇⁺³² kcal/mole, and a pre-exponential factor (D₀) = 2.2 × 10⁸ cm²/s.

Diffusion mechanisms are not a viable means of degassing He from the mantle. Olivine phenocrysts can be expected to retain previously trapped He, during cooling in extrusive basalts, provided the flow units are thinner than 50 m; xenoliths will retain mantle He signatures only if magma transport times are less than 50 years, or if the He fugacity in the magma is high enough to prevent xenolith degassing. The lower oceanic crust is probably substantially degassed of He. Trapped He will be qualitatively retained in quenched submarine basalt glass only if the cooling rate is faster than 5 × 10¹⁴°C/m.y.; glass at several centimeters depth in a basalt flow (near the spherulite zone) will have cooling rates lower than this, so He loss may be significant in many basalt glass samples.     

Andesitic magma degassing investigated through H2O vapour–melt partitioning of halogens at Soufrière Hills Volcano, Montserrat (Lesser Antilles)

Magma degassing at Soufrière Hills Volcano (SHV) is characterised by an almost permanent SO₂ flux and a HCl production rate which mainly depends on dome growth rate. Degassing processes have been studied through textural, H₂O and halogen analyses of clasts collected between 1995 and 2006 on the dome and in pyroclastic flows and vulcanian eruption deposits. Cl, Br and I are strongly depleted in melts during H₂O degassing with no significant Cl–Br–I fractionation, whereas F is almost unaffected. All magmas erupted at SHV have followed a multi-step degassing path from the magma chamber up to a shallow depth ( 1 km, P  20 MPa). From that depth, however, effusive and explosive paths are distinct; vulcanian eruptions are the result of closed system degassing (CSD), while effusive dome growth is the result of CSD up to a very shallow depth (≤ 200 m, P  5–2 MPa) followed by open system degassing (OSD). CSD is modelled using the H₂O solubility law, the perfect gas law and partition coefficients of halogens between a rhyolitic melt and H₂O vapour (d_v − lⁱ). Gas loss characteristic of OSD is modelled using a Rayleigh law. Degassing induced crystallisation is introduced through the ratio of crystallisation and degassing rates, which ranges from 150–500. d_v − l^Cl for OSD ranges between 50–300, increasing with melt Cl content. For CSD, the lower effective d_v − l^Cl ( 20) is attributed to kinetic effects.

Dome forming activity has a greater impact on atmospheric chemistry than vulcanian eruptions because OSD is much more efficient at extracting halogens. The model shows that HCl flux is a good proxy for the dome forming eruption rate. Comparison between model and measured gas compositions suggests a high HBr–BrO conversion rate (BrO/Total Br  1/3) in the SHV gas plume.

The degassing behaviour of Cl, Br and I implies similar Cl/Br ( 160) and Br/I ( 90) in initial melts, volcanic clasts and high temperature gases. The low Cl/Br at SHV compared to other island arcs ( 250–300) is attributed to a shallow, pre-eruptive Br enrichment. The almost permanent dome extrusion at SHV since 1995 has likely had a significant regional atmospheric impact because of the very efficient effusive degassing and the high conversion rate of halogens into reactive species within the gas plume.     

Water flow path interactions with soil hydraulic properties in till soil at Gårdsjön, Sweden

In 1989, in a hydrological research programme within a deacidification project in the Gårdsjön area in southwest Sweden, flow paths and residence times of soil water and groundwater in microcatchments were examined to support the interpretation of the hydrochemical changes. Saturated hydraulic conductivity and soil water retention were analysed on more than 100 cylinder samples. The catchments have shallow sandy-silty till soil with a mean depth in the main catchment of 43 cm. Porosity of the mineral soil in the main catchment was high and ranged from 38 to 85%. The samples from the B-horizon had generally higher porosity. Porosity and the content of organic matter were correlated. The soil water retention was relatively high at all tensions, likely owing to the high content of organic matter. Dissolved organic substances were most probably transported from the shallow soil on the steep sides of the catchment down to the valley where it precipitated. The high porosities could be a consequence of long-term weathering, provided that the organic substances present have increased the leaching of the weathering products. Measured values of saturated hydraulic conductivity were close to log-normally distributed with a mean for all samples of 3 × 10⁻⁵ m s⁻¹. There was a significant increase in conductivity toward the ground surface with the mean conductivity of the samples in the uppermost 10 cm of the mineral soil of 4 × 10⁻⁵ m s⁻¹, which was about 13 times higher than the conductivity of 3 × 10⁻⁶ m s⁻¹ at 1 m depth. From the relationship between runoff at the catchment outlet and groundwater levels, the conductivity was estimated to be 15–200 times higher in the upper soil layer than in the deeper ones. In one profile, 44–64% of the yearly lateral flow was estimated to occur above 30 cm depth. The conductivity was correlated with the content of drainable water, which indicated the importance of the largest pores for the saturated hydraulic conductivity.     

Formation of a salt plume in the Coastal Plain aquifer of Israel: the Be''er Toviyya region

The formation and development of a salt plume (salinity up to 800 mg Cl 1⁻¹) in the inner part of the Coastal Plain aquifer of Israel is analyzed. Massive groundwater exploitation during the 1950s caused a large drop in the water level and formation of a hydrologic depression in the Be'er Toviyya-Kefar Warburg area. The depression reached a maximal depth during the late 1960s; thereafter a reduction in the rate of pumpage led to restoration of water levels and shallowing of the depression, until its complete disappearance towards the end of the 1980s. A spot of high salinity first appeared in 1956, following a deep drawdown in the water levels. This saline plume has been continuously expanding with increasing salinity concentrations (200–800 mg Cl 1⁻¹) in its center. The average rate of radial expansion was about 50 m year⁻¹. The expansion and salinization did not cease as the depression disappeared. Rather, equalization of water levels in wells situated within the plume area with those of situated along its margins resulted in the salinization of the latter within a period of 1 year.

Mass balances for water and chloride contents were made for the period 1967–1990. Taking into consideration the storage change, pumpage, natural replenishment and artificial recharge, the lateral inflow to the depression is estimated as 60 × 10⁶ m³. Upon addition of the chloride balance, and taking into consideration the chloride concentrations of the surrounding fresh water and the apparent possible end-member of the saline source (based on geochemical considerations), the saline inflow is estimated as (40–60) × 10⁶ m³. These estimates indicate that a large amount of saline water penetrated into the aquifer, of about half of the natural replenishment of the study area, with an estimated salinity of 1900–2700 mg Cl 1⁻¹.

It is suggested that the salt plume was formed as a result of a drop in water level combined with a flow of underlying saline water bodies from deeper strata. The chemical composition of the groundwater points to the existence of two saline water bodies of Ca-chloride composition and a marine Br/Cl ratio: (1) saline water with low Na/Cl (0.6), So₄/Cl, and B/Cl ratio; (2) saline water with higher Na/Cl (> 0.6), So₄/Cl, and B/Cl ratios. These chemical compositions resemble Ca-chloride saline waters found in other locations in the Coastal Plain aquifer and in underlying formations. The saline water bodies may occur in either pockets at the bottom of the aquifer or lumachelle and sandstone layers of high hydraulic conductivity in underlying sediments.     

PSYCHIC – A process-based model of phosphorus and sediment transfers within agricultural catchments. Part 2. A preliminary evaluation

This paper describes the preliminary evaluation of the PSYCHIC catchment scale (Tier 1) model for predicting the mobilisation and delivery of phosphorus (P) and suspended sediment (SS) in the Hampshire Avon (1715 km²) and Herefordshire Wye (4017 km²) drainage basins, in the UK, using empirical data. Phosphorus and SS transfers to watercourses in the Wye were predicted to be greater than corresponding delivery in the Avon; SS, 249 vs 33 kg ha⁻¹ yr⁻¹; DP, 2.57 vs 1.26 kg ha⁻¹ yr⁻¹; PP, 2.20 vs 0.56 kg ha⁻¹ yr⁻¹. The spatial pattern of the predicted transfers was relatively uniform across the Wye drainage basin, whilst in the Avon, delivery to watercourses was largely confined to the river corridors and small areas of drained land. Statistical performance in relation to predicted exports of P and SS, using criteria for relative error (RE) and root mean square error (RMSE), reflected the potential shortcomings associated with using longer-term climate data for predicting shorter-term (2002–2004) catchment response and the need to refine calculations of point source contributions and to incorporate additional river basin processes such as channel bank erosion and in-stream geochemical processing. PSYCHIC is therefore best suited to characterising longer-term catchment response.     

Be andCl depth profiles in an Apollo 15 drill core

Cosmic-ray-produced¹⁰Be (t_1/2 = 1.6 × 10⁶ years) and³⁶Cl (t_1/2 = 3.0 × 10⁵ years) have been measured in the Apollo 15 long core for study of galactic cosmic ray production profiles using tandem accelerator mass spectrometry. From these experiments, the half-attenuation length for¹⁰Be production and³⁶Cl production were calculated to be 120 g/cm² and 132 g/cm² (150–400 g/cm² region). The measured half-attenuation length for¹⁰Be is slightly longer than that predicted by the Reedy-Arnold theoretical model. The flatter and somewhat deeper maximum seen in the³⁶Cl profile compared to the¹⁰Be,²⁶Al and⁵³Mn profiles can be explained by production from secondary thermal neutrons on³⁵Cl.     

Kinetic energy–rainfall intensity relationship for Central Cebu, Philippines for soil erosion studies

In the study of soil erosion, specifically on detachment of soil particles by raindrop impact, kinetic energy is a commonly suggested indicator of the raindrop's ability to detach soil particles from the soil mass. Since direct measurement of kinetic energy requires sophisticated and costly instruments, the alternative approach is to estimate it from rainfall intensity. The present study aims at establishing a relationship between rainfall intensity and kinetic energy for rainfalls in Central Cebu, Philippines as a preface of a wider regional investigation.

Drop size distributions of rainfalls were measured using the disdrometer RD-80. There are two forms of kinetic energy considered here. One is kinetic energy per unit area per unit time (KE_R, J m⁻² h⁻¹) and the other is kinetic energy per unit area per unit depth (KE, J m⁻² mm⁻¹). Relationships between kinetic energy per unit area per unit time (KE_R) and rainfall intensity (I) were obtained using linear and power relations. The exponential model and the logarithmic model were fitted to the KE–I data to obtain corresponding relationships between kinetic energy per unit area per unit depth of rainfall (KE) and rainfall intensity (I). The equation obtained from the exponential model produced smaller standard error of estimates than the logarithmic model.     

随钻地层测试压力响应数学模型及物理实验考察

基于随钻地层测试基本原理和渗流理论,建立了圆柱坐标系下压力下降和压力恢复阶段的压力扩散微分方程,确定出随钻地层压力测试定解问题,通过对定解问题的边界进行齐次化处理,并采用有限差分方法建立了四维空间下的显格式差分方程,以解析方法和有限差分方法分析了点源单相单测试室均质球形渗流问题,两种方法计算误差小于2.00%,检验了有限差分方法的正确性.采用有限差分方法定量分析了不同渗透率（0.1×10^-3μm²、1.0×10^-3μm²和10.0×10^-3μm²）地层下,表皮效应、储集效应、地层非均质性、抽吸速率、抽吸探头半径等因素对压力响应的影响,分析表明,这些因素对压力响应曲线的影响十分显著,主要表现在抽吸压降、压降速率、压力恢复速率、恢复时间等几个方面,而且不同渗透率地层的敏感程度差异比较大.建立了地层压力测试模拟实验平台,模拟了3种不同渗透率（108.81×10^-3μm²、16.10×10^-3μm²和4.79×10^-3μm²）地层的压力抽吸测试,并采用有限差分方法对实验进行了数值模拟.结果表明,数值模拟与实测压力响应一致性非常好,二者绝对误差小于1.60 MPa,相对误差低于5.00%（最大误差4.92%）,验证了有限差分数值方法的正确性和准确性.     

The physico-chemical conditions of Turkwel Gorge Reservoir, a new man made lake in Northern Kenya

Variations in some physical, chemical, and nutrient conditions were investigated at Turkwel Gorge Reservoir and its inflowing river, Suam between 1994 and 1995. Seasonal changes in inflow volume had the greatest impact on the reservoir and river conditions investigated. A wide fluctuation in inflow volume combined with a regulated outflow independent of season resulted in a draw down of over 10 m in each year. Flood inflows during the wet season resulted in the lowest values of Secchi depth (range, 0.09–2.16 m), electrical conductivity (EC, range = 140–200 mS cm⁻¹) and total alkalinity (TA, range = 75–111 mg l⁻¹) while the highest values were measured during the dry season. A functional relation between EC and TA (TA = 0.529 mg l⁻¹, EC: R² = 0.876) suggests a predominance of carbonates among the anions. Vertical profiles of temperature and dissolved oxygen (DO) revealed that the reservoir is monomictic with a wide variation in the depth of the daily mixed layer. High values of pH (range = 6.7–8.9) and DO (range = 4.9–9.2 mg l⁻¹) were associated with periods of peak phytoplankton photosynthesis while the lowest values followed reservoir mixing. Peak total nitrogen (TN, range = 119–526 μg l⁻¹) and total phosphorus (TP, range = 8.9–71.6 μg l⁻¹) levels during the wet season resulted from increased river loading. Values of dissolved reactive silica (DRS, range = 0.41–9.77 mg l⁻¹) showed a wet season decline which was related to diatom depletion during the wet season. Annual reservoir areal loading rates of 27.38, 10.90 and 408.5 mg m⁻² were computed for TN, TP and DRS respectively based on estimates of inflowing river loads in 1994.

At the inflowing river Suam, low levels of EC (range = 107–210 μS cm⁻¹) and TA (range = 62–125 mg l⁻¹) occurred during the wet season while the highest levels occurred shortly before the river dried up. The first flood water at the resumption of river inflow in March was characterized by very low levels of DO (range = 1.8–8.2 mg l⁻¹) and high levels of TN (range = 205–3354 μg l⁻¹) and TP (102–1259 μg l⁻¹). River pII (6.9–7.7) and DRs (range = 9.01–19.93 mg l⁻¹) varied irregularly throughout the year.