Focal mechanisms of recent earthquakes in the Southern Korean Peninsula

We evaluate the stress field in and around the southern Korean Peninsula with focal mechanism solutions, using the data collected from 71 earthquakes ( M_L = 1.9–5.2) between 1999 and 2004. For this, the hypocentres were relocated and well-constrained fault plane solutions were obtained from the data set of 1270 clear P -wave polarities and 46 SH / P amplitude ratios. The focal mechanism solutions indicate that the prevailing faulting types in South Korea are strike-slip-dominant-oblique-slip faultings with minor reverse-slip component. The maximum principal stresses (σ₁) estimated from fault-slip inversion analysis of the focal mechanism solutions show a similar orientation with E–W trend (269°–275°) and low-angle plunge (10°–25°) for all tectonic provinces in South Korea, consistent with the E–W trending maximum horizontal stress (σ_Hmax) of the Amurian microplate reported from in situ stress measurements and earthquake focal mechanisms. The directions of the intermediate (σ₂) and minimum (σ₃) principal stresses of the Gyeongsang Basin are, however, about 90 deg off from those of the other tectonic provinces on a common σ₂–σ₃ plane, suggesting a permutation of σ₂ and σ₃. Our results incorporated with those from the kinematic studies of the Quaternary faults imply that NNW- to NE-striking faults (dextral strike-slip or oblique-slip with a reverse-slip component) are highly likely to generate earthquakes in South Korea.     

An analytical solution to multi-component NAPL dissolution equations

This note presents a novel method for determining the changing composition of a multi-component NAPL body dissolving into moving groundwater, and the consequent changes in the aqueous phase solute concentrations in the surrounding pore water. A canonical system of coupled non-linear governing equations is derived which is suitable for representation of both pooled and residual configurations, and this is solved. Whereas previous authors have handled such problems numerically, it is shown that these governing equations succumb to analytical solution. By a suitable substitution, the equations become decoupled, and the problem collapses to a single first-order equation. The final result is expressed implicitly, with time as a function of the number of moles of the least soluble component, m₁. The number of moles of each other component is expressed explicitly in terms of m₁. It is shown that the time-m₁ relationship has a well behaved inverse. An example is given in which the analytic solution is verified against traditional finite difference analysis, and its computational efficiency is shown.     

Effective macroscopic transport parameters between parallel plates with constant concentration boundaries

A macroscopic transport model is developed, following the Taylor shear dispersion analysis procedure, for a 2D laminar shear flow between parallel plates possessing a constant specified concentration. This idealized geometry models flow with contaminant dissolution at pore-scale in a contaminant source zone and flow in a rock fracture with dissolving walls. We upscale a macroscopic transient transport model with effective transport coefficients of mean velocity, macroscopic dispersion, and first-order mass transfer rate. To validate the macroscopic model the mean concentration, covariance, and wall concentration gradient are compared to the results of numerical simulations of the advection–diffusion equation and the Graetz solution. Results indicate that in the presence of local-scale variations and constant concentration boundaries, the upscaled mean velocity and macrodispersion coefficient differ from those of the Taylor–Aris dispersion, and the mass transfer flux described by the first-order mass transfer model is larger than the diffusive mass flux from the constant wall. In addition, the upscaled first-order mass transfer coefficient in the macroscopic model depends only on the plate gap and diffusion coefficient. Therefore, the upscaled first-order mass transfer coefficient is independent of the mean velocity and travel distance, leading to a constant pore-scale Sherwood number of 12. By contrast, the effective Sherwood number determined by the diffusive mass flux is a function of the Peclet number for small Peclet number, and approaches a constant of 10.3 for large Peclet number.     

天津及邻区1999～2000年地震震源机制研究

文中使用的是天津市、河北省地震台的 P,Pn波初动方向资料 ,研究了 1 999～ 2 0 0 0年发生在天津及邻区的 4个地震的震源机制。其中 1 999年 8月 1 1日河北省大城县 ML 3.1地震发生在大城断裂上 ,震源机制解的节面 (走向 2 0 8°,倾角 54°,滑动角 1 2°)与大城断裂的走向一致。2 0 0 0年 1月 2 6日蓟县 ML 3.3地震发生在蓟县 -邦均断裂上 ,震源机制解的节面 (走向 2 85°,倾角 77°,滑动角 1 38°)与蓟县 -邦均断裂的走向一致。表明这两个断层在活动 ,这两个地震发生在多年来 (有感 )地震不多的地区 ,并且是发生在首都圈中部 ,应该引起注意     

Zinc-rich Pyrite from the TAG Active Mound, the TAG Hydrotherma Field, Mid-Atlantic Ridge

Abstract: Pyrite rich in Zn, up to 3.1 wt%, was found in the TAG active mound of the TAG hydrothermal field, the slow-spreading Mid-Atlantic Ridge at 26°08'N and 44°49'W. The Zn-rich pyrite is characterized by an optical homogeneity, a homogeneous distribution of Zn in the back-scattered electron images, both at a magnification of about 500, a negative correlation between Fe and Zn contents of the pyrite and a rather small unit cell edge (a₀ = 5.4117 ± 0.0008Å), strongly indicating that the detected Zn is present in the pyrite in solid solution. Such Zn concentrations are observed exclusively in dendritic pyrite, suggesting that the Znrich pyrite grew from hydrothermal fluids of a high degree of supersaturation due to quenching on the seafloor.     

Element budgets of pine stands on lignite and pyrite containing mine soils

To analyze the development of pine ecosystems on lignite and pyrite containing mine soils, four pine stands with ages of 3–35 years were investigated in a chronosequence approach. Bulk precipitation, throughfall and soil solution in depths of 20, 40, 70 and 100 cm were studied over a three-year period to determine element fluxes in these forest ecosystems on extreme acidic and saline soils.Element budgets are controlled by the processes induced by pyrite oxidation such as intensive weathering of primary minerals, precipitation and leaching of secondary phases.Despite low water fluxes, element outputs can reach extraordinary high values due to very high concentrations in the soil solution. Although element outputs decrease drastically with stand age, respectively, site age, they exceed those of comparable pine stands on non-mine sites. Nitrogen release from the lignite fraction represents a special characteristic of the examined substrates. Nitrogen losses beneath the rooted zone can be 34 kg N ha⁻¹ yr⁻¹ in 100 cm depth. Element input and output in the examined ecosystems are far from balance. Closed cycling of nutrients seems to be recognizable in the case of potassium.     

Numerical algorithm for the dynamic analysis of base-isolated structures with dry friction

This paper deals with the techniques for solving ordinary differential equations with essential nonlinearity arising from the representation of frictional force by the sign function of the relative velocity. This problem is connected with dynamic behavior studies of base-isolated structures with dry friction devices.Mathematical evidence is presented concerning the inaccuracies of the solution which are connected with the representation of the frictional force by the sign function. An approximation of this sign model is made on the basis of the appropriate linearization of a small part of the frictional force function. A numerical algorithm for the dynamic analysis of base-isolated structures with dry friction is proposed. Numerical tests are carried out for establishing the exactness of the proposed techniques. The obtained results prove the accuracy of the techniques.     

An approximate analytical solution for the deposition of heavy particles released from an elevated line source

This work re-examines and further develops an analytical solution for the deposition swath of heavy particles released in the atmosphere from an elevated source over uniform terrain, correcting the particle diffusivity for the crossing trajectory effect. The revised (approximate) analytical solution proves to be accurate within 20% over a wide range of micrometeorological conditions and particle size, despite its neglect of the turbulence component of the deposition flux. It compares very satisfactorily with experimental data and with the simulations of a Lagrangian stochastic model, provided the variable U(H)/w _g ≤7 (ratio of the mean horizontal wind speed at source height to the particle settling velocity). In this domain of validity, simple formulae relating the statistics of the deposition swath to U(H)/w _g are derived.     

Thermoelastic stress due to a rectangular heat source in a semi-infinite medium. Presentation of an analytical solution

The thermoelastic response due to a time-dependent rectangular heat source in a semi-infinite medium is analyzed. The problem originates from studies of nuclear waste repositories in rock. Canisters containing heat-emitting nuclear waste are deposited over a large rectangular area deep below the ground surface. The solution for a time-dependent heat source is obtained from the corresponding instantaneous heat source by superposition. The thermoelastic problem for the instantaneous rectangular heat source in a infinite surrounding is solved exactly. An important step is the introduction of so-called quadrantal heat sources. The solution for the rectangle is obtained from four quadrantal solutions. The solution for the quadrantal heat source depends on the three dimelasionless coordinates only. Time occurs in the scale factors only. The condition of zero normal and shear stresses at the ground surface is fulfilled by using a mirror heat source and a boundary solution. The boundary solution accounts for the residual normal stress at the ground surface. Using a Hertzian potential, a surprisingly simple solution is obtained. The final analytical solution is quite tractable considering the complexity of the initial problem. The solution may be used to test numerical models for coupled thermoelastic processes. It may also be used in more detailed numerical simulations of the process near the heat sources as boundary conditions to account for the three-dimensional global process.