Coal seams in the Waikato Coal Measures can vary significantly in thickness over distances of hundreds of meters to kilometers. Previously, the primary depositional controls on these variations have been inferred to be syn-depositional normal faulting and pre-depositional paleotopography. The data presented in support of these models are typically equivocal and which, if any, of these processes provide the principal control on the geometry and spatial distribution of coal seams in the Waikato Coal Region is uncertain. This study utilizes a large database of drill-logs, seismic-reflection lines and mine exposures in four areas (Huntly, Maramarua, North HuntlylWaikare and Rotowaro Coalfields) to test whether syn-depositional faulting and/or paleotopography influence coal seam architecture. These data were used to construct cross sections across faults and basement topography, which in turn, offer information on the relative timing of faulting and coal measure deposition, together with information on the spatial relations between seam thicknesses, faulting and paleotopography. Cross sections and isopach maps together with examination of spatial and temporal variations in fault displacements reveal that syn-depositional normal faulting had little or no impact on the deposition of the Waikato Coal Measures. Only in the Maramarua study area was any evidence found of fault control on coal measure deposition, with the Landing Fault accruing displacement between deposition of the Kupalrupa Seam and the end of coal measure sedimentation. The vast majority of faults in the Waikato Coalfield were, however, active following coal measure deposition. For example, the Foote, Kimihia and Pukekapia faults show evidence of displacement accrual, which commenced during deposition of the Mangakotuku Formation (37-35 Ma BP). The duration of this episode of faulting is difficult to determine, but may have ceased about 30 Ma ago. In addition, a number of faults (e.g. Foote Fault) display evidence oflate stage extension during the last 5 Ma. Given the lack of stratigraphic evidence for fault displacements during deposition of coal measures, it is suggested that the Mangakotuku and Waipuna basement scarps are erosional rather than tectonic features. Cross sections, together with structure contour and isopach maps in each of the four study areas examined, indicate that basement topography was the dominant control on the spatially variable accumulation of peat. These data show coal seams both thinning into, and away from, topographic lows. To account for this observation a model is proposed, in which peat accumulation is controlled by basement relief and sediment supply to parts of the depositional system. In the model it is postulated that the Waikato Coal Measures depositional system was a continuum between two end members. In one end member, with a high sediment supply, sediment is channeled into the lowest topographic areas and peat accumulates mainly on topographic highs. In the other end member, with little or no sediment supply, peat accumulates to its greatest thickness in areas of relatively low topography, in addition to on basement ridges. In the Rotowaro and North Huntly/Waikare study areas, the thickest peat developed on basement highs and the lows acted as a conduit for sedimentation. On basement highs, peat mires were largely sheltered from clastic sediment influx. In the Huntly East and Maramarua study areas, the thickest peat accumulated in basement lows, with comparable clastic sedimentation in highs and lows. The proposed model has application to other coalfields where peat accumulated on an undulating topographic surface and sediment supply was channelised. Prediction of coal seam thickness, as well as lithological types, is crucial in coal exploration and development. The methodology developed and employed in this study can be applied to other basins to access and model coal and clastic sediment distribution.

An integrated field and structural analysis of the Mitten Park fault-fold structure, northwestern Colorado and northeastern Utah, examines its structural origin. The Mitten Park structure is a modified fault-propagation-fold. This new model incorporates faulting, folding, and fracturing in one deformational event to produce the Mitten Park fault and associated monocline. The largest structure in the study area is the Mitten Park fault and associated monocline. The Mitten Park fault has approximately 127 meters (415 feet) of net slip, strikes S28°W and dips 55°WNW. In the footwall, net shortening was accommodated by reverse and normal faulting. Faulting was the result of northwest-southeast directed shortening. Reverse faulting accommodated the majority of the fault-related strain along the fault's trace and resulted in net shortening. However, normal faults in the overturned limb of the footwall of the Mitten Park fault also accommodated northwest-southeast directed shortening. Folds in the study area are asymmetrical and statistically cylindrical in both the footwall and the hanging wall. Folding facilitates northwest-southeast directed shortening. There is a direct correlation between changes in the strike and dip of the fault plane and changes in the trend and plunge of fold axis in the footwall. Fracture orientations show no significant variation in geometry from hanging wall to footwall. Fracture intensity increases with proximity to the Mitten Park fault. Balanced cross sections of the Mitten Park area use a modified fault-propagation- fold model and are also constrained by field observations and interlimb angles of folds. Total shortening in the study area is 13.5% and was accommodated by the hanging wall, the footwall, and the Mitten Park fault. The hanging wall accommodated 70.8% of total shortening, the footwall accommodated 14.9% of total shortening, and the Mitten Park fault accommodated 14.3% of total shortening. The significant amount of strain in the footwall of the fault is different from classical models of fault-propagation-folds, which depict a rigid undeformed footwall.

In this thesis, I use multiple techniques to investigate the active faulting and deformation of the Altay Mountains, Western Mongolia. The Altay are an intracontinental transpressional mountain range, which are deforming in the far-field of the India-Asia collision. An anastomosing network of dextral faults strikes NNW-SSE, and accommodates NE-SW oriented shortening by rotating anticlockwise about vertical axes. I begin by characterising the Altay faults, and add to what is already known about their surface expression with new observations of active faulting and three previously undescribed ancient earthquake ruptures. I use ¹⁰Be cosmogenic dating and uranium-series dating on pedogenic carbonate to estimate the average Quaternary rate of slip for two of the major fault zones in the Altay. The slip rate on the Ölgiy fault is constrained to 0.3-2.1 mm/yr^-1. Results from the Hovd fault are ambiguous, demonstrating the complications encountered with application of Quaternary dating techniques. I measure palaeomagnetic directions from Cretaceous to Pliocene-aged sediments in the eastern Altay to constrain the degree of anticlockwise rotation. Results from thermal demagnetisation of specimens indicate that the eastern Altay has not undergone significant rotation, in contrast with previous studies from the Siberian Altay that reveal almost 40 degrees of anticlockwise rotation. This suggests that the eastern-most Altay fault is too young to have experienced significant rotation, or is kinematically different from the Siberian Altay. I apply apatite fission track (AFT) dating and track length modeling to the central Altay. Results from AFT dating show rapid cooling in the late Cretaceous due to the distal assembly of Central Asia, suggesting that there was pre-existing topography at the start of the Late Cenozoic phase of deformation, the timing of which is constrained to have initiated at least 20 Myr ago. My work demonstrates that combining results from techniques that cover a variety of time scales quantifies the evolution of active faulting and deformation in the region.

Over the past few decades, remote sensing has emerged as a powerful tool for studying active faulting in continental regions. However, the commonly used remote sensing techniques, including radar interferometry, visual inspection of imagery, and image matching, cannot measure three-dimensional (3D) surface displacements in earthquakes, limiting our ability to investigate faulting. The improvement of very high-resolution (VHR) optical imaging systems (stereo in particular) in recent years has made it possible for earth scientists to measure 3D surface deformation remotely. In this thesis, I contribute to assessing the capability of VHR optical imagery, by determining earthquake deformation from four different types of earthquakes (different in sense of slip and date of the event). In the case of the 2010 M_w 7.2 El Mayor-Cucapah, Mexico earthquake, I show that digital elevation models (DEMs) derived from Pleiades stereo imagery are comparable to light detection and ranging (LiDAR) surveys, and differencing pre- and post-earthquake DEMs can measure 3D displacements, which will be very useful for studying future earthquakes. For the 2013 M_w 7.7 Balochistan, Pakistan earthquake, I determine the vertical motion from a post-earthquake Pleiades DEM and find constant fault kinematics throughout the Late Quaternary. This study has resolved a current controversy of the Balochistan earthquake, in which it has been argued that kinematics of the Hoshab fault switches between strike-slip and dip-slip. Applying historical aerial, KH-9 Hexagon spy satellite, SPOT-2 and modern SPOT-6 images to the 1978 M_w 7.3 Tabas-e-Golshan earthquake, I measure the coseismic and postseismic displacements, and show that the Tabas fold system in eastern Iran may exhibit characteristic slip behaviour. Combining Pleiades imagery, fieldwork and geological dating techniques, I determine slip in the 1556 Huaxian earthquake in China and the recurrence interval for similar events. These examples demonstrate the usefulness of high-resolution optical imagery in investigating past and present earthquakes.

Structures within very large displacement, mature fault zones, such as the North Branch San Gabriel Fault (NBSGF), are the product of a complex combination of processes. Off-fault damage within a damage zone and first-order geometric asperities, such as bends and steps, are thought to affect earthquake rupture propagation and energy radiation, but the effects are not completely understood. We hypothesize that the rate of accumulation of new damage decreases as fault maturity increases, and damage magnitude saturates in very large displacement faults. Nonetheless, geometric irregularities in the fault surface may modify damage zone characteristics. Accordingly, we seek to investigate the orientation, kinematics, and density of features at a range of scales within the damage zone adjacent to an abrupt 13 degree bend over 425 m in the NBSGF in order to constrain the relative role of the initiation of new damage versus the reactivation of preexisting damage adjacent to a bend. Field investigation and microstructural study focused on structural domains before, within, and after the fault bend on both sides of the fault. Subsidiary fault fabrics are similar in all domains outside the bend, which suggests a steady state fracture density and orientation distribution is established on the straight segments before and after the bend. The density of fractures within and outside the bend is similar; however, subsidiary fault orientations and kinematics are different within the bend relative to the straight segments. These observations are best explained by relatively low rates of damage generation relative to rates of fault reactivation during the later stages of faulting on the NBSGF, and that damage zone kinematics is reset as the host rock moves into the bend and again upon exiting the bend. Consequently, significant energy released during earthquake unloading can be dissipated by reactivation and slip on existing fractures in the damage zone, particularly adjacent to mesoscale faults. Thus, areas of heightened reactivation of damage, such as adjacent to geometric irregularities in the fault surface, could affect earthquake rupture dynamics.

Late Ordovician Red River strata of southeastern Saskatchewan were deposited in a broad epicontinental sea. In the lower part, the Yeoman and Herald formations comprise two cycles of carbonateevaporite sequences. Although these units possess an overall layer-cake aspect, thickness variations especially in the Herald Formation show that accumulation was affected by syndepositional flexure, differential subsidence and displacement of fault-bounded blocks. The mainly laminated dolomudstones and anhydrites of the Lake Alma and Coronach members of the Herald Formation were deposited under relatively tranquil conditions. These units host different kinds of synsedimentary deformation features, interpreted to have been induced by earthquakes generated because of movements along basement faults thought to have been oriented orthogonally NE−SW and NW−SE. The low-energy environmental setting was conducive to preserving these features, referred to as seismites.<p> The variety of seismites in the Herald Formation is related to the varying rheology of the carbonate or evaporite sediment, as well as shaking intensity. Brittle and quasi-brittle failure is represented by faults, microfaults, shear-vein arrays and pseudo-intraclastic breccias, mostly in dolomudstones which must have been stiff at the time of deformation. Plastic behaviour is recorded by soft-sediment deformation, comprising a family of features that includes loop bedding, folded laminae and convolute bedding. Indeed, these structures in enterolithic anhydrite are more reasonably interpreted as due to deformation than crystal growth, volume expansion and displacement, the more usual explanations. Sediment shrinkage and concomitant fluidization are recorded by dikelets containing injected carbonate mud or granular gypsum, the latter now preserved as anhydrite. Evidence for wholesale liquefaction, however, was not observed. These rheological differences were caused by the primary nature of the sediment plus modifications due to early diagenesis and burial confinement. Shaking intensity is difficult to gauge, but it is presumed that a minimum of VI on the modified Mercalli scale was required to produce these features. Consequently, shaking of lesser magnitude was probably not recorded.<p> The geographic distribution of seismites should reflect the location of basement faults presumed to have been active during deposition, and indeed there is a concentration adjacent to the known location of syndepositonal fault lineaments. In addition, the stratigraphic distribution of seismites records higher frequencies of activity of these same faults. These distributions show that earthquake-induced ground motion was common during deposition of the Lake Alma Member in southeastern Saskatchewan but less so during deposition of the Coronach Member.<p> Seismites serve as proxies for the activity of relatively nearby syndepositional faults making up the tectonic fabric of sedimentary basins. They also point to basement features that, if re-activated, can induce fracture porosity or influence subsurface fluid flow. Syndepositional tectonism undoubtedly had a much more profound influence on many successions than is presently accepted, and its effects are more widespread than currently appreciated.

The eastern California shear zone (ECSZ) and Walker Lane belt represent an important inland component of the Pacific-North America plate boundary. Current geodetic data indicate accumulation of transtensional shear at a rate of ~9.2 ± 0.3 mm/yr across the region, more than double the total geologic rate (<3.5 mm/yr) for faults in the northern ECSZ over the late Pleistocene [Bennett et al., 2003, Kirby et al., 2006, Lee et al., 2009, Frankel et al., 2007]. Unraveling the strain puzzle of the Walker Lane is therefore essential to understanding both how deformation is distributed through the lithosphere along this transtensional part of the Pacific-North America plate boundary and how the plate boundary is evolving through time. The observed mismatch between geodetic and geologic slip rates in the central Walker Lane is characteristic of other active tectonic settings, including the nearby Mojave segment of the ECSZ [Oskin et al., 2008] and the Altyn Tagh fault in China [Cowgill, 2007]. In each case, lack of fault slip data spanning multiple temporal and spatial scales hinders interpretation of fault interactions and their implications for lithospheric dynamics. The discrepancy between geodetic and geologic slip rates in the central Walker Lane indicates that if strain rates have remained constant since the late Pleistocene [e.g. Frankel et al., in press], then the "missing" strain is distributed on structures other than the two major dextral faults at this latitude (Death Valley-Fish Lake Valley fault and White Mountains fault). Otherwise the region could presently be experiencing a strain transient similar to that of the nearby Mojave section of the ECSZ [e.g., Oskin et al., 2008], or the rate of strain accumulation could actually increasing over the late Pleistocene [e.g. Reheis and Sawyer, 1997; Hoeft and Frankel, 2010]. The Silver Peak-Lone Mountain extensional complex (SPLM), to which the Clayton Valley faults belong, is the prime candidate to account for the "missing" strain. The down-to-the-northwest orientation of the SPLM faults makes them the most kinematically suitable structures to accommodate the regional pattern of NW-SE dextral shear. We use differential GPS to measure fault offset and terrestrial cosmogenic nuclide (TCN) geochronology to date offset landforms. Using these tools, we measure extension rates that are time-invariant, ranging from 0.1 ± 0.1 to 0.3 ± 0.1 mm/yr for fault dips of 30° and 60°. These rates are not high enough to account for the discrepancy between geologic and geodetic data in the ECSZ-Walker Lane transition zone. Based on geologic mapping and previously published geophysical data [Davis, 1981; Zampirro, 2005], deformation through Clayton Valley appears to be very widely-distributed. The diffuse nature of deformation leads to geologic slip rates that are underestimated due to the effects of off-fault deformation and unrecognized fault strands. Our results from Clayton Valley suggest that the discrepancy between geodetic and geologic strain rates at the latitude of the northern ECSZ is a result of long-term geologic rates that are underestimated. If the true geologic rates could be calculated, they would likely be significantly higher and therefore in closer agreement with geodetic data, as is the case everywhere else in the ECSZ north of the Garlock fault [Frankel et al., 2007a, in press; Kirby et al., 2008; Lee et al., 2009a].

Paleochannels are a major cause of roof failure in underground coal mines in southeastern Kentucky. Models that predict the location and geometry of paleochannels are essential to assist in mine planning and development. Data from approximately 506 coal exploration drill holes were subjected to second-order trend-surface analysis to identify stacking or offsetting relationships between sandstone bodies in adjacent stratigraphic intervals. The stacking of sandstone bodies within adjacent intervals suggests the presence of synsedimentary faulting. This model suggests that continued movement along the faults created topographic lows attracted paleodrainages and accommodated thick accumulations of sandstone in approximately the same areas through time. Trend-surface residuals analysis successfully located areas of potential synsedimentary faulting within the study area. An additional 7,189 elevation data points for the top of the Newman Limestone, interpreted from oil and gas records, were utilized to locate sub-Pennsylvanian System faults within the study area. The correlation between faults associated with the coal measures identified using second-order trend-surface analysis and faults affecting the Newman Limestone suggests Pennsylvanian synsedimentary faults were preceded by older Paleozoic fault movement. The greater availability of oil and gas subsurface data makes this relationship an important tool for predicting locations of fault-controlled coal measure paleochannels.

abstract: The tectonic significance of the physiographic transition from the low-relief Tibetan plateau to the high peaks, rugged topography and deep gorges of the Himalaya is the source of much controversy. Some workers have suggested the transition may be structurally controlled (e.g. Hodges et al., 2001), and indeed, the sharp change in geomorphic character across the transition strongly suggests differential uplift between the Himalayan realm and the southernmost Tibetan Plateau. Most Himalayan researchers credit the South Tibetan fault system (STFS), a family of predominantly east-west trending, low-angle normal faults with a known trace of over 2,000 km along the Himalayan crest (e.g. Burchfiel et al., 1992), with defining the southern margin of the Tibetan Plateau in the Early Miocene. Inasmuch as most mapped strands of the STFS have not been active since the Middle Miocene (e.g., Searle & Godin, 2003), modern-day control of the physiographic transition by this fault system seems unlikely. However, several workers have documented Quaternary slip on east-west striking, N-directed extensional faults, of a similar structural nature but typically at a different tectonostratigraphic level than the principal STFS strand, in several locations across the range (Nakata, 1989; Wu et al., 1998; Hurtado et al., 2001). In order to explore the nature of the physiographic transition and determine its relationship to potential Quaternary faulting, I examined three field sites: the Kali Gandaki valley in central Nepal (~28&#730;39'54"N; 83&#730;35'06"E), the Nyalam region of south-central Tibet (28°03'23.3"N, 86°03'54.08"E), and the Ama Drime Range in southernmost Tibet (87º15'-87º50'E; 27º45'-28º30'N). Research in each of these areas yielded evidence of young faulting on structures with normal-sense displacement in various forms: the structural truncation of lithostratigraphic units, distinctive fault scarps, or abrupt changes in bedrock cooling age patterns. These structures are accompanied by geomorphic changes implying structural control, particularly sharp knickpoints in rivers that drain from the Tibetan Plateau, across the range crest, and down through the southern flank of the Himalaya. Collectively, my structural, geomorphic, and thermochronometric studies confirm the existence of extensional structures near the physiographic transition that have been active more recently than 1.5 Ma in central Nepal, and over the last 3.5 Ma in south-central Tibet. The structural history of the Ama Drime Range is complex and new thermochronologic data suggest multiple phases of E-W extension from the Middle Miocene to the Holocene. Mapping in the accessible portions of the range did not yield evidence for young N-S extension, although my observations do not preclude such deformation on structures south of the study area. In contrast, the two other study areas yielded direct evidence that Quaternary faulting may be controlling the position and nature of the physiographic transition across the central Tibetan Plateau-Himalaya orogenic system. / Dissertation/Thesis / Ph.D. Geological Sciences 2012

East Matagorda Peninsula in southwestern Texas is characterized geologically by active, regional-scale and near-surface growth faulting. Decimeter scale (up to 0.42 m) vertical displacement was recorded at the study site over a period of four years, not believed to be associated with growth faulting. This research tested the hypotheses that fault slip rates were correlated with sediment accumulation rates, and that the observed vertical displacement was produced by shrink-and-swell clays in near surface sediments. To quantify sediment accumulation rates, a suite of radionuclides (7Be, 137Cs, and 210Pb) were used. To understand the effects of shrink-and-swell clays, analyses including particle size distribution, X-ray diffraction (XRD), and Fourier transform infrared spectroscopy (FTIR) were completed. Additionally, the free swell index test (FSI) was used to record the swelling potential of the sediment. Strong correlation (R2 = 0.99) indicates coupling between mean fault slip rates and mean sediment accumulation rates. Near surface sediment clay size fraction percentages ranged from 0.96 - 6.26% containing more than 90% smectite. Based on FSI results, maximum volume change in the top six cm was determined to be 208%. The presence and behavior of shrink-and-swell clay minerals in the region is an important contributor to the vertical displacement observed.