USGS Science⌕ Search

Geology topics

C.A. Williams

Publications and source records attributed to C.A. Williams.

4 recordsLinked to original sources

A domain decomposition approach to implementing fault slip in finite-element models of quasi-static and dynamic crustal deformation

We employ a domain decomposition approach with Lagrange multipliers to implement fault slip in a finite-element code, PyLith, for use in both quasi-static and dynamic crustal deformation applications. This integrated approach to solving both quasi-static and dynamic simulations leverages common finite-element data structures and implementations of various boundary conditions, discretization schemes, and bulk and fault rheologies. We have developed a custom preconditioner for the Lagrange multiplier portion of the system of equations that provides excellent scalability with problem size compared to conventional additive Schwarz methods. We demonstrate application of this approach using benchmarks for both quasi-static viscoelastic deformation and dynamic spontaneous rupture propagation that verify the numerical implementation in PyLith.

Journal of Geophysical Research B: Solid Earth↗

Fault locking, block rotation and crustal deformation in the Pacific Northwest

We interpret Global Positioning System (GPS) measurements in the northwestern United States and adjacent parts of western Canada to describe relative motions of crustal blocks, locking on faults and permanent deformation associated with convergence between the Juan de Fuca and North American plates. To estimate angular velocities of the oceanic Juan de Fuca and Explorer plates and several continental crustal blocks, we invert the GPS velocities together with seafloor spreading rates, earthquake slip vector azimuths and fault slip azimuths and rates. We also determine the degree to which faults are either creeping aseismically or, alternatively, locked on the block-bounding faults. The Cascadia subduction thrust is locked mainly offshore, except in central Oregon, where locking extends inland. Most of Oregon and southwest Washington rotate clockwise relative to North America at rates of 0.4–1.0 ° Myr −1 . No shear or extension along the Cascades volcanic arc has occurred at the mm/yr level during the past decade, suggesting that the shear deformation extending northward from the Walker Lane and eastern California shear zone south of Oregon is largely accommodated by block rotation in Oregon. The general agreement of vertical axis rotation rates derived from GPS velocities with those estimated from palaeomagnetic declination anomalies suggests that the rotations have been relatively steady for 10–15 Ma. Additional permanent dextral shear is indicated within the Oregon Coast Range near the coast. Block rotations in the Pacific Northwest do not result in net westward flux of crustal material—the crust is simply spinning and not escaping. On Vancouver Island, where the convergence obliquity is less than in Oregon and Washington, the contractional strain at the coast is more aligned with Juan de Fuca—North America motion. GPS velocities are fit significantly better when Vancouver Island and the southern Coast Mountains move relative to North America in a block-like fashion. The relative motions of the Oregon, western Washington and Vancouver Island crustal blocks indicate that the rate of permanent shortening, the type that causes upper plate earthquakes, across the Puget Sound region is 4.4 ± 0.3 mm yr −1 . This shortening is likely distributed over several faults but GPS data alone cannot determine the partitioning of slip on them. The transition from predominantly shear deformation within the continent south of the Mendocino Triple Junction to predominantly block rotations north of it is similar to changes in tectonic style at other transitions from shear to subduction. This similarity suggests that crustal block rotations are enhanced in the vicinity of subduction zones possibly due to lower resisting stress.

Geophysical Journal International↗

Ice-affected streamflow records using tracer-dilution discharge methods

Accurate ice-affected streamflow records are difficult to obtain for several reasons. Problems measuring stage, variable backwater conditions, access limitations in wintertime, and problems measuring flowing water under ice cover all contribute to make ice-affected streamflow records less accurate than open-channel streamflow records. The inaccuracy of ice-affected streamflow records is particularly troublesome for small streams where Instream-Flow water rights exist. The Colorado Water Conservation Board uses these water rights to protect in-stream aquatic communities. In January and February 2002, the U.S. Geological Survey, in cooperation with the Colorado Water Conservation Board, conducted an experiment using a sodium chloride tracer to determine streamflow under ice cover. The purpose of this study is to determine the usefulness and accuracy of ice-affected streamflow records using a sodium chloride tracer that was automatically injected into the stream. The tracer was injected at two gaging stations once per day for up to 25 days. Multiple-parameter water-quality sensors at the two gaging stations monitored background and peak tracer concentrations and conductance. These data were used to determine discharge at each site. A comparison of current-meter measurements to tracer-dilution discharge measurements shows an underestimation of discharge due to inaccuracy of current-meter measurements with ice cover and inconsistent tracer-pump rates caused by partial freezing of the tracer solution in the injection lines.

Conference Paper↗