Multidisciplinary inferences on a newly recognized active east‐dipping extensional system in Central Italy

We use a multidisciplinary approach to gather preliminary evidence for a Quaternary east‐dipping extensional detachment in Central Italy. This structure crops out in the Sabini‐Eastern Simbruini (SES) and would be hidden at mid‐crustal depths beneath the L'Aquila 2009 (M_w6.3) epicentral area. The SES geometry is reconstructed through geological mapping, structural analysis and seismic line interpretation. The geometry of the mid‐crustal segment, referred to as the Ocre Segment (OS), is interpreted through seismological analyses of the largest aftershock (M_w5.4) of the L'Aquila 2009 sequence. The kinematic compatibility between the SES and the OS under a common SW–NE tensional field is tested through stress inversion of both geological and seismological data. The reliability of OS activation is tested through slip tendency analysis. Like other Italian cases, the SES and the OS are preliminarily interpreted as expressions at different depths of the same unknown east‐dipping extensional detachment, characterized by a ramp–flat–ramp geometry.     

利用卫星长波长地球物理位场对东亚地区深部结构的研究

有研究表明从中国北部经黄海、朝鲜半岛、日本海至日本本州岛存在一个大型断裂带,被称为“北纬40°断裂带”.为了验证该断裂带的存在性,本文通过GRACE和CHAMP卫星的位场数据中的大尺度区域重力和磁力异常来研究该区的大型构造特征.区域性研究范围为24°N~56°N,90°E~150°E,详细研究范围为32°N~42°N,122°E~132°E.2008年Taylor等利用CHAMP卫星2005年6月至12月收集的轨道数据,从测量数据中消除地核和外部产生磁场效应,经化磁极之后得到磁异常图,从中来追踪构造特征,并根据地质和地球物理资料建立了数学模型来进行解释.本文在此基础上,增加了GRACE卫星2003年10月份测量的重力异常场数据.实测差分位场数据利用高斯-勒让德求积法转换为垂向重力异常值,通过减去最新地球重力场模型（EGM96）的场值来去除波长超过1100 km的长波长异常,所得的重力异常数据和磁场数据进行比较.为了能从中得到研究区的共同构造特征,建立一个球面棱柱模型进行正演模拟、并进行重、磁异常场的波数相关性分析.研究发现在重、磁异常中均发现一个呈负相关的纬向构造特征,东亚地区“北纬40°断裂带”横穿了整个朝鲜半岛.     

Elasticity of antigorite,seismic detection of serpentinites,and anisotropy in subduction zones

Serpentinization of the mantle wedge is an important process that influences the seismic and mechanical properties in subduction zones. Seismic detection of serpentines relies on the knowledge of elastic properties of serpentinites, which thus far has not been possible in the absence of single-crystal elastic properties of antigorite. The elastic constants of antigorite, the dominant serpentine at high-pressure in subduction zones, were measured using Brillouin spectroscopy under ambient conditions. In addition, antigorite lattice preferred orientations (LPO) were determined using an electron back-scattering diffraction (EBSD) technique. Isotropic aggregate velocities are significantly lower than those of peridotites to allow seismic detection of serpentinites from tomography. The isotropic V_P/V_S ratio is 1.76 in the Voigt–Reuss–Hill average, not very different from that of 1.73 in peridotite, but may vary between 1.70 and 1.86 between the Voigt and Reuss bonds. Antigorite and deformed serpentinites have a very high seismic anisotropy and remarkably low velocities along particular directions. V_P varies between 8.9 km s^? 1 and 5.6 km s^? 1 (46% anisotropy), and 8.3 km s^? 1 and 5.8 km s^? 1 (37%), and V_S between 5.1 km s^? 1 and 2.5 km s^? 1 (66%), and 4.7 km s^? 1 and 2.9 km s^? 1 (50%) for the single-crystal and aggregate, respectively. The V_P/V_S ratio and shear wave splitting also vary with orientation between 1.2 and 3.4, and 1.3 and 2.8 for the single-crystal and aggregate, respectively. Thus deformed serpentinites can present seismic velocities similar to peridotites for wave propagation parallel to the foliation or lower than crustal rocks for wave propagation perpendicular to the foliation. These properties can be used to detect serpentinite, quantify the amount of serpentinization, and to discuss relationships between seismic anisotropy and deformation in the mantle wedge. Regions of high V_P/V_S ratios and extremely low velocities in the mantle wedge of subduction zones (down to about 6 and 3 km.s^?1 for V_P and V_S, respectively) are difficult to explain without strong preferred orientation of serpentine. Local variations of anisotropy may result from kilometer-scale folding of serpentinites. Shear wave splittings up to 1–1.5 s can be explained with moderately thick (10–20 km) serpentinite bodies.     

Daily morphological changes determined by high‐energy events on an embayed beach: a qualitative model

This paper examines the daily morphological responses of Sununga Beach, an embayed beach located on the south‐eastern Brazilian coast, to storms in the South Atlantic Ocean. The main mechanisms and timing of beach erosion and accretion, the relationship between wave height and direction, and beach volume changes are considered, to establish a qualitative model for short‐term embayed beach morphological changes. The methodology consisted of daily topographic surveys during the month of May in 2001, 2002, and 2003, using an RTK‐GPS (real‐time kinematics – global positioning system). Weather and wave model results were used to correlate hydrodynamics and beach morphology. The results indicate that the morphodynamics of Sununga Beach are characterized by a process of beach rotation, which occurred more or less clearly during all three surveys. Unlike what has been commonly described in the literature for longer time intervals and alternations of fair and stormy weather, the beach rotation processes on Sununga Beach occurred under conditions of moderate‐to‐high wave energy change (wave heights greater than 2 m). An integrated evaluation of the behaviour of the meteorological aspects, together with beach morphology, enabled us to recognize that extra‐tropical cyclones were the most important agent in remobilizing the beach planform, whether in beach rotation or in cross‐shore erosion. Copyright © 2010 John Wiley & Sons, Ltd.     

Hydrothermal calderas

The standard model of caldera formation is related to the emptying of a magma chamber and ensuing roof collapse during large eruptions or subsurface withdrawal. Although this model works well for numerous volcanoes, it is inappropriate for many basaltic volcanoes (with the notable exception of Hawaii), as these have eruptions that involve volumes of magma that are small compared to the collapse. Many arc volcanoes also have similar oversized depressions, such as Poas (Costa Rica) and Aoba (Vanuatu). In this article, we propose an alternative caldera model based on deep hydrothermal alteration of volcanic rocks in the central part of the edifice. Under certain conditions, the clay-rich altered and pressurized core may flow under its own weight, spread laterally, and trigger very large caldera-like collapse. Several specific mechanisms can generate the formation of such hydrothermal calderas. Among them, we identify two principal modes: mode 1: ripening with summit loading and flank spreading and mode II: unbuttressing with flank subsidence and flank sliding. Processes such as summit loading or flank subsidence may act simultaneously in hybrid mechanisms. Natural examples are shown to illustrate the different modes of formation. For ripening, we give Aoba (Vanuatu) as an example of probable summit loading, while Casita (Nicaragua) is the type example of flank spreading. For unbuttressing, Nuku Hiva Island (Marquesas) is our example for flank subsidence and Piton de la Fournaise (La Réunion) is our example of flank sliding. The whole process is slow and probably needs (a) at least a few tens of thousands of years to deeply alter the edifice and reach conditions suitable for ductile flow and (b) a few hundred years to achieve the caldera collapse. The size and the shape of the caldera strictly mimic that of the underlying weak core. Thus, the size of the caldera is not controlled by the dimensions of the underlying magma reservoir. A collapsing hydrothermal caldera could generate significant phreatic activity and trigger major eruptions from a coexisting magmatic complex. As the buildup to collapse is slow, such caldera-forming events could be detected long before their onset.     

Definition and interests of reciprocity and reciprocity gap principles for groundwater flow problems

We introduce the reciprocity and reciprocity gap principles for flow problems in hydrogeology and illustrate their interest in addressing identification problems. The reciprocity principle is derived from mechanics and establishes for flow problems a relationship between different sets of forcing terms, including sources, sinks and boundary conditions, and the resulting head fields. The reciprocity gap principle compares different head fields resulting from the same forcing terms applied to different structures. We give general 2D expressions of the reciprocity and reciprocity gap principles for transient flow problems and give two examples of applications for the identification of transmissivity values and interfaces between different transmissivities. Identification capacities of the reciprocity and reciprocity gap principles yielding direct inversion methods could be used as initial guesses for more advanced inverse problem methodologies.     

Statistical characteristics of flow as indicators of channeling in heterogeneous porous and fractured media

We introduce two new channeling indicators D_ic and D_cc based on the Lagrangian distribution of flow rates. On the basis of the participation ratio, these indicators characterize the extremes of both the flow-tube width distribution and the flow rate variation along flow lines. The participation ratio is an indicator biased toward the larger values of a distribution and is equal to the normalized ratio of the square of the first-order moment to the second-order moment. Compared with other existing indicators, they advantageously provide additional information on the flow channel geometry, are consistently applicable to both porous and fractured media, and are generally less variable for media generated using the same parameters than other indicators. Based on their computation for a broad range of porous and fracture permeability fields, we show that they consistently characterize two different geometric properties of channels. D_ic gives a characteristic scale of low-flow zones in porous media and a characteristic distance between effectively flowing structures in fractured cases. D_cc gives a characteristic scale of the extension of high-flow zones in porous media and a characteristic channel length in fractured media. D_ic is mostly determined by channel density and permeability variability. D_cc is, however, more affected by the nature of the correlation structure like the presence of permeability channels or fractures in porous media and the length distribution in fracture networks.     

Sky luminance distribution in Pamplona (Spain) during the summer period

In this work the outdoor daylight conditions in Pamplona (South Europe) during the summer period have been studied. The selected sky type (from fifteen standards) at a given moment is the one exhibiting the lowest RMSD when comparing the theoretical and experimental luminance distributions in the sky hemisphere. Two year data of luminance distribution registered every 10 min in 145 positions of the sky hemisphere have been used for selecting the sky type. The most frequent sky type in Pamplona is V.5 (cloudless polluted with a broad solar corona), with an occurrence of 29.5%. This result coincides with the one observed in a previous study in Athens.Six types of sky (V.5, IV.4, III.4, III.3, V.4 y II.2) out of the fifteen standards become practically the 80% of all the studied ones. Regarding a possible use in daylight climate studies, the frequency of occurrence of the fifteen types of sky for fourteen solar elevation intervals has been included.     

The roles of surface heat flux and ocean heat transport convergence in determining Atlantic Ocean temperature variability

The temperature variability of the Atlantic Ocean is investigated using an eddy-permitting (1/4°) global ocean model (ORCA-025) forced with historical surface meteorological fields from 1958 to 2001. The simulation of volume-averaged temperature and the vertical structure of the zonally averaged temperature trends are compared with those from observations. In regions with a high number of observations, in particular above a depth of 500 m and between 22° N and 65° N, the model simulation and the dataset are in good agreement. The relative contribution of variability in ocean heat transport (OHT) convergence and net surface heat flux to changes in ocean heat content is investigated with a focus on three regions: the subpolar and subtropical gyres and the tropics. The surface heat flux plays a relatively minor role in year-to-year changes in the subpolar and subtropical regions, but in the tropical North Atlantic, its role is of similar significance to the ocean heat transport convergence. The strongest signal during the study period is a cooling of the subpolar gyre between 1970 and 1990, which subsequently reversed as the mid-latitude OHT convergence transitioned from an anomalously weak to an anomalously strong state. We also explore whether model OHT anomalies can be linked to surface flux anomalies through a Hovmöller analysis of the Atlantic sector. At low latitudes, increased ocean heat gain coincides with anomalously strong northward transport, whereas at mid-high latitudes, reduced ocean heat loss is associated with anomalously weak heat transport.     

Capturing the residence time boundary layer—application to the Scheldt Estuary

At high Peclet number, the residence time exhibits a boundary layer adjacent to incoming open boundaries. In a Eulerian model, not resolving this boundary layer can generate spurious oscillations that can propagate into the area of interest. However, resolving this boundary layer would require an unacceptably high spatial resolution. Therefore, alternative methods are needed in which no grid refinement is required to capture the key aspects of the physics of the residence time boundary layer. An extended finite element method representation and a boundary layer parameterisation are presented and tested herein. It is also explained how to preserve local consistency in reversed time simulations so as to avoid the generation of spurious residence time extrema. Finally, the boundary layer parameterisation is applied to the computation of the residence time in the Scheldt Estuary (Belgium/The Netherlands). This timescale is simulated by means of a depth-integrated, finite element, unstructured mesh model, with a high space–time resolution. It is seen that the residence time temporal variations are mainly affected by the semi-diurnal tides. However, the spring–neap variability also impacts the residence time, particularly in the sandbank and shallow areas. Seasonal variability is also observed, which is induced by the fluctuations over the year of the upstream flows. In general, the residence time is an increasing function of the distance to the mouth of the estuary. However, smaller-scale fluctuations are also present: they are caused by local bathymetric features and their impact on the hydrodynamics.