— To elucidate the spatial complexity of damage and evolution of localized failure in the transitional regime from brittle faulting to cataclastic ductile flow in a porous sandstone, we performed a series of triaxial compression experiments on Rothbach sandstone (20% porosity). Quantitative microstructural analysis and X-ray computed tomography (CT) imaging were conducted on deformed samples. Localized failure was observed in samples at effective pressures ranging from 5 MPa to 130 MPa. In the brittle faulting regime, dilating shear bands were observed. The CT images and stereological measurements reveal the geometric complexity and spatial heterogeneity of damage in the failed samples. In the transitional regime (at effective pressures between 45 MPa and 130 MPa), compacting shear bands at high angles and compaction bands perpendicular to the maximum compression direction were observed. The laboratory results suggest that these complex localized features can be pervasive in sandstone formations, not just limited to the very porous aeolian sandstone in which they were first documented. The microstructural observations are in qualitative agreement with theoretical predictions of bifurcation analyses, except for the occurrence of compaction bands in the sample deformed at effective pressure of 130 MPa. The bifurcation analysis with the constitutive model used in this paper is nonadequate to predict compaction band formation, may be due to the neglect of bedding anisotropy of the rock and multiple yield mechanisms in the constitutive model. 相似文献
It gradually becomes a common work using large seismic wave data to obtain source parameters, such as seismic moment, break radius, stress drop, with completingof digital seismic network in China (Hough, et al, 1999; Bindi, et al, 2001). These parameters are useful on earthquake prediction and seismic hazard analysis.Although the computation methods of source parameters are simple in principle and the many research works have been done, it is not easy to obtain the parameters accurately. There are two factors affecting the stability of computation results. The first one is the effect of spread path and site respond on signal. According to the research results, there are different geometrical spreading coefficients on different epicenter distance. The better method is to introduce trilinear geometrical spreading model (Atkinson, Mereu, 1992; Atkinson, Boore, 1995; WONG, et al, 2002). In addition, traditional site respond is estimated by comparing with rock station, such as linear inversion method (Andrews, 1982), but the comparative estimation will introduce some errors when selecting different stations. Some recent research results show that site respond is not flat for rock station (Moya, et al, 2000; ZHANG,. et al, 2001; JIN, et al, 2000; Dutta, et al, 2001). The second factor is to obtain low-frequency level and corner frequency fromdisplacement spectrum. Because the source spectrum model is nonlinear function,these values are obtained by eye. The subjectivity is strong. The small change of corner frequency will affect significantly the result of stress drop. 相似文献
— The Papua New Guinea (PNG) tsunami of 1998 is re-examined through a detailed review of the field survey as well as numerous numerical computations. The discussion of the field survey explores a number of possible misinterpretations of the recorded data. The survey data are then employed by a numerical model as a validation tool. A Boussinesq model and a nonlinear shallow water wave (NLSW) model are compared in order to quantify the effect of frequency dispersion on the landslide-generated tsunami. The numerical comparisons indicate that the NLSW model is a poor estimator of offshore wave heights. However, due to what appears to be depth-limited breaking seaward of Sissano spit, both numerical models are in agreement in the prediction of maximum water elevations at the overtopped spit. By comparing three different hot-start initial profiles of the tsunami wave, it is shown that the initial shape and orientation of the tsunami wave is secondary to the initial displaced water mass in regard to prediction of water elevations on the spit. These numerical results indicate that agreement between numerical prediction of runup values with field recorded values at PNG cannot be used to validate either a NLSW tsunami propagation model or a specific landslide tsunami hot-start initial condition. Finally, with the use of traditional tsunami codes, a new interpretation of the PNG runup measurements is presented. 相似文献
The thermal structure of Archean and Proterozoic lithospheric terranes in southern Africa during the Mesozoic was evaluated by thermobarometry of mantle peridotite xenoliths erupted in alkaline magmas between 180 and 60 Ma. For cratonic xenoliths, the presence of a 150–200 °C isobaric temperature range at 5–6 GPa confirms original interpretations of a conductive geotherm, which is perturbed at depth, and therefore does not record steady state lithospheric mantle structure.
Xenoliths from both Archean and Proterozoic terranes record conductive limb temperatures characteristic of a “cratonic” geotherm (40 mW m−2), indicating cooling of Proterozoic mantle following the last major tectonothermal event in the region at 1 Ga and the probability of thick off-craton lithosphere capable of hosting diamond. This inference is supported by U–Pb thermochronology of lower crustal xenoliths [Schmitz and Bowring, 2003. Contrib. Mineral. Petrol. 144, 592–618].
The entire region then suffered a protracted regional heating event in the Mesozoic, affecting both mantle and lower crust. In the mantle, the event is recorded at 150 Ma to the southeast of the craton, propagating to the west by 108–74 Ma, the craton interior by 85–90 Ma and the far southwest and northwest by 65–70 Ma. The heating penetrated to shallower levels in the off-craton areas than on the craton, and is more apparent on the southern margin of the craton than in its western interior. The focus and spatial progression mimic inferred patterns of plume activity and supercontinent breakup 30–100 Ma earlier and are probably connected.
Contrasting thermal profiles from Archean and Proterozoic mantle result from penetration to shallower levels of the Proterozoic lithosphere by heat transporting magmas. Extent of penetration is related not to original lithospheric thickness, but to its more fertile character and the presence of structurally weak zones of old tectonism. The present day distribution of surface heat flow in southern Africa is related to this dynamic event and is not a direct reflection of the pre-existing lithospheric architecture. 相似文献
The interactions between climatic and volcanic forcing on diatom communities contained in a 50,000-year sedimentary sequence from Lake Massoko, Tanzania, were examined. At the century scale, 19 discrete tephra inputs to the lake isolated the sedimentary nutrient supply and shifted the diatom communities to those tolerant of low phosphorus levels, whereas at the millennial scale, diatom-inferred shifts in precipitation–evaporation based on conductivity optima and diatom life-form ratios were broadly similar to lake-level reconstructions from Lake Rukwa, Lake Malawi, and others in the region. Some fluctuations of Lake Massoko are consistent with the precession-driven changes in insolation, but the major climate shifts do not relate directly to orbital forcing of summer insolation south of the equator and show more consistency with records from the equatorial and northern tropics that receive rainfall from the passing of the intertropical convergence zone. Sea surface temperatures are strongly correlated to multimillennial-scale climate patterns over this region of Africa. 相似文献
To speed up the calculation of the field Jacobian for 2-D magnetoteliuric inversion using finite elements, the principle of electromagnetic reciprocity is applied. The governing relationship for the Jacobian of the field along strike is obtained by differentiating the Helmholtz equation with respect to the resistivity of each region in the finite-element mesh. The result is a similar Helmholtz equation for the Jacobian, with new sources distributed over all nodes within the parameter medium. However, according to the principle of electromagnetic reciprocity, the roles of sources and receivers are interchangeable. Utilizing reciprocity, the field values obtained from the original forward problem and for new unit sources imposed at the receivers are then utilized in the calculation of the Jacobian by simple multiplication and summation with finite-element terms at each rectangle in the mesh. For the auxiliary (across-strike) fields, the Jacobian terms are obtained by solving source vectors loaded with parabola coefficients used in the approximation to Maxwell's equations. Jacobian terms for the apparent resistivity ( p a), the impedance phase (φ) and the vertical magnetic field ( K zy) are then calculated utilizing the parallel- and auxiliary-field Jacobians. Comparison of Jacobian values obtained from reciprocity calculations and by differencing two forward solutions show that the reciprocity method is accurate and can be used to decrease the number of calculations required to obtain sensitivities by one to two orders of magnitude. 相似文献