Finite-Element Simulations of Earthquakes

This thesis discusses simulations of earthquake ground motions using prescribed ruptures and dynamic failure. Introducing sliding degrees of freedom led to an innovative technique for numerical modeling of earthquake sources. This technique allows efficient implementation of both prescribed ruptures and dynamic failure on an arbitrarily oriented fault surface. Off the fault surface the solution of the three-dimensional, dynamic elasticity equation uses well known finite-element techniques. We employ parallel processing to efficiently compute the ground motions in domains containing millions of degrees of freedom.

Using prescribed ruptures we study the sensitivity of long-period near-source ground motions to five earthquake source parameters for hypothetical events on a strike-slip fault (M_w 7.0 to 7.1) and a thrust fault (M_w 6.6 to 7.0). The directivity of the ruptures creates large displacement and velocity pulses in the ground motions in the forward direction. We found a good match between the severity of the shaking and the shape of the near-source factor from the 1997 Uniform Building Code for strike-slip faults and thrust faults with surface rupture. However, for blind thrust faults the peak displacement and velocities occur up-dip from the region with the peak near-source factor. We assert that a simple modification to the formulation of the near-source factor improves the match between the severity of the ground motion and the shape of the near-source factor.

For simulations with dynamic failure on a strike-slip fault or a thrust fault, we examine what constraints must be imposed on the coefficient of friction to produce realistic ruptures under the application of reasonable shear and normal stress distributions with depth. We found that variation of the coefficient of friction with the shear modulus and the depth produces realistic rupture behavior in both homogeneous and layered half-spaces. Furthermore, we observed a dependence of the rupture speed on the direction of propagation and fluctuations in the rupture speed and slip rate as the rupture encountered changes in the stress field. Including such behavior in prescribed ruptures would yield more realistic ground motions.

Understanding Micro- and Macro-Mechanics of Soil Liquefaction: A Necessary Step for Field-Scale Assessment

Liquefaction is a devastating instability associated with saturated, loose, and cohesionless soils. It poses a significant risk to distributed infrastructure systems that are vital for the security, economy, safety, health, and welfare of societies. In order to make our cities resilient to the effects of liquefaction, it is important to be able to identify areas that are most susceptible. Some of the prevalent methodologies employed to identify susceptible areas include conventional slope stability analysis and the use of so-called liquefaction charts. However, these methodologies have some limitations, which motivate our research objectives. In this dissertation, we investigate the mechanics of origin of liquefaction in a laboratory test using grain-scale simulations, which helps (i) understand why certain soils liquefy under certain conditions, and (ii) identify a necessary precursor for onset of flow liquefaction. Furthermore, we investigate the mechanics of liquefaction charts using a continuum plasticity model; this can help in modeling the surface hazards of liquefaction following an earthquake. Finally, we also investigate the microscopic definition of soil shear wave velocity, a soil property that is used as an index to quantify liquefaction resistance of soil. We show that anisotropy in fabric, or grain arrangement can be correlated with anisotropy in shear wave velocity. This has the potential to quantify the effects of sample disturbance when a soil specimen is extracted from the field. In conclusion, by developing a more fundamental understanding of soil liquefaction, this dissertation takes necessary steps for a more physical assessment of liquefaction susceptibility at the field-scale.

Charles Henry Gilbert (1859-1928), Naturalist-in-Charge: the 1906 North Pacific Expedition of the Steamer Albatross

Fishery science pioneers often faced challenges in their field work that are mostly unknown to modern biologists. Some of the travails faced by ichthyologist and, later, fishery biologist Charles Henry Gilbert (1859-1928) during his service as Naturalist-in-Charge of the North Pacific cruise ofthe U.S. Bureau of Fisheries Steamer Albatross in 1906, are described here, as are accomplishments of the cruise. The vessel left San Francisco, Calif., on 3 May 1906, just after the great San Francisco earthquake, for scientific exploration of waters of the Aleutian islands, Bering Sea, Kamchatka, Sakhalin, and Japan, returning to San Francisco in December. Because the expedition occurred just after the war between Japan and Russia of 1904-05 floating derelict mines in Japanese waters were often a menace. Major storms caused havoc in the region, and the captain of the Albatross, Lieutenant Commander LeRoy Mason Garrett (1857-1906), U.S.N., was lost at sea, apparently thrown from the vessel during a sudden storm on the return leg of the cruise. Despite such obstacles, Gilbert and the Albatross successfully completed their assigned chores. They occupied 339 dredging and 48 hydrographic stations, and discovered over 180 new species of fishes and many new species of invertebrates. The expedition's extensive biological collections spawned over 30 descriptive publications, some of which remain today as standards of knowledge.