Historical 1942 Ecuador and 1942 Peru subduction earthquakes and earthquake cycles along Colombia-Ecuador and Peru subduction segments

Two large shallow earthquakes occurred in 1942 along the South American subduction zone inclose proximity to subducting oceanic ridges: The 14 May event occurred near the subducting Carnegie ridge off the coast of Ecuador, and the 24 August event occurred off the coast of southwestern Peru near the southern flank of the subducting Nazca ridge. Source parameters for these for these two historic events have been determined using long-periodP waveforms,P-wave first motions, intensities and local tsunami data.We have analyzed theP waves for these two earthquakes to constrain the focal mechanism, depth, source complexity and seismic moment. Modeling of theP waveform for both events yields a range of acceptable focal mechanisms and depths, all of which are consistent with underthrusting of the Nazca plate beneath the South American plate. The source time function for the 1942 Ecuador event has one simple pulse of moment release with a duration of 22 suconds, suggesting that most of the moment release occurred near the epicenter. The seismic moment determined from theP waves is 6–8×10²⁰N·m, corresponding ot a moment magnitude of 7.8–7.9. The reported location of the maximum intensities (IX) for this event is south of the main shock epicenter. The relocated aftershcks are in an area that is approximately 200 km by 90 km (elongated parallel to the trench) with the majority of aftershocks north of the epicenter. In contrast, the 1942 Peru event has a much longer duration and higher degree of complexity than the Ecuador earthquake, suggesting a heterogeneous rupture. Seismic moment is released in three distinct pulses over approximately 74 seconds; the largest moment release occurs 32 seconds after rupture initiation. the seismic moment as determined from theP waves for the 1942 Peru event is 10–25×10²⁰N·m, corresponding to a moment magnitude of 7.9–8.2. Aftershock locations reported by the ISS occur over a broad area surrounding the main shock. The reported locations of the maximum intensities (IX) are concentrated south of the epicenter, suggesting that at least part of the rupture was to the south.We have also examined great historic earthquakes along the Colombia-Ecuador and Peru segments of the South American subduction zone. We find that the size and rupture length of the underthrusting earthquakes vary between successive earthquake cycles. This suggests that the segmentation of the plate boundary as defined by earthquakes this century is not constant.     

尕斯库勒油田N1—N2^1油藏砂西区Ⅱ油组储层非均质特征评价

储层非均质特征评价是储层特征研究很重要的一部分内容，储层非均质特征是描述储层的参数变化及非均质特性，非均质性的强弱直接影响着开发方案的设计和实施。本文以尕斯库勒油田N1—N2^1油藏砂西区Ⅱ油组为例，研究储层非均质性特征。用来描述储层非均质性的特征参数主要有：连通系数、重叠系数、夹层厚度、夹层频数以及均值、最大值、最小值、变异系数、级差、非均质系数。     

A new surface runoff parameterization with subgrid-scale soil heterogeneity for land surface models

Soil heterogeneity plays an important role in determining surface runoff generation mechanisms. At the spatial scales represented by land surface models used in regional climate model and/or global general circulation models (GCMs) for numerical weather prediction and climate studies, both infiltration excess (Horton) and saturation excess (Dunne) runoff may be present within a studied area or a model grid cell. Proper modeling of surface runoff is essential to a reasonable representation of feedbacks in the land–atmosphere system. In this paper, a new surface runoff parameterization that dynamically represents both Horton and Dunne runoff generation mechanisms within a model grid cell is presented. The new parameterization takes into account of effects of soil heterogeneity on Horton and Dunne runoff. A series of numerical experiments are conducted to study the effects of soil heterogeneity on Horton and Dunne runoff and on soil moisture storage under different soil and precipitation conditions. The new parameterization is implemented into the current version of the hydrologically based variable infiltration capacity (VIC) land surface model and tested over three watersheds in Pennsylvania. Results show that the new parameterization plays a very important role in partitioning the water budget between surface runoff and soil moisture in the atmosphere–land coupling system. Significant underestimation of the surface runoff and overestimation of subsurface runoff and soil moisture could be resulted if the Horton runoff mechanism were not taken into account. Also, the results show that the Horton runoff mechanism should be considered within the context of subgrid-scale spatial variability of soil properties and precipitation. An assumption of time-invariant spatial distribution of potential infiltration rate may result in large errors in surface runoff and soil moisture. In addition, the total surface runoff from the new parameterization is less sensitive to the choice of the soil moisture shape parameter of the distribution.     

Some critical factors for engineering and environmental microgravity investigations

The gravity method is one of the geophysical tools used for engineering and environmental investigations where the detection of cavities, karst phenomena, subsoil irregularities, or landfills is essential. In many cases, deep or small-scale heterogeneities generating low-amplitude anomalies have to be detected and the reliability of further interpretation requires highly accurate measurements, carefully corrected for any quantifiable disturbing effects. The purpose of this study is to investigate the factors likely to limit measurement quality and how to make improvements.Calibrations of a Scintrex gravimeter were made between French relative and absolute base stations, and the relative uncertainties on the calibration factors were estimated for these links. Ranging from 10⁻³, for calibration on an old gravity net, to 10⁻⁴, for a high amplitude absolute base line, this accuracy will be generally sufficient for microgravity surveys.Continuous gravity recordings of Scintrex gravimeters, installed at the same stable site, enabled the estimation of the stability and accuracy of the instruments and revealed that some of the time variations of g measurements, such as instrumental drift, tidal effects and seismic noise, are not entirely removed by standard processing procedures. The accuracy of corrected gravity measurements is mainly limited by inadequate corrections of tidal effects and by a poor estimation of ocean loading effects. In comparison with residual defaults in tidal corrections, instrumental and seismic noises are taken more properly into account by statistical data processing.In field operation, residual tidal effects are generally integrated into an experimental terrain drift estimated on the basis of frequent repeated measurements. A differential gravity approach, based on a fixed gravimeter reference whose recordings are used to correct measurements made with a mobile gravimeter, has also been investigated at a test site. Compared to standard processing, this method can help improve repeatability of gravity measurements.Microgravity surveys in the urban environment require effective and accurate consideration of the effects of infrastructures, nearby buildings and basements, as well as those of topography, in the vicinity of a gravity station. Correction procedures, applied at the same experimental site, where gravity points are located close to buildings, walls and basement slope, appear to have almost totally eliminated these disturbances.     

Nd-isotope systematics of ∼2.7 Ga adakites, magnesian andesites, and arc basalts, Superior Province: evidence for shallow crustal recycling at Archean subduction zones

An association of adakite, magnesian andesite (MA), and Nb-enriched basalt (NEB) volcanic flows, which erupted within ‘normal’ intra-oceanic arc tholeiitic to calc-alkaline basalts, has recently been documented in ∼2.7 Ga Wawa greenstone belts. Large, positive initial ?_Nd values (+1.95 to +2.45) of the adakites signify that their basaltic precursors, with a short crustal residence, were derived from a long-term depleted mantle source. It is likely that the adakites represent the melts of subducted late Archean oceanic crust. Initial ?_Nd values in the MA (+0.14 to +1.68), Nb-enriched basalts and andesites (NEBA) (+1.11 to +2.05), and ‘normal’ intra-oceanic arc tholeiitic to calc-alkaline basalts and andesites (+1.44 to +2.44) overlap with, but extend to lower values than, the adakites. Large, tightly clustered ?_Nd values of the adakites, together with Th/Ce and Ce/Yb systematics of the arc basalts that rule out sediment melting, place the enriched source in the sub-arc mantle. Accordingly, isotopic data for the MA, NEBA, and ‘normal’ arc basalts can be explained by melting of an isotopically heterogeneous sub-arc mantle that had been variably enriched by recycling of continental material into the shallow mantle in late Archean subduction zones up to 200 Ma prior to the 2.7 Ga arc. If the late Archean Wawa adakites, MA, and basalts were generated by similar geodynamic processes as their counterparts in Cenozoic arcs, involving subduction of young and/or hot ocean lithosphere, then it is likely that late Archean oceanic crust, and arc crust, were also created and destroyed by modern plate tectonic-like geodynamic processes. This study suggests that crustal recycling through subduction zone processes played an important role for the generation of heterogeneity in the Archean upper mantle. In addition, the results of this study indicate that the Nd-isotope compositions of Archean arc- and plume-derived volcanic rocks are not very distinct, whereas Phanerozoic plumes and intra-oceanic arcs tend to have different Nd-isotopic compositions.     

下地幔侧向不均匀及绝对板块运动对大地水准面的影响

本文提出并研究了板块上大地水准面受下地幔侧向不均匀性及板块绝对运动的影响而随时间变化的问题。利用已有的下地幔侧向不均匀性模型、板块绝对运动速度模型进行了计算,并分核幔边界不存在起伏和存在起伏两种情况进行了讨论。结果表明,在相当长时间内,其垂直变化绝对幅度可达 200米。确定古大地水准面应考虑这一变化。     

Seismic Modelling of Lateral Heterogeneity At the Base of the Mantle

Summary. Results from several recent studies suggest that there are lateral heterogeneities of up to a few per cent in the lowermost 150–200 km of the mantle (Bullen's D " region). Inferred anomaly sizes span the range from less than 50 km to greater than 1000 km.
In this study differences in the velocity structure among regions at the base of the mantle were inferred from an analysis of amplitude ratios of PKPAB and PKPDF for given earthquake-station pairs at distances greater than 155° (Sacks, Snoke & Beach). We distinguish two kinds of regions: A (anomalous) regions in which the mean, median and spread in AB/DF amplitude ratios are significantly higher (> 50 per cent) than for a reference radial earth model and N (normal) regions in which the distribution of the amplitude ratios is as expected.
The AB branch has near-grazing incidence to the core and therefore maximum sensitivity to velocity structure compared to the near-normal incident DF phases. Using an iterative, forward-modelling approach, we have determined general characteristics of the velocity structure for regions at the base of the mantle which can produce amplitude-ratio distributions similar to those for an A region. Agreement between model and data is obtained over the period range from 0.5 s to greater than 10 s using a laterally heterogeneous model for the D " region. the model consists of cells which are 200 km in lateral extent with velocity variations of up to ±1 per cent. This structure is modulated by a region-wide (1000km) perturbation which increases smoothly from zero at the edges of the region to a negative 1 per cent at the centre. Small cells (∼40 km) cannot produce anomalously large amplitude, long-period AB arrivals, and larger cells (∼1000km) cannot match the observed scatter. the ∼200 km scale anomalies could be small-scale convection cells confined to the D " region.     

Calculation of Well Index for Nonconventional Wells on Arbitrary Grids

Wells are seldom modeled explicitly in large scale finite difference reservoir simulations. Instead, the well is coupled to the reservoir through the use of a well index, which relates wellbore flow rate and pressure to grid block quantities. The use of an accurate well index is essential for the detailed modeling of nonconventional wells; i.e., wells with an arbitrary trajectory or multiple branches. The determination of a well index for such problems is complicated, particularly when the simulation grid is irregular or unstructured. In this work, a general framework for the calculation of accurate well indices for general nonconventional wells on arbitrary grids is presented and applied. The method entails the use of an accurate semianalytical well model based on Green's functions as a reference single phase flow solution. This result is coupled with a finite difference calculation to provide an accurate well index for each grid block containing a well segment. The method is demonstrated on a number of homogeneous example cases involving deviated, horizontal and multilateral wells oriented skew to the grid. Both Cartesian and globally unstructured multiblock grids are considered. In all these cases, the method is shown to provide results that are considerably more accurate compared to results using standard procedures. The method is also applied to heterogeneous problems involving horizontal wells, where it is shown to be capable of approximating the effects of subgrid heterogeneity in coarse finite difference models.