In historical, large-magnitude earthquakes, lateral spreading has been a very damaging type of ground failure. When a subsurface soil deposit liquefies, intact blocks of surficial soil can move downslope, or toward a vertical free face, even when the ground surface is nearly level. A lateral spread is defined as the mostly horizontal movement of gently sloping ground (less than 5% surface slope) due to elevated pore pressures or liquefaction in undelying, saturated soils. Here, lateral spreading is defined specifically to exclude liquefaction failures of steeper embankments and retaining walls, which can also produce lateral surface deformations. Lateral spreads commonly occur at waterfront sites underlain by saturated, recent sediments and are particularly threatening to buried utilities and transportation networks. While the occurrence of soil liquefaction and lateral spreading can be predicted at a given site, methods are needed to estimate the magnitude of the resulting deformations. In this research effort, an empirical model was developed for predicting horizontal and vertical surface displacements due to liquefaction-induced lateral spreading. The resulting model is called "EPOLLS" for Empirical Prediction Of Liquefaction-induced Lateral Spreading. Multiple linear regression analyses were used to develop model equations from a compiled database of historical lateral spreads. The complete EPOLLS model is comprised of four components: (1) Regional-EPOLLS for predicting horizontal displacements based on the seismic source and local severity of shaking, (2) Site-EPOLLS for improved predictions with the addition of data on the site topography, (3) Geotechnical-EPOLLS using additional data from soil borings at the site, and (4) Vertical-EPOLLS for predicting vertical displacements. The EPOLLS model is useful in phased liquefaction risk studies: starting with regional risk assessments and minimal site information, more precise predictions of displacements can be made with the addition of detailed site-specific data. In each component of the EPOLLS model, equations are given for predicting the average and standard deviation of displacements. Maximum displacements can be estimated using probabilities and the gamma distribution for horizontal displacements or the normal distribution for vertical displacements. / Ph. D.

There are numerous examples in the literature of time-dependent changes in the proper-ties of sands, or aging effects. Most of these aging effects are of increases in the cone penetration resistance. Time-dependent increases in penetration resistance have been measured in hydraulically placed fills and freshly densified deposits, with the largest in-creases following the use of ground modification techniques such as vibrocompaction, dynamic compaction, and blast densification. It is not known what causes these increases in penetration resistance to occur. The objective of this research was to gain an understanding of the possible mechanisms responsible for aging effects in sands. Current hypotheses to explain what causes aging effects in sands include increased interlocking of particles, internal stress arching, and precipitation of silica or carbonate minerals at the contacts between grains. To date, no unambiguous evidence has been presented to support these hypotheses. A laboratory testing program was developed to study the influence of different variables on the pres-ence and magnitude of aging effects. Three different sands were tested in rigid wall cells and buckets. Samples were aged under different effective stresses, densities, tempera-tures, and pore fluids. In every rigid wall cell, three independent measurements were made to monitor property changes during the aging process: small strain shear modulus using bender elements, electrical conductivity, and mini-cone penetration resistance. At the end of each test, detailed mineralogical tests were performed to assess changes in the chemistry of the samples and pore fluids. A total of 22 tests in rigid wall cells were per-formed with periods of aging ranging from 30 to 118 days. Mini-cone penetration resis-tances were measured in the buckets before and at various times during the aging process. Increases in the small strain shear modulus were measured with time. It was found that sand type and pore fluid composition greatly influenced the amount of increase in small strain shear modulus. Density was also found to influence the amount of increase in small strain shear modulus. Temperature was found to have little influence on the in-crease in small strain shear modulus with time. Changes in the chemistry of the samples were also measured with time. The dissolution and precipitation of minerals in solution was monitored with electrical conductivity measurements. In most of the tests, there was continual dissolution of minerals with time. Mineralogical studies and conductivity measurements indicated precipitation of carbonates and silica in two of the tests; however, scanning electron micrographs showed no visible evidence of precipitation. Despite the measured increases in small strain shear modulus and evidence of mineral precipitation, there were no increases in the mini-cone penetration resistance with time. This finding is significant and suggests that small-scale laboratory experiments do not capture the mechanism(s) that are responsible for time-dependent increases in penetration resistance in the field. / Ph. D.

The feasibility of using ground improvement at existing highway bridges to mitigate the risk of earthquake-induced liquefaction damage has been studied. The factors and phenomena governing the performance of the improved ground were identified and clarified. Potential analytical methods for predicting the treated ground performance were investigated and tested. Key factors affecting improved ground performance are the type, size, and location of the treated ground. The improved ground behavior is influenced by excess pore water pressure migration, ground motion amplification, inertial force phasing, dynamic component of liquefied soil pressure, presence of a supported structure, and lateral spreading forces. Simplified, uncoupled analytical methods were unable to predict the final performance of an improved ground zone and supported structure, but provided useful insights. Pseudostatic stability and deformation analyses can not successfully predict the final performance because of their inability to adequately account for the transient response. Equivalent-linear dynamic response analyses indicate that significant shear strains, pore water pressures and accelerations will develop in the improved ground when the treated-untreated soil system approaches resonance during shaking. Transient seepage analyses indicate that evaluating pore pressure migration into a three-dimensional improved zone using two-dimensional analyses can underestimate the pore pressures in the zone. More comprehensive, partially-coupled analyses performed using the finite difference computer program FLAC provided better predictions of treated ground performance. These two-dimensional, dynamic analyses based on effective stresses incorporated pore pressure generation, non-linear stress-strain behavior, strength reduction, and groundwater flow. Permanent movements of structures and improved soil zones were predicted within a factor of approximately two. Predictions of ground accelerations and pore water pressures were less accurate. Dynamic analyses were performed with FLAC for an example bridge pier and stub abutment on an approach embankment supported on shallow foundations and underlain by thick, liquefiable soils with and without improved ground zones. Ground improvement that restricted movements of the pier and stub abutment to tolerable levels included improved zones of limited size extending completely through the underlying liquefiable soils and formed through densification by compaction grouting or cementation by chemical grouting or jet grouting. A buttress fill at the abutment was unsuccessful. / Ph. D.

Many problems are involved with inspecting and evaluating the condition of bridges in the United States. Concrete bridge deck inspection and evaluation presents one of the largest problems. The deterioration of these concrete decks progresses more rapidly than any other bridge component, which leads to early concrete deck replacements that must be done before the bridge superstructure needs to be replaced. The primary cause of deterioration in these concrete bridge decks is corrosion-induced concrete cracking, which frequently results in delaminations. Delamination distress increases the life cycle cost of maintaining a concrete bridge deck, particularly when it is not detected early on. Early detection of delamination distress can facilitate economical repair and rehabilitation work, but bridge engineers must recommend deck replacement if repairs are delayed too long or inspection tools cannot detect delaminations early enough. The Federal Highway Administration has responded to the need for a better bridge deck inspection tool by contracting Lawrence Livermore National Laboratory to develop two new prototype ground penetrating radar systems. These two systems generate three-dimensional data that provide a representation of features that lie below the bridge deck surface. Both of these systems produce large amounts of data for an individual bridge deck, which makes automated data processing very desirable. The primary goal of the automated processing is to characterize bridge deck distress represented in the data. This study presents data collected from sample bridge deck sections using one of the prototype systems. It also describes the development and implementation of appropriate methods for automating data processing. The automated data processing is accomplished using image processing and pattern recognition algorithms developed in the study. / Ph. D.

Ground penetrating radar (GPR) technology has been used for the past 20 years for a variety of applications to assess transportation infrastructure. However, the main issue after all these years remains: "How well does GPR work and under what conditions?" Results show that GPR works well for some situations, but is not an appropriate tool for other situations. It is not used currently on a routine basis by the US Departments of Transportation (DOTs) due mainly to difficulties encountered with data interpretation. Data interpretation difficulties are mainly attributed to the fact that images obtained from the reflected signals are not photographs of the features that are beneath the surface being investigated. The images show the amplitude of the radar-reflected signals from the interfaces with different dielectric properties. Therefore, a considerable amount of experience and operator skill may be required to correctly interpret sub-surface radar results. To better understand reflected GPR signals, this research was conducted with the following objectives: to determine the dielectric properties of concrete over the used GPR frequency range; to synthesize the reflected air-coupled radar signals and compare them with measured waveforms; to model and study the effects of simulated defects in concrete on the reflected air-coupled and ground-coupled radar signals; and to validate the research results in the field by predicting layer thicknesses of flexible pavements and detecting moisture in flexible pavement systems. Several concrete slabs, 1.5x1.5 m, were constructed with known thicknesses, simulated defects, and different reinforcement configurations. The concrete mixes included four different bridge deck mixes and one concrete pavement mix used in the State of Virginia. Results have shown that the dielectric constant of concrete is frequency and mix dependent. However, modeling the reflected signals using an average complex dielectric constant over the entire radar frequency range led to modeled waveforms comparable to the measured waveforms. Although air- and water-filled voids did distort the reflected waveforms, a model was developed to predict the reflected waveforms from the simulated defects. Reinforcement was found to affect the reflected waveforms only when it was oriented in a direction perpendicular to the GPR antennas. A model was also developed to predict the GPR waveforms obtained from flexible pavements. This model could be used in a procedure to measure layer thicknesses more accurately by including losses that occur inside the pavement materials. Two different case studies, where a ground-coupled GPR system was used to locate moisture at different layers, have led to the conclusion that the ground-coupled GPR is a feasible tool to detect moisture inside pavements. / Ph. D.

Since its development in the 1920s, bolting has become the most dominant support method in underground construction. However, because of the geological environment, the design process for roof bolt systems is an art rather than a science. To quantify the selection of bolting systems a MSBT (minimum solid beam thickness) approach was developed. The ultimate goal of this bolt design paradigm was achieved by optimizing bolt length, bolt density, and bolt pretension during installation. The impact of the number of strata layers within bolting range and pretension applied to bolts upon the stability of an opening was investigated using FLAC model. Four statistical models for predicting optimum bolt supports using a minimum solid beam thickness were established, and based on these results, a design criterion was proposed. To meet support needs in various geological and geotechnical settings, a variety of bolt types have been developed. The installation of such bolt-based support systems is often complex and specialized, and thus imposes a challenge for engineers to identify the specific cause and to take appropriate remedial measures once problems arise. To solve these problems, a knowledge-based expert system (KBES) has been developed. The knowledge base includes the data accumulated from years of laboratory and field investigations conducted by the Mine Safety and Health Administration of the US Department of Labor. A user-friendly Windows-based program was implemented using KAPPA environment. After identifying the problem, the KBES searches its knowledge base and reasons out the most likely, secondary, and other potential causes, then provides solutions according to users' input. The results of this research are validated and demonstrated using case studies. / Ph. D.

Estimating local site effects and modifying the uncertainty in ground motion predictions are two indispensable parts of partially non-ergodic site-specific PSHA. Local site effects can be estimated using site response simulations or recorded ground motions at the site. When such predictions are available, the aleatory variability of ground motions used in PSHA can be changed to the single station sigma value. However, in these cases, the epistemic uncertainty in predicting site effects must be incorporated into the hazard analyses. This research focuses on the challenges specific to conducting partially non-ergodic site-specific PSHA using recorded ground motions or site response analysis. The main challenge in estimating local site effects using recorded data is whether ground motions collected in a relatively short time can be used to estimate site effects for long return period events. We first develop a database for recorded ground motions at the KiK-net array to investigate this question and use this database to develop a predictive model for the Fourier Amplitude Spectra of ground motions. The ground motion model (GMM) residuals are used to investigate the stability of site terms across different tectonic regimes. We observe that empirical site terms are stable across different tectonic regimes. This observation allows the use of ground motions from any tectonic regime (whether they belong to the tectonic regime that controls the hazard or not) to estimate local site effects. Moreover, in Fourier amplitude, site effects are not dependent on event magnitude and source to site distance; therefore, estimates of site effects from low magnitude events can be easily extrapolated to larger events. The Fourier amplitude GMM developed in this study adds to the library of Fourier amplitude models to be used in future partially non-ergodic site-specific PSHAs. In practice, one of the most common tools for simulating wave propagation is 1-D site response analysis. Two central assumptions in 1-D site response analysis are that the soil profile is comprised of horizontal soil layers of infinite extent and that the vertically propagating SH-waves control the horizontal component of ground motion. SH-waves tend to propagate vertically near the surface because as earthquake waves hit softer layers traveling from the source to the site, they refract until the path becomes steeply inclined. The validity of both assumptions in 1-D site response depends on the geological setting at the site and the geology between the earthquake source and the site, raising the question of which sites are suitable for 1-D site response analysis and what the model error in 1-D site response analysis is. We use the GMM developed for FAS to estimate observed and empirical site terms. The empirical site effects are then compared with the theoretical site effects to determine whether sites are amenable to 1-D site response analyses, and to quantify the model error in the analyses. / Doctor of Philosophy / It is impossible to predict future earthquake-induced ground motions due to randomness in the process and a lack of knowledge. In fact, there are significant uncertainties not only in predicting the location, time, and magnitude of a future earthquake but also in predicting the intensity of ground motion induced by a given future earthquake. Therefore, assessing the safety of the human environment against earthquake hazards requires a method that considers all sources of uncertainties. To this end, Earthquake Engineers have developed Probabilistic Seismic Hazard Analysis(PSHA) framework. Structural engineers use the results of PSHA to design a new structure or assess the safety of an existing building. The accuracy of PSHA estimations leads to designs that are both safe and cost-efficient. The distribution of possible ground motions induced by a given earthquake scenario significantly controls the result of PSHA. This distribution should consider the effect of source, source to site path, and local site effects. This research focuses on improving PSHA results by estimating local site effects using recorded ground motions or simulating wave propagation in the site. In estimating local site effects using recorded data, the local site effect observed in ground motions collected in a relatively short time window is used to estimate hazards from all scenarios. However, the collected ground motions usually belong to frequent low magnitude events that are different from large magnitude events that control the hazard. This difference requires either using a measure of local site effect that is independent of the magnitude and distance of the earthquake or considering the effect of magnitude and distance on the local site effect estimate. Moreover, since frequent events sample different sources and paths than large events, we need to make sure the local site effect is consistent across different sources and paths. This research develops Ground Motion Models(GMMs) for Fourier amplitude, a linear function of ground motion times series, using Japanese ground motions. The ratio of Fourier amplitude at the surface over bedrock is a measure of local site effect that is not dependant on magnitude and distance. The model is then used to see if the trade-off between source and site effect and path and site effect is significant or not. In practice, one of the most common tools for simulating wave propagation is 1-D site response analysis. Two central assumptions in 1-D site response analysis are that the soil profile comprises horizontal soil layers of infinite extent and that the vertically propagating horizontal shear waves (SH-waves) control the horizontal component of ground motion. SH-waves tend to propagate vertically near the surface because as earthquake waves hit softer layers traveling from the source to the site, they refract until the path becomes vertically inclined. The validity of both assumptions in 1-D site response depends on the geological setting at the site and the geology between the earthquake source and the site, raising the question of which sites are suitable for 1-D site response analysis and what the model error in 1-D site response analysis is. We use the GMM developed for FAS to estimate empirical local site effects. The empirical site effects are then compared with the theoretical site effects to determine whether sites are amenable to 1-D site response analyses and quantify the model error in the analyses.

An analytical modal response tool was developed to investigate the characteristics of and to estimate static and dynamic launch vehicle responses to ground wind loads (GWL). The motivation of this study was to estimate the magnitude of response of the Ares I-X launch vehicle to ground winds and wind-induced oscillation (WIO) during roll-out and on the pad. This method can be extended to other launch vehicle designs or structures that possess a nearly cylindrical cross-section. Presented in this thesis is an overview of the theory used, a comparison of the theory with wind tunnel data, further investigation of the data to support the assumptions used within the analysis, and a prediction of the full-scale Ares I-X response. Additionally, an analytical investigation is presented that estimates the effect of atmospheric turbulence on WIO response. Most of the wind tunnel data presented in this report is taken from the GWL Checkout Model tested in the NASA Langley Transonic Dynamics Tunnel (TDT) in April 2007. The objective of the GWL Checkout Model was to reestablish and evaluate the capability of the facility to conduct GWL testing and to operate the associated equipment. This wind tunnel test was not necessarily intended to predict the full scale Ares vehicle response to GWL; however, it can be used to help validate the newly developed analytical method described in this thesis. A detailed GWL test incorporating updated vehicle designs and launch pad configurations of the Ares I-X flight test vehicle was also conducted in the TDT during the fall of 2008. This test provides more accurate predictions of the second bending mode response of the Ares I-X, and it models effects of the nearby tower and support structures. The proposed analytical method is also compared to select data from the Ares I-X GWL test; however, it is presented as normalized values to protect the sensitivity of the data. Results of the proposed analytical method show reasonable correlation to wind tunnel data. Also, this method was the first to determine that second bending mode WIO response was not only possible for the Ares I-X, but will also produce the most critical loads. Finally, an explanation is offered in this thesis regarding discrepancies between wind tunnel and full-scale WIO response data. / Master of Science

This work focuses on the operational and safety issues associated with karst voids in large opening underground mines. Issues include water inrush, structural instability, and engineering uncertainty in these environments. Coupled with the fracturing prevalent in folded sedimentary rocks, karsts are complex and challenging ground control risks. Traditional methods of predicting karst void locations such as probe-drilling are impeded by the inconsistent spatial distribution and variable sizes of the features. Ground penetrating radar (GPR) is a geophysical technique that transmits radio waves into a medium and subsequently detects reflected waves via a receiver. The travel time and energy of received signals are then processed and interpreted. The difference in material properties between limestone and open karst voids causes strong reflections. This work summarizes a series of 2D and 3D GPR surveys for karst void mapping within a mine pillar and within sill pillars between mine levels in a large opening underground limestone mine. In this case study mine, karst voids are hazardous ground control risks that interact with geologic discontinuities, creating free blocks within the rock mass. As tunnels are advanced via blasting, unknown karst voids may be exposed and pose risks to mining personnel. The karst voids also form a hydrogeological network of water reservoirs with spatial locations throughout the rock mass that are difficult to predict with traditional methods such as drilling. While GPR has been utilized throughout several industries for anomaly detection, mapping, and validating other geophysical data sets, this technique has not seen the same proliferation within the mining industry. Regarding published literature, there is a lack of works that detail the applicability of GPR in underground mining scenarios. The aim of this work is to expand on previous methodologies establishing GPR as a useful tool in underground mining applications, and to discuss the benefits and limitations GPR data in such scenarios. / Master of Science / This work focuses on the operational and safety issues associated with karst voids in large opening underground mines. Typical issues include water flooding into the tunnels and rocks falling out from the roof and walls, among other things. Sedimentary rock structures sometimes are geologically complex, and karst voids only add to that complexity. Engineers usually predict karst void locations with drilling or statistics, but this is often challenging as karst voids have various shapes and orientations. Ground penetrating radar (GPR) is a geophysical technique that sends electric signals into the rock; these signals can reflect off of karst voids and other anomalies. The travel time and energy of signals that come back to the antennas are then processed and interpreted. The difference in material properties between limestone and open karst voids causes strong reflections. This work shows a series of 2D and 3D GPR surveys for karst void mapping within a pillar in a stone mine and also below the floor of mine tunnels. In this mine, karst voids are very dangerous and the miners spend significant time and resources to ensure the tunnels walls are stable. As tunnels are blasted, hidden karst voids may be exposed and pose unpredicted risks to miners. The karst voids are also connected by cracks and discontinuities, providing a path for water to travel along. While GPR has been used in various ways among the construction, civil engineering, and tunneling industries, there is not enough literature pertaining to its benefits for mines. The goal of this work is to grow the available literature on GPR in mining and to talk about the best practices for GPR use as a means of improving health and safety for miners underground.

The key goal in this work is to enable a flexible ground station that is not constrained to a particular mission or set of hardware. In addition, with the concepts and software produced in this thesis, it will play a significant role in educating engineers and students by providing critical infrastructure and a sandbox for ground station operations. Key pieces of software were developed in this work to create a flexible and robust software-defined ground station. Several digital transmission modes were developed in order to allow communication between the ground station and common amateur radio CubeSats and SmallSats. In order to handle distributed tasks and process at a ground station with multiple servers and controllers, a specialized actor framework was written in Python for ease of use. Actors have the ability to send messages to one another over a network, and they maintain their own memory in order to avoid synchronization problems that come with sharing memory. In addition to the software developed in this work, a novel Peer-to-Peer (P2P) protocol for a network of ground stations is proposed in order to increase coverage and access to spacecraft without requiring centralized server infrastructure. This protocol provides the method to scale the developed software architecture beyond a single ground station. Since the Virginia Tech Ground Station (VTGS) will have many concurrent processes running across multiple servers, it was necessary to apply the actor model in order to simplify the design of the system. The purpose of this thesis is to describe the developed software for the VTGS as well as the P2P protocol for a larger network of ground stations. There are three primary repositories: planck-dsp, gr-vtgs, and pystation. The planck-dsp library and gr-vtgs Out-of-tree (OOT) make up the primary digital signal processing and communications toolboxes, where GNU Radio serves as the scheduler for signal processing blocks used in flow graphs. The pystation module is the extensible software actor framework that connects various systems both locally and remotely. It is also responsible for scheduling and handling ground station requests. While the software was primarily created for the VTGS, it is general enough to apply to other ground station implementations. / Master of Science