Why the disagreement about seismic risk at New Madrid?
The New Madrid Seismic Zone presents a difficult problem. We know that large earthquakes have happened in the past. If earthquakes of that magnitude happened today, the damage and recovery would be difficult. Here is the problem: how big were those historical earthquakes actually? How likely are they to happen again? How should the cost of retrofitting be weighed against the predicted cost of a large earthquake? Scientists and policymakers have different training. Scientists are trained to assess the recurrence interval and estimate the ground motion of hypothetical events, while policymakers are trained to assess normative problems (i.e. given a seismic risk at some level, what should we do about it?)
In the data analyses you just completed, you became familiar with earthquake catalogs, including their strengths and limitations. You practiced looking at frequency-magnitude diagrams and you used this data to estimate the recurrence interval for earthquakes of various sizes. In fact, seismological data is just one of the tools scientists use to estimate earthquake recurrence interval. In the reading activity on the next page, you will break up into groups to investigate other methods of studying the NMSZ.
Over the past ten or fifteen years, global positioning system satellite data has become an invaluable tool for measuring plate motion and strain accumulation across faults. This data is gathered by installing geodetic markers in the ground. Scientists then use GPS receivers at the sites of the markers to find out their exact locations from satellites. Over time, the position of some markers may shift relative to each other; for example, markers on opposite sides of a fault may move closer together or further apart or be offset laterally as the years go by. This motion can be used to infer the strain rate in the crust. In the case of the New Madrid Seismic Zone, the faults are buried, so GPS data can help to find out exactly where the faults are and to determine the direction and extent of motion along them.
After several years of repeated measurements, the motion of the markers over the measurement time period is assessed. At active plate boundaries, such as along the San Andreas Fault on the West Coast of the United States, geodetic surveys have been used in concert with detailed records of seismicity to estimate stress buildup on faults and to predict seismic hazard. For example, a suite of geodetic markers may be placed around a fault of interest. After many measurements, the motion of the markers relative to each other can confirm the sense of motion on the fault, how fast the plates on either side of the fault are moving, and whether the fault itself is creeping or locked.
There have been several GPS campaigns over the last decade whose purpose has been to discover how much strain is building up at the New Madrid Seismic Zone. This work has been tricky because the faults involved are not well mapped, so the decision about where to place the markers hasn't been straightforward. The debate is still ongoing concerning whether the strain rates are high, thus posing a great seismic risk, or whether the strain rates are low, thus posing a lesser seismic risk to the area.
The map below shows current GPS stations operating in the USA.
Some faults can be excavated and mapped geologically in order to find out about the recurrence interval for large earthquakes. This sort of work is often done by digging a big trench with a backhoe and then trying to date any large offsets that are found. This technique is useful because the largest possible earthquakes of even quite active faults usually happen several hundred years apart. (Recall the ballpark range of recurrence intervals you estimated in your data analysis exercise.) We simply don’t have seismicity records that go back that far in this country. Dates for prehistoric earthquakes can be estimated by using the dates of the sediments that have been interrupted by an earthquake or some bit of organic material, such as charcoal, in an adjacent layer that can be dated. In the New Madrid Seismic Zone, stream offsets and evidence of liquefaction (sand blows and dikes) caused by strong shaking are also clues to past earthquakes. Paleoseismologists use all these clues to try to put together a timeline of recurrence interval and the approximate earthquake magnitude for a particular fault. These data can be linked with seismicity catalogs and geodetic surveys to get a fuller picture of seismic hazard.
To see excavation and mapping in action, check out this short video from Teachers' Domain and NOVA Online about the work of Kerry Sieh, a paleoseismologist at The Earth Observatory of Singapore (He was a prof at Caltech when they made this video).
Video: NOVA: The Work of Kerry Sieh (3:19)
The Earth produces heat from the decay of radioactive elements in its interior. This heat drives mantle convection and therefore the movement of tectonic plates. Heat flow is routinely measured in boreholes around the planet. These measurements are compiled to produce a map of heat flow for the Earth's surface. Some degree of estimation and smoothing must be applied to the measurements because the boreholes are not evenly spaced and some are on continents while other measurements are taken in oceanic crust. The map below shows global heat flow.
This map shows color-coded contours of the global distribution of heat flow at the surface of the Earth's crust. Major plate boundaries and continent outlines are also shown. The fundamental data embodied in this map are the more than 24,000 field measurements in both continental and oceanic terrains, supplemented by estimates of the heat flow in the unsurveyed regions. The estimates are based on empirically determined characteristic values for the heat flux in various geological and tectonic settings. Observations of the oceanic heat flux have been corrected for heat loss by hydrothermal circulation through the oceanic crust. The global data set so assembled was then subjected to a spherical harmonic analysis. The map is a representation of the heat flow to spherical harmonic degree and order 12.
What does this have to do with the New Madrid Seismic Zone? By looking at the map above, you can see that the amount of heat flowing out of the Earth is not uniform over the surface of the planet. Some areas have much higher heat flow than others and these areas are usually associated with tectonic activity such as volcanism and plate boundaries. For example, the boundaries of the North American plate, the Mid-Atlantic Ridge and the San Andreas fault system, both show up as "warm" places on this map. Heat flow measurements have been made in the New Madrid Seismic Zone to see whether this is a high heat flow area compared to what would be expected for the interior of a continent. (Conventional geophysical wisdom holds that the interior of continents should be old, cold, and stable.) If heat flow is higher than expected, this would be evidence for why earthquakes happen in this area. This remains a point of scientific contention. Past surveys concluded that heat flow was high in the NMSZ, but the most recent studies disagree with those earlier findings.
The map below comes from the Global Heat Flow Database, at the University of North Dakota. They have heat flow maps and data files for all different parts of the world. This particular map shows borehole data for the United States. Warm colors denote higher heat flow than cool colors (see the legend, which shows milliwatts per meter squared values color-coded). Notice that this map looks different than the global map above. It looks different because it shows exact borehole measurements as opposed to smoothed values that have been interpolated over the whole map. Where is the heat flow highest? Where is it lowest? Compare this map with the map above to see whether they are consistent for the US.