Permeability, overburden stress and friction coefficients of ODP Hole 190-1173A and 190-1174B

Permeability measured on three samples in a triaxial cell under effective confining pressure from 0.2 to 2.5 MPa ranges from 10**-18 to 10**-19 m**2. Overall, results indicate that permeability decreases with effective confining pressure up to 1.5 MPa; however, measurements at low effective pressure are too dispersed to yield a precise general relationship between permeability and pressure. When the effective pressure is increased from 1.5 to 2.5 MPa, permeability is roughly constant (~1-4 x 10**-19 m**2). Samples deformed in the triaxial cell developed slickenlined fractures, and permeability measurements were performed before and after failure. A permeability increase is observed when the sample fails under low effective confining pressure (0.2 MPa), but not under effective pressure corresponding to the overburden stress. Under isotropic stress conditions, permeability decrease related to fracture closure occurs at a relatively high effective pressure of ~1.5 MPa. Coefficients of friction on the fractures formed in the triaxial cell are ~0.4.

Measurements of velocity-dependent and slip-dependent frictional strength in laboratory friction experiments on natural samples from IODP Expedition 316 and ODP Hole 190-1174B

Slowslip forms part of the spectrum of fault behaviour between stable creep and destructive earthquakes. Slow slip occurs near the boundaries of large earthquake rupture zones and may sometimes trigger fast earthquakes. It is thought to occur in faults comprised of rocks that strengthen under fast slip rates, preventing rupture as a normal earthquake, or on faults that have elevated pore-fluid pressures. However, the processes that control slow rupture and the relationship between slow and normal earthquakes are enigmatic. Here we use laboratory experiments to simulate faulting in natural rock samples taken from shallow parts of the Nankai subduction zone, Japan, where very low-frequency earthquakes - a form of slow slip - have been observed.We find that the fault rocks exhibit decreasing strength over millimetre-scale slip distances rather than weakening due to increasing velocity. However, the sizes of the slip nucleation patches in our laboratory simulations are similar to those expected for the very lowfrequency earthquakes observed in Nankai. We therefore suggest that this type of fault-weakening behaviour may generate slow earthquakes. Owing to the similarity between the expected behaviour of slow earthquakes based on our data, and that of normal earthquakes during nucleation, we suggest that some types of slow slip may represent prematurely arrested earthquakes.

Code and readme of a MATLAB-based graphical user interface program for computing functionals of the geopotential

We present a novel graphical user interface program GrafLab (GRAvity Field LABoratory) for spherical harmonic synthesis (SHS) created in MATLAB®. This program allows to comfortably compute 38 various functionals of the geopotential up to ultra-high degrees and orders of spherical harmonic expansion. For the most difficult part of the SHS, namely the evaluation of the fully normalized associated Legendre functions (fnALFs), we used three different approaches according to required maximum degree: (i) the standard forward column method (up to maximum degree 1800, in some cases up to degree 2190); (ii) the modified forward column method combined with Horner's scheme (up to maximum degree 2700); (iii) the extended-range arithmetic (up to an arbitrary maximum degree). For the maximum degree 2190, the SHS with fnALFs evaluated using the extended-range arithmetic approach takes only approximately 2-3 times longer than its standard arithmetic counterpart, i.e. the standard forward column method. In the GrafLab, the functionals of the geopotential can be evaluated on a regular grid or point-wise, while the input coordinates can either be read from a data file or entered manually. For the computation on a regular grid we decided to apply the lumped coefficients approach due to significant time-efficiency of this method. Furthermore, if a full variance-covariance matrix of spherical harmonic coefficients is available, it is possible to compute the commission errors of the functionals. When computing on a regular grid, the output functionals or their commission errors may be depicted on a map using automatically selected cartographic projection.