Population dynamics of earthquakes and mathematical modeling     

Toru Ouchi 《Pure and Applied Geophysics》1993,140(1):15-28

We present a mathematical model that describes temporal variations of earthquakes. This model is represented as $$dn(t)/dt = n(t)\left[ {\alpha - \beta n(t) - \int_{ - \infty }^t {n(s)h(t - s)ds} } \right].$$ Heren(t) shows the numberof earthquakes per unit time in a certain region. α and β are constants. The functionh(t) denotes the hysteresis effect of the earthquake occurrences and can take the following forms depending on the physical conditions of the crusts; (A)h(t)=0: the equation represents a logistic type increase or decrease and approaches a stationary state asymptotically. This describes aftershock series of large earthquakes and earthquake swarms of large scale such as the Wakayama and Matsushiro swarms in Japan; (B)h(t)=constant (β=0): frequencyn(t) increases initially and then decreases gradually and shows some kind of volcanic swarms; (C)h(t) = κ · {exp(?γ₁ t) ? exp(γ₂ t)}, (γ₂ > γ₁): this denotes time delay effects and the model shows periodic patterns of bursts or “rhythms” of earthquakes, which are observed in earthquake swarms. When external effects are taken into consideration, the model is further generalized and can describe various seismic patterns. These effects represent various influences of the circumstances like the earth tide and fluctuations of plate motions, etc. Whenh(t) takes type (A) and the external effect is random, the equation displays repetitive random patterns with bursts. Particularly interesting cases may be those whenh(t) is type (C) and the external force is periodic like the earth tide. Various nonperiodic as well as periodic patterns of earthquakes appear. These are the phenomena of “chaos” and “entrainment”, etc. and can be commonly observed. Varieties of actual earthquake patterns seem to be, at least partly, explained by the nonlinear coupling between the tidal forces and autonomous rhythms of earthquakes.  相似文献   

Estimation of lead sources in a Japanese cedar ecosystem using stable isotope analysis     

Yuko Itoh  Kyotaro NoguchiMasamichi Takahashi  Toru OkamotoShuichiro Yoshinaga 《Applied Geochemistry》2007

Anthropogenic Pb affects the environment worldwide. To understand its effect on forest ecosystem, Pb isotope ratios were determined in precipitation, various components of vegetation, the forest floor, soil and parent material in a Japanese cedar (Cryptomeria japonica D. Don) forest stand. The average ²⁰⁶Pb/²⁰⁷Pb ratio in bulk precipitation was 1.14 ± 0.01 (mean ± SD), whereas that in the subsoil (20–130 cm) was 1.18 ± 0.01. Intermediate ratios ranging from 1.15 to 1.16 were observed in the vegetation, the forest floor, and the surface soil (0–10 cm). Using the ²⁰⁶Pb/²⁰⁷Pb ratios, the contribution of anthropogenic sources to Pb accumulated in the forest were estimated by the simple binary mixing model. Sixty-two percent of the Pb in the forest floor, 71% in the vegetation, and 55% in the surface soil (0–10 cm) originated from anthropogenic sources, but only 16% in the sub-surface soil (10–20 cm) was anthropogenic. These results suggest that internal Pb cycling occurs mainly between surface soil and vegetation in a Japanese cedar ecosystem, and that anthropogenic Pb strongly influences Pb cycling. Although the Japanese cedar ecosystem has a shallow forest floor, very little atmospherically derived Pb migrated downward over 10 cm in depth.  相似文献   

Earthquake-induced residual deformation of Ariake clay deposits with leaching   总被引：1，自引：0，他引：1  

H. Nagase  K. Shimizu  A. Hiro-oka  Y. Tanoue  Y. Saitoh 《Soil Dynamics and Earthquake Engineering》2006,26(2-4):209-220

It is widely known that the sensitivity ratio of Ariake marine clay increases, if the clay is leached. Therefore, it is assumed that the mechanical properties of such clay considerably change and the shear strength decreases, if it is subjected to a cyclic load as a disturbance force, like earthquakes. In the present study, dynamic deformation and strength tests, as well as consolidation tests were conducted on Ariake clay that had leached. Static shear tests were also carried out on the clay after it had been subjected to a cyclic load. The stability of man-made earth structures on Ariake clay deposits that had leached during earthquakes was investigated. Undisturbed samples of the clay sedimented in the sea area were used in all the tests. Furthermore, an earthquake-induced residual deformation analysis was also conducted on the basis of the test results described above. Consequently, the following behaviors were observed in the study. (1) When Ariake marine clay was leached, the dynamic strength decreased and the shear modulus after subjecting it to a cyclic load remarkably decreased. (2) The effect of the leaching on the behavior of Ariake clay deposits was considerable during earthquakes.  相似文献   

Magnetotelluric resistivity modeling for 3D characterization of geothermal reservoirs in the Western side of Mt. Aso, SW Japan   总被引：1，自引：1，他引：1  

Hisafumi Asaue  Katsuaki Koike  Toru Yoshinaga  Shinichi Takakura   《Journal of Applied Geophysics》2006,58(4):296

Geothermal reservoirs are usually located at a depth range of 2 to 5 km, so to efficiently utilize such resources an advanced prospecting method is needed to detect these deep geologic structures. This study aimed to three-dimensionally characterize geothermal reservoirs by a combination of Magnetotelluric (MT) survey, inversion analysis of apparent resistivity, and interpolation of the resistivity data obtained. The western side of Mt. Aso crater, southwest Japan, was chosen as the case study area. Three hot springs exist there and a fault is assumed to go in the direction connecting them. A MT survey was carried out at 26 sites and the data processed by a remote reference method to reduce artificial noises. Based on skewness and Mohr circle analyses of the impedance tensor, the local geologic structure at each site could be approximated as horizontally layered and therefore, a one-dimensional inversion analysis was applied to the MT raw data. The resultant resistivity column data were then interpolated by the three-dimensional optimization principle method. The resistivity distributions obtained clarified continuous conductors with especially low resistivity (less than 10 Ω·m) at the hot springs along the fault. Because the resistivity decreases largely with an abundance of clay minerals, the conductors could be considered to correspond with the cap rocks. Thus, two geothermal reservoirs, whose shapes were estimated to be pillar, were detected under the cap rocks in an elevation range of − 1000 to − 3000 m. By comparing the resistivity distributions with the temperature distributions based on fluid-flow calculations at a steady state using FEM, the validity of the location and dimension of the estimated reservoirs were confirmed.  相似文献   

Concentrations of 222 Rn, Its Short-Lived Daughters And 212 Pb And Their Ratios Under Complex Atmospheric Conditions And Topography     

Toshio Kataoka  Eiji Yunoki  Mitsuo Shimizu  Tadashige Mori  Osamu Tsukamoto  Satoshi Takahashi  Hironori Fudeyasu  Yukitaka Ohashi  Ken Sahashi  Toshihiko Maitani  Koh'ichi Miyashita  Toru Iwata  Takayuki Sasaki  Yoko Fujikawa  Akira Kudo  Roger H. Shaw 《Boundary-Layer Meteorology》2003,107(1):219-249

Atmospheric activity concentrations of ²¹²Pb and short-lived ²²²Rndaughters, together with meteorological elements, have been observed continuously atthree sites at Kamisaibara Village in Japan. In addition, atmospheric activity concentrationof ²²²Rn, equilibrium-equivalent concentration of ²²²Rn and conditionsof the lower atmosphere were observed for three intensive observation periods at Akawase,one of the three sites in Kamisaibara Village. The equilibrium-equivalent concentration of²²²Rn is almost the same as the atmospheric activity concentration of short-lived²²²Rn daughters.The activity concentrations of ²¹²Pb and the short-lived ²²²Rn daughtersand their ratio were low in the daytime owing to convective mixing, and high at nightowing to the surface-based inversion during periods of no precipitation. Their variationshave several patterns corresponding to the scale of the drainage wind or weak mixing.Mechanical mixing due to strong winds through both day and night during the first andsecond observation periods made the atmospheric activity concentrations of ²¹²Pb and the short-lived ²²²Rn daughters continuously low. However, their ratios werecontinuously high during the first period yet continuously low during the second period.This difference can be explained by the effect ofextraction of ²²⁰Rn and ²²²Rndue to strong winds and snow cover. There were also cases in which the ratio of theatmospheric activity concentration of ²¹²Pb to that of the short-lived ²²²Rndaughters at night was equal to or less than the ratio in the daytime. Thisinverse trend, asin the periods of no precipitation mentioned above, is considered to be due to near-neutralconditions on these nights.We find a difference in the ratio of the equilibrium-equivalent concentration of²²²Rn (the activity concentration of short-lived ²²²Rn daughters) tothe activity concentration of ²²²Rn during the first observation period and thatduring the second. The difference can be explained by snow cover on the ground. Wealso find differences among the ratios of the activity concentration of the short-lived²²²Rn daughters to that of ²²²Rn during the three observation periods.These differences can be explained by the submergence of paddy fields.  相似文献   

High Resolution, Deep-Tow Seismic Survey to Investigate the Methane Hydrate Stability Zone in the Nankai Trough     

Peter Ward    Eiichi Asakawa  Shoshiro Shimizu 《Resource Geology》2004,54(1):115-124

Abstract. A high frequency deep-tow seismic survey was carried out in the Nankai Trough area in 1996. The objective of the survey was to obtain high resolution seismic sections and velocity profiles of the methane hydrate zone, inferred from the strong BSR events seen on conventional seismic data in the area. A special feature of the survey is that both the source and the streamer cable are towed close to the seabed. This special acquisition geometry requires special data processing to handle the varying source and receiver depths. A CMP floating datum processing sequence was designed which led to high quality sections of the shallow geology. A key step in the processing was devising a residual statics technique to compensate for errors in the measured depths.
The processing sequence was applied to a number of lines, totaling 200 km. The final data quality was highly variable. Some lines produced high quality sections and others, much poorer sections with few interpretable events. Conventional seismic data in the area also shows variation in the data quality so part of the reason is a variation in the sub-sea geology, but the deep-tow data is much more sensitive to change in conditions than conventional data. With the current acquired data and processing sequence the deep-tow system offers most advantages when 1) the water depth is around 1 km or greater, 2) the seabed and underlying geology is not too complex, and 3) the acquisition proceeds smoothly with regular shotpoints, slowly varying depths, and with accurate positioning.  相似文献   

Formation of the Yakuno ophiolite; accretionary subduction under medium-pressure-type metamorphic conditions     

Soichi Osozawa  Hiroshi Takeuchi  Toru Koitabashi 《Tectonophysics》2004,393(1-4):197

The notion that the Yakuno ophiolite and overlying Maizuru Group represents an accretionary prism formed during the Permian evolution of Japan on the Yakuno eruptive sequence, association of hemipelagic mudstone with silicic tuff, exotic fossiliferous limestones derived from previously accreted sea-mounts, upward coarsening of sequences terrigenous sandstone and conglomerate, and mildly deformed Permian and Triassic forearc basin formations. The most important indicator, however, is the seaward imbrication and repetition observed in both the Maizuru Group and the ophiolite itself. D1 deformation structures include axial–planar foliations (pressure-solution cleavage for the Maizuru Group and granulite–amphibolite metamorphic layering in the ophiolite), flattening type strain, symmetric pressure shadows and fringes, and isoclinal folds showing axial–planar foliations and thrust faulting at their overturned limb. The exceptional asymmetry observed indicates seaward-directed shearing near the thrust, while D1 structures in the Maizuru zone are explained by off-scraping, above the basal decollement. The later Jurassic D2 kink fold structure includes a first-order asymmetric kink with a brittle thrust at its overturned limb, more-or-less coeval with M2 retrograde metamorphism. Medium-pressure M1 prograde metamorphism in the Yakuno ophiolite produced layering of granulite and amphibolite, and in the Maizuru Group, formation of illite along pressure-solution cleavage of mudstones. The metamorphic grade is controlled by the stratigraphic relationships and appears typical of that in ocean floor regions. However, there was only one episode of M1 prograde metamorphism which occurred contemporaneously with D1 off-scraping. Given that subduction zones are normally characterized by high P/T metamorphic regimes, the observed P/T history appears to reflect relatively unusual conditions. Such high thermal gradients may plausibly reflect the approach of a young, hot oceanic plate which continued subducting beneath the Japanese arc. Accordingly, the Yakuno ophiolite was probably formed at the trench–trench–ridge triple junction.  相似文献   

Constraints on the pre-metamorphic trace element composition of Barberton komatiites from ion probe analyses of preserved clinopyroxene     

Stephen W. Parman  Nobumichu Shimizu  Timothy L. Grove  Jesse C. Dann 《Contributions to Mineralogy and Petrology》2003,144(4):383-396

  相似文献   

中国东部富钾埃达克岩成因的实验约束   总被引：27，自引：21，他引：27  

RobertP.PAPP  肖龙 《岩石学报》2002,18(3):293-302

Adakite在地球化学上具明显特征的火山岩和深成花岗岩类岩石，见于洋内岛孤环境和大陆孤，如安底斯孤。在洋内岛孤，由热的消减的大洋岩石圈熔融形成（叫做“板片熔融”），而在大陆孤，熔融曾发生在构造或岩浆加厚的下地壳底（叫做“下地壳熔融”）。在这两种产状环境中，adakite的鲜明地球化学特征被认为是起因子，一种不同程度含水的变质基性原岩在足够深度上的部分熔融，这里的足够深度是指可使石榴子石在残余结晶组合（即石榴角闪石和/或榴辉石的残余）中保持稳定的深度。“原始”或“母”adakite熔体一旦形成，便可能在其向上运移和侵位上地壳期间受到同化作用（或是地幔，或是大陆物质）和结晶分异作用的改造。中国东部晚中生代（早中白垩世，160-110Ma）的adakite,与见于同一地区和其它地方的钠质adakite相比，通常富含钾(K2O)和其它大离子亲石元素（如Ba,Th,U），有较低的Na2O/K2O比值（-1.0-1.1）,类似于玄武岩在石榴角闪岩-榴辉岩相含水熔融实验中所产生adakite熔体，要么是由洋壳板片熔融所形成，要么是由不同成分的玄武质下地壳原岩部分熔融所形成。尽管有些成分差异，它们的总体化学特征仍然可将中国东部的富钾花岗岩类岩石定均adakite。我们把这些富钾的adakite的独特化学行特征，归因于成分来源的特殊性，或adakite母岩浆遭受了同化混染和结晶分异(AFC)作用的改造。虽然中国东部与消减带环境明显不同这一点表明，那里的adakite可由板块底部侵位的（岩浆加厚的）镁铁质下地壳部分熔融所形成，但燕山运动期间中国东部存在“平坦”俯冲的地球动力学环境是可能被排除的。  相似文献   

Mantle heterogeneity beneath the southern Mid-Atlantic Ridge: trace element evidence for contamination of ambient asthenospheric mantle   总被引：1，自引：0，他引：1  

P.J le Roux  A.P le Roex  J.-G SchillingN Shimizu  W.W PerkinsN.J.G Pearce 《Earth and Planetary Science Letters》2002,203(1):479-498

We report new trace element data for an extensive suite of quench basalt glasses dredged from the southern Mid-Atlantic Ridge (MAR) between 40°S and 52.5°S. Ratios between highly incompatible trace elements are strongly correlated and indicate a systematic distribution of incompatible element enriched mid-ocean ridge basalt (MORB) (E-type: Zr/Nb=5.9-19, Y/Nb=0.9-8.4, (La/Sm)_n=1.0-2.9) and incompatible element depleted MORB (N-type: Zr/Nb=30-69, Y/Nb=11-29, (La/Sm)_n=0.48-0.79) along this section of the southern MAR. A notable feature of N-type MORB from the region is the higher than usual Ba/Nb (4-9), La/Nb (1.2-2.4) and primitive mantle normalised K/Nb ratios (>1). Ba/Nb ratios in E-type MORB samples from 47.5 to 49°S are especially elevated (>10). The occurrence and geographic distribution of E-type MORB along this section of the southern MAR can be correlated with the ridge-centred Shona and off-axis Discovery mantle plumes. In conjunction with published isotope data for a subset of the same sample suite [Douglass et al., J. Geophys. Res. 104 (1999) 2941], a model is developed whereby prior to the breakup of Gondwana and the opening of the South Atlantic Ocean, the underlying asthenospheric mantle was locally contaminated by fluids/melts rising from the major Mesozoic subduction zone along the south-southwest boundary of Gondwana, leaving a subduction zone geochemical imprint (elevated (K/Nb)_n and ⁸⁷Sr/⁸⁶Sr ratios, decreased ¹⁴³Nd/¹⁴⁴Nd ratios). Subsequent impingement of three major mantle plume heads (Tristan/Gough, Discovery, Shona) resulted in heating and thermal erosion of the lowermost subcontinental lithosphere and dispersal into the convecting asthenospheric mantle. With the opening of the ocean basin, continued plume upwelling led to plume-ridge interactions and mixing between geochemically enriched mantle derived from the Shona and Discovery mantle plumes, material derived from delamination of the subcontinental lithosphere, and mildly subduction zone contaminated depleted asthenospheric mantle.  相似文献