Frictional properties of faults : from observation on the Longitudinal Valley Fault, Taiwan, to dynamic simulations

Faults can slip either aseismically or through episodic seismic ruptures, but we still do not understand the factors which determine the partitioning between these two modes of slip. This challenge can now be addressed thanks to the dense set of geodetic and seismological networks that have been deployed in various areas with active tectonics. The data from such networks, as well as modern remote sensing techniques, indeed allow documenting of the spatial and temporal variability of slip mode and give some insight. This is the approach taken in this study, which is focused on the Longitudinal Valley Fault (LVF) in Eastern Taiwan. This fault is particularly appropriate since the very fast slip rate (about 5 cm/yr) is accommodated by both seismic and aseismic slip. Deformation of anthropogenic features shows that aseismic creep accounts for a significant fraction of fault slip near the surface, but this fault also released energy seismically, since it has produced five M_w>6.8 earthquakes in 1951 and 2003. Moreover, owing to the thrust component of slip, the fault zone is exhumed which allows investigation of deformation mechanisms. In order to put constraint on the factors that control the mode of slip, we apply a multidisciplinary approach that combines modeling of geodetic observations, structural analysis and numerical simulation of the "seismic cycle". Analyzing a dense set of geodetic and seismological data across the Longitudinal Valley, including campaign-mode GPS, continuous GPS (cGPS), leveling, accelerometric, and InSAR data, we document the partitioning between seismic and aseismic slip on the fault. For the time period 1992 to 2011, we found that about 80-90% of slip on the LVF in the 0-26 km seismogenic depth range is actually aseismic. The clay-rich Lichi M\'elange is identified as the key factor promoting creep at shallow depth. Microstructural investigations show that deformation within the fault zone must have resulted from a combination of frictional sliding at grain boundaries, cataclasis and pressure solution creep. Numerical modeling of earthquake sequences have been performed to investigate the possibility of reproducing the results from the kinematic inversion of geodetic and seismological data on the LVF. We first investigate the different modeling strategy that was developed to explore the role and relative importance of different factors on the manner in which slip accumulates on faults. We compare the results of quasi dynamic simulations and fully dynamic ones, and we conclude that ignoring the transient wave-mediated stress transfers would be inappropriate. We therefore carry on fully dynamic simulations and succeed in qualitatively reproducing the wide range of observations for the southern segment of the LVF. We conclude that the spatio-temporal evolution of fault slip on the Longitudinal Valley Fault over 1997-2011 is consistent to first order with prediction from a simple model in which a velocity-weakening patch is embedded in a velocity-strengthening area.

Examination of the flow rheological and textural properties of polymer gels composed of poly(methylvinylether-co-maleic anhydride) and poly(vinylpyrrolidone): rheological and mathematical interpretation of textural parameters.

The purpose of this study was to mathematically characterize the effects of defined experimental parameters (probe speed and the ratio of the probe diameter to the diameter of sample container) on the textural/mechanical properties of model gel systems. In addition, this study examined the applicability of dimensional analysis for the rheological interpretation of textural data in terms of shear stress and rate of shear. Aqueous gels (pH 7) were prepared containing 15% w/w poly(methylvinylether-co-maleic anhydride) and poly(vinylpyrrolidone) (PVP) (0, 3, 6, or 9% w/w). Texture profile analysis (TPA) was performed using a Stable Micro Systems texture analyzer (model TA-XT 2; Surrey, UK) in which an analytical probe was twice compressed into each formulation to a defined depth (15 mm) and at defined rates (1, 3, 5, 8, and 10 mm s-1), allowing a delay period (15 s) between the end of the first and beginning of the second compressions. Flow rheograms were performed using a Carri-Med CSL2-100 rheometer (TA Instruments, Surrey, UK) with parallel plate geometry under controlled shearing stresses at 20.0°?±?0.1°C. All formulations exhibited pseudoplastic flow with no thixotropy. Increasing concentrations of PVP significantly increased formulation hardness, compressibility, adhesiveness, and consistency. Increased hardness, compressibility, and consistency were ascribed to enhanced polymeric entanglements, thereby increasing the resistance to deformation. Increasing probe speed increased formulation hardness in a linear manner, because of the effects of probe speed on probe displacement and surface area. The relationship between formulation hardness and probe displacement was linear and was dependent on probe speed. Furthermore, the proportionality constant (gel strength) increased as a function of PVP concentration. The relationship between formulation hardness and diameter ratio was biphasic and was statistically defined by two linear relationships relating to diameter ratios from 0 to 0.4 and from 0.4 to 0.563. The dramatically increased hardness, associated with diameter ratios in excess of 0.4, was accredited to boundary effects, that is, the effect of the container wall on product flow. Using dimensional analysis, the hardness and probe displacement in TPA were mathematically transformed into corresponding rheological parameters, namely shearing stress and rate of shear, thereby allowing the application of the power law (??=?k?n) to textural data. Importantly, the consistencies (k) of the formulations, calculated using transformed textural data, were statistically similar to those obtained using flow rheometry. In conclusion, this study has, firstly, characterized the relationships between textural data and two key instrumental parameters in TPA and, secondly, described a method by which rheological information may be derived using this technique. This will enable a greater application of TPA for the rheological characterization of pharmaceutical gels and, in addition, will enable efficient interpretation of textural data under different experimental parameters.

Using space geodesy to constrain variations in seismogenic behavior on subduction megathrusts

The concept of seismogenic asperities and aseismic barriers has become a useful paradigm within which to understand the seismogenic behavior of major faults. Since asperities and barriers can be thought of as defining the potential rupture area of large megathrust earthquakes, it is thus important to identify their respective spatial extents, constrain their temporal longevity, and to develop a physical understanding for their behavior. Space geodesy is making critical contributions to the identification of slip asperities and barriers but progress in many geographical regions depends on improving the accuracy and precision of the basic measurements. This thesis begins with technical developments aimed at improving satellite radar interferometric measurements of ground deformation whereby we introduce an empirical correction algorithm for unwanted effects due to interferometric path delays that are due to spatially and temporally variable radar wave propagation speeds in the atmosphere. In chapter 2, I combine geodetic datasets with complementary spatio-temporal resolutions to improve our understanding of the spatial distribution of crustal deformation sources and their associated temporal evolution – here we use observations from Long Valley Caldera (California) as our test bed. In the third chapter I apply the tools developed in the first two chapters to analyze postseismic deformation associated with the 2010 Mw=8.8 Maule (Chile) earthquake. The result delimits patches where afterslip occurs, explores their relationship to coseismic rupture, quantifies frictional properties associated with inferred patches of afterslip, and discusses the relationship of asperities and barriers to long-term topography. The final chapter investigates interseismic deformation of the eastern Makran subduction zone by using satellite radar interferometry only, and demonstrates that with state-of-art techniques it is possible to quantify tectonic signals with small amplitude and long wavelength. Portions of the eastern Makran for which we estimate low fault coupling correspond to areas where bathymetric features on the downgoing plate are presently subducting, whereas the region of the 1945 M=8.1 earthquake appears to be more highly coupled.

Characterization of diverse megathrust fault behavior related to seismic supercycles, Mentawai Islands, Sumatra

Long paleoseismic histories are necessary for understanding the full range of behavior of faults, as the most destructive events often have recurrence intervals longer than local recorded history. The Sunda megathrust, the interface along which the Australian plate subducts beneath Southeast Asia, provides an ideal natural laboratory for determining a detailed paleoseismic history over many seismic cycles. The outer-arc islands above the seismogenic portion of the megathrust cyclically rise and subside in response to processes on the underlying megathrust, providing uncommonly good illumination of megathrust behavior. Furthermore, the growth histories of coral microatolls, which record tectonic uplift and subsidence via relative sea level, can be used to investigate the detailed coseismic and interseismic deformation patterns. One particularly interesting area is the Mentawai segment of the megathrust, which has been shown to characteristically fail in a series of ruptures over decades, rather than a single end-to-end rupture. This behavior has been termed a seismic “supercycle.” Prior to the current rupture sequence, which began in 2007, the segment previously ruptured during the 14th century, the late 16th to late 17th century, and most recently during historical earthquakes in 1797 and 1833. In this study, we examine each of these previous supercycles in turn.

First, we expand upon previous analysis of the 1797–1833 rupture sequence with a comprehensive review of previously published coral microatoll data and the addition of a significant amount of new data. We present detailed maps of coseismic uplift during the two great earthquakes and of interseismic deformation during the periods 1755–1833 and 1950–1997 and models of the corresponding slip and coupling on the underlying megathrust. We derive magnitudes of Mw 8.7–9.0 for the two historical earthquakes, and determine that the 1797 earthquake fundamentally changed the state of coupling on the fault for decades afterward. We conclude that while major earthquakes generally do not involve rupture of the entire Mentawai segment, they undoubtedly influence the progression of subsequent ruptures, even beyond their own rupture area. This concept is of vital importance for monitoring and forecasting the progression of the modern rupture sequence.

Turning our attention to the 14th century, we present evidence of a shallow slip event in approximately A.D. 1314, which preceded the “conventional” megathrust rupture sequence. We calculate a suite of slip models, slightly deeper and/or larger than the 2010 Pagai Islands earthquake, that are consistent with the large amount of subsidence recorded at our study site. Sea-level records from older coral microatolls suggest that these events occur at least once every millennium, but likely far less frequently than their great downdip neighbors. The revelation that shallow slip events are important contributors to the seismic cycle of the Mentawai segment further complicates our understanding of this subduction megathrust and our assessment of the region’s exposure to seismic and tsunami hazards.

Finally, we present an outline of the complex intervening rupture sequence that took place in the 16th and 17th centuries, which involved at least five distinct uplift events. We conclude that each of the supercycles had unique features, and all of the types of fault behavior we observe are consistent with highly heterogeneous frictional properties of the megathrust beneath the south-central Mentawai Islands. We conclude that the heterogeneous distribution of asperities produces terminations and overlap zones between fault ruptures, resulting in the seismic “supercycle” phenomenon.

Loosely packed self-assembled monolayer of N-hexadecyl-3,6-di(p-mercaptophenylacetylene)carbazole on gold and its application in biomimetic membrane research

Dithiols of N-hexadecyl-3,6-di(p-mercaptophenylacetylene)carbazole (HDMC) have been synthesized and employed to form self-assembled monolayers (SAMs) on gold. One characteristic of the HDMC molecule is its peculiar molecular structure consisting of a large and rigid headgroup and a small and flexible alkyl-chain tail. HDMC adsorbates can attach to gold substrates by a strong Au-S bond with weak van der Waals interactions between the alkyl-chain tails, leading to a loosely packed hydrophobic SAM. In this way we can couple hybrid bilayer membranes (HBMs) to gold surfaces with more likeness to a cell bilayer than the conventional HBMs based on densely packed long-chain alkanethiol SAMs. The insulating properties and stability of the HDMC monolayer as well as the HDMC/lipid bilayer on gold have been investigated by electrochemical techniques including cyclic voltammetry and impedance spectroscopy. To test whether the quality of the bilayer is sufficiently high for biomimetic research, we incorporated the pore-forming protein a-hemolysin) and the horseradish peroxidase into the bilayers, respectively.

Self-lubricating silicone elastomer biomaterials

Silicone has a relatively high coefficient of friction and silicone medical devices therefore lack inherent lubricity, leading to pain on device insertion and potential tissue trauma. In this study, higher molecular weight tetra(alkoxy) silanes, particularly tetra(oleyloxy) silane, have been used as crosslinkers in the condensation cure of a hydroxy end-functionalised linear poly(dimethylsiloxane). The resulting elastomers displayed a persistent lubricous surface of oleyl alcohol, and coefficients of friction (static and dynamic) approaching zero. Chemical structures of the synthesised silanes and surface alcohol exudate were confirmed by nuclear magnetic resonance spectroscopy. Mechanical properties of the elastomers, which were chemically identical to conventionally cured systems, suggested that an 80/20 mixture of tetra(oleyloxy) silane and tetra(propoxysilane) gave the best compromise between desirable mechanical and frictional properties.

Formulation and characterisation of tetracycline-containing bioadhesive polymer networks designed for the treatment of periodontal disease.

This study described the drug release, rheological (dynamic and flow) and textural/mechanical properties of a series of formulations composed of 15% w/w polymethylvinylether-co-maleic anhydride (PMVE-MA), 0-9% w/w polyvinylpyrrolidone (PVP) and containing 1-5% w/w tetracycline hydrochloride, designed for the treatment of periodontal disease. All formulations exhibited pseudoplastic flow with minimal thixotropy. Increasing the concentration of PVP sequentially increased the zero-rate viscosity (derived from the Cross model) and the hardness and compressibility of the formulations (derived from texture profile analysis). These affects may be accredited to increased polymer entanglement and, in light of the observed synergy between the two polymers with respect to their textural and rheological properties, interaction between PVP and PMVE-MA. Increasing the concentration of PVP increased the storage and loss moduli yet decreased the loss tangent of all formulations, indicative of increased elastic behaviour. Synergy between the two polymers with respect to their viscoelastic properties was observed. Increased adhesiveness, associated with increased concentrations of PVP was ascribed to the increasing bioadhesion and tack of the formulations. The effect of increasing drug concentration on the rheological and textural properties was dependent on PVP concentration. At lower concentrations (0, 3% w/w) no effect was observed whereas, in the presence of 9% w/w PVP, increasing drug concentration increased formulation elasticity, zero rate viscosity, hardness and compressibility. These observations were ascribed to the greater mass of suspended drug in formulations containing the highest concentration of PVP. Drug release from formulations containing 6 and 9% PVP (and 5% w/w drug) was prolonged and swelling/diffusion controlled. Based on the drug release, rheological and textural properties, it is suggested that the formulation containing 15% w/w PMVE-MA, 6% w/w PVP and tetracycline hydrochloride (5% w/w) may be useful for the treatment of periodontal disease.

Physicochemical characterisation of injectable polymeric gels for use in the treatment of periodontal disease

The physicochemical characteristics of the injectable polymeric gels for use in the treatment of periodontal disease were investigated. The hardness, compressibility, and mucoadhesive properties of the gel were determined using a TA-XT2 Texture analyzer. The effect of of polymer concentration on the various viscoelastic, textural, bioadhesive properties and drug release were also analyzed using multifactorial analysis of variance. It was found that increased polymer concentration significantly increased gel structure, reduced polymer chain mobility and subsequently decreased the swelling/erosion and diffusion properties of the gel networks.

Application of fibre assemblies as damping elements in the automotive industry

This investigation aims to characterise the damping properties of the nonwoven materials with potential applications in automotive and aerospace industry. Nonwovens are a popular choice for many applications due to their relatively low manufacturing cost and unique properties. It is known that nonwovens are efficient energy dispersers for certain applications such as acoustic damping and ballistic impact. It is anticipated that these energy absorption properties could eventually be used to provide damping for mechanical vibrations. However the behaviour of nonwovens under dynamic load and vibration has not been investigated before. Therefore we intend to highlight these aspects of the behaviour of the nonwovens through this research. In order to obtain an insight to the energy absorption properties of the nonwoven fabrics, a range of tests has been performed. Forced vibration of the cantilever beam is used to explore damping over a range of resonance modes and input amplitudes. The tests are conducted on aramid, glass fibre and polyester fabrics with a range of area densities and various coatings. The tests clarified the general dynamic behaviour of the fabrics tested and the possible response in more real application condition as well. The energy absorption in both thickness and plane of the fabric is tested. The effects of the area density on the results are identified. The main absorption mechanism is known to be the friction. The frictional properties are improved by using a smaller fibre denier and increasing fibre length, this is a result of increasing contact surface between fibres. It is expected the increased friction result in improving damping. The results indicate different mechanism of damping for fiber glass fabrics compared to the aramid fabrics. The frequency of maximum efficiency of damping is identified for the fabrics tested. These can be used to recommend potential applications.

Use of fused airborne scanning laser altimetry and digital map data for urban flood modelling

Flood modelling of urban areas is still at an early stage, partly because until recently topographic data of sufficiently high resolution and accuracy have been lacking in urban areas. However, Digital Surface Models (DSMs) generated from airborne scanning laser altimetry (LiDAR) having sub-metre spatial resolution have now become available, and these are able to represent the complexities of urban topography. The paper describes the development of a LiDAR post-processor for urban flood modelling based on the fusion of LiDAR and digital map data. The map data are used in conjunction with LiDAR data to identify different object types in urban areas, though pattern recognition techniques are also employed. Post-processing produces a Digital Terrain Model (DTM) for use as model bathymetry, and also a friction parameter map for use in estimating spatially-distributed friction coefficients. In vegetated areas, friction is estimated from LiDAR-derived vegetation height, and (unlike most vegetation removal software) the method copes with short vegetation less than ~1m high, which may occupy a substantial fraction of even an urban floodplain. The DTM and friction parameter map may also be used to help to generate an unstructured mesh of a vegetated urban floodplain for use by a 2D finite element model. The mesh is decomposed to reflect floodplain features having different frictional properties to their surroundings, including urban features such as buildings and roads as well as taller vegetation features such as trees and hedges. This allows a more accurate estimation of local friction. The method produces a substantial node density due to the small dimensions of many urban features.

Uses of airborne scanning laser altimetry for river flood modelling and coastal geomorphology

Two ongoing projects at ESSC that involve the development of new techniques for extracting information from airborne LiDAR data and combining this information with environmental models will be discussed. The first project in conjunction with Bristol University is aiming to improve 2-D river flood flow models by using remote sensing to provide distributed data for model calibration and validation. Airborne LiDAR can provide such models with a dense and accurate floodplain topography together with vegetation heights for parameterisation of model friction. The vegetation height data can be used to specify a friction factor at each node of a model’s finite element mesh. A LiDAR range image segmenter has been developed which converts a LiDAR image into separate raster maps of surface topography and vegetation height for use in the model. Satellite and airborne SAR data have been used to measure flood extent remotely in order to validate the modelled flood extent. Methods have also been developed for improving the models by decomposing the model’s finite element mesh to reflect floodplain features such as hedges and trees having different frictional properties to their surroundings. Originally developed for rural floodplains, the segmenter is currently being extended to provide DEMs and friction parameter maps for urban floods, by fusing the LiDAR data with digital map data. The second project is concerned with the extraction of tidal channel networks from LiDAR. These networks are important features of the inter-tidal zone, and play a key role in tidal propagation and in the evolution of salt-marshes and tidal flats. The study of their morphology is currently an active area of research, and a number of theories related to networks have been developed which require validation using dense and extensive observations of network forms and cross-sections. The conventional method of measuring networks is cumbersome and subjective, involving manual digitisation of aerial photographs in conjunction with field measurement of channel depths and widths for selected parts of the network. A semi-automatic technique has been developed to extract networks from LiDAR data of the inter-tidal zone. A multi-level knowledge-based approach has been implemented, whereby low level algorithms first extract channel fragments based mainly on image properties then a high level processing stage improves the network using domain knowledge. The approach adopted at low level uses multi-scale edge detection to detect channel edges, then associates adjacent anti-parallel edges together to form channels. The higher level processing includes a channel repair mechanism.

A comparison of pilling results tested with the ICI pill box and a new perspex™ pill box

Identical single jersey wool knit fabrics produced inconsistent results when tested on the same type of ICI pill box in 4 independence testing laboratories.  The hypothesis formulated to investigate this event was that continued pilling testing in the ICI pill box changes the frictional properties of the cork liners, resulting in the difference pilling ratings.  The aim of this study was to decrease the inconsistencies in the pilling results by replacing the conventional cork lined pill box with a Perspex^TM pill box.  The smae or similar pilling results found for several of the knit facots tested in the two difference pill boxes, while for others there were some significant differences found.  In addition, while using the Perspex^TM pill box observations of the movement of the pilling tubes during testing provided evidence of greater fabric-to-fabric abrasion, rather the expected high degree of fabric-to-cork abrasion (or fabric-to-Perspex^TM) during pilling testing.

Lateral strength of zero bond masonry walls subjected to wind loads

Masonry walls are usually laid with the individual masonry units along a course overlapping units in the course below. Commonly, the perpend joints in the course occur above the mid-points of the units below to form a ‘half-bond’ or above a third point to form a ‘third-bond’. The amount of this overlap has a profound influence on the strength of a wall supported on three or four sides, where lateral pressures from wind cause combined vertical and horizontal flexure. Where masonry units are laid with mortar joints, the torsional shear bond resistance between the mortar and overlapping units largely determines the horizontal flexural strength. If there is zero bond strength between units, then the horizontal flexural strength is derived from the frictional resistance to torsion on the overlapping bed-faces of the units. This thesis reports a theoretical and experimental investigation into the frictional properties of overlapping units when subjected to combinations of vertical and horizontal moments and vertical axial compression. These basic properties were used to develop a theory to predict the lateral strength of walls supported on two, three or four sides. A plastic theory of behaviour was confirmed by experiment. The theory was then used to determine maximum unbraced panel sizes for particular boundary conditions. Design charts were developed to determine temporary bracing requirements for panels during construction.