Model Overview
Model overview. The regional structure is constructed by extending the physical properties of most of the well-defined units to less well-constrained areas in the Mississippi Embayment and the CUS: Paleozoic, Crystalline Upper Crust, Lower Crust, Modified Lower Crust and Upper Mantle. The Mesozoic to Cenozoic materials are represented by several units within the Mississippi Embayment and a Vs30 based layer outside the embayment.

CUSVM Version 1.3 Model and Documentation

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. A research paper on this modeling effort (4.7MB PDF) was published in the August 2015 issue of the Bulletin of the Seismological Society of America.

Topographic map showing the extent of the CUSVM
Topographic map showing the extent of the CUSVM. Red circles indicate earthquakes greater than magnitude 2.5 that occurred after 1972 (largest earthquake in this time period was the Mw 5.4 Mt. Carmel earthquake in southeastern Illinois). The Mississippi Embayment and the Reelfoot Rift (RFR) boundaries are depicted with red dashed and black continuous lines, respectively. The plausible geometry of three major faults (Northern, Reelfoot and Axial or Cottonwood Grove) is displayed in white discontinuous line. Some of the geologic structures of the region are also indicated, OU=Ozark Uplift, ND=Nashville Dome, CA=Cincinnati Arch, IB=Illinois Basin and RCG=Rough Creek Graben.

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.2 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.

Computer Simulation of a Magnitude 7.7 Earthquake in the New Madrid Seismic Zone

This animation shows the simulated ground motion of the Earth’s surface that could occur in the central U.S. region for a magnitude 7.7 strike-slip earthquake on the southern section of the New Madrid seismic zone. The simulation is based on a detailed geologic model of the central U.S. and is intended to show the general character of the long-period ground shaking expected for a large earthquake in the New Madrid seismic zone. It demonstrates the profound focusing effect the Reelfoot rift has on ground shaking produced by earthquakes. Further, the relatively unconsolidated material in the Mississippi embayment amplifies the ground shaking and also traps seismic energy to prolong the duration of ground shaking.

Simulation of M7.7
View simulation of magnitude 7.7 earthquake in the New Madrid Seismic Zone.

The solid straight line in the middle of the New Madrid seismic zone is the surface projection of the modeled fault, which ruptures in the simulation. The colors are keyed to the peak intensity of ground velocity at the surface. As the seismic waves propagate away from the hypocenter in all directions, we see successive waves of strong ground shaking that begins moving along the Reelfoot rift and tends to be focused northeast toward Paducah, Kentucky, and southwest toward Little Rock, Arkansas. These higher amplitude waves are the more damaging surface waves. These seismic waves, trapped by the softer sedimentary rock in the Rift and Embayment and amplified by the softer sediments above, cause the duration of long-period shaking to be about 30–45 seconds in some areas, including Memphis, Little Rock, and Paducah, Kentucky.

Technical Summary

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 geophysical 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 below).

Data and Resources

Most sources used to construct the velocity model are available upon request.

Scientific Staff