Central US Seismic Velocity Model (CUSVM) Version 1
Compared to what we know about earthquake hazards in California, less is known in the Central US. The main reason for this difference is the relatively few numbers of large earthquakes in the Central US compared to California. The more frequent occurrence of small and larger earthquakes in California gives scientists data that can be used to infer the effects of historical earthquakes in the region and estimate the effects of future large earthquakes. A lack of large earthquakes in the Central US means a lack of data on stronger ground motions from big earthquakes. When there is little observational data, scientists use models to generate estimated data. A model for estimating the effects of an earthquake includes an earthquake source and the earth through which the seismic waves travel.
Scientists can use recorded small earthquakes or mathematically simulated earthquakes as a model for a large earthquake source, but what about the earth part? The subsurface geophysical structure, that is in this case, a combination of how fast sound waves travel through different parts of the crust and the density of these materials, can be a proxy for crustal structure, which can provide the needed crustal model. The measurements of sound wave speeds and density in the earth come from earthquake recordings, man-made seismic sources, and drill-hole data.
The goals of this research project were to 1) build a model of the crustal structure in a 650,000 km2 area including several urban areas, such as Little Rock, AR; Evansville, IN; Memphis, TN; Nashville, TN; Louisville, KY; and St. Louis, MO, that 2) can be used to generate simulations of shaking from earthquake sources. To test the model, the source of the M5.4 Mt Carmel, Illinois, earthquake on April 18, 2008 was used with the model to simulate data, and the simulated data were then compared to the real data recorded at many seismograph stations in the central U.S.for that earthquake. The two data sets had a good match, indicating the crustal structure model is a good approximation of actual geologic structure.
The next step will be to use this crustal structure model to simulate the shaking effects from the 1811 and 1812 earthquakes, both to estimate what happened then and be able to constrain forecasts about what could happen in a similar large earthquake in the Central US in the future.
The CUS model that we present covers a region of approximately 650,000 km2 and is home to more than 17 million people. The model includes several urban areas, such as Little Rock, AR; Evansville, IN; Memphis, TN; Nashville, TN; Louisville, KY; and St. Louis, MO. The CUSVM is parameterized by extending the physical properties of most of the well defined units, to less well-constrained areas, in the Mississippi Embayment and the CUS: stiff to unconsolidated Mesozoic and Cenozoic sediments, Paleozoic, crystalline upper crust, lower Crust, modified Lower Crust and Upper Mantle (e.g. Mooney et al., 1983; Ginzburg et al., 1983, Stewart, 1968; Prodehl et al., 1984).
The post-Paleozoic structure in the Mississippi Embayment is subdivided into units whose tops are defined by prominent reflectors in available refraction/reflection lines and geologic unit boundaries detected using geophysical and geologic logs. In addition, details of the model are locally enhanced by embedding, with small variations, the St. Louis, MO and Memphis, TN metropolitan models used in the respective hazard mapping projects (Chung and Rogers, 2010; Gomberg et al., 2003; Cramer et al. 2004). Outside the embayment and the urban areas mentioned above, where shallow Vs measurements are sparse, we relied on Vs30 information derived from the correlation between Vs30 and topographic slope (Allen and Wald, 2007, 2009).
The units are bounded with surfaces computed using a Kriging algorithm (see e.g. Isaaks and Srivastava, 1989; see Data and Resources), and in some cases corrected manually to remove spurious artifacts. We note that using units with well-defined boundaries is only an approximation of the complexity of the crust. It has been reported however, that the very sharp vertical contrasts, which potentially define our boundaries, in the geophysical properties in certain regions of the CUS seen in the abundant crustal imaging data described above, may not be present in less-well studied areas such as central Tennessee, or at least are not evident (Prodehl et al., 1984).
Depending on the availability of information, the properties within each unit have vertical variability in all cases and in some regions horizontal variability as well. In the absence of information, except where it is mentioned in the text, we employ Brocher (2005b) Vp to Vs and mass density ()-Vp relationships. These empirical relationships, valid within the rage 1.5<Vp<8.5 km/s, were derived based on a compilation of Vp, Vs and densities for a wide variety of common lithologies and have been used in other velocity models (e.g. Brocher 2005a; Olsen et al., 2008). The model is available and distributed as a C program where Vp, Vs and density can be queried as a function of depth from the surface, henceforth depth, or elevation with mean sea level as datum (see Data and Resources).
Data and Resources
- The Kriging Matlab toolbox software used in this research was developed by Dezhang Chu with funding from the National Science Foundation through the U.S. GLOBEC Georges Bank Project's Program Service and Data Management Office.
- Preliminary Integrated Geologic Map Databases for the United States.
- Digital Surfaces and thickness of the hydro-geologic of the Mississippi Embayment units.
- USGS North America compilation of seismic properties.
- The elevation data used is derived using the 1-arc second USGS Digital Elevation Model of the region.
- The Vs30 outside the Mississippi Embayment and St. Louis, MO were obtained from the USGS Global Vs30 map server.
- The Illinois depth to bedrock was obtained from the Illinois State Geological Survey.
- The Indiana depth to bedrock was compiled by the Illinois Geological Survey.
- The Tennessee depth to bedrock was compiled by the USGS.
- The ground motion data was obtained from the St. Louis University Earthquake Center and Incorporated Research Institute for Seismology (IRIS) sites.
Most sources used to construct the velocity model are available upon request.