2008年汶川大地震震源机制的时空变化

本文提出了一种基于恒定破裂速度和固定子事件震源时间函数的假定、利用远场地震波形资料获取大地震震源机制的时空变化图像的线性反演方法,并利用这种方法及全球范围内48个台站的长周期波形资料反演建立了2008年汶川M_S8.0地震的震源机制随时间和空间变化的图像.根据这个图像可知,汶川大地震断层的西南端震源机制接近于逆冲,随着破裂向东北方向延伸,震源机制的走滑分量逐渐增大,走滑分量超过逆冲分量的转折点在震中东北大约190 km的位置.为了检验反演方法的有效性和反演结果的可靠性,我们特别设计了一个数值试验对反演结果进行了检验.检验结果表明,我们在本文中提出的反演方法是有效的,关于汶川大地震的反演结果也是可靠的（除长周期信号较弱的一段外).通过比较发现,反演结果与震后野外考察的结果也相当吻合.     

利用区域地震台网、振幅比资料测定小震震源参数

本文提出一种利用区域地震台网垂直向位移记录中和最大振幅的比值,由计算层状介质中点源位错的、综合地震图测定小震震源参数的方法。用人工数据检验该测定方法的可行性;并用此方法测定了海城地震（1975年2月4日,M_s=7.3）前震序列震源参数,与其他方法得到的结果作了对比分析。结果表明,本文提出的方法具有需用记录资料少的优点。     

The influence of urbanization on karst rivers based on nutrient concentration and nitrate dual isotopes: an example from southwestern China

China is experiencing rapid urbanization that has changed the water quality of rivers, especially nutrient loads. In this study, a typical urban river located in a karst area, Chengguan River, was chosen to explore the influence of urbanization on river ecosystems based on nutrient concentration and nitrate isotopes. The results show monthly variability of water chemistry and nutrient concentration. Nutrient concentration in two tributaries and the mainstem showed significant spatial variability, with heavy N and P pollution in one tributary near a suburban area,indicating a response to different levels of urbanization.Measurements of nitrate dual isotopes suggest thatvolatilization, assimilation, nitrification, and denitrification all occur in the polluted river. Water chemistry and nitrate isotopes show that major nitrogen sources included domestic waste and agricultural input, such as chemical fertilizer and manure. The results suggest that urbanization increases nutrient concentrations and accelerates the riverine nitrogen dynamic, and point to the need to manage point sources of sewage effluents to improve the water quality of urban rivers in southwestern China.     

A Screening Procedure for Heavy Metals in Soil Extracts by Alcohol Dehydrogenase and Urease Inhibition Eine Screening-Methode für Schwermetalle in Extrakten aus Bodenproben,basierend auf der Inhibierung der Enzyme Alkoholdehydrogenase und Urease

A screening method for heavy metals in aqueous extracts of soil is presented which is based on inhibition of the enzymes urease and alcohol dehydrogenase. The method is suitable to detect cupric and mercury ions in concentrations below 0.01 mg/L and several other heavy metal ions in 1000 fold higher concentration. It is shown that the test may be used for screening of mercury ion concentrations exceeding 0.03 mg/L in aqueous solution when copper chelators are added to the test system. The usefulness of the presented tests to detect heavy metals eluted from soil was verified with samples from ore mining waste. The concentration of copper, lead, and zinc eluted from these samples to different amount was determined by atomic absorption spectrometry and was in good agreement with the enzyme inhibition data obtained with these samples.     

Identification of contaminant point source in surface waters based on backward location probability density function method

A backward location probability density function (BL-PDF) method capable of identifying location of point sources in surface waters is presented in this paper. The relation of forward location probability density function (FL-PDF) and backward location probability density, based on adjoint analysis, is validated using depth-averaged free-surface flow and mass transport models and several surface water test cases. The solutions of the backward location PDF transport equation agreed well to the forward location PDF computed using the pollutant concentration at the monitoring points. Using this relation and the distribution of the concentration detected at the monitoring points, an effective point source identification method is established. The numerical error of the backward location PDF simulation is found to be sensitive to the irregularity of the computational meshes, diffusivity, and velocity gradients. The performance of identification method is evaluated regarding the random error and number of observed values. In addition to hypothetical cases, a real case was studied to identify the source location where a dye tracer was instantaneously injected into a stream. The study indicated the proposed source identification method is effective, robust, and quite efficient in surface waters; the number of advection–diffusion equations needed to solve is equal to the number of observations.     

The Effect of Fatty Acid Diffusion in Leachate on the Propagation of Concentration Waves in the Process of Municipal Solid Waste Decomposition

A distributed model of municipal solid waste decomposition is proposed. The main model parameters are concentrations of municipal solid waste, volatile fatty acids, and the biomass of methane-producing microorganisms. Numerical analysis of the model is made. Diffusion of fatty acids in leachate, which facilitates the formation of anaerobic conditions, is shown to cause the formation of concentration chemical waves propagating in the space. The area of initiation methane production is found to expand.     

Histological Changes in a Few Tissues of Notopterus notopterus Induced by Sublethal and Lethal Concentrations of Pulp,Paper and Mixed Pulp and Paper Wastes

In the present work, histological changes in the gill, liver, kidney and intestine of Notopterus notopterus, induced by pulp waste, paper waste and mixed pulp and paper waste, have been studied both at sublethal and lethal concentrations. The pulp waste produces much more histological changes both at sublethal and lethal concentrations in the gill, liver, kidney and intestine as compared to paper waste and mixed pulp and paper waste. The changes are much more remarkable in the case of sublethal concentration, as compared to a lethal concentration. In a sublethal concentration, the localized necrosis, atrophy and partial and total liquefaction of epithelial cells, pilaster cells, mucous cells of the gill, hepatic cells of the liver, mesenchymatous cells, renal tubule cells and Bowman's capsule of the kidney, mucosal and submucosal cells and blood vessels of intestine have been observed. In lethal concentration, the changes are much more marked in the gill, as compared to the liver, kidney and intestine. The histological changes, produced in the different types of cells and in the tissue as a whole, can be correlated with the different chemical components present in the pulp waste, paper waste and mixed pulp and paper waste.     

基于解析时间波场外推与波场分解的逆时偏移方法研究

双程波方程逆时深度偏移是复杂介质高精度成像的有效技术,但其结果中通常包含成像方法引起的噪音和假象,一般的滤波方法会破坏成像剖面上的振幅,其中的假象也会给后续地质解释带来困扰.将波场进行方向分解然后实现入射波与反射波的相关成像能够有效地消除这类成像噪音,并提高逆时偏移成像质量.波传播方向的分解通常在频率波数域实现,它会占用大量的存储和计算资源,不便于在沿时间外推的逆时深度偏移中应用.本文提出解析时间波场外推方法,可以在时间外推的每个时间片上实现波传播方向的显式分解,逆时深度偏移中利用分解后的炮检波场进行对应的相关运算,实现成像噪音和成像信号的分离.在模型和实际数据上的测试表明,相比于常规互相关逆时偏移成像结果,本文方法能够有效地消除低频成像噪音和特殊地质构造导致的成像假象.     

Water column monitoring near oil installations in the North Sea 2001-2004

Fisheries have been vital to coastal communities around the North Sea for centuries, but this semi-enclosed sea also receives large amounts of waste. It is therefore important to monitor and control inputs of contaminants into the North Sea. Inputs of effluents from offshore oil and gas production platforms (produced water) in the Norwegian sector have been monitored through an integrated chemical and biological effects programme since 2001. The programme has used caged Atlantic cod and blue mussels. PAH tissue residues in blue mussels and PAH bile metabolites in cod have confirmed exposure to effluents, but there was variation between years. Results for a range of biological effects methods reflected exposure gradients and indicated that exposure levels were low and caused minor environmental impact at the deployment locations. There is a need to develop methods that are sufficiently sensitive to components in produced water at levels found in marine ecosystems.     

Evaluation of shear wave velocity profile using SPT based uphole method

The uphole method is a field seismic test which uses receivers on the ground surface and an underground source. A modified form of the uphole method is introduced in order to obtain efficiently the shear wave velocity (V_S) profile of a site. This method is called the standard penetration test (SPT)-uphole method because it uses the impact energy of the split spoon sampler in the SPT test as a source. Since the SPT-uphole method can be performed simultaneously with the SPT test it is economical and not labor intensive compared to the original uphole methods which use small explosives or a mechanical source. Field testing and interpretation procedures for the proposed method are described. To obtain reliable travel time information of the shear wave, the first peak point of the shear wave using two component geophones is recommended. Through a numerical study using the finite element method (FEM), the procedure of the proposed method was verified. Finally, the SPT-uphole method was performed at several sites, and the field applicability of the proposed method was verified by comparing the V_S profiles determined by the SPT-uphole method with the profiles determined by the downhole, the spectral analysis of surfaces waves (SASW) method and from the SPT-N values.     

Examining the Use of USEPA's Generic Attenuation Factor in Determining Groundwater Screening Levels for Vapor Intrusion

A value of 0.001 is recommended by the United States Environmental Protection Agency (USEPA) for its groundwater‐to‐indoor air Generic Attenuation Factor (GAFG), used in assessing potential vapor intrusion (VI) impacts to indoor air, given measured groundwater concentrations of volatile chemicals of concern (e.g., chlorinated solvents). The GAFG can, in turn, be used for developing groundwater screening levels for VI given target indoor air quality screening levels. In this study, we examine the validity and applicability of the GAFG both for predicting indoor air impacts and for determining groundwater screening levels. This is done using both analysis of published data and screening model calculations. Among the 774 total paired groundwater‐indoor air measurements in the USEPA's VI database (which were used by that agency to generate the GAFG) we found that there are 427 pairs for which a single groundwater measurement or interpolated value was applied to multiple buildings. In one case, up to 73 buildings were associated with a single interpolated groundwater value and in another case up to 15 buildings were associated with a single groundwater measurement (i.e., that the indoor air contaminant concentrations in all of the associated buildings were influenced by the concentration determined at a single point). In more than 70% of the cases (390 of 536 paired measurements in which horizontal building‐monitoring well distance was recorded) the monitoring wells were located more than 30 m (and one up to over 200 m) from the associated buildings. In a few cases, the measurements in the database even improbably implied that soil gas contaminant concentrations increased, rather than decreased, in an upward direction from a contaminant source to a foundation slab. Such observations indicate problematic source characterization within the data set used to generate the GAFG, and some indicate the possibility of a significant influence of a preferential contaminant pathway. While the inherent value of the USEPA database itself is not being questioned here, the above facts raise the very real possibility that the recommended groundwater attenuation factors are being influenced by variables or conditions that have not thus far been fully accounted for. In addition, the predicted groundwater attenuation factors often fall far beyond the upper limits of predictions from mathematical models of VI, ranging from screening models to detailed computational fluid dynamic models. All these models are based on the same fundamental conceptual site model, involving a vadose zone vapor transport pathway starting at an underlying uniform groundwater source and leading to the foundation of a building of concern. According to the analysis presented here, we believe that for scenarios for which such a “traditional” VI pathway is appropriate, 10^?4 is a more appropriately conservative generic groundwater to indoor air attenuation factor than is the EPA‐recommended 10^?3. This is based both on the statistical analysis of USEPA's VI database, as well as the traditional mathematical models of VI. This result has been validated by comparison with results from some well‐documented field studies.     

A nonlinear mixed model for the assessment of traffic noise levels in urban areas

Quantification of traffic noise levels in urban areas can be especially useful for environmental and public health policies. A statistical model is proposed here for this quantification in small urban areas. The model is nonlinear and mixed, with two components: a fixed component for traffic noise emission (source effect) and a random component for assessing traffic noise propagation. The proposed model requires only sample data and has been used to estimate basic noise indicators established in European standards. A composite indicator, the Lden indicator, has been computed and used to draw up a noise ranking of urban areas. The proposed approach is applied to the city of Leganés (Spain).     

Probabilistic-fuzzy health risk modeling

Health risk analysis of multi-pathway exposure to contaminated water involves the use of mechanistic models that include many uncertain and highly variable parameters. Currently, the uncertainties in these models are treated using statistical approaches. However, not all uncertainties in data or model parameters are due to randomness. Other sources of imprecision that may lead to uncertainty include scarce or incomplete data, measurement error, data obtained from expert judgment, or subjective interpretation of available information. These kinds of uncertainties and also the non-random uncertainty cannot be treated solely by statistical methods. In this paper we propose the use of fuzzy set theory together with probability theory to incorporate uncertainties into the health risk analysis. We identify this approach as probabilistic-fuzzy risk assessment (PFRA). Based on the form of available information, fuzzy set theory, probability theory, or a combination of both can be used to incorporate parameter uncertainty and variability into mechanistic risk assessment models. In this study, tap water concentration is used as the source of contamination in the human exposure model. Ingestion, inhalation and dermal contact are considered as multiple exposure pathways. The tap water concentration of the contaminant and cancer potency factors for ingestion, inhalation and dermal contact are treated as fuzzy variables while the remaining model parameters are treated using probability density functions. Combined utilization of fuzzy and random variables produces membership functions of risk to individuals at different fractiles of risk as well as probability distributions of risk for various alpha-cut levels of the membership function. The proposed method provides a robust approach in evaluating human health risk to exposure when there is both uncertainty and variability in model parameters. PFRA allows utilization of certain types of information which have not been used directly in existing risk assessment methods.     

外部扰动重力场的频谱响应质点模型

本文提出一种新的虚拟质点建模方式，旨在实现地球外部重力场的快速精密赋值．其基本思想在于：根据内部场源与外部场的频谱响应关系，视质点系为若干子系组成，让其中每个子系具有与外部场的相应频段分量响应的谱频特性．质点系的频谱响应性质将会大大削弱通常质点模型产生的所谓且Ｒｕｎｇｅ现象：泄漏与增频，从而显著提高逼近精度．文中还详细讨论了频谱响应质点模型的原理和构制方法，并对模型精度进行了分析和估算．     

Evaluation of longitudinal and transverse dispersivities/distance ratios for tracer test in a radially convergent flow field with scale-dependent dispersion

This study presents a novel mathematical model for analysis of non-axisymmetrical solute transport in a radially convergent flow field with scale-dependent dispersion. A two-dimensional, scale-dependent advection–dispersion equation in cylindrical coordinates is derived based on assuming that the longitudinal and transverse dispersivities increase linearly with the distance of the solute transported from its injected source. The Laplace transform finite difference technique is applied to solve the two-dimensional, scale-dependent advection–dispersion equation with variable-dependent coefficients. Concentration contours for different times, breakthrough curves of average concentration over concentric circles with a fixed radial distance, and breakthrough curves of concentration at a fixed observation point obtained using the scale-dependent dispersivity model are compared with those from the constant dispersivity model. The salient features of scale-dependent dispersion are illustrated during the non-axisymmetrical transport from the injection well into extraction well in a convergent flow field. Numerical tests show that the scale-dependent dispersivity model predicts smaller spreading than the constant-dispersivity model near the source. The results also show that the constant dispersivity model can produce breakthrough curves of averaged concentration over concentric circles with the same shape as those from the proposed scale-dependent dispersivity model at observation point near the extraction well. Far from the extracting well, the two models predict concentration contours with significantly different shapes. The breakthrough curves at observation point near the injection well from constant dispersivity model always produce lesser overall transverse dispersion than those from scale-dependent dispersivity model. Erroneous dimensionless transverse/longitudinal dispersivity ratio may result from parametric techniques which assume a constant dispersivity if the dispersion process is characterized by a distance-dependent dispersivity relationship. A curve-fitting method with an example is proposed to evaluate longitudinal and transverse scale-proportional factors of a field with scale-dependent dispersion.     

弹性半空间位错内源的数值解

本文首先将一般剪切位错点源表示成四种基本点位错的组合,断层面倾角和滑动系数仅与组合系数有关.利用将三维动力学问题化为两维的处理方法,将每一基本点源问题化成求解位移函数的轴对称问题,最后将人工透射边界的处理方法同有限差分结合起来,提出一种简便、节省的计算理论地震图的方法.利用这种方法获得了四种基本点位错源的典型近场理论地震图,并初步讨论了近场地面运动的若干特征,表明浅源地震近场地面位移中震相不能分离,且点位错源在近场产生大小相当的水平与垂直位移.本文的方法可用于研究复杂的近场地震波传播问题.      

Probability Theory as Logic: Data Assimilation for Multiple Source Reconstruction

Probability theory as logic (or Bayesian probability theory) is a rational inferential methodology that provides a natural and logically consistent framework for source reconstruction. This methodology fully utilizes the information provided by a limited number of noisy concentration data obtained from a network of sensors and combines it in a consistent manner with the available prior knowledge (mathematical representation of relevant physical laws), hence providing a rigorous basis for the assimilation of this data into models of atmospheric dispersion for the purpose of contaminant source reconstruction. This paper addresses the application of this framework to the reconstruction of contaminant source distributions consisting of an unknown number of localized sources, using concentration measurements obtained from a sensor array. To this purpose, Bayesian probability theory is used to formulate the full joint posterior probability density function for the parameters of the unknown source distribution. A simulated annealing algorithm, applied in conjunction with a reversible-jump Markov chain Monte Carlo technique, is used to draw random samples of source distribution models from the posterior probability density function. The methodology is validated against a real (full-scale) atmospheric dispersion experiment involving a multiple point source release.     

The Antarctic meteorite Yamato 74123 — a new ureilite

The rare gases He, Ne, Ar, Kr and Xe were measured in bulk samples of Yamato 74123. The ³He and ²¹Ne exposure ages are found to be 5.50 Ma and 2.83 Ma, respectively. In addition to the cosmogenic component the samples contain primordial rare gases of the fractionated type in amounts typical of ureilites. In a three-isotope plot neon turns out to be a mixture of planetary neon and cosmogenic neon.The elements Na, Mg, Al, Si, P, S, K, Ca, Cr, Mn, Fe, Co, and Ni have been determined by spark source mass spectrometry in Yamato 74123 and for comparison in the ureilites Haveröand Kenna. The chemical composition as well as the noble gas abundance pattern identify Yamato 74123 as an ureilite.     

Exposure reconstruction using a physiologically based toxicokinetic model with cumulative amount of metabolite in urine: a case study of trichloroethylene inhalation

The use of a physiologically based toxicokinetic (PBTK) model to reconstruct chemical exposure using human biomonitoring data, urinary metabolites in particular, has not been fully explored. In this paper, the trichloroethylene (TCE) exposure dataset by Fisher et al. (Toxicol Appl Pharm 152:339–359, 1998) was reanalyzed to investigate this new approach. By treating exterior chemical exposure as an unknown model parameter, a PBTK model was used to estimate exposure and model parameters by measuring the cumulative amount of trichloroethanol glucuronide (TCOG), a metabolite of TCE, in voided urine and a single blood sample of the study subjects by Markov chain Monte Carlo (MCMC) simulations. An estimated exterior exposure of 0.532 mg/l successfully reconstructed the true inhalation concentration of 0.538 mg/l with a 95% CI (0.441–0.645) mg/l. Based on the simulation results, a feasible urine sample collection period would be 12–16 h after TCE exposure, with blood sampling at the end of the exposure period. Given the known metabolic pathway and exposure duration, the proposed computational procedure provides a simple and reliable method for environmental (occupational) exposure and PBTK model parameter estimation, which is more feasible than repeated blood sampling.     

Robust Regulation of Alkalinity in Highly Uncertain Continuous Anaerobic Digestion Processes

In this paper, a multiple‐input multiple‐output (MIMO) model‐based robust control scheme is proposed for the indirect control of both total alkalinity and the ratio (intermediate alkalinity)/(total alkalinity) by regulating volatile fatty acid concentrations and strong ions concentration, while guaranteeing the so‐called operational stability in anaerobic digestion (AD) processes. The proposed MIMO regulator is an adaptive controller derived from an AD model which incorporates the physicochemical equilibrium of the system as well as the use of a robust interval observer to estimate key process bounds that are used in the computation of the control efforts. Numerical simulations were carried out for a number of operating conditions under the most uncertain scenarios. Results showed that the proposed multivariable control law is able to recover the system stability around a pre‐determined set point in the face of parameter uncertainty and load disturbances.