Plate generation and two-phase damage theory in a model of mantle convection

The formation of narrow, rapidly deforming plate boundaries separating strong plate interiors are integral components of the generation of plate tectonics from mantle convection. The development of narrow plate boundaries requires the interaction of a non-linear rheology and convection. One such non-linear rheology is two-phase damage theory which employs a non-equilibrium relation between interfacial surface energy, pressure and viscous deformation, thereby forming a theoretical model for void generation. Two-phase damage theory was recently extended to allow for deformational work to increase the fineness (reduce the grain size) of the matrix phase. We present results testing two-phase damage theory in a 2-D convectively driven system where we allow for (1) pure void-generating damage, (2) pure fineness-generating damage and (3) combined void- and fineness-generating damage. Pure void-generating damage is found to be unsuccessful at producing plate-like features. Fineness-generating damage is successful at inducing plate-like behaviour in certain circumstances, including increasing viscosity sensitivity to fineness and certain regimes of damage input and healing rate. Cases with combined void- and fineness-generating damage produce significantly more localization than the end-members due to the apparent increase of deformational work input into fineness generation. The interaction of microcracks and grain size reduction in two-phase damage theory suggests a rheological model for shear localization necessary for the formation of plate tectonic boundaries.     

A multidisciplinary study of a slow-slipping fault for seismic hazard assessment: the example of the Middle Durance Fault (SE France)

Assessing seismic hazard in continental interiors is difficult because these regions are characterized by low strain rates and may be struck by infrequent destructive earthquakes. In this paper, we provide an example showing that interpretations of seismic cross sections combined with other kinds of studies such as analysis of microseismicity allow the whole seismogenic source area to be imaged in this type of region. The Middle Durance Fault (MDF) is an 80-km-long fault system located southeastern France that has a moderate but regular seismicity and some palaeoseismic evidence for larger events. It behaves as an oblique ramp with a left-lateral-reverse fault slip and has a low strain rate. MDF is one of the rare slow active fault system monitored by a dedicated dense velocimetric short period network. This study showed a fault system segmented in map and cross section views which consists of staircase basement faults topped by listric faults ramping off Triassic evaporitic beds. Seismic sections allowed the construction of a 3-D structural model used for accurate location of microseismicity. Southern part of MDF is mainly active in the sedimentary cover. In its northern part and in Alpine foreland, seismicity deeper than 8 km was also recorded meaning active faults within the crust cannot be excluded. Seismogenic potential of MDF was roughly assessed. Resulting source sizes and estimated slip rates imply that the magnitude upper limit ranges from 6.0 to 6.5 with a return period of a few thousand years. The present study shows that the coupling between 3-D fault geometry imaging and accurate location of microseismicity provides a robust approach to analyse active fault sources and consequently a more refined seismic hazard assessment.     

基于SD和CA的城镇土地扩展模拟模型——以江苏省南通地区为例

充分利用系统动力学模型(System Dynamics,SD)在情景模拟和宏观因素反映上的优势和元胞自动机模型(Cellular Automata,CA)在微观土地利用空间格局反映上的优势,构建一个耦合SD和CA的城镇土地扩展模拟模型,并以江苏省南通地区为例,对模型的实证应用做了进一步的验证。结果表明,这种耦合模型不仅能够对研究区域未来城镇土地扩展数量给予一个比较好的预测,而且还对其空间分布效果做了一定精度上的模拟,这使得城市规划在土地利用预测方面有一个相对科学的依据。     

Gravity constraints on the dynamics of the crust–mantle system during Calabrian subduction

The thermomechanic evolution of the lithosphere–upper mantle system during Calabrian subduction is analysed using a 2-D finite element approach, in which the lithosphere is compositionally stratified into crust and mantle. Gravity and topography predictions are cross-checked with observed gravity and topography patterns of the Calabrian region. Modelling results indicate that the gravity pattern in the arc-trench region is shaped by the sinking of light material, belonging to both the overriding and subduction plates. The sinking of light crustal material, up to depths of the order of 100–150 km is the ultimate responsible for the peculiar gravity signature of subduction, characterized by a minimum of gravity anomaly located at the trench, bounded by two highs located on the overriding and subducting plates, with a variation in magnitude of the order of 200 mGal along a wavelength of 200 km, in agreement with the isostatically compensated component of gravity anomaly observed along a transect crossing the Calabrian Arc, from the Tyrrhenian to the Ionian Seas. The striking agreement between the geodetic retrieved profiles and the modelled ones in the trench region confirms the crucial role of compositional stratification of the lithosphere in the subduction process and the correctness of the kinematic hypotheses considered in our modelling, that the present-day configuration of crust–mantle system below the Calabrian arc results from trench's retreat at a rate of about 3 cm yr⁻¹, followed by gravitational sinking of the subducted slab in the last 5 Myr.     

Numerical modelling of rise and fall of a dense layer in salt diapirs

Numerical models are used to study the entrainment of a dense anhydrite layer by a diapir. The anhydrite layer is initially horizontally embedded within a viscous salt layer. The diapir is down-built by aggradation of non-Newtonian sediments ( n = 4, constant temperature) placed on the top of the salt layer. Several parameters (sedimentation rate, salt viscosity, perturbation width and stratigraphic position of the anhydrite layer) are studied systematically to understand their role in governing the entrainment of the anhydrite layer. High sedimentation rates during the early stages of the diapir evolution bury the initial perturbation and, thus, no diapir forms. The anhydrite layer sinks within the buried salt layer. For the same sedimentation rate, increasing viscosity of the salt layer decreases the rise rate of the diapir and reduces the amount (volume) of the anhydrite layer transported into the diapir. Model results show that viscous salt is capable of carrying separate blocks of the anhydrite layer to relatively higher stratigraphic levels. Varying the width of the initial perturbation (in our calculations 400–800 m), from which a diapir triggers, shows that wider diapirs can more easily entrain an embedded anhydrite layer than the narrower diapirs. The anhydrite layer is entrained as long as rise rate of the diapir exceeds the descent rate of the denser anhydrite layer. We conclude that the four parameters mentioned above govern the ability of a salt diapir to entrain an embedded dense layer. However, the model results show that the entrained blocks inevitably sink back if the rise rate of the diapir is less than the rate of descent of the anhydrite layer or the diapir is permanently covered by a stiff overburden in case of high sedimentation rates.     

Quasi-static analysis of strike fault growth in layered media

We study the effects of structural inhomogeneity on the quasi-static growth of strike-slip faults. A layered medium is considered, made up of an upper layer bounded by a free surface and welded to a lower half-space with different elastic property. Mode III crack is employed as a mathematical model of strike-slip fault, which is nucleated in the lower half-space and then propagates towards the interface. We adopt FEM-β, newly proposed analysis method for failure, to simulate the quasi-statistic crack growth governed by the stress distribution in layered media. Our results show that along planar traces across interfaces a compliant upper layer has significant effects on promoting/suppressing crack growth before/after its extension into the layer and vice versa for a rigid one. This proposes a possibility that surface breaks due to strike-slip faulting could be arrested by deposit layers at the topmost part of the Earth's crust.     

金字塔沙丘地表气流场及其动力学过程研究

沙丘沉积构造指示及风洞实验结果表明：莫高窟顶金字塔沙丘形成于３组风向条件下（ＮＥ、ＳＷ及ＮＷ）,是链状沙丘的一种变形、分离形式。从沙丘形态特征分析,三面体金字塔沙丘及其它复杂沙丘类型（线形、格状、多面体金字塔沙丘）之间存在某种内在本质的联系。即多信风条件下产生的风向组合,通过主、副脊（梁）的不同形式发育组合过程,形成适应于地表气流的地貌形态。同时认为：沙丘沉积构造及风洞实验相结合的综合研究比以往单一过程研究更为全面、系统、深刻。     

Note on rain-triggered earthquakes and their dependence on karst geology

Recently reported rain-triggered seismicity from three separate storms occurred exclusively in karst geology. In this paper, I discuss how the hydrogeology of karst controls rain-triggered seismicity by channeling of the watershed after intense rainfall directly into the karst network. Such channeling results in very large increases in hydraulic head, and more importantly, substantially increases the vertical stress acting on the underlying pore-elastic media. Rapid loading upon a pore-elastic media induces seismicity by increasing pore pressure at depth in a manner similar to that observed from reservoir impounding. Using a simple 1-D model of a pore-elastic medium, it is shown that the instantaneous fluid pressure increase at depth is a substantial fraction of the pressure step applied at the boundary, followed by time-dependent pore pressure increases associated with the typical linear diffusion problem. These results have implications for the change in fluid pressure necessary to trigger earthquakes, and leads to the following hypothesis to be tested: Unambiguous rain-triggered seismicity will only occur in karst regions.     

Long-term evolution of the aerosol debris cloud produced by the 2009 impact on Jupiter

We present a study of the long-term evolution of the cloud of aerosols produced in the atmosphere of Jupiter by the impact of an object on 19 July 2009 (Sánchez-Lavega, A. et al. [2010]. Astrophys. J. 715, L155-L159). The work is based on images obtained during 5 months from the impact to 31 December 2009 taken in visible continuum wavelengths and from 20 July 2009 to 28 May 2010 taken in near-infrared deep hydrogen-methane absorption bands at 2.1-2.3 μm. The impact cloud expanded zonally from ∼5000 km (July 19) to 225,000 km (29 October, about 180° in longitude), remaining meridionally localized within a latitude band from 53.5°S to 61.5°S planetographic latitude. During the first two months after its formation the site showed heterogeneous structure with 500-1000 km sized embedded spots. Later the reflectivity of the debris field became more homogeneous due to clump mergers. The cloud was mainly dispersed in longitude by the dominant zonal winds and their meridional shear, during the initial stages, localized motions may have been induced by thermal perturbation caused by the impact’s energy deposition. The tracking of individual spots within the impact cloud shows that the westward jet at 56.5°S latitude increases its eastward velocity with altitude above the tropopause by 5-10 m s⁻¹. The corresponding vertical wind shear is low, about 1 m s⁻¹ per scale height in agreement with previous thermal wind estimations. We found evidence for discrete localized meridional motions with speeds of 1-2 m s⁻¹. Two numerical models are used to simulate the observed cloud dispersion. One is a pure advection of the aerosols by the winds and their shears. The other uses the EPIC code, a nonlinear calculation of the evolution of the potential vorticity field generated by a heat pulse that simulates the impact. Both models reproduce the observed global structure of the cloud and the dominant zonal dispersion of the aerosols, but not the details of the cloud morphology. The reflectivity of the impact cloud decreased exponentially with a characteristic timescale of 15 days; we can explain this behavior with a radiative transfer model of the cloud optical depth coupled to an advection model of the cloud dispersion by the wind shears. The expected sedimentation time in the stratosphere (altitude levels 5-100 mbar) for the small aerosol particles forming the cloud is 45-200 days, thus aerosols were removed vertically over the long term following their zonal dispersion. No evidence of the cloud was detected 10 months after the impact.     

The reflection seismic survey of project TIPTEQ—the inventory of the Chilean subduction zone at 38.2° S

We describe results of an active-source seismology experiment across the Chilean subduction zone at 38.2°S. The seismic sections clearly show the subducted Nazca plate with varying reflectivity. Below the coast the plate interface occurs at 25 km depth as the sharp lower boundary of a 2–5 km thick, highly reflective region, which we interpret as the subduction channel, that is, a zone of subducted material with a velocity gradient with respect to the upper and lower plate. Further downdip along the seismogenic coupling zone the reflectivity decreases in the area of the presumed 1960 Valdivia hypocentre. The plate interface itself can be traced further down to depths of 50–60 km below the Central Valley. We observe strong reflectivity at the plate interface as well as in the continental mantle wedge. The sections also show a segmented forearc crust in the overriding South American plate. Major features in the accretionary wedge, such as the Lanalhue fault zone, can be identified. At the eastern end of the profile a bright west-dipping reflector lies perpendicular to the plate interface and may be linked to the volcanic arc.