Conditioning mean steady state flow on hydraulic head and conductivity through geostatistical inversion

Nonlocal moment equations allow one to render deterministically optimum predictions of flow in randomly heterogeneous media and to assess predictive uncertainty conditional on measured values of medium properties. We present a geostatistical inverse algorithm for steady-state flow that makes it possible to further condition such predictions and assessments on measured values of hydraulic head (and/or flux). Our algorithm is based on recursive finite-element approximations of exact first and second conditional moment equations. Hydraulic conductivity is parameterized via universal kriging based on unknown values at pilot points and (optionally) measured values at other discrete locations. Optimum unbiased inverse estimates of natural log hydraulic conductivity, head and flux are obtained by minimizing a residual criterion using the Levenberg-Marquardt algorithm. We illustrate the method for superimposed mean uniform and convergent flows in a bounded two-dimensional domain. Our examples illustrate how conductivity and head data act separately or jointly to reduce parameter estimation errors and model predictive uncertainty.This work is supported in part by NSF/ITR Grant EAR-0110289. The first author was additionally supported by scholarships from CONACYT and Instituto de Investigaciones Electricas of Mexico. Additional support was provided by the European Commission under Contract EVK1-CT-1999-00041 (W-SAHaRA-Stochastic Analysis of Well Head Protection and Risk Assessment).     

Random domain decomposition for flow in heterogeneous stratified aquifers

We study two-dimensional flow in a layered heterogeneous medium composed of two materials whose hydraulic properties and spatial distribution are known statistically but are otherwise uncertain. Our analysis relies on the composite media theory, which employs random domain decomposition in the context of groundwater flow moment equations to explicitly account for the separate effects of material and geometric uncertainty on ensemble moments of head and flux. Flow parallel and perpendicular to the layering in a two-material composite layered medium is considered. The hydraulic conductivity of each material is log-normally distributed with a much higher mean in one material than in the other. The hydraulic conductivities of points within different materials are uncorrelated. The location of the internal boundary between the two contrasting materials is random and normally distributed with given mean and variance. We solve the equations for (ensemble) moments of hydraulic head and flux and analyze the impact of unknown geometry of materials on statistical moments of head and flux. We compare the composite media approach to approximations that replace statistically inhomogeneous conductivity fields with pseudo-homogeneous random fields. This work was performed under the auspices of the US Department of Energy (DOE): DOE/BES (Bureau of Energy Sciences) Program in the Applied Mathematical Sciences contract KC-07–01–01 and Los Alamos National Laboratory under LDRD 98604. This work made use of STC shared experimental facilities supported by the National Science Foundation under Agreement No. EAR-9876800. This work was supported in part by the European Commission under Contract No. EVK1-CT-1999–00041 (W-SAHaRA).     

Foreword

The main objective of this special issue of the Bulletin of Earthquake Engineering is to bring to the earthquake engineering community the outcomes of the research project QUAKER (Quantification and Reduction of Seismic Risk through the Application of Advanced Geotechnical Engineering Techniques) financed by European Commission under contract EVG1-CT-2002-00064.     

Simplified approach for design of raft foundations against fault rupture. Part Ⅱ： soil-structure interaction

This is the second paper of two, which describe the results of an integrated research effort to develop a four-step simplified approach for design of raft foundations against dip-slip （normal and thrust） fault rupture. The first two steps dealing with fault rupture propagation in the free-field were presented in the companion paper. This paper develops an approximate analytical method to analyze soil-foundation-structure interaction （SFSI）, involving two additional phenomena： （i） fault rupture diversion （Step 3）; and （ii） modification of the vertical displacement profile （Step 4）. For the first phenomenon （Step 3）, an approximate energy-based approach is developed to estimate the diversion of a fault rupture due to presence of a raft foundation. The normalized critical load for complete diversion is shown to be a function of soil strength, coefficient of earth pressure at rest, bedrock depth, and the horizontal position of the foundation relative to the outcropping fault rupture. For the second phenomenon （Step 4）, a heuristic approach is proposed, which ＂scans＂ through possible equilibrium positions to detect the one that best satisfies force and moment equilibrium. Thus, we account for the strong geometric nonlinearities that govern this interaction, such as uplifting and second order （P-△） effects. Comparisons with centrifuge-validated finite element analyses demonstrate the efficacy of the method. Its simplicity makes possible its utilization for preliminary design.     

Simplified approach for design of raft foundations against fault rupture. Part Ⅰ： free-field

Over the past few decades, earthquake engineering research mainly focused on the effects of strong seismic shaking. After the 1999 earthquakes in Turkey and Taiwan, and thanks to numerous cases where fault rupture caused substantial damage to structures, the importance of faulting-induced deformation has re-emerged. This paper, along with its companion （Part Ⅱ）, exploits parametric results of finite element analyses and centrifuge model testing in developing a four-step semi-analytical approach for analysis of dip-slip （normal and thrust） fault rupture propagation through sand, its emergence on the ground surface, and its interaction with raft foundations. The present paper （Part Ⅰ） focuses on the effects of faulting in the absence of a structure （i.e., in the free-field）. The semi-analytical approach comprises two-steps： the first deals with the rupture path and the estimation of the location of fault outcropping, and the second with the tectonically- induced displacement profile at the ground surface. In both cases, simple mechanical analogues are used to derive simplified semi-analytical expressions. Centrifuge model test data, in combination with parametric results from nonlinear finite element analyses, are utilized for model calibration. The derived semi-analytical expressions are shown to compare reasonably well with more rigorous experimental and theoretical data, thus providing a useful tool for a first estimation of near-fault seismic hazard.     

The Barents Sea impact study (BASIS): methodology and first results

Global and climate changes are subject to scientific, societal and political debates. Recent observational evidence and results of global climate models have identified the circumpolar North as a region particularly susceptible to future climate change. To understand and assess the consequences of these changes for environmental and societal components of the European Arctic, the Barents Sea Impact Study (BASIS), an EU-funded integrated regional impact study (IRIS) has been carried out (ENV4-CT-97-0637). In common with global integrated assessments (IAs), IRISs also take a holistic view on climate change and its impact. Contrary to IAs, however, IRISs adopt a regional to sub-regional spatial scale. BASIS was carried out by an interdisciplinary team of specialists from 13 institutions in 6 countries. Major results pertain to impacts of possible climate change on marine and terrestrial ecosystems, freshwater hydrology, marine trace gas budgets, forestry and fishery. However, in this paper we focus on the major methodological aspects of an IRIS in general and on methods applied in BASIS in particular.     

Assessing morphological degradation in running waters using Blackfly communities (Diptera, Simuliidae): Can habitat quality be predicted from land use?

Based on benthic macroinvertebrate samples from 57 sites of streams in the European Central Highlands (Ecoregion 9; Illies (1978). Limnofauna Europaea, Stuttgart), the composition and the abundance of Simuliid species were analysed in relation to hydromorphological and land use parameters. Sampling sites were located at two stream types differing in catchment geology, stream morphology and channel width. Land use data were taken from the official German information system for cartography and topography (ATKIS) and the German River Habitat Survey. Land use categories suitable to describe the sampling sites were ‘agricultural land’, ‘urban’ and near-natural areas. Hydromorphological parameters of the sites were recorded using the AQEM site protocol (AQEM Consortium (2002). Manual for the application of the AQEM system. A comprehensive method to assess European streams using benthic macroinvertebrates, developed for the purpose of the Water Framework Directive. EVK1-CT1999-00027, Version 1.0. Available via the Internet from www.aqem.de). Relevant parameters controlling Simuliid distribution in streams were the mineral substrates and the biotic microhabitats. Blackfly larvae and pupae were sampled at each site for 15 min. Statistical analysis was performed by CANOCO^® (Ter Braak & Smilauer (1997). CANOCO Version 4.5. Biometrics Plant Research International, Wageningen, The Netherlands) using redundancy analysis (RDA).

Our results show a stream type-specific composition of the Blackfly fauna. The analyses reveal sensitivity of the Simuliid species to morphological degradation, which is indicated by the shift in the longitudinal zonation of the Simuliid communities. Especially, Prosimulium hirtipes (Fries, 1824) and Simulium argyreatum Meigen, 1838 are typical representatives of headwater streams. While they seem to indicate undisturbed conditions of this stream type and react sensitively to the degradation of stream habitats, Simulium ornatum Meigen, 1818 and Simulium equinum (Linnaeus, 1758) are more tolerant to stream degradation.

On catchment scale, ‘% natural forest’ and ‘% agricultural land use’ illustrate the degree of degradation of the two selected stream types. ‘Average stream depth’ and ‘relation riffles/pools’ account for hydromorphological degradation reflected by Simuliid species on the smaller reach scale. The analysis of habitat quality revealed that ‘% woody debris’ represents an important parameter of morphological degradation reflected by the Blackfly community.     

Eutrophication and some European waters of restricted exchange

Regions of Restricted Exchange (RREs) are an important feature of the European coastline. They are historically preferred sites for human settlement and aquaculture and their ecosystems, and consequent human use, may be at risk from eutrophication. The OAERRE project (EVK3-CT1999-0002) concerns ‘Oceanographic Applications to Eutrophication in Regions of Restricted Exchange’. It began in July 2000, and studies six sites. Four of these sites are fjords: Kongsfjorden (west coast of Spitzbergen); Gullmaren (Skagerrak coast of Sweden); Himmerfjärden (Baltic coast of Sweden); and the Firth of Clyde (west coast of Scotland). Two are bays sheltered by sand bars: Golfe de Fos (French Mediterranean); and Ria Formosa (Portuguese Algarve). Together they exemplify a range of hydrographic and enrichment conditions. The project aims to understand the physical, biogeochemical and biological processes, and their interactions, that determine the trophic status of these coastal marine RRE through the development of simple screening models to define, predict and assess eutrophication. This paper introduces the sites and describes the component parts of a basic screening model and its application to each site using historical data. The model forms the starting point for the OAERRE project and views an RRE as a well-mixed box, exchanging with the sea at a daily rate E determined by physical processes, and converting nutrient to phytoplankton chlorophyll at a fixed yield q. It thus uses nutrient levels to estimate maximum biomass; these preliminary results are discussed in relation to objective criteria used to assess trophic status. The influence of factors such as grazing and vertical mixing on key parameters in the screening model are further studied using simulations of a complex ‘research’ model for the Firth of Clyde. The future development of screening models in general and within OAERRE in particular is discussed. In addition, the paper looks ahead with a broad discussion of progress in the scientific understanding of eutrophication and the legal and socio-economic issues that need to be taken into account in managing the trophic status of RREs.     

A new Eocene catostomid (Teleostei: Cypriniformes) from northeastern China and early divergence of Catostomidae

A new Eocene catostomid genus and species, Plesiomyxocyprinus arratiae, is described from Huadian, Jilin Province, northeastern China. The materials include a nearly complete skeleton, dozens of pharyngeal bones with teeth, and a number of disarticulated bones. The new articulated specimen is large-sized and deep-bodied, with an estimated standard length of ca. 300 mm and body depth of 156 mm or about half of its standard length. The assignment of the fish to the Catostomidae is based on its falciform pharyngeal bone with one row of numerous (more than 50) compressed teeth, and the bone is much smaller than in cyprinids, especially given the size of the fish. The new fish distinguishes itself from all known catostomids (both extinct and extant) in its long anal fin with four unbranched and 17–18 branched rays, and its extremely short caudal peduncle that is only about one fourth of its depth. Plesiomyxocyprinus arratiae resembles the Eocene-Oligocene transpacific-distributed Amyzon in many general skeletal characters. However, it shows a few characters uniquely shared with the Recent catostomid Myxocyprinus asiaticus. Those include a very long dorsal fin with about 50 branched fin rays, the end of dorsal fin rays being close to the caudal fin base, and anal rays stretching posteriorly beyond the base of caudal fin. It is the first fossil catostomid that shows a close relationship to the endemic Myxocyprinus now living in the Yangtze River and Minjiang River, China. The discovery of Plesiomyxocyprinus arratiae, along with two previously described possible catostomid genera Jianghanichthys and Vasnetzovia, may indicate that the divergence of the Catostomidae started much earlier, in the middle Eocene or earlier, on the western side of the Pacific than on its eastern side. Supported by National Natural Science Foundation of China (Grant No. 40432003, J0630965), Cypriniformes Tree of Life under the U.S. National Science Foundation to R. Mayden (Grant No. EF0431326), and Ecocarp Project (European Commission, INCO-DEV Programme) (Grant No. ICA4-CT-2001- 10024)     

Total mercury in terrestrial systems (air-soil-plant-water) at the mining region of San Joaquín,Queretaro, Mexico

Deficient management of cinnabar mining left the San Joaquín region with high concentrations of mercury in its soils (2.4 – 4164 mg kg^-1). Numerous cinnabar mines have contributed to the dispersion of mercury into agricultural (0.5 –314 mg kg^-1) and forest (0.2 – 69 mg kg^-1) soils. Sediments are a natural means of transportation for mercury, causing its spreading, especially in areas near mine entrances (0.6 – 687 mg kg^-1). The nearness of maize crops to mines favors mercury accumulation in the different plant structures, such as roots, stems, leaves, and grain (0.04 – 8.2 mg kg^-1); these being related to mercury volatilization and accumulation in soils. Mercury vapor present in the settlements could indicate a constant volatilization from lands and soils (22 – 153 ng m^-3). The mercury levels found in the soils, in maize grain, and in the air resulted greater than the standards reported by the Official Mexican Norm (NOM) and the World Health Organization (WHO). Mercury in rainwater is due mainly to the presence of suspended atmospheric particles, later deposited on the surface (1.5 – 339 μg |^-1). Mercury dissolution was found in the drinking water (10 – 170 ng |^-1), with concentrations below those established by the NOM and the WHO. The contamination existing in the San Joaquín region does not reach the levels of the world’s greatest mercury producers: Almaden (Spain) and Idrija (Slovenia). It is, however, like that found in other important second degree world producers such as Guizhou (China). The population of San Joaquín, as well as its surrounding environment, are constantly exposed to mercury contamination, thus making a long term monitoring necessary to determine its effects, especially to people.     

Mercury distribution in Douro estuary (Portugal)

Determinations of dissolved reactive and total dissolved mercury, particulate and sedimentary mercury, dissolved organic carbon (DOC), particulate organic carbon (POC) and suspended particulate matter (SPM) have been made in the estuary of river Douro, in northern Portugal. The estuary was stratified by salinity along most of its length, it had low concentrations of SPM, typically <20 mg dm(-3), and concentrations of DOC in the range <1.0-1.8 mg dm(-3). The surface waters had a maximum dissolved concentration of reactive mercury of about 10 ng dm(-3), whereas for the more saline bottom waters it was about 65 ng dm(-3). The surface waters had maximum concentrations of total suspended particulate mercury of approximately 7 microg g(-1) and the bottom waters were always <1 microg g(-1). Concentrations of mercury in sediments was low and in the range from 0.06 to 0.18 microg g(-1). The transport of mercury in surface waters was mainly associated with organic-rich particulate matter, while in bottom waters the dissolved phase transport of mercury is more important. Lower particulate organic matter, formation of chlorocomplexes in more saline waters and eventually the presence of colloids appear to explain the difference of mercury partitioning in Douro estuarine waters.     

Toxicities of antifouling biocide Irgarol 1051 and its major degraded product to marine primary producers

Irgarol 1051 (2-methythiol-4-tert-butylamino-6-cyclopropylamino-s-triazine) is an algaecide commonly used in antifouling paints. It undergoes photodegradation which yields M1 (2-methylthio-4-tert-butylamino-6-amino-s-triazine) as its major and most stable degradant. Elevated levels of both Irgarol and M1 have been detected in coastal waters worldwide; however, ecotoxicity effects of M1 to various marine autotrophs such as cyanobacteria are still largely unknown. This study firstly examined and compared the 96 h toxicities of Irgarol and M1 to the cyanobacterium Chroococcus minor and two marine diatom species, Skeletonema costatum and Thalassiosira pseudonana. Our results suggested that Irgarol was consistently more toxic to all of the three species than M1 (96 h EC50 values: C. minor, 7.71 microug L(-1) Irgarol vs. > 200 microg L(-1) M1; S. costatum, 0.29 microg L(-1) Irgarol vs. 11.32 microg L(-1)M1; and T. pseudonana, 0.41 microg L(-1) Irgarol vs. 16.50 microg L(-1)M1). Secondly, we conducted a meta-analysis of currently available data on toxicities of Irgarol and M1 to both freshwater and marine primary producers based on species sensitivity distributions (SSDs). Interestingly, freshwater autotrophs are more sensitive to Irgarol than their marine counterparts. For marine autotrophs, microalgae are generally more sensitive to Irgarol than macroalgae and cyanobacteria. With very limited available data on M1 (i.e. five species), M1 might be less toxic than Irgarol; nonetheless this finding warrants further confirmation with additional data on other autotrophic species.     

Respondent uncertainty and reliability in contingent valuation—A case of the Dalai Lake protected area

This paper focuses on examining the validity and reliability of the contingent valuation method (CVM) in the case of evaluation of ecotourism resources in Dalai Lake Protected Area (DLPA). The CVM is a primary method to evaluate resource value, however, it is subject to criticism because of its bias. Many factors result in CVM evaluation bias, which decreases validity of the method. An important factor is that the CVM assumes that all of the respondents are certain about their answers. In fact, many respondents are not very certain about their answer, therefore, identifying the respondents’ certainty level and exploring a proper method to address it are critical issues. There are three objectives in the present study: first, to identify the uncertainty treatment that can more significantly reduce the hypothetical bias, resulting in a greater closeness to the Willingness to Pay (WTP) in a real setting; second, to identify the determinants of the respondents’ uncertainty; and third, to test the reliability of the respondents’ payment in a real setting and a hypothetical setting by comparison with our study in 2010. One thousand three hundred fifty questionnaires about the DLPA were distributed in face-to-face interviews targeted at Hulunbeier inhabitants, and 1142 questionnaires were collected. The primary main conclusions are summarized as follows. First, the treatment of a real setting reduced the hypothetical bias, increased the WTP accuracy and had better reliability compared with the hypothetical setting. Second, income, impact (the impact of resource development on DLPA) and environmental attitudes had a significant positive influence on the respondents’ certainty in hypothetical and real settings, and income and impact of resource development on DLPA had more significant influences in the real setting than in the hypothetical setting. Third, the computation of the mean WTP showed that the exposure model and the asymmetric uncertainty model YES6 fell within the 95% interval around the point estimated from the real model. The WTP of the SUM (symmetric uncertainty model) was much higher than that of the real model and the hypothetical model.     

Water/rock interaction in the Kızılcahamam Geothermal Field, Galatian Volcanic Province (Turkey): a modelling study of a geothermal system for reinjection well locations

The Kızılcahamam geothermal field is emplaced in Tertiary-aged volcanic units 70 km NW of Ankara (Turkey). Data for this low-temperature (74–86°C) geothermal field regarding the fracture zone system were obtained from surface manifestations (hot springs, alteration zones), five exploration wells (MTA-2, -3, -4, -5, -6) and two production wells (KHD-1, MTA-1). The Kızılcahamam reservoir developed along the Kızılcahamam fault zone and so the production wells (180–1556 m) effect each other due to their limited separation. Meteoric water enters from a recharge area NE of Kızılcahamam, circulates and gains heat through the fault zone, and discharges to the surface.Detailed petrographical studies have been carried out with samples taken from surface rocks, cores, and cuttings from three wells (KHD-1, MTA-2, -3). X-ray diffraction techniques were also used in the present study. Kaolinite and montmorillonite zones were identified at outcrop samples. Chloritization, clay mineralization, sericitization and carbonization were determined in the ground mass of samples from wells. The observed alteration mineral assemblages indicate that Kızılcahamam geothermal system has been cooling since the alteration minerals formed.The exploration well MTA-3 seems to be more suitable for a reinjection well than the other wells (MTA-2, -6), even if the cost of surface piping to transport the waste water to MTA-3 is higher than to another well (MTA-6).     

Recharge estimates for various land uses in the Guarani Aquifer System outcrop area

ABSTRACT

The Guarani Aquifer System (GAS) is a subsurface reservoir that contains the largest volume of fresh groundwater in South America. Despite the relevance of the GAS, a lack of attention has been paid to land use effects on its recharge. We present the most detailed long-term (2004–2011) results of land-use effects on recharge in an outcrop area of the GAS. Water table fluctuations (WTFs) were measured at 11 monitoring wells, which are distributed between different land uses (i.e. eucalyptus, sugarcane, citrus and grassland). Recharge was estimated using a point-scale method (WTF) for each monitored well. The annual recharge estimates for different land uses are eucalyptus forest (135 mm year^-1), sugarcane (248 mm year^-1), citrus areas (296 mm year^-1) and grassland (401 mm year^-1). The results indicate that the evapotranspiration seems to be a key parameter in the assessment of recharge in the study area.     

Metallothionein in intestine of red mullet, Mullus barbatus as a biomarker of copper exposure in the coastal marine areas

Dietary metal exposure in fish is studied in native red mullet, Mullus barbatus, involving metallothionein (MT) response in the intestine, a primary site for metal uptake. Intestinal MTs (microg mg(-1) proteins), total proteins (mgml(-1) S50), Cu, Zn, Fe, Mn and Cd (microg ml(-1) S50) are considered as indicators of metal exposure in specimens from near-shore and off-shore locations of the Eastern Adriatic Sea. Metallothionein and Cu are statistically significantly higher in 68 specimens from the Kastela Bay, a near-shore area (30.0+/-6.71microg MT mg(-1) proteins; 0.30+/-0.08microg Cu ml(-1) S50), in comparison to 51 specimens from off-shore area, Solta Island (25.9+/-4.21microg MT mg(-1) proteins; 0.26+/-0.09microg Cu ml(-1) S50). Cytosolic Cu, a MT inducer, is significantly correlated (p<0.05, r=0.44) with MT in specimens from the Kastela Bay, independent on fish biometry. Our study indicates that cytosolic intestinal Cu and MT represent a measurable biochemical response in red mullet specimens, dwelling above sediments contaminated with Cu.     

Organochlorine compounds and aliphatic hydrocarbons in Pacific walrus blubber

Blubber samples were collected from 8 male and 19 female Pacific walrus (Odobenus rosmarus divergens) taken during a 1991 joint USA/USSR cruise traveling widely through the Bering Sea. Dieldrin was found at a level similar to that reported 10 years earlier; oxychlordane was found at a slightly higher concentration than reported previously (Taylor et aL, 1989). Heptachlor epoxide was detected for the first time and found at a low concentration. An initial testing for alpha-, beta- and gamma-HCH detected concentrations similar to those in other Bering Sea pinnipeds. Mean summation of PCB was 0.45 microg g(-1) wet weight in males and 0.16 microg g(-1) in females; only one sample was > 1 microg g(-1). Traces of aliphatic hydrocarbons were detected in all sampled animals, only pristane (x = 0.48 microg g(-1)) was found in concentrations > 1 microg g(-1). Small sample sizes, a lack of samples from immature animals, and uniformly low concentrations of contaminants precluded meaningful analysis of age-related effects and regional differences.     

拟建哈—佳铁路工程地震安评抗断参数研究——以依兰-伊通断裂为例

在对依兰-伊通断裂通河段晚第四纪活动参数获取的基础上,对该断裂的几何学特征和晚第四纪活动特征及有关拟建铁路地震安全的设防参数进行了探讨.研究结果表明,依兰-伊通断裂通河段总体走向N30-40°E,由3条断层组成,中支最新活动时代为全新世.该断裂未来发生7级地震时的水平位移量可能达2.2m左右,垂直位错量约为1.1m左右.该断裂全新世以来的活动性质表现为右旋走滑为主,兼具逆冲活动.根据探槽结果和野外地震地质调查得出断裂垂直位错约（1.0±0.2）m,右旋位错量约为（2.7±0.1）m,（1730-30）aB.P.以来的垂直滑动速率和水平滑动速率分别约为（0.57±0.11）mrn/a和（1.57±0.06）mm/a.未来100年内若遭遇地震,其最大水平位错量约2.87m,垂直位错值为1.04m.断层影响带宽度约为8m.该研究结果为拟建哈—佳铁路工程可能遭受的断层影响和抗震设防提供了一定的数据基础,也为地震安评中线状工程的抗断评价问题提供了一定的参考.     

Trace metals in fish and invertebrates of three California coastal Wetlands

Metal concentrations were measured in selected fish and invertebrate species from Mugu Lagoon, Malibu Lagoon and Ballona Wetlands in southern California in order to assess the extent of metal contamination in these three wetlands. Ranges of element concentrations (in microgram/g) found in biota were: Zn 12-650; Cu 1.9-440; Ni < 1-37; Cr < 1-55; Pb < 0.5-6.8; As < 1-8.5; Se < 1-3.8; Cd < 0.2-0.90; and Ag < 0.3-5.9. Relative to previous studies of California biota, the highest metal concentrations found were for chromium and nickel. The highest levels were in one of the two bottom-dwelling fish (juvenile Leptocottus armatus) (55 micrograms/g) and the two water-column fish sampled (Fundulus parvipinnis and Atherinops affinis) (30 and 24 micrograms/g). At Ballona Lagoon, elevated levels of copper and silver were found in the bivalve Tagelus californianus (440 and 5.9 micrograms/g). Chromium and nickel appeared to be most persistent in fish from Mugu (4.6-55 and 2.6-37 micrograms/g), the most northern site and an active military base, and Ballona (< 1-30 and < 1-16 micrograms/g), believed to be the most metal-contaminated site. Compared to previously measured metal concentrations in species of California coastal waters, these regions revealed higher levels of chromium, nickel, silver, arsenic, zinc, copper and, to a lesser extent, cadmium and selenium. Chromium and silver were present at high enough levels at all three sites to be considered environmental health hazards.     

GNSS影像及其时空特征初探

地壳垂向形变在连续空间和时间域内呈现显著特征,探求其时空变化特征有助于理解地球物理过程,为研究地球内部相互作用机制提供支持.本文使用美国西部地区PBO与中国大陆CMONOC两个GNSS网测站的坐标时间序列,通过基于中位数并顾及年际差异的非参数方法(MIDAS方法)估计测站的速度与不确定性;建立空间结构函数(SSF)并确定区域内各测站间的相对关系;以此为基础,构建顾及空间结构的滤波器(MSF)以剔除粗差,增强区域共性;最后,基于MSF与图像处理技术对速度场进行空间加密,生成了研究区域空间内连续的地壳垂向形变图,即区域GNSS影像.随后,针对两个研究区域,分别利用MSF验证实验与棋盘格检测验证了GNSS成像方法的合理性及生成GNSS影像的可靠性;并通过对比使用顾及空间结构滤波前后的各测站速度与不确定性生成的GNSS影像,验证了顾及空间结构的滤波方法在GNSS影像生成中的必要性,并分析了其中存在过度平滑与规则圆弧状突变边缘的问题,讨论了可能的解决方案.最终,将两区域GNSS影像结果与已有的大地测量学及地球动力学结果进行了对比,发现美国西部地区的GNSS影像正确反映出了海岸山脉以峰值速度为2mm·a~(-1),内华达山脉以峰值速度为3mm·a~(-1),以及赫布根湖地区以峰值速度为1.5mm·a~(-1)隆升;洛杉矶地区(峰值速度为-2.5mm·a~(-1)),维多利亚河及其河谷地区(速度为-1.5mm·a~(-1)),以及蛇河平原东部、蒙大拿州西南部(速度为-1mm·a~(-1)左右)的沉降运动;中国大陆的GNSS影像同样反映出喜马拉雅山脉与青藏高原南部(速度呈现为1.0mm·a~(-1)),华北地块与天山地块(速度为1.5mm·a~(-1)与0.3～0.6mm·a~(-1))等典型区域的隆升;长江下游地区以苏锡常地区(速度为-2.1mm·a~(-1))为中心,向外速度逐渐减小的沉降运动,以及华南地区(速度呈现为-0.6～-1.5mm·a~(-1))、东北地区(速度呈现为-0.6～-1.5mm·a~(-1))、塔里木盆地(速度呈现为-1.2mm·a~(-1))等区域的沉降运动.因此,本文认为GNSS影像具有合理性与正确性,有助于地壳垂向形变的整体时空分布特征研究.