Explicit, approximate expressions for the resolution and a posteriori covariance of massive tomographic systems

We present an approximate method to estimate the resolution, covariance and correlation matrix for linear tomographic systems Ax = b that are too large to be solved by singular value decomposition. An explicit expression for the approximate inverse matrix A ⁻ is found using one-step backprojections on the Penrose condition AA ⁻≈ I , from which we calculate the statistical properties of the solution. The computation of A ⁻ can easily be parallelized, each column being constructed independently.
The method is validated on small systems for which the exact covariance can still be computed with singular value decomposition. Though A ⁻ is not accurate enough to actually compute the solution x , the qualitative agreement obtained for resolution and covariance is sufficient for many purposes, such as rough assessment of model precision or the reparametrization of the model by the grouping of correlating parameters. We present an example for the computation of the complete covariance matrix of a very large (69 043 × 9610) system with 5.9 × 10⁶ non-zero elements in A . Computation time is proportional to the number of non-zero elements in A . If the correlation matrix is computed for the purpose of reparametrization by combining highly correlating unknowns x _i , a further gain in efficiency can be obtained by neglecting the small elements in A , but a more accurate estimation of the correlation requires a full treatment of even the smaller A _ij . We finally develop a formalism to compute a damped version of A ⁻.     

Long-period geomagnetic variations and mantle conductivity: an inversion using Bailey's method

Summary. We have implemented an algorithm which is based on Bailey's solution of the inverse problem of electromagnetic induction in the Earth. The study was motivated by recent determinations of very long period data and also benefited from recent redeterminations of high frequency data. The algorithm has been successfully tested to provide reliable estimates of conductivity down to a depth of 2000 km, using synthetic data in the period range from 4 days to 11 years. Smooth data sets, which are required for the inversion, were constructed from various sources. At a given depth, the range of inverted models is less than one order of magnitude. Due to the lack of high frequency data, the conductivity of the upper 600 km of the mantle, which is found to be of the order of 10⁻¹Ω⁻¹ m⁻¹, may be overestimated. The algorithm performs well in the middle mantle, where conductivity rises steadily from 1 to 50 Ω⁻¹ m⁻¹. The lack of very low frequency data and limitations of the algorithm prevent one from obtaining meaningful estimates in the lower mantle. However, the study of the propagation of the late 1960s secular variation acceleration provides an estimate of the mean conductivity of the whole mantle. Thus, a complete mantle profile can be constructed. It is found that deep mantle conductivity probability does not exceed a few hundred Ω⁻¹ m⁻¹.     

New approach in the kinematic k−2 source model for generating physical slip velocity functions

In an attempt to improve the ground motion modelling, the characteristics of the slip velocity functions (SVF) generated using the kinematic k ⁻² source are investigated and compared to the dynamic solutions proposed in the literature. Several numerical simulations were performed to test the influence of the model parameters on the SVF modelling. Overall, the shapes of SVF are very complex and exhibit a large variability in time and space. However, we found out that the mean SVF is a simple boxcar with duration equal to the largest rise time value. In the areas of weak slip, the SVFs are characterized by the existence of negative values, whereas in large slip areas, the SVF is more impulsive. Overall, on the examples investigated, the SVFs modelled with this k ⁻² source model are different from a typical Kostrov's solution. The critical analysis of the kinematic k ⁻² source led us to identify the Fourier decomposition of the slip to be responsible for these difficulties, and to propose a new recombination scheme. It consists of adding a positive correction to the Fourier slip components. The slip is described as the sum of positive contributions at various scales. The SVFs modelled using this new scheme are greatly improved. Moreover, through several parametrical analyses performed to qualify this new approach, we show that the SVF are corrected while preserving the essential quality of the k ⁻² modelling, that is, the ω² spectral shape and C _d apparent directivity of the synthetic accelerograms. Strong ground motion modelling in the near-fault region was made and numerical ground motion parameters were compared to the empirical relationships. We show that predicted peak ground motion is consistent with near-source attenuation laws.     

Deglaciation effects on the rotation of the Earth

Summary. The viscoelastic response of the Earth to the mass displacements caused by late Pleistocene deglaciation and concomitant sea level changes is shown to be capable of producing the secular motion of the Earth's rotation pole as deduced from astronomical observations. The calculations for a viscoelastic Earth yield a secular motion in the direction of 72° W meridian which is in excellent agreement with observed values. The average Newtonian viscosity and the relaxation time obtained from polar motion data are about (1.1 ± 0.6)10²³ poise (P) and 10⁴ (1 ± 0.5) yr. The non-tidal secular acceleration of the Earth can also be attributed to the viscoelastic response to deglaciation and results in an independent viscosity estimate of 1.6 × 10²³ P with upper and lower limits of 1.1 × 10²³ and 2.8 × 10²³ P. These values are in agreement with those based on the polar drift analysis and indicate an average mantle viscosity of 1–2 × 10²³ P.     

Role of pore pressure generation in sediment transport within half-grabens

Sediment transport and overpressure generation are coupled primary through the impact of effective stress on subsidence and compaction. Here, we use mathematical modeling to explore the interactions between groundwater flow and diffusion-controlled sediment transport within alluvial basins. Because of lateral variation in permeability, proximal basin facies will have pore pressure close to hydrostatic levels while distal fine-grained facies can reach near lithostatic levels. Lateral variation in pore pressure leads to differential compaction, which deforms basins in several ways. Differential compaction reduces basin size, bends isochron surfaces across the sand–clay interface, restricts basinward progradation of sand facies, and reduces the amplitude of oscillation in the lateral position of the sand–clay interface especially in the deepest part of the section even when temporal sediment supply are held constant. Overpressure generation was found to be sensitive to change in sediment supply in permeable basins (at least 10⁻¹⁷ m² in our model). We found that during basin evolution, temporal variations in overpressure and sediment supply fluctuations are not necessarily in phase with each other, especially in tight (low permeability) basins (<10⁻¹⁷ m² in our model).     

Transient and impulse responses of a one-dimensional linearly attenuating medium—II. A parametric study

Summary. We investigate one-dimensional waves in a standard linear solid for geophysically relevant ranges of the parameters. The critical parameters are shown to be T*= t_u/Q_m where t _u is the travel time and Q_m the quality factor in the absorption band, and τ⁻¹ _m , the high-frequency cut-off of the relaxation spectrum. The visual onset time, rise time, peak time, and peak amplitude are studied as functions of T* and τ _m. For very small τ _m , this model is shown to be very similar to previously proposed attenuation models. As τ _m grows past a critical value which depends on T* , the character of the attenuated pulse changes. Seismological implications of this model may be inferred by comparing body wave travel times with a'one second'earth model derived from long-period observations and corrected for attenuation effects assuming a frequency independent Q over the seismic band. From such a comparison we speculate that there may be a gap in the relaxation spectrum of the Earth's mantle for relaxation times shorter than about one second. However, observational constraints from the attenuation of body waves suggest that such a gap might in fact occur at higher frequencies. Such a hypothesis would imply a frequency dependence of Q in the Earth's mantle for short-period body waves.     

The elastodynamic near field

Summary. The elastodynamic fields of point forces and shear dislocations of finite source duration are analysed with the aim of establishing the frequency and time-domain characteristics of the field in the near-source region. Criteria are obtained for amplitude dominance in regions where the source–sensor distance is much smaller than the wavelength.
It is shown that in the frequency domain , the Green's tensor (and hence the displacement field of a single point force) attenuates like R ⁻¹ in the near-source region and there exists no region in which the 'near-field' term becomes dominant such that the 'far-field' term can be neglected. Hence, there is no real 'near-field' term for the elastodynamic Green's tensor. The near-field terms of the displacements, velocities and accelerations excited by a shear dislocation attenuate like R ⁻², since the R ⁻³ and R ⁻⁴ terms tend to be eliminated due to mutual cancellation of P and S motions in the near-source region.
In the time domain , the corresponding near field of the displacement field is defined for the steady amplitude interval (away from transients) R /β < t < R /α+ T by the condition R ≤βT where β is the shear velocity and T is the source's duration. The relative strengths of all other arrivals will depend on the particular time window under consideration.
The particle motion patterns due to a single force in the near-source region are shown to be similar to rotating hyperbolas with an axis along the force direction, which are quite different from the 'smoke ring' motion patterns of the so-called 'near-field' term itself.     

Dynamic behaviour of a phase boundary under non-uniform surface loads

Summary. A method is outlined to determine the dynamic behaviour of a phase boundary in the Earth when non-uniform time-varying pressure and temperature conditions are assumed at the Earth's surface. An integral equation describing the phase boundary motion is derived and it is solved under a linearizing assumption. The solution is obtained in the form of a double integral transform. Short and long time-expansions of the solution can be obtained from series expansion and integration of the Laplace transform along a branch cut. The method is illustrated by considering a stepwise change in surface pressure conditions.
For short times, the solution exhibits the same type of time dependence (i.e. the first-order term is in t ^1/2) as the one obtained in the one-dimensional case (i.e. uniform pressure perturbation at the Earth's surface).
For long times, it is shown that the time dependence of the phase boundary motion is almost identical to the one derived for the one- dimensional case if the wavenumber k _L of the surface excitation is such that κ k ²_Lτ≤ 1 (where τ is the relaxation time associated with the one-dimensional phase boundary motion and κ is the thermal diffusivity). If κ k ²_Lτ > 1, then the relaxation time for the phase boundary motion in two dimensions is of the order of κ⁻¹ k ⁻²_L.
When considering parameters that would be appropriate for a basalt to eclogite phase transition at Moho depth, the latter situation is met only when the load wavelength is smaller than 35 km.     

A new method of palaeofield magnitude correction for thermally altered samples and its application to Lower Carboniferous lavas

Summary. Some lavas of Early Carboniferous age from Scotland and Derbyshire have been analysed for palaeofield magnitude by the application of a new correction to data obtained by the Shaw method. Its application has yielded field magnitude values from data which had previously been rejected. A modern pottery specimen and recent lavas from Sicily and Westman Island have also been analysed.
An average virtual dipole moment (VDM, Smith) of 2.5 ± 0.4 × 10²² A m² was determined for the Early Carboniferous. This is only 30 per cent of the present-day VDM.     

Oxygen consumption in Macrotrachela musculosa and Trichotria truncata (Rotatoria) from the High Arctic

Oxygen consumption by rotifers Macrotrachela musculosa and Trichotria truncata from Spitsbergen tundra (77°N) was measured using the method of Cartesian divers. The metabolic rate of M. musculosa was: 0.205 10⁻³mm³ 0₂ per g 10⁻⁶ per hour at 2°C, 0.201 10⁻⁶mm³ at 6°C and 0.616 10⁻³mm³ 0₂ per g 10⁻⁶ per hour at 10°C. The metabolic rate of Trichotria truncata at 6° was 0.103 10⁻³mm³ per g 10⁻⁶ per hour. The relation between body weight and oxygen consumption by M. musculosa at 2°C is expressed with the equation R = 0.18W^0.67, with R – oxygen consumption in mm³10⁻³ per individual per hour and W – wet weight of an animal in g 10⁻⁶.     

Palaeosecular variation in Central Mexico over the last 30 000 years: the record from lavas

13 lava flows of known age (ages from ¹⁴C dating), which have been erupted in the last 30 000 years, have been studied to determine the palaeosecular variation of the geomagnetic field in Central Mexico. Samples were taken from two different monogenetic volcanic fields: the Michoacan-Guanajuato volcanic field (six sites) and the Chichinautzin Formation (seven sites), both part of the Transmexican Volcanic Belt. The lavas were studied in detail using rock magnetic methods (magnetic susceptibility at room temperature, low-temperature susceptibility behaviour, hysteresis loops, Curie temperatures), combined with reflected light microscopy, in order to deduce their magnetic mineralogy and the domain states of the magnetic minerals. The magnetic carriers are titanomagnetites, which show differing degrees of high-temperature deuteric oxidation, and seem to be predominantly pseudo-single domain (PSD), though in many cases are probably a mixture of domain states. Mean palaeomagnetic directions and palaeointensity values using Shaw and Thellier techniques were obtained using several specimens from each flow. Our data seem to indicate a sharp easterly swing in declination about 5000 years ago, which is also observed in lake sediments from Central Mexico. The calculated values of the virtual dipole moment (VDM) range from 3.1 to 14.9 × 10²² A m². Our data indicate that the virtual dipole moment seems to have increased gradually in magnitude over the last 30 kyr, with a peak at about 9000 years BP. These are features that have been observed in other parts of the globe and are probably caused by variations in the dipole part of the geomagnetic field.     

On the consistency of earthquake moment release and space geodetic strain rates: Europe

In this paper, approximately 100 VLBI/SLR/GPS velocities map European strain rates from <0.09 × 10⁻⁸ to >9.0 × 10⁻⁸ yr⁻¹ with regional uncertainties of 20 to 40 per cent. Kostrov's formula translates these strain-rate values into regional geodetic moment rates M¯˙ _geodetic . Two other moment rates, M¯˙ _seismic , extracted from a 100-year historical catalogue and M¯˙ _plate , taken from plate-tectonic models, contrast the geodetic rates. In Mediterranean Europe, the ratios of M¯˙ _seismic to M¯˙ _geodetic are between 0.50 and 0.71. In Turkey the ratio falls to 0.22. Although aseismic deformation may contribute to the earthquake deficit ( M¯˙ _seismic values less than M¯˙ _geodetic ), the evidence is not compelling because the magnitudes of the observed shortfalls coincide with the random variations expected in a 100-year catalogue. If the lack of aseismic deformation inferred from the 100-year catalogue holds true for longer periods, then much of Europe's strain budget would have to be accommodated by more frequent or larger earthquakes than have been experienced this century to raise the ratios of M¯˙ _seismic to M¯˙ _geodetic to unity. Improved geological fault data bases, longer historical earthquake catalogues, and densification of the continent's space geodetic network will clarify the roles of aseismic deformation versus statistical quiescence.     

Coseismic and other short-term strain changes recorded with Sacks-Evertson strainmeters in a deep mine, South Africa

Summary. During 1977 March and April, three Sacks-Evertson borehole dilatometers were installed at the ends of boreholes drilled into the sidewall of an experimental tunnel at a depth of 3.1 km in the ERPM gold mine near Johannesburg. In the following year coseismic strain changes ranging from 5 ± 10⁻¹⁰ to values exceeding 5 ± 10⁻⁶ were recorded for hundreds of mine tremors in the magnitude range -1 to 3.7 and at hypocentral distances of 50 m to about 2 km. Hypocentral coordinates and magnitudes were determined from seismograms recorded from an underground array of geophones. Amplitudes and polarities of the coseismic strain steps are generally in excellent agreement with theoretical expectations based on point-source dislocation theory; specifically, the strain steps are proportional to the seismic moment divided by the cube of hypocentral distance. At a strain level of 5 ± 10⁻⁹ or greater the tremors do not appear to be preceded by any short-term indications of instability even for tremors producing coseismic steps greater than 5 ± 10⁻⁶ and for which the strainmeters were within a source radius of the hypocentre. Continuous strain changes observed at the times when the mine excavation, at a distance of about 100 m, is extended are in good agreement with calculated changes based on the theory of elasticity. A similar calculation is consistent with post-seismic strain changes observed to follow some of the closer tremors. These post-seismic strains show a logarithmic dependence on time following the tremor and appear to be due to the interaction of a tremor with the adjacent mine excavation rather than to deformation within the actual seismic source region.     

Numerical simulation of the propagation of P waves in fractured media

We study the propagation of P waves through media containing open fractures by performing numerical simulations. The important parameter in such problems is the ratio between crack length and incident wavelength. When the wavelength of the incident wavefield is close to or shorter than the crack length, the scattered waves are efficiently excited and the attenuation of the primary waves can be observed on synthetic seismograms. On the other hand, when the incident wavelength is greater than the crack length, we can simulate the anisotropic behaviour of fractured media resulting from the scattering of seismic waves by the cracks through the time delay of the arrival of the transmitted wave. The method of calculation used is a boundary element method in which the Green's functions are computed by the discrete wavenumber method. For simplicity, the 2-D elastodynamic diffraction problem is considered. The rock matrix is supposed to be elastic, isotropic and homogeneous, while the cracks are all empty and have the same length and strike direction. An iterative method of calculation of the diffracted wavefield is developed in the case where a large number of cracks are present in order to reduce the computation time. The attenuation factor Q ⁻¹ of the direct waves passing through a fractured zone is measured in several frequency bands. We observe that the attenuation factor Q ⁻¹ of the direct P wave peaks around kd = 2, where k is the incident wavenumber and d the crack length, and decreases proportionally to ( kd ) ⁻¹ in the high-wavenumber range. In the long-wavelength domain, the velocity of the direct P wave measured for two different crack realizations is very close to the value predicted by Hudson's theory on the overall elastic properties of fractured materials.     

Stress drops, wave spectra and recurrence intervals of great earthquakes — implications of the Etorofu earthquake of 1958 November 6

Summary . The great Etorofu earthquake of 1958 November 6 is characterized by a relatively small aftershock area (70 × 150 km²) and an extremely large felt area. The felt area is more extensive than those of any other large earthquakes which have occurred in the southern Kurile to northern Japan arc since the beginning of this century. The mechanism is a pure thrust fault typical of most great earthquakes in island arcs. A body wave magnitude of m _b = 8.2 is obtained at periods around 6 s using more than 40 observations, although an m _b value of only 7.6–7.7 would be expected empirically from the observed surface wave magnitude of M _s= 8.1–8.2. Both an unusually large felt area and a high m _b indicate a dominance of high-frequency components in the seismic waves. A seismic moment of M _o= 4.4 × 10²⁸ dyne cm is determined from long-period surface waves from which a high stress drop of Δσ = 78 bar is obtained using a relatively small aftershock area. Historic data indicate an anomalously long time interval between the 1958 event and any earlier great earthquake from the same source region. The observed high stress drop can be interpreted as a consequence of this long intervening period through which strain built up. The dominance of the high-frequency seismic waves can then be interpreted as a result of this high stress drop. Stress drops, seismic wave spectra and recurrence intervals of great earthquakes are in this way closely related to each other. The 1958 event may represent a high strength extreme of stochastic fluctuation of fracture strength relevant to great earthquakes.     

A near-regional explosion source model for tuff

Summary. A variety of near-regional (300 km) data, including spectral amplitudes of Pg , surface-wave forms, and close-in (5–10 km) accelerograms have been used to build an elastic seismic source model for a 1-Mton explosion in tuff at near-regional distances. The model consists of: (1) a pressure pulse which injects 3 × 10¹² cm³ of volume into the medium, (2) a vertical, upward force impulse that imparts 10¹⁸ dyn-s of momentum to the medium, each source component having a time duration of 0.6 s and a depth of 1.3 km. The force impulse appears to be required by two considerations: (a) the striking similarity, apart from sign, of explosion surface waves with those of their cavity collapses, (b) the observation of considerable SV energy leaving the source of the 1-Mton explosions JORUM and HANDLEY . Scaling curves have been constructed which fit the proposed source model. These scaling curves employ: very slow decrease, as (yield)^−0.10 of the primary corner frequency; decay as (frequency)⁴ or (frequency)³ to high frequency. While these scaling curves are unconventional, they appear to be the only ones which can satisfy the near-regional data. The slow scaling with yield of the spectral carner frequency suggests that it is caused by something other than the equivalent elastic radius, e.g. the time duration of motion at the source. The results, at odds with similar studies at teleseismic distances, suggest that significantly different equivalent elastic sources are required at near-regional (as compared with teleseismic) distances; therefore, the effect of the upward impulse might not be seen at teleseismic distances. Consequently, these results probably do not pertain to the seismic discrimination problem at teleseismic distances.     

生物模块沙障对流动沙丘土壤微生物分布的影响

生物模块沙障是一种新型沙袋沙障,同时具备防风固沙、局部土壤改良和植被建植功能.区别于其他沙袋沙障,生物模块沙障外包装为黄麻纤维,内部除填充风沙土外,还添加了玉米芯和羊粪.为了解生物模块沙障的局部土壤改良功能,从土壤微生物角度出发,测定和研究了生物模块沙障内部及周围土壤空间内细菌、真菌、放线菌数量分布情况.结果表明:生物模块沙障内部土壤细菌、真菌、放线菌数量及三者总数量分别达到了2.41×10⁷ cfu·g^-1、4.51×10⁷ cfu·g^-1、1.19×10⁷ cfu·g^-1及8.11×10⁷ cfu·g^-1,而对照仅分别为1.56×10⁶ cfu·g^-1、3.01×10⁶ cfu·g^-1、1.26×10⁶ cfu·g^-1及5.83×10⁶ cfu·g^-1;生物模块沙障内部3种微生物及三者总数量均高于其覆盖的下方土层,下方各土层3种微生物数量均显著高于未设置沙障区;在生物模块沙障迎风侧50 cm至背风侧50 cm区域,3种微生物总数量大小顺序为沙障放置处 > 沙障背风侧20 cm处 > 沙障迎风侧20 cm处 > 沙障背风侧50 cm处 > 沙障迎风侧50 cm处,在此防护区内土壤细菌、真菌及放线菌数量平均值分别是流动沙丘的19.56、24.47、8.39倍.生物模块沙障设置5年后,其内部3种微生物数量显著高于周围土壤,说明生物模块沙障对周围土壤微生物的影响依然存在潜力.     

Geodetic observations of ice flow velocities over the southern part of subglacial Lake Vostok, Antarctica, and their glaciological implications

In the austral summer seasons 2001/02 and 2002/03, Global Positioning System (GPS) data were collected in the vicinity of Vostok Station to determine ice flow velocities over Lake Vostok. Ten GPS sites are located within a radius of 30 km around Vostok Station on floating ice as well as on grounded ice to the east and to the west of the lake. Additionally, a local deformation network around the ice core drilling site 5G-1 was installed.
The derived ice flow velocity for Vostok Station is  2.00 m a⁻¹± 0.01 m a⁻¹  . Along the flowline of Vostok Station an extension rate of about 10⁻⁵ a⁻¹ (equivalent to 1 cm km⁻¹ a⁻¹) was determined. This significant velocity gradient results in a new estimate of 28 700 years for the transit time of an ice particle along the Vostok flowline from the bedrock ridge in the southwest of the lake to the eastern shoreline. With these lower velocities compared to earlier studies and, hence, larger transit times the basal accretion rate is estimated to be 4 mm a⁻¹ along a portion of the Vostok flowline. An assessment of the local accretion rate at Vostok Station using the observed geodetic quantities yields an accretion rate in the same order of magnitude. Furthermore, the comparison of our geodetic observations with results inferred from ice-penetrating radar data indicates that the ice flow may not have changed significantly for several thousand years.     

Palaeointensity results from Ethiopian basalts: implications for the Oligocene geomagnetic field strength

From a large collection of Ethiopian flood basalts (~30  Myr old) sampled for magnetostratigraphy, ⁴⁰Ar/³⁹Ar geochronology and geochemical investigations, 47 samples were selected in order to test their suitability for Thellier palaeointensity experiments. Only 17 samples from eight individual flows yielded reliable palaeointensity estimates, with flow-mean virtual dipole moments ranging from 3.0 to 10.5 × 10²²  A  m² .
  A critical review of the Oligocene palaeointensity data set, including these new Ethiopian data, indicates an Oligocene mean virtual dipole moment of 5.1 ± 2.5 × 10²²  A  m² for the complete data set. After applying mild selection criteria, the reduced data set yields a mean value of only 4.6 ± 1.9 × 10²²  A  m² . This value is significantly lower than the present-day field strength but is higher than the Mesozoic dipole low mean field. This low Oligocene field might be in agreement with the high palaeosecular variation and rather high non-dipole field invoked around 30  Ma. However, the Oligocene data set is largely dependent on the palaeointensity determinations from Armenia, obtained mainly from baked contacts, which show amazingly low dispersion at both flow and between-flow levels. More data are needed to reduce the weight of these determinations on the mean value and avoid a possible bias.     

Temporal and spatial variations of post-seismic deformation following the 1999 Chi-Chi, Taiwan earthquake

We use GPS displacements collected in the 15 months after the 1999 Chi-Chi, Taiwan earthquake  ( M _w 7.6)  to evaluate whether post-seismic deformation is better explained by afterslip or viscoelastic relaxation of the lower crust and upper mantle. We find that all viscoelastic models tested fail to fit the general features in the post-seismic GPS displacements, in contrast to the satisfactory fit obtained with afterslip models. We conclude that afterslip is the dominant mechanism in the 15-month period, and invert for the space–time distribution of afterslip, using the Extended Network Inversion Filter. Our results show high slip rates surrounding the region of greatest coseismic slip. The slip-rate distribution remains roughly stationary over the 15-month period. In contrast to the limited coseismic slip on the décollement, afterslip is prominent there. Maximum afterslip of 0.57 m occurs downdip and to the east of the hypocentral region. Afterslip at hypocentral depths is limited to the southern part of the main shock rupture, with little or no slip on the northern section where coseismic slip was greatest. Whether this results from along strike variations in frictional properties or dynamic conditions that locally favour stable sliding is not clear. In general, afterslip surrounds the area of greatest coseismic slip, consistent with post-seismic slip driven by the main shock stress change. The total accumulated geodetic afterslip moment is  3.8 × 10¹⁹ N m  , significantly more than the seismic moment released by aftershocks,  6.6 × 10¹⁸ N m  . Afterslip and aftershocks appear to have different temporal evolutions and some spatial correlations, suggesting that aftershock rates may not be completely controlled by the rate of afterslip.