Polycyclic aromatic hydrocarbons in coastal sediments of southwest Taiwan: An appraisal of diagnostic ratios in source recognition

Fifty-seven surface sediment samples were collected from the coast of southwest Taiwan and analyzed for polycyclic aromatic hydrocarbons (PAHs). Concentrations of total PAHs (28 PAH compounds) ranged from 15 to 907 ng g⁻¹ dry weight. Diagnostic ratios showed that PAHs in the sediments of the Gaoping estuary were predominantly of petroleum origin, whereas sediments from the Kaohsiung coast contained principally combustion-derived PAHs. Principal component analysis indicated that emissions from automobiles and coal burning were the main sources of combustion-derived PAHs. The relatively high ratios of perylene/penta-aromatic PAH isomers in sediments from the Tainan coast and some off-shore stations on the Kaohsiung coast suggest a significant diagenetic PAH contribution. The study shows that certain diagnostic ratios are useful and sensitive in delineating the distribution of PAHs from specific sources in southwest Taiwan. The phenanthrene/anthracene ratio is a better indicator than the methylphenanthrenes/phenanthrene ratio for tracing petrogenic PAHs, and the benzo(a)anthracene/chrysene and indeno(1,2,3-c,d)pyrene/benzo(g,h,i)perylene ratios are more specific than the benzo(a)pyrene/benzo(e)pyrene and benzo(b)fluoranthcene/benzo(k)fluoranthcene ratios in distinguishing PAHs from various pyrogenic sources.     

Comparison of alternative upscaled model formulations for simulating DNAPL source dissolution and biodecay

This paper compares the performance of analytical and numerical approaches for modeling DNAPL dissolution with biodecay. A solution derived from a 1-D advective transport formulation (“Parker” model) is shown to agree very closely with high resolution numerical solutions. A simple lumped source mass balance solution in which with decay is assumed proportional to DNAPL mass (“Falta1” model) over- or underpredicts aqueous phase biodecay depending on the magnitude of the exponential factor governing the relationship between dissolution rate and DNAPL mass. A modification of the Falta model that assumes decay proportional to the source exit concentration is capable of accurately simulating source behavior with strong aqueous phase biodecay if model parameters are appropriately selected or calibrated (“Falta2” model). However, parameters in the lumped models exhibit complex interdependencies that cannot be quantified without consideration of transport processes within the source zone. Combining the Falta2 solution with relationships derived from the Parker model was found to resolve these limitations and track the numerical model results. A method is presented to generalize the analytical solutions to enable simulation of partial mass removal with changes in source parameters over time due to various remedial actions. The algorithm is verified by comparison with numerical simulation results. An example application is presented that demonstrates the interactions of partial mass removal, enhanced biodecay, enhanced mass transfer and source zone flow reduction applied at various time periods on contaminant flux reduction. Increasing errors that arise in numerical solutions with coarse discretization and high decay rates are shown to be controlled by using an adjusted decay coefficient derived from the Parker analytical solution.     

A seismological overview of long-period ground motion

Long-period ground motion has become an increasingly important consideration because of the recent rapid increase in the number of large-scale structures, such as high-rise buildings and oil storage tanks. Large subduction-zone earthquakes and moderate to large crustal earthquakes can generate far-source long-period ground motions in distant sedimentary basins with the help of path effects. Near-fault long-period ground motions are generated, for the most part, by the source effects of forward rupture directivity. Far-source long-period ground motions consist primarily of surface waves with longer durations than near-fault long-period ground motions. They were first recognized in the seismograms of the 1968 Tokachi-oki and 1966 Parkfield earthquakes, and their identification has been applied to the 1964 Niigata earthquake and earlier earthquakes. Even if there is no seismogram, we can identify far-source long-period ground motions through the investigation of tank damage by liquid sloshing.     

汶川8.0级地震前紫坪铺水库小震活动及震源参数研究

利用区域和水库地震台网记录的数字地震波资料,研究了汶川8.0级地震前的2000年1月1日～2008年5月11日紫坪铺水库区的地震活动及震源参数.结果表明:汶川8.0级地震前该区域小震活动持续,地震活动水平为M_L2～3,自2005年9月紫坪铺水库蓄水以来,特别是2008年2月地震活动明显增加,但从长期的序列来看,地震活动频次和强度仍在区域地震活动的正常起伏范围.汶川8.0级地震发生前近3个月内,中小地震的视应力差值△σ_(app)普遍大于0,计算得到的视应力总大于其拟合值.汶川8.0级地震震中位于低应力区,而相对较高的视应力差值分布在震中以东地区,局部出现视应力增加的现象.     

Analysis of tracer tomography using temporal moments of tracer breakthrough curves

Hydraulic/partitioning tracer tomography (HPTT) was recently developed by Yeh and Zhu [Yeh T-CJ, Zhu J. Hydraulic/partitioning tracer tomography for characterization of dense nonaqueous phase liquid source zones, Water Resour Res 2007;43:W06435. doi:10.1029/2006WR004877.] for estimating spatial distribution of dense nonaqueous phase liquids (DNAPLs) in the subsurface. Since discrete tracer concentration data are directly utilized for the estimation of DNAPLs, this approach solves the hyperbolic convection–dispersion equation. Solution to the convection–dispersion equation however demands fine temporal and spatial discretization, resulting in high computational cost for an HPTT analysis. In this work, we use temporal moments of tracer breakthrough curves instead of discrete concentration data to estimate DNAPL distribution. This approach solves time independent partial differential equations of the temporal moments, and therefore avoids solving the convection–dispersion equation using a time marching scheme, resulting in a dramatic reduction of computational cost. To reduce numerical oscillations associated with convection dominated transport problems such as in inter-well tracer tests, the approach uses a finite element solver adopting the streamline upwind Petrov–Galerkin method to calculate moments and sensitivities. We test the temporal moment approach through numerical simulations. Comparing the computational costs between utilizing moments and discrete concentrations, we find that temporal moments significantly reduce the computation time. We also find that tracer moment data collected through a tomographic survey alone are able to yield reasonable estimates of hydraulic conductivity, as indicated by a correlation of 0.588 between estimated and true hydraulic conductivity fields in the synthetic case study.     

The influence of dimensionality on simulations of mass recovery from nonuniform dense non-aqueous phase liquid (DNAPL) source zones

The influence of model dimensionality on predictions of mass recovery from dense non-aqueous phase liquid (DNAPL) source zones in nonuniform permeability fields was investigated using a modified version of the modular three-dimensional transport simulator (MT3DMS). Thirty-two initial two- (2D) and three-dimensional (3D) tetrachloroethene–DNAPL source zone architectures, taken from a recent modeling study, were used as initial conditions for this analysis. Commonly employed source zone metrics were analyzed to determine differences between 2D and 3D predictions: (i) down-gradient flux-averaged contaminant concentration, (ii) reductions in contaminant mass flux through a down-gradient boundary, (iii) source zone ganglia-to-pool (GTP) ratio, and (iv) time required to achieve a remediation objective. 3D flux-averaged contaminant concentrations were approximately 3.5 times lower than concentrations simulated in 2D. This difference was attributed to dilution of the contaminant concentrations down gradient of the source zone. Contaminant flux reduction predictions for a given mass recovery were generally 5% higher in 3D simulations than in 2D simulations. The GTP ratio declined over time as mass was recovered in both 2D and 3D simulations. Although the source longevity (i.e., time required to achieve 99.99% mass recovery) differed between individual 2D and 3D realizations, the mean source longevity for the 2D and 3D simulation ensembles was within 2%. 2D simulations tended to over-predict the time required to achieve lower mass recovery levels (e.g. 50% mass recovery) due to a smaller contaminated area exposed to uncontaminated water. These findings suggest that ensemble averages of 2D numerical simulations of DNAPL migration, entrapment, dissolution, and mass recovery in statistically homogenous, nonuniform media may provide reasonable approximations to average behavior obtained using simulations conducted in fully three-dimensional domains.     

Ground motions observed during the 15 August 2007 Pisco, Peru, earthquake

A M_w 7.9 earthquake event occurred on 15 August 2007 off the coast of central Peru, 60 km west of the city of Pisco. This event is associated with subduction processes at the interface of the Nazca and South American plates, and was characterised by a complex source mechanism involving rupture on two main asperities, with unilateral rupture propagation to the southeast. The rupture process is clearly reflected in the ground motions recorded during this event, which include two separate episodes of strong shaking. The event triggered 18 accelerographic stations; the recordings are examined in terms of their characteristics and compared to the predictions of ground-motion prediction equations for subduction environments, using the maximum-likelihood-based method of Scherbaum et al. (Bull Seismol Soc Am 94(6):2164–2185, 2004). Additionally, macroseismic observations and damage patterns are examined and discussed in the light of local construction practices, drawing on field observations gathered during the post-earthquake reconnaissance missions.     

Tmax of asphaltenes: a parameter for oil maturity assessment

A new maturity parameter determined on both oil and bitumen samples, the asphaltene T_max, is proposed and discussed. This parameter could be very useful to address the maturity of the source rock. The asphaltene T_max is measured by programmed Rock-Eval pyrolysis, using a modified temperature program. Some phases of the experimental procedure, such as the asphaltene preparation and the Rock-Eval measurement substratum choice, are crucial in order to achieve reliable data. Laboratory simulations were carried out in order to assess the possible effects of both primary and secondary migration on asphaltene T_maxin the expelled oil: the original value of the asphaltene T_max in the bitumen is not substantially modified and it is very close to that measured on kerogen. Examples of the determination of asphaltene T_max on many samples, collected from different areas and with different organic matter composition, are given. Results show that T_max values from oil asphaltenes are reasonable indicators of source rock maturity.     

Simulation of near source ground motion and its characteristics

It is fact that the severe ground motions of shear waves have a strong effect on the dynamic behavior of buildings and civil structures. We simulate near source strong motions of a pure shear wave and synthesize small motions, using the parameters based on the recorded accelerograms at the site that is regarded as a base rock in the Osaka basin, Japan. By making use of a stochastic technique, we can easily introduce higher frequency contents in the motions and apply the technique to the synthesis of small waves regarding as a green function. We also introduce to the analysis the useful relationships among the time duration T_d, the seismic moment M₀, the corner frequency f_c and the high cutoff frequency f_max which were regressed by a simple representation scheme. Considering two active faults that may affect severe damage on buildings and civil structures, we try to predict strong ground motions in Osaka basin and show the characteristics of them.     

震源硬化模型的理论、实验及观测事实依据

从前兆异常共性特征、有关理论分析及岩石实验结果等方面,论证了震源硬化模型的科学性和合理性。