Shear Wave Velocity Profiles at Sites Liquefied by the 1999 Kocaeli, Turkey Earthquake

This paper presents shear wave velocity profiles for 15 sites liquefied by the 1999 Kocaeli, Turkey earthquake. These profiles are used in order to evaluate each liquefaction site by the simplified shear wave velocity procedure. This procedure allowed for the identification of a potentially liquefiable region within the subsurface at each site. Locating this region at each site allowed for the separation of soils that were too stiff to liquefy from soils that were soft enough to liquefy. Once these soft regions had been identified, they were evaluated to separate granular soils expected to liquefy, from fine-grained soils expected not to liquefy. At sites where actual soil samples were available, this was accomplished by using the Chinese Criteria and the Andrews and Martin Criteria. At sites where only CPT data were available, this was accomplished by developing profiles of soil behavior type index (lc). Granular soil layers were located within the liquefiable region at 11 of the liquefaction test sites. It is assumed that these layers were the ones responsible for the observed liquefaction. The depth and thickness of each of these layers have been identified. However, at four of the liquefaction sites, only soils predicted as not susceptible to liquefaction were encountered. In these cases, the layer coming closest to fulfilling the Chinese Criteria and the Andrews and Martin Criteria was chosen as the one most likely to have liquefied. At each of these four sites, this layer appeared to be primarily made up of non-plastic silts having 2 μm clay contents ranging from 15 - 25

Using explainability to design physics-aware CNNs for solving subsurface inverse problems

We present a novel method of using explainability techniques to design physics-aware neural networks. We demonstrate our approach by developing a convolutional neural network (CNN) for solving an inverse problem for shallow subsurface imaging. Although CNNs have gained popularity in recent years across many fields, the development of CNNs remains an art, as there are no clear guidelines regarding the selection of hyperparameters that will yield the best network. While optimization algorithms may be used to select hyperparameters automatically, these methods focus on developing networks with high predictive accuracy while disregarding model explainability (descriptive accuracy). However, the field of Explainable Artificial Intelligence (XAI) addresses the absence of model explainability by providing tools that allow developers to evaluate the internal logic of neural networks. In this study, we use the explainability methods Score-CAM and Deep SHAP to select hyperparameters, such as kernel sizes and network depth, to develop a physics-aware CNN for shallow subsurface imaging. We begin with a relatively deep Encoder-Decoder network, which uses surface wave dispersion images as inputs and generates 2D shear wave velocity subsurface images as outputs. Through model explanations, we ultimately find that a shallow CNN using two convolutional layers with an atypical kernel size of 3x1 yields comparable predictive accuracy but with increased descriptive accuracy. We also show that explainability methods can be used to evaluate the network's complexity and decision-making. We believe this method can be used to develop neural networks with high predictive accuracy while also providing inherent explainability.Comment: 26 pages, 14 figures, 4 table

Using Convolutional Neural Networks to Develop Starting Models for 2D Full Waveform Inversion

Non-invasive subsurface imaging using full waveform inversion (FWI) has the potential to fundamentally change engineering site characterization by enabling the recovery of high resolution 2D/3D maps of subsurface stiffness. Yet, the accuracy of FWI remains quite sensitive to the choice of the initial starting model due to the complexity and non-uniqueness of the inverse problem. In response, we present the novel application of convolutional neural networks (CNNs) to transform an experimental seismic wavefield acquired using a linear array of surface sensors directly into a robust starting model for 2D FWI. We begin by describing three key steps used for developing the CNN, which include: selection of a network architecture, development of a suitable training set, and performance of network training. The ability of the trained CNN to predict a suitable starting model for 2D FWI was compared against other commonly used starting models for a classic near-surface imaging problem; the identification of an undulating, two-layer, soil-bedrock interface. The CNN developed during this study was able to predict complex 2D subsurface images of the testing set directly from their seismic wavefields with an average mean absolute percent error of 6%. When compared to other common approaches, the CNN approach was able to produce starting models with smaller seismic image and waveform misfits, both before and after FWI. The ability of the CNN to generalize to subsurface models which were dissimilar to the ones upon which it was trained was assessed using a more complex, three-layered model. While the predictive ability of the CNN was slightly reduced, it was still able to achieve seismic image and waveform misfits comparable to the other commonly used starting models. This study demonstrates that CNNs have great potential as a tool for developing good starting models for FWI ...Comment: 29 pages, 16 figures, submitted to Geophysic

Shear-Wave Velocity Characterization of the USGS Hawaiian Strong-Motion Network on the Island of Hawaii and Development of an NEHRP Site-Class Map

To assess the level and nature of ground shaking in Hawaii for the purposes of earthquake hazard mitigation and seismic design, empirical ground-motion prediction models are desired. To develop such empirical relationships, knowledge of the subsurface site conditions beneath strong-motion stations is critical. Thus, as a first step to develop ground-motion prediction models for Hawaii, wspectral-analysis-of-surface-waves (SASW) profiling was performed at the 22 free-field U.S. Geological Survey (USGS) strong-motion sites on the Big Island to obtain shear-wave velocity (V(S)) data. Nineteen of these stations recorded the 2006 Kiholo Bay moment magnitude (M) 6.7 earthquake, and 17 stations recorded the triggered M 6.0 Mahukona earthquake. V(S) profiling was performed to reach depths of more than 100 ft. Most of the USGS stations are situated on sites underlain by basalt, based on surficial geologic maps. However, the sites have varying degrees of weathering and soil development. The remaining strong-motion stations are located on alluvium or volcanic ash. V(S30) (average V(S) in the top 30 m) values for the stations on basalt ranged from 906 to 1908 ft/s [National Earthquake Hazards Reduction Program (NEHRP) site classes C and D], because most sites were covered with soil of variable thickness. Based on these data, an NEHRP site-class map was developed for the Big Island. These new V(S) data will be a significant input into an update of the USGS statewide hazard maps and to the operation of ShakeMap on the island of Hawaii.George E. Brown, Jr. Network for Earthquake Engineering Simulation (NEES) under NSF CMS-0086605FEMA HSFEHQ-06-D-0162, HSFEHQ-04-D-0733U.S. Geological Survey, Department of the Interior 08HQGR0036Geotechnical Engineering Cente

hapConstructor: automatic construction and testing of haplotypes in a Monte Carlo framework

Summary: Haplotypes carry important information that can direct investigators towards underlying susceptibility variants, and hence multiple tagging single nucleotide polymorphisms (tSNPs) are usually studied in candidate gene association studies. However, it is often unknown which SNPs should be included in haplotype analyses, or which tests should be performed for maximum power. We have developed a program, hapConstructor, which automatically builds multi-locus SNP sets to test for association in a case-control framework. The multi-SNP sets considered need not be contiguous; they are built based on significance. An important feature is that the missing data imputation is carried out based on the full data, for maximal information and consistency. HapConstructor is implemented in a Monte Carlo framework and naturally extends to allow for significance testing and false discovery rates that account for the construction process and to related individuals. HapConstructor is a useful tool for exploring multi-locus associations in candidate genes and regions

Deep Downhole Seismic Testing at the Waste Treatment Plant Site, Hanford, WA,Volume IV.S-Wave Measurements in Borehole C4993 Seismic Records, Wave-Arrival Identifications and Interpreted S-Wave Velocity Profile.

In this volume (IV), all S-wave measurements are presented that were performed in Borehole C4993 at the Waste Treatment Plant (WTP) with T-Rex as the seismic source and the Lawrence Berkeley National Laboratory (LBNL) 3-D wireline geophone as the at-depth borehole receiver. S-wave measurements were performed over the depth range of 370 to 1300 ft, typically in 10-ft intervals. However, in some interbeds, 5-ft depth intervals were used, while below about 1200 ft, depth intervals of 20 ft were used. Shear (S) waves were generated by moving the base plate of T-Rex for a given number of cycles at a fixed frequency as discussed in Section 2. This process was repeated so that signal averaging in the time domain was performed using 3 to about 15 averages, with 5 averages typically used. In addition, a second average shear wave record was recorded by reversing the polarity of the motion of the T-Rex base plate. In this sense, all the signals recorded in the field were averaged signals. In all cases, the base plate was moving perpendicular to a radial line between the base plate and the borehole which is in and out of the plane of the figure shown in Figure 1.1. The definition of “in-line”, “cross-line”, “forward”, and “reversed” directions in items 2 and 3 of Section 2 was based on the moving direction of the base plate. In addition to the LBNL 3-D geophone, called the lower receiver herein, a 3-D geophone from Redpath Geophysics was fixed at a depth of 22 ft in Borehole C4993, and a 3-D geophone from the University of Texas (UT) was embedded near the borehole at about 1.5 ft below the ground surface. The Redpath geophone and the UT geophone were properly aligned so that one of the horizontal components in each geophone was aligned with the direction of horizontal shaking of the T-Rex base plate. This volume is organized into 12 sections as follows. Section 1: Introduction, Section 2: Explanation of Terminology, Section 3: Vs Profile at Borehole C4993, Sections 4 to 6: Unfiltered S-wave records of lower horizontal receiver, reaction mass, and reference receiver, respectively, Sections 7 to 9: Filtered S-wave signals of lower horizontal receiver, reaction mass and reference receiver, respectively, Section 10: Expanded and filtered S-wave signals of lower horizontal receiver, and Sections 11 and 12: Waterfall plots of unfiltered and filtered lower horizontal receiver signals, respectively