Histopathological effects of petroleum hydrocarbons and heavy metals on the American oyster (Crassostrea virginica) from Tabasco, Mexico

Oyster tissues may be affected by the concentration of ions in the water (i.e. salinity) and by contaminants such as petroleum hydrocarbons and heavy metals. Oyster populations from three coastal lagoons (Mecoacan, Camen and Machona) in the Mexican state of Tabasco, in the southern Gulf of Mexico, were sampled for pollution studies during June, September and November 1992 and May 1993. No statistically significant relationships were found between the concentration of contaminants in the body tissues of organisms and their shell length, soft tissue weight and particulate matter; however, a significant correlation was found between condition index and salinity (r = −0.72). Generally, the lowest weights, shell lengths and indices of condition were found in Mecoacan. The concentrations of cadmium and zinc were inversely related to salinity (r = −0.52 and r = −0.32, respectively), so a riverine input is suspected. On average, 63% of the individual oysters showed histopathological lesions, which can be related to salinity and to the concentrations of cadmium and the unresolved fraction of hydrocarbons (UCM). The response of each of the tissues analysed was different. The percentage of individuals with damage in the digestive diverticulum increased linearly with UCM (r = 0.71), but in a saturation-response fashion (Y = B₀ − B₀/(1 + X^B1) r = 0.66) with cadmium. The percentage of individuals with damage in the gills increased linearly with weight (r = 0.68), cadmium (r = 0.60) and UCM (r = 0.60). The lesions in the connective tissue decreased linearly (r = −0.82) with salinity, but increased in a saturation-response way (r = 0.83) with cadmium. Finally, the percentage of individuals with lesions in the digestive tube decreased linearly with salinity (r = −0.59). Only the damage to the gills and digestive diverticula were dependent on gonadal maturity, while damage to the connective tissue was dependent on the sex of the individual.     

Volcanic lake systematics I. Physical constraints

 Volcanic lakes have a wide range of characteristics, and we make an attempt to delineate the limiting physical conditions for several lake classes. The ratio between heat input and heat dissipation capacity of a lake constrains the temperature for perfectly mixed steady-state volcanic lakes. Poorly mixed lakes are also conditioned by this ratio, but their temperature structure is also strongly influenced by the hydrodynamics resulting from different mechanisms of heat transfer. The steady-state temperatures of volcanic lakes are largely determined by the magnitude of the volcanic heat influx relative to the surface area of the lake. Small lakes have only a small capacity for heat dissipation and their temperature rises quickly with only small heat inputs; large lakes are buffered against variations in heat input. Both the heat dissipation and meteoric water input into a lake are functions of lake surface area and therefore each lake water temperature demands a certain precipitation rate for mass conservation, independent of lake size. The results of energy/mass-balance modeling shows that under common atmospheric conditions, most steady-state volcanic lakes are unlikely to maintain a temperature in excess of 45–50  °C. Validation of the volcanic lake model was performed using published data from Yugama Lake (Japan) and the Keli Mutu lakes (Indonesia). Also, the model was applied to 24 natural systems to provide a baseline assessment of energy fluxes under the model assumptions so future work on those systems can identify nuances in individual systems that deviate from the simple model conditions. We recommend the model for use in assessing temperature variations and volcanic lake stability in settings with known physical and atmospheric conditions. Application of the energy/mass balance calculations of model lakes provides a genetic classification scheme largely based on physical process parameters. Received: 28 February 1996 / Accepted: 23 November 1996     

Using high hydraulic conductivity nodes to simulate seepage lakes

In a typical ground water flow model, lakes are represented by specified head nodes requiring that lake levels be known a priori. To remove this limitation, previous researchers assigned high hydraulic conductivity (K) values to nodes that represent a lake, under the assumption that the simulated head at the nodes in the high-K zone accurately reflects lake level. The solution should also produce a constant water level across the lake. We developed a model of a simple hypothetical ground water/lake system to test whether solutions using high-K lake nodes are sensitive to the value of K selected to represent the lake. Results show that the larger the contrast between the K of the aquifer and the K of the lake nodes, the smaller the error tolerance required for the solution to converge. For our test problem, a contrast of three orders of magnitude produced a head difference across the lake of 0.005 m under a regional gradient of the order of 10(-3) m/m, while a contrast of four orders of magnitude produced a head difference of 0.001 m. The high-K method was then used to simulate lake levels in Pretty Lake, Wisconsin. Results for both the hypothetical system and the application to Pretty Lake compared favorably with results using a lake package developed for MODFLOW (Merritt and Konikow 2000). While our results demonstrate that the high-K method accurately simulates lake levels, this method has more cumbersome postprocessing and longer run times than the same problem simulated using the lake package.     

Imposex in Morula Granulata as Bioindicator of Tributyltin (TBT) Contamination in the Dampier Archipelago, Western Australia

This study was undertaken to examine the use of imposex in the tropical neogastropod whelk Morula granulata as bioindicator of tributyltin (TBT) contamination. Samples were collected from sites throughout the Dampier Archipelago; both impact sites close to intensive shipping activity and remote control sites. TBT in this area originates mainly from commercial vessels using the Port of Dampier that are >25 m in length and exempt from legislation controlling the use of TBT in antifouling paints. A field survey was undertaken in July 1997; 100 M. granulata were taken at each of 18 sites with varying vessel activity, and examined for imposex. The level of TBT contamination in the Archipelago was low; ranging from <0.3 to 25 ng Sn l⁻¹ in water,<0.3 to 80 ng Sn g−¹ in oyster tissue (Saccostrea cuccullata – prey of M. granulata), and <0.3 to 33 ng Sn g⁻¹ in whelks. Percentage of imposex in M. granulata ranged from 0% to 57%, and was correlated with distance to vessel activity (r=0.646, p<0.05). The only significant relationship of imposex with butyltin concentrations was with TBT in oysters (r=0.826, p<0.05), and dibutyltin (DBT) in whelks (r=0.783, p<0.05). It is postulated that this is because of interspecies variation in ability to metabolize TBT to DBT and MBT (monobutyltin), as well as induction of imposex in earlier life stages, so that imposex in adults could be a reflection of prior contamination. M. granulata is a suitable bioindicator of TBT contamination in the Dampier Archipelago, which could be used throughout the Indo-Pacific given its abundance in the region.     

Simulating the response of a closed‐basin lake to recent climate changes in tropical West Africa (Lake Bosumtwi,Ghana)

Historical changes in the level of Lake Bosumtwi, Ghana, have been simulated using a catchment‐scale hydrological model in order to assess the importance of changes in climate and land use on lake water balance on a monthly basis for the period 1939–2004. Several commonly used models for computing evaporation in data‐sparse regions are compared, including the Penman, the energy budget, and the Priestley–Taylor methods. Based on a comparison with recorded lake level variations, the model with the energy‐budget evaporation model subcomponent is most effective at reproducing observed lake level variations using regional climate records. A sensitivity analysis using this model indicates that Lake Bosumtwi is highly sensitive to changes in precipitation, cloudiness and temperature. However, the model is also sensitive to changes in runoff related to vegetation, and this factor needs to be considered in simulating lake level variations. Both interannual and longer‐term changes in lake level over the last 65 years appear to have been caused primarily by changes in precipitation, though the model also suggests that the drop in lake level over the last few decades has been moderated by changes in cloudiness and temperature over that time. Based on its effectiveness at simulating the magnitude and rate of lake level response to changing climate over the historical record, this model offers a potential future opportunity to examine the palaeoclimatic factors causing past lake level fluctuations preserved in the geological record at Lake Bosumtwi. Copyright © 2006 John Wiley & Sons, Ltd.     

A water flow model of the active crater lake at Aso volcano,Japan: fluctuations of magmatic gas and groundwater fluxes from the underlying hydrothermal system

The first crater of Nakadake, peak of Aso volcano, Japan, contains a hot water lake that shows interesting variations in water level and temperature. These variations were discovered by precise, continuous observations of the lake independent of precipitation. We developed a numerical model of a hot crater lake and compared with observational data for the period from July 2006 to January 2009. The numerical model revealed seasonal changes in mass flux (75–132 kg/s) and enthalpy (1,840–3,030 kJ/kg) for the fluid supplied to the lake. The relation between the enthalpy and mass flux indicates that the bottom input fluid is a mixture of high- and low-temperature fluids. Assuming a mixture of high-temperature steam at 800°C and liquid water at 100°C, we evaluated the liquid and steam fluxes. The liquid water flux shows a seasonal increase lagging behind the rainy season by 2 months, suggesting that the liquid water is predominantly groundwater. The fluctuation pattern in the flux of the high-temperature steam shows a relation with the amplitude of volcanic tremor, suggesting that heating of the hydrothermal system drives the tremor. Consequently, precise observations of a hot crater lake represent a potential method of monitoring volcanic hydrothermal systems in the shallow parts of the volcanoes.     

Analysis of the hydrological response of a tropical terminal lake,Lake Abiyata (Main Ethiopian Rift Valley) to changes in climate and human activities

Integrated dynamic water and chloride balance models with a catchment‐scale hydrological model (PRMS) are used to investigate the response of a terminal tropical lake, Lake Abiyata, to climate variability and water use practices in its catchment. The hydrological model is used to investigate the response of the catchment to different climate and land‐use change scenarios that are incorporated into the lake model. Lake depth–area–volume relationships were established from lake bathymetries. Missing data in the time series were filled using statistical regression techniques. Based on mean monthly data, the lake water balance model produced a good agreement between the simulated and observed levels of Lake Abiyata for the period 1968–83. From 1984 onwards the simulated lake level is overestimated with respect to the observed one, while the chloride concentration is largely underestimated. This discrepancy is attributed to human use of water from the influent rivers or directly from the lake. The simulated lake level and chloride concentration are in better agreement with observed values (r² = 0·96) when human water use for irrigation and salt exploitation are included in the model. A comparison of the simulation with and without human consumption indicates that climate variability controls the interannual fluctuations and that the human water use affects the equilibrium of the system by strongly reducing the lake level. Sensitivity analysis based on a mean climatic year showed that, after prolonged mean climatic conditions, Lake Abiyata reacts more rapidly to an abrupt shift to wetter conditions than to dry conditions. This study shows the significant sensitivity of the level and salinity of the terminal Lake Abiyata to small changes in climate or land use, making it a very good ‘recorder’ of environmental changes that may occur in the catchment at different time scales. Copyright © 2004 John Wiley & Sons, Ltd.     

Inter‐annual changes of alpine inland lake water storage on the Tibetan Plateau: Detection and analysis by integrating satellite altimetry and optical imagery

The effects of climate change have a substantial influence on the extremely vulnerable hydrologic environment of the Tibetan Plateau. The estimation of alpine inland lake water storage variations is essential to modeling the alpine hydrologic process and evaluating water resources. Due to a lack of historical hydrologic observations in this remote and inaccessible region, such estimations also fill a gap in studies on the continuous inter‐annual and seasonal changes in the inland lake water budget. Using Lake Siling Co as a case study, we derived a time‐series of lake surface extents from MODIS imagery, and scarce lake water level data from the satellite altimetry of two sensors (ICESat/GLAS and ENVISAT RA‐2) between 2001 and 2011. Then, based on the fact that the rise in lake water levels is tightly dependent on the expansion of the lake extent, we established an empirical model to simulate a continuous lake water level dataset corresponding to the lake area data during the lake's unfreezing period. Consequently, from three dimensions, the lake surface area, water level and water storage variations consistently revealed that Lake Siling Co exhibited a dramatic trend to expand, particularly from 2001 to 2006. Based on the statistical model and lake area measurements from Landsat images since 1972, the extrapolated lake water level and water storage indicate that the lake has maintained a continual expansion process and that the cumulative water storage variations during 1999–2011 account for 66.84% of the total lake water budget (26.87 km³) from 1972 to 2011. Copyright © 2013 John Wiley & Sons, Ltd.     

山体遮挡对滇池风生流的影响初探

用二维风生流数值模型模拟滇池湖流运动。滇池在主导风向西南风作用下,假定湖面风场是均匀的,数值模拟的湖流流态与实测湖流结果相差很大。而考虑山体遮挡影响,根据实测湖流期间现有的风情资料,在湖面上构造一非均匀风场,数值模拟结果与实测值基本一致。山体遮挡对滇池风生流的影响是不容忽视的。建议进一步进行湖流和湖面风向、风速监测,并建立过山气流数学模型,深入研究山体遮挡对湖泊风生流的影响。     

Mercury,organic-mercury and selenium in small cetaceans in Taiwanese waters

Total Hg (∑Hg), organic-Hg (O-Hg) and Se bioaccumulations in small cetaceans distributed in Taiwanese waters of the Taiwan Strait and the southwestern Pacific have been investigated for the first time. The results could represent the baseline metal concentrations of marine mammals in the southwestern Pacific, where volcanic activities are possibly the major source of mercury to the environments. Muscle samples of four species of small cetaceans were collected from animals accidentally caught by tuna-longline fisheries from 1994 through 1995. In total, 53 pantropical spotted dolphins, Stenella attenuata, nine spinner dolphins, S. longirostris, five bottlenose dolphins, Turiops truncatus and four Risso's dolphins, Grampus giseus were analyzed. In addition, two stranded pantropical spotted dolphins were investigated. Cold vapour AAS and ICP-MS were used in the analysis of Hg and Se, respectively. Significant species difference was found in the four species of small cetaceans. Among them, the pantropical spotted dolphin showed the highest mean concentration (mg/kg wet wt.) of both ∑Hg (3.64±2.19) and O-Hg (2.79±1.23), whereas the Risso's dolphin had the highest mean concentrations of Se (1.77±1.29). There was no significant sex difference with respect to metal bioaccumulation in the samples of S. attenuata. Significant correlations between body length (BL) and ∑Hg, as well as O-Hg concentrations were observed in pantropical spotted (Sa) and spinner dolphins (Sl). The linear relationships were Sa: ∑Hg=−8.290+0.066BL, r=0.421; Sl: ∑Hg=−2.735+0.025BL, r=0.875; Sa: O-Hg=−3.723+0.036BL, r=0.408; and Sl: O-Hg=−3.017+0.025BL, r=0.870. However, a demethylation phenomenon that decreasing the percentage of O-Hg coupled with increasing levels of Se was observed when the ∑Hg concentrations in the muscle tissues of dolphins reached 4 mg/kg wet wt.     

鄱阳湖湖泊流域系统水文水动力联合模拟

本文以鄱阳湖湖泊流域系统为研究对象,鉴于该湖泊流域系统尺度较大,下垫面自然属性呈现高度空间异质性且具有流域-平原区-湖泊不同机制的水文水动力过程,为了真实描述湖泊流域间的水文水动力联系及反映不同过程间的作用机制,构建了鄱阳湖湖泊流域联合模拟模型.该模型基于自主研发的流域分布式水文模型WATLAC和湖滨平原区产流模型以及水动力模型MIKE 21 3个不同功能子模型的连接来实现该复杂系统的模拟.模型的联合采用输入-输出驱动及子模型的顺序执行进程,即将五大子流域与平原区入湖径流量作为输入条件来驱动湖泊水动力模型,模拟湖泊水位对流域入湖径流量的响应.以2000-2005年鄱阳湖流域6个水文站点的河道径流量、流域基流指数以及湖泊4个站点的水位资料来率定模型,其中各站点日径流量拟合的纳希效率系数E_ns为0.71~0.84,确定性系数R²介于0.70~0.88之间,而湖泊各站点水位拟合的纳希效率系数E_ns变化为0.88~0.98,确定性系数R²为0.96~0.98,均取得令人满意的率定结果.本文提出的鄱阳湖湖泊流域系统水文水动力联合模拟模型能较为理想再现湖泊水位对流域降雨-径流过程的响应.水位模拟结果进一步表明,该联合模型能用来获取重要的水动力空间变化特征.该模型可作为有效工具定量揭示湖泊流域系统水文水动力过程对气候变化和流域人类活动的响应.     

Evaporation estimation on Lake Titicaca: a synthesis review and modelling

The aim of this study was to validate evaporation models that can be used for palaeo‐reconstructions of large lake water levels. Lake Titicaca, located in a high‐altitude semi‐arid tropical area in the northern Andean Altiplano, was the object of this case study. As annual evaporation is about 90% of lake output, the lake water balance depends heavily on the yearly and monthly evaporation flux. At the interannual scale, evaporation estimation presents great variability, ranging from 1350 to 1900 mm year^?1. It has been found that evaporation is closely related to lake rainfall by a decreasing relationship integrating the implicit effect of nebulosity and humidity. At the seasonal scale, two monthly evaporation data sets were used: pan observations and estimations derived from the lake energy budget. Comparison between these data sets shows that (i) there is one maximum per year for pan evaporation and two maxima per year for lake evaporation, and (ii) pan evaporation is greater than lake evaporation by about 100 mm year^?1. These differences, mainly due to a water depth scale factor, have been simulated with a simple thermal model θ_w(h, t) of a free‐surface water column. This shows that pan evaporation (h = 0·20 m) is strongly correlated with direct solar radiation, whereas the additional maximum of lake evaporation (h = 40 m) is related to the heat restitution towards the atmosphere from the water body at the end of summer. Finally, five monthly evaporation models were tested in order to obtain the optimal efficiency/complexity ratio. When the forcing variables are limited to those that are most readily available in the past, i.e. air temperature and solar radiation, the best results are obtained with the radiative Abtew model (r = 0·70) and with the Makkink radiative/air temperature model (r = 0·67). Copyright © 2007 John Wiley & Sons, Ltd.     

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.     

Application of a spreadsheet hydrological model for computing the long-term water balance of Lake Awassa,Ethiopia

Abstract

The water balance of the closed freshwater Lake Awassa was estimated using a spreadsheet hydrological model based on long-term monthly hydrometeorological data. The model uses monthly evaporation, river discharge and precipitation data as input. The net groundwater flux is obtained from model simulation as a residual of other water balance components. The result revealed that evaporation, precipitation, and runoff constitute 131, 106 and 83 × 10⁶ m³ of the annual water balance of the lake, respectively. The annual net groundwater outflow from the lake to adjacent basins is 58 × 10⁶ m³. The simulated and recorded lake levels fit well for much of the simulation period (1981–1999). However, for recent years, the simulated and recorded levels do not fit well. This may be explained in terms of the combined effects of land-use change and neotectonism, which have affected the long-term average water balance. With detailed long-term hydrogeological and meteorological data, investigation of the subsurface hydrodynamics, and including the effect of land-use change and tectonism on surface water and groundwater fluxes, the water balance model can be used efficiently for water management practice. The result of this study is expected to play a positive role in future sustainable use of water resources in the catchment.     

A multiscale investigation of ground water flow at Clear Lake, Iowa

Ground water flow was investigated at Clear Lake, a 1468-ha glacial lake in north-central Iowa, as part of a comprehensive water quality study. A multiscale approach, consisting of seepage meters (and a potentiomanometer), Darcy's law, and an analytic element (AE) model, was used to estimate ground water inflow to and outflow from the lake. Estimates from the three methods disagreed. Seepage meters recorded a median-specific discharge of 0.25 mum/s, which produced a lake inflow rate between 90,750 and 138,200 m3/d, but no detectable outflow. A wave-induced Bernoulli effect probably compromised both inflow and outflow measurements. Darcy's law was applied to 11 zones around the lake, producing inflow and outflow values of 10,500 and 5000 m3/d, respectively. The AE model, GFLOW, coupled with the parameter estimation model, UCODE, simulated ground water flow in a 700-km2 region using 31 hydraulic head and base flow measurements as calibration targets. The model produced ground water inflow and outflow rates of 14,300 and 9200 m3/d, respectively. Although not a substitute for field data, the model's ability to simulate ground water flow to the lake and the region, estimate uncertainty for model parameters, and calculate a lake stage and associated lake water balance makes it a powerful tool for water quality management and an attractive alternative to the traditional methods of ground water/lake investigation.     

Numerical modeling of EM fields over local anomalies with a vertical axis of symmetry

A method for numerical modeling of quasi-stationary electromagnetic fields in axially symmetric media is proposed. It is based on the direct finite-element method and the use of special basis functions. Assuming cylindrical coordinates r, φ, z, the three-dimensional solution is presented as a superposition of fields with the azimuthal dependence exp(inφ). For each case this results in a system of two equations of elliptical type in two scalar functions in the (r, z) plane.

The discretization leads to the conservative nine-points difference scheme. The system of linear equations is solved by means of the LU-decomposition technique, the band structure of the matrix being taken into account.

The program is tested using analytical results (DC asymptote) for a near-surface inhomogeneity. Comparison also is made with 2D results (H-polarization) for the model of a local well conducting inclusion in a three-layered Earth.     

An improved methodology to measure evaporation from bare soil based on comparison of surface temperature with a dry soil surface

Accurate estimation of the resistances to water vapor movement is a major difficulty in evaluating evaporation from soil. By including the temperature of a dry soil surface (the temperature of the surface of a dry soil column buried in the field), a method for estimating evaporation from soil is proposed. The necessary input variables for the suggested method are temperature, net radiation, and soil heat flux. There are three advantages of the proposed method over the conventional methods. First, soil surface resistance and aerodynamic resistance are not required. Second, the variables included are fewer. Third, measurement and analysis of the parameters involved are relatively easy. Sensitivity analysis shows that the suggested method is sensitive to temperatures. Test experiments were conducted in a sandy field, where a weighing lysimeter was installed. Evaporation from soil, together with the variables specified above, were measured. For temperatures measured by thermocouples, experimental results showed that the mean absolute error (MAE) for the daily evaporation over 22 days was 0.17 mm day⁻¹. The regression between calculated and measured evaporation was highly significant (r²=0.89). Moreover, the intercept and slope of the regression equation were not significantly different from zero and unity, respectively, at the 0.05 probability level. Furthermore, by using the temperatures measured by infrared thermometers, the MAE between measured evaporation and estimated evaporation was 0.15 mm day⁻¹. The regression between them was highly significant (r²=0.94). In addition, the intercept and slope of the regression equation were not significantly different from zero and unity, respectively, at the 0.05 probability level. These results show that evaporation calculated using the proposed method is in good agreement with lysimeter measured values. By comparing with the temperature difference method, it was shown that the suggested method estimated soil evaporation more accurately than the temperature difference method. Therefore, it is concluded that the proposed method is not only a simple way for application, but also an accurate way to estimate soil evaporation.     

南水北调东线调度对南四湖水质的影响

南水北调东线一期工程2013年通水,必将对南四湖水质产生重大影响.借鉴既有成功理论方法,建立南四湖三维水力调配和水质模型,并经实测值验证.利用该三维模型,模拟南水北调东线一个完整调水期中的南四湖常规、干旱调水工况调度运行过程,分析南四湖上级湖、下级湖的水力和水质变化规律,提出加强和改善南水北调东线调度管理的措施和建议,既可充分发挥南水北调工程效益,又可改善南四湖水质.     

Extreme value theory based on the r largest annual events: a robust approach

In Smith (1986, J. Hydrol. 86, 27–43), a family of statistical distributions and estimators for extreme values based on a fixed number r > = 1 of the largest annual events are presented. The method of estimation was numerical maximum likelihood. In this paper, we consider the robust estimation of parameters in such families of distributions. The estimation technique, which is based on optimal B-robust estimates, will assign weights to each observation and give estimates of the parameters based on the data which are well modeled by the distribution. Thus, observations which are not consistent with the proposed distribution can be identified and the validity of the model can be assessed. The method is illustrated on Venice sea level data.