Counting chronology and climate records with about 1000 annual layers of a Holocene stalagmite from the Water Cave in Liaoning Province, China

One active stalagmite from the Water Cave in Liaoning Province contains growth layers of three sizes. Based on thermal ionization mass spectrometry 230Th dating, we found that middle size layers are annual layers, with each middle layer consisting of one narrow dark layer and a wide bright layer. The small layers within middle layers are sub-annual layers and the large layers are multi-year layers. Based on the layer-counting method, we established a high-resolution time scale for layer thickness. Our results reveal two dramatic century-scale climate cycles over the past 1000 years in this region.     

The use of modern statistical theories in the assessment of earthquake hazard, with application to quiet regions of eastern North America

It is generally accepted that both deterministic and statistical approaches are useful for the characterization of earthquake hazard. Although the most reliable estimates of seismic hazard can only be based on an improved understanding of the earthquake mechanism, efficient utilization of the appropriate methods provided by recent statistical theories is also important in seismic risk analysis. This is especially true in regions where the connection between seismicity and geologic structure is tenuous at best. We are particularly interested in developing better statistical treatments of data for regions with little known seismic activity. To this end, we have applied three statistical methods to the historical record of seismicity in relatively quiet regions of eastern North America. These are: (1) the threshold method for tail inference, a new theory for modeling earthquakes with sizes above a given threshold, (2) the ‘bootstrap’ technique in which the characteristics of an unknown population are simulated by replacing the true population by an estimated one, and (3) a technique to estimate the number of earthquakes below a given size, in order to compensate for the under-reporting of small earthquakes in most catalogs. A combination of these techniques has been used to estimate the probabilities of future large earthquakes for the regions studied. Because of limitations imposed by existing catalogs, the size estimate used has been maximum intensity.     

中朝地台东北缘地区的地震层析成像

根据中朝地台东北缘地区 (东经 1 1 7°0 0′— 1 2 6°0 0′ ,北纬 36°0 0′— 44°0 0′) 1 980— 1 997年的 380 0 0余条P波走时数据 ,利用正交投影法重建了该区地壳和上地幔的三维速度结构 .通过分析及同人工地震测深剖面的详细对比 ,证明了成像结果的可靠性 .结果表明 :中朝地台东北缘地区地壳上地幔介质存在显著的横向不均匀性 ,直至 1 2 0km深度处依然明显 ;地壳上部的速度图像清楚地反映了不同岩石单元的分布与该区不同性质的基岩分布基本吻合 ;从上、中地壳的速度图像中发现了研究区存在海城、朝阳、义县、丹东南、唐山等几个低速异常区(即速度逆反层区 ) ,其中海城、唐山、朝阳等地区的壳内低速层已由深地震测深资料所证实 ;研究区陆地发生的几次强震均发生于壳内低速层上方的高速脆性介质内 ,而渤海发生的强震 ,此现象不明显 ,但都发生于横向介质速度显著突变的位置 ;在地壳不同深度上发现了普兰店至山海关横跨渤海的北西向低速异常带 ;地震层析二维速度图像与深地震测深资料的对比表明 ,研究区利用地震层析成像技术 ,在一定条件下可以获得与人工地震测深相似的效果     

非规则测区三维地震偏移的边界吸收层法

对不规则测区的地震资料直接进行三维偏移不仅可减少计算量,而且可提高偏移剖面的质量. 本文在单程波方程中增加了一项简单的阻尼因子,其作用如同在实际不规则测区的外部有一个薄吸收层. 因为波穿越这个薄层时,波场值迅速衰减,所以偏移时把吸收层外边界处的波场值取为零所产生的反射非常微弱. 又因为不规则测区零边界条件偏移程序极易编制,所以就可以实现对不规则测区的地震资料直接进行三维偏移. 通过对实际地震资料的处理,说明了该方法的有效性.     

Bedrock structure in Adapazari, Turkey—a possible cause of severe damage by the 1999 Kociaeli earthquake

The 1999 Kocaeli earthquake brought serious damage to downtown of Adapazari. To study why strong motions were generated at the town, a bedrock structure was investigated on the basis of Bouguer gravity anomaly, and SPAC and H/V analyses of microseisms. It was revealed that, the basin consists of three narrow depressions of bedrock with very steep edges, extending in E–W or NE–SW directions along the North Anatolia faults, and the depth to bedrock reaches 1000 m or more. Downtown of Adapazari is located 1–2 km apart from the basin-edge. It is considered that, the specific configuration of bedrock amplifies ground motions at the downtown area by focusing of seismic waves and/or interference between incident S-waves and surface-waves secondarily generated at the basin-edge. Studying 3D bedrock structure is an urgent issue for microzoning an urban area in a sedimentary basin.     

Effective stress analyses of liquefaction-induced deformation in river dikes

A series of effective stress analyses is carried out on the seismic performance of river dikes based on the case histories during the 1993 Hokkaido-Nansei-oki and 1995 Hyogoken-Nambu earthquakes in Japan. Seven case histories selected for the analyses involve a crest settlement ranging from none to 2.7 m in the dikes 3–6 m high with evidence of liquefaction at foundation soil. The effective stress model used is based on a multiple shear mechanism and was developed by one of the authors. The soil parameters are evaluated based on the site investigation and laboratory test results. The results of the analyses are basically consistent with the observed performance of the river dikes. In particular, the effective stress model shows a reasonable capability to reproduce the varying degree of settlements depending on the geotechnical conditions of foundation soils beneath the dikes. The analyses also indicate that the effect of a cohesive soil layer mixed with the liquefiable sand layers beneath the dikes can be a primary factor for reducing the liquefaction-induced deformation of dikes.     

Efficient methods for solving water flow in variably saturated soils under prescribed flux infiltration

The non-linear solvers in numerical solutions of water flow in variably saturated soils are prone to convergence difficulties. Many aspects can give rise to such difficulties and in this paper we address the gravity term and the prescribed-flux boundary in the Picard iteration. The problem of the gravity term in the Picard iteration is iteration-to-iteration oscillation as the gravity term is treated, by analogy with the time-step advance technique, ‘explicitly’ in the iteration. The proposed method for the gravity term is an improvement of the ‘implicit’ approach of Zhang and Ewen [Water Resour. Res. 36 (2000) 2777] by extending it to heterogeneous soil and approximating the inter-nodal hydraulic conductivity in the diffusive term and the gravity term with the same scheme. The prescribed-flux boundary in traditional methods also gives rise to iteration-to-iteration oscillation because there is no feedback to the flux in the solution at the new iteration. To reduce such oscillation, a new method is proposed to provide such a feedback to the flux. Comparison with traditional Picard and Newton iteration methods for a wide range of problems show that a combination of these two proposed methods greatly improves the stability and consequently the computational efficiency, making the use of small time step and/or under-relaxation solely for convergence unnecessary.     

Stable isotope values of Costa Rican surface waters

Stable isotope data of surface waters from the humid tropics in general, and Costa Rica in particular, are scarce. To improve our understanding of the spatial distribution of stable isotopes in surface waters, we measured δ¹⁸O and δD in river and lake (n=63) and precipitation (n=3) samples from Costa Rica. We also present data from the IAEA/WMO isotopes in precipitation network as context for our study. Surface water isotope values do not strongly correlate with elevation, stream head elevation, stream length, distance from Caribbean Sea, or estimated mean annual precipitation for the country as a whole. However, the data show distinct regional trends. The δ¹⁸O and δD values downwind of mountain ranges are inversely related to the altitude of the ranges the air masses traverse. In the lee of the high Talamanca Range, δ¹⁸O values are 6–8‰ lower, while in the lee of the lower Tilarán Range δ¹⁸O values are 2–3‰ lower than upwind sites along the Caribbean Slope. An altitude effect of −1.4‰ δ¹⁸O/km is present on the Pacific slope of southern Costa Rica, equivalent to a temperature effect of −0.3‰/°C. The Nicoya and Osa Peninsulas have higher values than upwind sites, suggesting input of Pacific-sourced moisture, evaporative enrichment, or decreased condensation temperatures. Elevated and increasing d-excess values inland along the Nicaragua Trough suggest a recycled component may be an important contributor to the water budget. These data provide preliminary stable isotope information for Costa Rica, and will benefit paleoclimatic research in the region. More detailed studies would be beneficial to our understanding of the controls on stable isotope composition of tropical waters.     

Interpretation of charging on fracture or frictional slip surface of rocks

Fracturing and frictional sliding of quartz and granite under dry condition generates fractoluminescence, charged particle emission and electromagnetic radiation. Various kinds of experiments indicate that surface charge density on fracture or frictional slip surface of quartz and granite is 10⁻⁴ to 10⁻² C/m² which is larger than bound charges induced by the disappearance of piezoelectricity due to the release of stress. Hole and electron trapping centers, which is found in semiconductor devices with the Si–SiO₂ system, are causes of surface charging on fracture or frictional slip surface of quartz crystal. The quantity of the surface charge is enough to cause corona discharge that can generate earthquake lights. The mechanism considering the hole and electron trapping centers has a probability to explain why non-piezoelectric minerals or rocks generate electromagnetic phenomena. It can be one of origins of seismo-electromagnetic phenomena (SEP).     

Improved treatments for evaluating horizontal magnetic components through the 3-D FDM in E -polarization induction problems

To improve the accuracy of the numerical evaluation through the 3-D finite difference method, the surface boundary conditions are added to modify the old program. The author has tested the new program by making calculations for the model constructed by Wanamaker, et al (1984). The comparison between the numerical results obtained from this paper and those by Wannamaker, et al (1984) indicates that a pronounced improvement is realized in the evaluation of the horizontal magnetic components. Moreover, better calculations for the vertical magnetic components are also obtainable by using the new program.