Causes and Effects of a Late 20th Century Fault Slip Event in Coastal Louisiana

Wetlands loss in the delta plain of coastal Louisiana is best understood as an episode that peaked in the late 20th century, and has been generally declining for the past three decades.  Losses were concentrated in the Barataria and Terrebonne hydrologic basins, as well as some intense hot spots in the Mississippi River Delta.  This land loss episode had a strong temporal correlation with an episode of higher subsidence rate estimates derived from tide gauge records (Penland et al, 1988) and elevation benchmark leveling surveys (Shinkle and Dokka, 2004).  Many of the major hot spots of wetlands loss associated with this episode are located along the surface traces of major faults, and it appears that a late 20th century fault slip event was a principal cause of increased rates of subsidence.  Wetlands loss was primarily a function of the submergence of marshes on the hanging walls of the faults.

The implications for coastal sustainability are striking.  The hot spots of wetlands loss are for the most part hot spots of higher subsidence rates associated with faults.  The objective of most coastal sustainability projects is to build and maintain the elevation of the marsh surface.  There is a much lower probability of maintaining elevation in areas with higher subsidence rates.  Conversely, stable areas outside of the subsidence hot spots have been effectively maintaining elevation for decades through natural sediment accretion.  Sustainability projects that seek to work with nature by supplementing natural accretion on stable marsh platforms will be the most effective.

Picture1.png

THE TERREBONNE LINKED TECTONIC SYSTEM

The fault slip event occurred mainly on the major faults in the “Terrebonne Linked Tectonic System”.  The faults and salt domes in this interconnected system form the subsurface structural framework of the delta plain.  There appears to be a strong genetic relationship among the individual elements of the system.  Salt and faults appear to have moved interactively throughout the geological history of the area. 

Picture2.png
Picture3.png

Most of the salt domes in the system have a rounded conical shape in the subsurface above depths of about 10,000 feet.  For much of the 20th century geologists assumed that most domes rose vertically from a sedimentary layer called the Louann Salt.  The advent of 3-D seismic made it clear that most domes have a distinctly tilted orientation at greater depths.  Salt appears to have been squeezed out laterally from the source layer under the weight of an advancing sedimentary load during much of the Cenozoic Era.  The inset profile across the Lake Washington Dome shows this configuration, and the relationship between the dome and the Magnolia and Barataria faults.  It is likely that there was a feedback loop between sediment accumulation, movement of the salt, and movement on the faults throughout their geologic history.  The combined effects of fault and salt movement have created accommodation capacity for deltas throughout much of the Cenozoic Era, and the weight of the accumulated sediment has induced further movement.

The faults are arranged as connected linear segments within the system.  Faults along the north side of the system tend to be down to the south.  Their surface traces, which are coincident with the 0-depth contour of the fault plane map, are shown in light blue.  Faults along the south side of the system tend to be down to the north, and their surface traces are shown in red.  This arrangement tends to group the faults into conjugate pairs of opposing faults that create large structural grabens between them.  Sediments of the Late Miocene Epoch within the major grabens are noticeably thicker than they are outside of the system.  This network of grabens that runs along the coastline is commonly called the Terrebonne Trough.

Picture4.png
Picture5.png

The regional profile across southern Louisiana from Nelson et al, 2000 shows the dramatic thickening of the Late Miocene sediments into the Terrebonne Trough.  It also indicates that faults of the Terrebonne Linked Tectonic System merge into a larger detachment surface in which sediments of the late Mesozoic and Cenozoic Eras have moved laterally to the south above deeper sediments.  A horizontal vector of motion toward the south is measured at the surface today, and it probably is an expression of this deeply-rooted lateral movement (Karegar et al, 2017).  Salt also appears to have moved laterally along this same detachment surface before following tilted upward paths to form diapirs that rise toward the surface.

Picture9.png

The grabens of the Terrebonne Linked Tectonic System provided accommodation capacity for deltas throughout the Miocene Epoch.  The black and purple contours of the faults planes and salt domes show the configuration of the linked tectonic system.  The outlines of the Late Miocene Regressive Phase deltas numbered for cycles XI – XIII (from Curtis, 1970) are overlain on the contours and a colored grid representing Late Miocene sediment thickness (from Wu & Galloway, 2006).  It appears that Late Miocene deltas “took up residence” in the grabens of the structural system where maximum sediment accumulation occurred.  More detailed studies of faults within this system (Bullock et al, 2018, Levesh et al, 2019, McLindon, 2017) have shown a strong relationship between sediment loading and inferred rates of fault slip.  Large-scale, long-term fault slip events are indicated by the expansion indices measured across these faults.

Picture6.png

The accumulation of thick deltaic sands within the grabens of the linked tectonic system and the continual movement of the faults in response to sedimentary loading created a very efficient system for removing water from the sedimentary column.  This allowed the sand and shales to dewater and compact.  Faults within the system have been shown to be good conduits for fluid movement, and the thick clean sand layers formed aquifers that were also very efficient at moving water.  On either side of the central grabens the movement of water in the subsurface was not as efficient.  Fluids became trapped in the pore spaces, which prevented the mineral grains of the sediment from compacting with increasing overburden.  The onset of geopressure in these sediments occured when the mineral grains of the sediment matrix could not increase grain to grain contact to support the increased overburden, and so the pore fluids took on a portion of the overburden stress resulting in increased fluid pressure. 

Within the major grabens pore fluid escape and compaction continued to greater depths than on either side of the graben.  This has resulted in a thick column of normally-pressured sediments with the grabens that is more compacted and denser than the sedimentary columns on either side of the graben.  This density differential between sediments downthrown to the graben faults and the sediments upthrown to those faults has been a continual impetus for fault movement since the Late Miocene.  The major delta systems of the Pliocene and Pleistocene Epochs continued to build out basinward and prograding the continental shelf.  The major deltas from these epochs primarily interacted with faults on the modern shelf, but Fisk, 1960 showed that Pliocene and Pleistocene sediments were also thicker within the area of the Terrebonne Linked Tectonic System.  This indicates that fault movement continued into the Pleistocene.  There is evidence of Holocene movement on the Bastian Bay, Barataria, and Ironton faults in the Pleistocene.

Picture7.png

The Holocene Epoch represents the end of last glacial cycle of the Pleistocene.  It is a reasonable expectation that episodic movement, which had been continual on the faults and domes of the Terrebonne Linked Tectonic System throughout the Miocene, Pliocene and Pleistocene Epochs would have continued into the Holocene.  One of the best indications of Holocene fault movement is the relationship among the configuration of the deltas, the thickness of the sediments and the linked tectonic system. The dashed colored outlines of the the most recent delta lobes (from Frazier, 1967) shows that each lobe appears to have had a distinct relationship with one of the conjugate pairs of faults that define the Terrebonne Trough.  The thickness of the “topstratum” of the Holocene (from Kulp, 2000) represents the entire sedimentary accumulation of the Holocene deltas.  Sediment thickness is greater within the grabens of the linked tectonic system than it is outside of them. These grabens provided accommodation capacity to the deltas, just as they had done with the Miocene deltas.  Each of the most recent delta lobes is underlain by the sediments of at least one previous Holocene delta lobes.  Frazier defined sixteen historical delta lobes of the late Holocene.  Most of them have submerged below the surface, which further evidences subsidence throughout the Holocene.

The dramatic increase in Holocene thickness in the far eastern graben associated with the Balize delta may be evidence of a sedimentation anomaly associated with the onset of agriculture in the Mississippi Basin the the 18th and 19th centuries.  This will be discussed in more detail in a later section on possible causes for the late 20th century fault slip event.

Picture8.png

THE BASTIAN BAY AND LAKE BOUDREAUX FAULTS

Research on the surface traces of the faults in the Terrbonne Linked Tectonic System has expanded significantly in the past decade.  Most of the traces shown here are documented by academic research.  Two of the major faults will be considered in more detail to examine characteristics that are shared all of the faults. 

Picture10.png
Picture11.png

The fault plane contours of the Bastian Bay fault have been mapped in the subsurface by the interpretation of the petrophysical logs from oil and gas wells and seismic data.  The trace of the fault where it crosses a Late Miocene sand layer at a depth of about 10,000 feet is parallel to the surface trace.  The fault creates a “rollover” anticlinal structure on the sand layer, which allows for the trapping of natural gas in the sand at the top of the anticline. Gas accumulations in several Late Miocene sands form Bastian Bay field.  Most of the major faults in the Terrebonne Linked Tectonic System have a similar relationship with oil and gas fields.  There is a causal relationship between the location of the surface traces of the faults and the locations of large oil and gas fields. 

Picture12.png

These block diagrams show the Bastian Bay fault in context with other elements of the Terrebonne Linked Tectonic System, and in isolation to reconstruct its history.  The fault probably formed as a low-relief slump feature on the continental slope some time before the start of the Miocene Epoch.  Sediment loading across the fault from the Miocene deltas caused episodic fault movement.  Fluids migrating along the fault plane during episodes of movement included hydrocarbons which had been cooked out of deeper sources rock.  The migrating hydrocarbons accumulated in Late Miocene sand layers forming the gas reservoirs of Bastian Bay field.  Episodic fault movement continued throughout the Pliocene and Pleistocene Epochs and into the Holocene and recent periods. 

Studies have shown that oil and gas reservoirs which have experienced significant depletion in reservoir pressure during extraction may have induced compaction of the reservoir sand resulting in subsidence at the surface (Chan and Zoback, 2007).  None of the reservoirs in the Bastian Bay field experienced adequate depletion in reservoir pressure to have induced compaction and subsidence.  Hydrocarbon extraction is not considered to be a viable cause for the subsidence measured at the surface in this area.  It is likely that the subsidence estimates made with data from the Grand Isle tide gauge reflect some contribution of recent movement on the Bastian Bay fault.

In the period between 1973 and 1975 a large area across the hanging wall of the Bastian Bay fault submerged and converted from mostly saline marsh to open water.  Gagliano et al (2003) estimated that the maximum total displacement on the fault at the surface was 3.0 to 3.5 feet based on his boring profiles and anecdotal evidence from oyster fishermen and camp owners.  Martin (2006) also quoted anecdotal evidence from Greg Linscombe of the Louisiana Department of Wildlife and Fisheries that the marsh break associated with the fault escarpment of the Adams Bay fault just to the north formed between 1976 and 1977.  The magnitude and pattern of these two fault slip events are similar, but they appear to be separated in time by 3 or 4 years.  The relationship between the Bastian Bay and Adams Bay fault slip events may provide insights into the nature of the connection between elements in the linked tectonic system.  It may be possible that the Bastian Bay fault slip event changed the stress fields in the near-surface in a way that affected and contributed to the cause of the Adams Bay fault slip event.

Picture13.png

The dramatic submergence of marshes on the hanging wall of the Bastian Bay fault can be seen in a comparison between 1973 and 1975 reconstructions of the marsh surface (from ProPublica, 2015).  This appears to have been the peak of the fault slip event, which is estimate of have extend across a period of about 18 years between 1970 and 1988 based on the records of the Grand Isle tide gauge.

Picture14.png
Picture15.png
Picture16.png

The Lake Boudreaux fault has a distinct northwest-southeast orientation within the linked tectonic system but its relationship with a rollover structure oil and gas field and a surface trace are very similar to the Bastian Bay fault.  The Lapeyrouse field was studied by Chan and Zoback, 2007 to investigate the relationships between oil and gas extraction and subsidence at the surface.  They determined that reservoirs in the Lapeyrouse field had experienced reductions in reservoir pressure due to extraction that were adequate to have induced compaction in the subsurface. They found that this could have caused subsidence at the surface.  Their models showed that subsidence associated with these processes would have been expressed as bowl-shaped depressions at the surface immediately above the outline of the depleted reservoirs.  The Madison Bay hot spot of wetlands loss lies along the surface trace of the Lake Boudreaux, and not above the depleted reservoirs.  There was no oil and gas production beneath the Madison Bay hot spot, and it is unlikely that hydrocarbon extraction contributed to subsidence in this area.  Subsidence in this area is more likely to have been caused by a slip event on the Lake Boudreaux fault.  Patterns of wetlands loss suggest that fault movement may have been happening nearly simultaneously on the Lake Boudreaux, Montegut and Isle de Jean Charles faults. 

Picture17.png

Morton et al, 2005 used boring profiles to study patterns of subsidence at the Madison Bay hot spot.  Their investigation shows that a layer of muddy sediment was deposited on top of the marsh surface beginning about 150 years ago.  The original thickness of the mud layer can be seen in the un-subsided portions of the profile.  Below this layer there is a layer of peat that indicates that the marshes maintained elevation through natural accretion and organic growth of marsh plants for a period from about 700 years ago until 150 years ago.  The magnitude of subsidence in the Madison Bay hot spot can be seen in the current elevation of boundary between the mud layer and the peat layer relative to the elevation of the original marsh surface.  The approximately 80 centimeters of subsidence inferred by this method is equivalent to about 2.6 feet.  This is roughly comparable to the value estimated by Gagliano at the Bastian Bay and Adams Bay faults, and it appears to have happened within the same time frame.  As will discussed in the section on possible causes of the fault slip event, the deposition of layer of muddy sediment on the marsh surface over the past 150 years may be associated with the same sedimentation anomaly responsible for the thick Holocene sediments at the eastern end of the Terrebonne Trough.  Sediment loading across the fault may have contributed to the cause of the fault slip event.

LATE 20TH CENTURY SUBSIDENCE EPISODE

Evidence for the patterns of recent movement on the Lake Boudreaux and Bastian Bay faults can be seen on faults throughout the Terrebonne Linked Tectonic System.  There is good evidence that the elements of system have moved interactively throughout the Cenozoic Era into the Holocene and recent.  It appears that as the fault slip event propagated throughout the system it created an episode of increased subsidence rates across the delta plain that was recorded in historical tide gauge records (Penland et al, 1988) and the re-surveying of elevation benchmark (Shinkle & Dokka, 2004).

Picture18.png
Picture19.png

Estimates of current subsidence rates were made using static GPS measurements by Byrnes et al, 2021.  These are direct measurements of the vertical velocity of the surface of the earth at the sites of the benchmarks.  The total value of subsidence measured at the surface is the cumulative result of the contribution of multiple processes.  Byrnes suggested that the very strong relationship between the magnitude of the subsidence values and the thickness of the Holocene topstratum at the site of the benchmark indicate that the compaction of the Holocene sediment is the primary driver of current subsidence rates.  It is likely that other processes including deep lithospheric flexure of the mantle, continuing compaction of Pleistocene through Miocene sediments and Glacial Isostatic Adjustment are contributing to the total value of subsidence being measured.  It is also likely that the magnitude of subsidence attributable to these processes is small, and that it would not vary significantly across the area.

None of the processes listed above provides a viable mechanism to account for the episode of increased subsidence rates measured in the late 20th century.  A sampling of the values of subsidence rates estimated from Penland et al, 1988 and Shinkle & Dokka, 2004 are shown in blue and green, respectively.  In some cases the elevated rates of subsidence are two to three times the value of current subsidence rates in the same area.  The most viable explanation for the episode of increased subsidence rates in the late 20th century is a fault slip event in the Terrebonne Linked Tectonic System.

Picture20.png

The temporal relationship between rates of subsidence and rates of wetlands loss over the past century is striking.  It appears likely that natural subsidence has been the primary driver of wetlands loss over that time period.  This would be consistent with the inferred relationships between subsidence and wetlands loss in the delta lobes delineated by Frazier, 1967 that are now submerged 20 to 30 feet below the surface.  The strength of the relationship between subsidence and wetlands loss also implies that patterns of wetlands loss may be used to add detail to the values of subsidence that have been measured at discrete locations over the past few centuries.

THE NATURE OF FAULT SLIP EVENTS

Very little is known about the nature of fault slip events in the Louisiana delta plain.  The measurement of accurate current subsidence rates has only been taking place for about two decades.  The density of measurement stations is much too low to have detected variations on subsidence rate that may be attributable to fault slip.  Earthquakes are periodically recorded in the state, and they very likely indicate fault slip events of some kind.  The only direct observation of a fault slip event that ruptured the surface occurred at Vacherie between April 12th and 15th in 1943.  The magnitude of that slip event was measured at 8 inches.  Shen et al, 2017 did precise dating of sedimentary layers and accurate elevation measurements of those layers on either side of faults in the Baton Rouge fault system.  They were able accurately estimate the average rate of fault slip over several millennia.  The total cumulative offset of one Pleistocene layer due to fault slip was estimated to be about 1100 millimeters over 32,000 years or an average rate of fault slip of 0.034 mm/year.  The offset of a Holocene layer on a different fault was measured by Shen et al at about 600 mm over 2780 years or an average slip rate of 0.216 mm/year.  The difference between these two values can be attributed to an inverse relationship between the magnitude of average rate of fault slip and the time span over which it was measured (Sadler, 1981).

Picture23.png

A semi-quantitative distribution of fault slip events may be constructed by combining the data from Shen et al, 2017 with Gagliano’s observation of 3 feet of displacement at the Bastian Bay fault, which is spread over a span of 18 years based on the span of increased subsidence rates measured by tide gauge and benchmark data, and the displacement of the Vacherie fault over 3 days in 1943.  What is more important than the derivation of any actual data from this very sparse population is the characterization of the nature of fault slip events.  The values estimated by Shen et al, are averages of a total population of fault slip events that happened on those faults over periods of 32,000 and 2,780 years respectively.  The average values of fault slip rate estimated by Shen et al are also the mean value of the population. The total population of fault slip events responsible for the total offset probably has a size distribution similar to the semi-quantitative power-law distribution in the figure.  In other words the total displacement of the Pleistocene layer across the fault measured by Shen et al, was the result of many individual fault slip events over the span of 32,000 years.  If slip events on this fault happened once a year there would be 32,000 individual events in this population.  If slip events happened once a century, there would be 320 in the population.  It is likely that rather than being a population of events of equal magnitude, that there would be a distribution of different magnitudes.  If this distribution followed a power law, it would consist of many smaller events and very few larger events.  Large slip events would be very rare, and small slip events would be more common.  The magnitude and frequency of the small slip events may be so low as to be manifest as a slow creeping motion of the fault rather than a distinct slip event.  If the pattern of accelerated subsidence in the late 20th century as seen in tide gauge data in the delta plain is a reflection of the nature of fault slip that caused it, then it is likely that rate of fault slip increased somewhat gradually over time rather than there being distinct slip events.  In other words the movement on faults in the Terrebonne Linked Tectonic System may have been a gradual creeping motion, and the rate of creep increased for a period of several years.

 

POSSIBLE CAUSE OF THE FAULT SLIP EVENT

A possible cause of the late 20th century fault slip event in the Louisiana delta plain is suggested by the sites of the subsidence hot spots relative to the distributary channels of the last few historical deltas and the surface traces of the major faults.  The timing of the fault slip event suggests a slow build-of surface stresses due to abnormal sedimentary loading in the 18th and 19th centuries and a release of that stress caused by a triggering event in the last three decades of the 20th century. 

Picture24.png

The map of the 1874 flood of the Mississippi River (Hardee, 1874) shows that without artificial levees major floods covered the entire delta plain.  Sediment-laden flood waters were carried throughout the entire distributary channel network.  It also shows that the modern Balize delta lobe has been built over the last two centuries.  Tweel and Turner, 2012 showed that agricultural expansion across the Mississippi River basin caused a dramatic increase in the suspended sediment lobe of the river at the beginning of the 19th century.  The anomalously thick accumulation of sediment in the Balize delta lobe seen on the Holocene topstratum thickness map is probably a direct result of the increased sediment load of the river.  The Hardee map suggests that sedimentary loading would have been spread across the delta plain.  It is likely that the layer of muddy sediment seen in the boring profile at the Madison Bay hot spot is a function of increased rates of sedimentation over the last 150 years.

Picture25.png
Picture26.png

The most logical sites for increased sediment accumulation in the delta plain during the period of higher suspended sediment load of the Mississippi River would have been off the flanks of the natural levees of the distributary channel networks where those networks branch out.  Because of the genetic relationship between the Holocene delta lobes and the faults of the Terrebonne Linked Tectonic System, the distributary channel networks of each delta lobe tend to branch out where they cross the major down to the south faults.  Sediment-laden flood waters would have fanned out across these branching channel networks and deposited their load on the flanks of the natural levees and across the hanging walls of the faults.  There is a very strong correlation between the sites where branching distributary channel networks cross the major faults and the hot spots of wetlands loss in the late 20th century.  It is likely that the accumulated sediment in these areas, such as the muddy layer seen in the Madison Bay borings, would have caused stress differentials on either side of the major faults.  It appears that these stress differentials were released in the fault slip event brought on by a trigger mechanism.

Historical tide gauge data has been used around the world to measure “relative sea level rise”, which is the combination of global “eustatic” sea level rise and subsidence.  By knowing the rate of global sea level rise it is possible to estimate a history of the rate of subsidence at the site of the tide gauge (Boon et al, 2010).  Blum and Roberts, 2012 compared water levels from the Pensacola and Grand Isle tide gauges to the established global mean sea level.

Picture22.png

The Pensacola gauge tracks very closely to the global mean sea level, but the Grand Isle gauge appears to indicate that sea level was rising nearly four times faster than it was at Pensacola.  This is because the rate of subsidence a Pensacola is near zero, and so the gauge is recording actual global sea level.  Conversely, the Grand Isle gauge began subsiding at an accelerated rate in the latter part of the 20th century, and the apparent higher water levels are due to subsidence. 

The level of the Grand Isle gauge relative to the station datum is shown in the figure below.  Increasing water level anomalies indicate increasing rates of subsidence.  These anomalies peak in the 1970s and 1980s, then begin to decline.  The relationship between the subsidence anomaly and the episode of land loss in the Terrebonne and Barataria basins is obvious, but there also appears to be a relationship between the discharge of the Mississippi River and the rate of subsidence measured by the tide gauge.

Picture21.png

This relationship may suggest a trigger mechanism for the fault slip event.  Increases in river discharge coincide with the elevation of the river stage within the channel.  These are major flood events on the Mississippi River.  Elevated river stages are also accompanied by the exchange of groundwater in near surface aquifers and the river channel.  Kolker et al, 2013 showed that groundwater fluxes during high river stages extended across the delta plain in the near-surface aquifers.  Because near-surface aquifers are the natural levee deposits of the last few historical deltas, the flux of increased groundwater flow during periods of high river discharge would have extended throughout the distributary channel network.  Increased flow of groundwater in the natural levee sands would have tended to dilate the aquifers and affect the stress fields in the area.  Changes in stress would have been concentrated in the distributary channel networks at the exact sites where sediment loading has caused stress differentials across the major faults.  In this manner the periods of high river discharge could have served as a trigger mechanism for a fault slip event.  In simplest terms the dilation of near surface aquifers may have served to “lubricate” the faults allowing for increases in the rate of slip movement. The increased rates of fault slip across the Terrebonne Linked Tectonic System would have been measured as increased rates of subsidence in the tide gauge and benchmark leveling surveys.

The geological history of the Terrebonne Linked Tectonic System is one of continual episodic movement in response to sediment loading from the deltas of the Mississippi River.  Throughout the Miocene Epoch and into the late Holocene there is evidence of fault movement in response to sediment loading.  It is logical that an increased suspended sediment load of the river during the 18th and 19th century would have resulted in loading at the sites where the distributary channels crossed the major faults.  A possible trigger mechanism for a fault slip event may have been caused by changes in the stress fields around the sites of sediment loading caused by the dilation of near-surface aquifers during periods of high river discharge.  This set of processes would accurately account for the location of the hot spots of a late 20th century episode of subsidence.  Recognizing the relationship between this episode of subsidence and wetlands loss in the delta plain has significant implications for coastal sustainability planning.

 

REFERENCES

 

Akintomide, A. O. and Dawers, N. H., 2019, Spatial and Temporal Variation of Fault Activity in the Terrebonne Salt Withdrawal Basin, Southeastern Louisiana: Response to Salt Evacuation and Sediment Loading, presentation, AAPG Annual Meeting, San Antonio, Texas.

Armstrong, C., Mohrig, D., Hess, T., George, T., Straub, K.M.,  2014, Influence of growth faults on coastal fluvial systems: Examples from the late Miocene to Recent Mississippi River Delta, Sedimentary Geology, v. 301, p. 120-132

Blum, M.D. and Roberts, H.H., 2012, The Mississippi Delta Region: Past, Present, and Future, Annual Review of Earth Planetary Science, v. 40, p. 655-83

Boon, J.D., Brubaker, J.M., Forrest, D.R., 2010, Chesapeake Bay Land Subsidence and Sea Level Change, Virginia Institute of Marine Science, Special Report No. 425, 81 p.

Bullock, J. S., Kulp, M. A., McLindon, C. D., 2018, Evaluation of the Magnolia growth fault, Plaquemines Parish, southeastern Louisiana, poster session, GSA Annual Meeting, Indianapolis, Indiana.

Burke, L., Kinney, S.A., Dubiel, R.F., and Pitman, J.K., 2013, Geopressure gadient maps of Southeastern Louisiana, state and vicinity, U.S.G.S.

Byrnes, M. R., Britsch, L. D., Burlinghoff, J. L., Johnson, R., Khalil, S., 2019a, Recent Subsidence rates for Barataria Basin, Louisiana, Geo-Marine Letters, 14 p. doi.org/10.1007/s00367-019-00573-3.

Byrnes, M. R., Britsch, L. D., Burlinghoff, J. L., Johnson, R., Khalil, S., 2019b. Determining Recent Subsidence Rates for Breton Sound and Eastern Pontchartrain Basins, Louisiana: Implications for Engineering and Design of Coastal Restoration Projects. Final Report prepared for Louisiana Coastal Protection and Restoration Authority. Contract 4400009020, Task 8, 58 p.

Chan, A.W. & Zoback, M.D., 2007. The role of hydrocarbon production on land subsidence and fault reactivation in the Louisiana coastal zone. Journal of Coastal Research, 23(3), 771–786

Culpepper, D.B., McDade, E. C., Dawers, N. H., Kulp, M. A., Zhang, R., 2019, Synthesis of Fault Traces in SE Louisiana Relative to Infrastructure,  TranSET Project No. 17GTLSU12

Curtis, D. M., 1970, Miocene Deltaic Sedimentation, Louisiana Gulf Coast, S.E.M.P. Special Publication 15, Deltaic Sedimentation, Modern and Ancient, p. 293-308

Dawers, N. 2018, Improving our Understanding of Faults and Salt Dome Interactions in Southeastern Louisiana’s Subsurface, State of the Coast, New Orleans, LA

Dokka, R. K., G. Sella, and T. H. Dixon, 2006, Tectonic control of subsidence and southward displacement of southeast Louisiana with respect to stable North America, Geophys. Res. Lett., 33, L23308, 5 p.

Dokka, R.K., 2011, The role of deep processes in late 20th century subsidence of New Orleans and coastal areas of southern Louisiana and Mississippi, Journal of Geophysical Research, v. 16, 25 pgs.

Fisk, H.N., 1944, Geological Investigation of the Alluvial Valley of the Lower Mississippi River, Mississippi River Commission.

Frank, J. P., 2017, Evidence of fault movement during the Holocene in Southern Louisiana: integrating 3-D seismic data with shallow high resolution seismic data, MS Thesis, University of New Orleans, 91 p.

Frazier, D.E., 1967, Recent deltaic deposits of the Mississippi River: their development and chronology, Trans. G.C.A.G.S., v. 17, p. 287‐315

Gagliano, S.M., Kemp III, E. B., Wicker, K. M., Wiltenmuth, K. S., 2003. Active Geological Faults and Land Change in Southeastern Louisiana. Prepared for U.S. Army Corps of Engineers, New Orleans District, Contract No. DACW 29-00-C-0034.

Hiatt, M., Snedden, G., Day, J.W., Rohli, R.V., Nyman, J.A., Lane, R., and Sharp, L.A., 2019, Drivers and impacts of water level fluctuations in the Mississippi River delta: Implications for delta restorations, Estuarine, Coastal and Shelf Science, v. 224, p. 117-37

Hopkins, M., Lopez, J., Songy, A. 2018, Subsidence rates from faulting determined by real-time kinematic (RTK) elevation surveys of bridges in Lake Pontchartrain,  presentation, State of the Coast Conference 2018, New Orleans, Louisiana.

Johnston, A., Zhang, R., Gottardi, R., Dawers, N. H., 2017, Investigating the relationships between tectonics and land loss near Golden Meadow, Louisiana by utilizing 3-D seismic and well log data, poster session, GSA Annual Meeting, Seattle, Washington.

Karegar, M.A., Dixon, T.H., Malservisi, R., 2015, A three-dimensional surface velocity field for the Mississippi River Delta: Implications for coastal restoration and flood potential, Geology, v. 43, p. 519-522

Kolker, A.S., Allison, M.A., and Hameed, S., 2011, An evaluation of subsidence rates and sea-level variability in the northern Gulf of Mexico, Geophysical Research Letters, v. 38, 6 p.

Kolker, A.S., Cable, J.E., Johannesson, K.H., and Inniss, L.V., 2013, Pathways and processes associated with the transport of groundwater in deltaic systems, Journal of Hydrology, v. 498, p. 319-34

Kuecher, G.J., Roberts, H.H., Thompson, M.D. and Matthews, L., 2001, Evidence for Active Growth Faulting in the Terrebonne Delta Plain, South Louisiana: Implications for Wetland Loss and the Vertical Migration of Petroleum., Environmental Geosciences, v. 8, p. 77‐94

Kulp M. 2000. Holocene stratigraphy, history, and subsidence of the Mississippi River delta region, north-central Gulf of Mexico. PhD thesis. Univ. Kentucky, Lexington. 283 pp.

Levesh, J. L,, Kulp, M. A., McLindon, C. D., 2019, Fault-slip history of the Delacroix Island fault system and its effect on Holocene salt marshes of the Mississippi River delta plain,  presentation, GSA Annual Meeting, Charleston, South Carolina

Lopez, J., Penland, S., Williams, J., 1997, Confirmation of Active Geologic Faults in Lake Pontchartrain in Southeast Louisiana, GCAGS Transactions, v. 47, p. 299-303

Martin, E., 2006, Fault Induced Subsidence Near Empire and Bastian Bay, Louisiana, MS Thesis, Tulane University, 163 p.

McLindon, C.D., 2017, History of fault slip and interaction with deltaic deposition from the middle Miocene to the Present – Barataria fault, coastal Louisiana, poster session, American Geophysical Union, annual meeting, New Orleans, Louisiana.

McLindon, C.D., Dawers, N.A., Culpepper, D., Kulp, M.A., and McDade, E., 2017, Comments to the U.S. Army Corps of Engineers, New Orleans Division in reference to the The Environmental Impact Statement for the Mid-Barataria Sediment Diversion.  11 pg.

Morton, R.A., Bernier, J.C., Barras, J.A, and Ferina, N.F., 2005, Rapid Subsidence and Historical Wetland Loss in the Mississippi Delta Plain: Likely Causes and Future Implications U.S. Geological Survey Open-File Report, 2005–1216, 124 p.

Penland, S.,, Ramsey, K.E., McBride, R.A., Mestayer, J.T. and Westphal, K.A., 1988, Relative Sea Level Rise and Delta-Plain Development in the Terrebonne Parish Region, LA Geological Survey Coastal Geology Technical Report No. 4, 140 p.

ProPublica, 2015, Losing Ground

Shinkle, K. D., and R. K. Dokka, 2004, Rates of vertical displacement at benchmarks in the Lower Mississippi Valley and the northern Gulf Coast, NOAA Tech. Rep. 50

Stadler, P.M., 1981, Sediment Accumulation Rates and the Completeness of Stratigraphic Sections, The Journal of Geology, Vol. 89 p. 569-584

Thorne, C., Harmar, O., Watson, C., Clifford, N., Biedenharn, D., and Measures, R., 2008, Current and historical sediment loads in the Mississippi River, Report to the European Research Office of the U.S. Army, 170 p.

Tweel, A.W. and Turner, R.E., 2012, Watershed land use and river engineering drive wetland formation and loss in the Mississippi River birdfoot delta, Limnology and Oceanography, v. 57, p. 18‐28.

Zou, L., Kent, J., Lam, N. S.-N., Cai, H., Qjang, Y., and Li, K., 2016, Evaluating Land Subsidence Rates and Their Implications for Land Loss in the Lower Mississippi River Basin, Water, v. 8, 15 p.

 

A Geological Evaluation of the vicinity of the Mid-Barataria Sediment Diversion

This evaluation was submitted as comments on the Draft Environmental Impact Statement for the Mid-Barataria Sediment Diversion. They propose that a thorough geological evaluation should be conducted at the site to determine the potential for sediment loading by the project to cause a major fault slip event.

These comments are being submitted on the Draft Restoration Plan and Environmental Impact Statement: Mid-Barataria Sediment Diversion (the “Draft EIS”).  Comments submitted prior to drafting this Statement (McLindon et al, 2017) recommended: ”that a thorough subsurface geological evaluation of the vicinity of MBSD [Mid-Barataria Sediment Diversion] be conducted to attempt to determine the location of geological faults, the recent history of fault movement and the effects of active faults on subsidence rates and variations in the thickness of highly compactible soils.”  This recommendation included evaluation including the use of oil and gas industry seismic data, the acquisition of high resolution seismic data, the acquisition of sediment cores, and the development of subsidence measuring capabilities.  The Draft EIS concluded in response to this recommendation that “there is insufficient information on which to evaluate the impact of faulting on the proposed Project or the impact of the proposed Project on the future fault movement.”  It was the intention of the McLindon et al (2017) recommendation to collect the information necessary to make such an evaluation.

The comments submitted here are in the form of a geological evaluation of the vicinity of the MBSD.  There are two principal objectives to this evaluation:

1.       To provide a framework within which to develop the ability to forecast a probabilistic distribution for the frequency and magnitude of fault slip events in the vicinity of MBSD and other large infrastructure projects on the coastal plain, and

2.       To provide insights into the historical relationships between sediment loading, fault slip events, and fault-induced subsidence that may be incorporated into predictive models for land elevation and land area gain in the MBSD project area.

This evaluation has been constructed with the best available data, and is in no way intended to replace or substitute for the thorough geological evaluation recommended in McLindon et al (2017).  The subsurface geological maps used here represent a compilation of geological interpretations over many decades.  Inputs to the interpretation have included well logs, biostratigraphic data from micropaleontology, seismic data, gravity data, published subsurface geological interpretations from peer-reviewed technical literature and atlases constructed by the New Orleans and Lafayette Geological Societies, and interpretations submitted to the Louisiana Office of Conservation in support of oil and gas unitization

The Ironton fault

The Ironton fault is the most likely geological feature to have a direct impact on the MBSD Project.The fault plane has been mapped in the subsurface, and it is expressed on subsurface structural contour maps as a fault trace delineating the intersection of the fault plane and the mapped stratigraphic horizon.The fault can also be seen in cross sections constructed with oil and gas well logs, as well as 2-D and 3-D seismic data.Gagliano et al (2003) described this fault and other faults and salt domes in the area as being part of a “linked tectonic system”.This implies that individual elements within the system may be impacted by activity on other elements such as episodic fault slip events, diapiric salt movement (halokenesis), or salt dissolution.It also implies that faults within the system are likely to share similar characteristics.There is not adequate data available on the Ironton fault to derive the objective framework for modeling fault slip activity and its impacts, so data from other faults in the system and within the basin will be used here to construct the framework.

Picture1.png

This evaluation will consider aspects of the ten faults labeled here to derive a set of characteristics that may be applied to the Ironton fault.  Taken together this set of characteristics can be used to make reasonable estimates of a probabilistic distribution of the frequency and magnitude of episodic slip events on the Ironton fault. The map of surface fault traces shown in Fig.1 is a compilation of numerous individual interpretations including university theses and dissertations and peer-reviewed publications (Akintomide & Dawers (2019), Armstrong et al (2014), Bullock et al (2018), Culpepper et al (2019), Frank (2017), Johnston et al (2017), Levesh et al (2019), McLindon (2017)).  A portion of this map is available in a GIS application on the Louisiana Department of Transportation and Development website. 

A subsurface contour map of the Ironton fault plane is shown in Fig. 2 along with subsurface contours on the top of salt at the Lafitte and Lake Hermitage salt domes.A fault plane map is constructed by integrating the interpretation of biostratigraphy, well log correlation and seismic data.The values in blue in Fig. 2 are the depths at which the fault can be interpreted on a well log by the “missing section” in the log relative to surrounding well logs.A cross section of six well logs across the fault plane is shown in Fig. 3.  Wells 2, 3

Picture2.png
Picture3.png

and 4 exhibit missing stratigraphic section due to the fault at depths of -2175, -9,195', and -10,512’ respectively.  These values are integrated with those from other wells that intersect the fault plane along with seismic data interpretation to construct the fault plane map.  The tops of the biostratigraphic intervals are shown as colored lines on the cross section.  Each of these interval tops is correlated between the well logs and could be a potential mapping horizon.  Mapped horizons at the top of the Pliocene and top of the Mid-Miocene Epoch are indicated on the cross section.  The subsurface structure maps for these horizons are shown in Figs. 6 and 7, respectively.  The vertical change in elevation across the fault, or fault throw, is shown for 5 horizons.  The throw of the fault is 35’ and 840’ for the top of Pliocene and top of Mid-Miocene horizons respectively.  These values could be measured from the cross section or by comparing the values of the subsurface elevation contours on either side of the fault on the subsurface structure maps.  The throw of the fault continually increases with depth in the classic form of a Gulf Coast “growth fault”.  Increasing throw with depth indicates that the fault has been continually, if episodically, moving throughout the geological timespan measured by biostratigraphic control.

Picture4.png

The interval thickness between each horizon is shown for a well on the upthrown, or foot wall side of the fault, and one on the downthrown or hanging wall side.  The ratio of the hanging wall thickness to the footwall thickness for any interval is called the “expansion index”, and it is used gain insight into historical fault movement.  The deepest correlated interval is 485’ thick on the foot wall side and 704’ thick on the hanging wall side.  The expansion index across the fault at this interval is 1.45, or the hanging wall side is 145% thicker than the foot wall side.  The cross section also shows that the Ironton fault has a classic listric shape, which flattens with depth.  The magnitude of the horizontal component of fault displacement increases with depth relative to the vertical component as the fault plane flattens.  This illustrates the significant lateral movement that is associated with fault slip over time.  It is likely that the horizontal vectors of movement measured at CORS stations in Fig. 33 are due to fault slip.

Fig. 4 is a portion of an oil and gas industry 2-D seismic profile across the Ironton fault.  Seismic data can be integrated with subsurface well log interpretations in the construction of the fault plane map and subsurface structure maps within the interval of biostratigraphic control.  It can also be used to map the fault plane and subsurface structure below this interval, allowing for the projection of fault plane contours to depths of up to 20,000 feet.  Oil and gas industry seismic data is acquired and processed to optimize imaging of the subsurface between depths of about 2,000 and 15,000 feet below the surface, where most oil and gas is found.  Seismic profiles such as that shown in Fig. 4 are less than optimal for imaging the near-surface extent of a fault, and interpretations generally have to be extrapolated to the surface.  The acquisition of new high resolution seismic data across the Ironton fault recommended in McLindon et al (2017), and reiterated here, is intended provide the necessary imaging to evaluate the near-surface extent of the fault.

Picture5.png

Fig. 5 is a graph of the variation of the expansion index across the Ironton fault through the span of geological time between the Mid-Miocene and Pliocene mapping horizons.  The expansion index is generally above the value of 1 throughout most the time span indicating continual movement of the fault.  As will be discussed in more detail in a subsequent section, the higher values of expansion index on the fault can be related to periods of increased sedimentary loading associated with an active delta system.  Graphs like this can be used to assign a generalized relative quantitative value to the distribution of fault slip events over geologic time.  It is likely that the history of fault movement has consisted of a distribution of individual events over time whose variation in magnitude is similar to other natural phenomenon such as earthquakes, floods and hurricanes.  Unlike tectonic fault movement associated with earthquakes, the movement of listric faults in a passive margin setting is generally aseismic.  The relative distribution of the magnitude of slip events over time may, however, be similar.  In other words, fault slip on these faults may consist primarily of low magnitude slip events that manifest as a creeping motion on the fault, but larger events may occur sporadically over time.  The primary intention of a thorough geological evaluation is to attempt to put ranges on the magnitude and frequency of slip events on the Ironton fault.

Picture6.png

The Ironton fault plane is represented as a colored grid overlain on the subsurface structure map for the top of Mid-Miocene mapped horizon in Fig. 6.  The subsurface trace of the fault is the brown polygon on the structure map based on the lateral component of displacement.  Black squares indicate the down-dropped side of the fault.  Well 3 from the cross section falls within the fault trace polygon for the Ironton fault on this horizon meaning that the horizon is “faulted out” of the well.  It can be seen on the cross

Picture7.png
Picture8.png

section that the Mid-Miocene mapped horizon is not present in the well, as it has been displaced vertically and laterally by the fault.  A comparison of the Mid-Miocene map in Fig. 6 and the Pliocene map in Fig. 7 shows how the trace of the fault on the map reflects the decrease in throw with depth.  Subsurface geological structures generally become simpler and flatter at shallower depths.  At the depth of the Pliocene mapped horizon only the major faults are still active.  These are the same faults that appear to reach the surface.  The 0’-depth contour of the Ironton fault plane is coincident with the surface trace of the fault shown in Fig. 1.  The most significant impacts of a potential future fault slip event should be expected along this trace, and future geological evaluation of the fault should be focused on more fully delineating the near-surface fault plane with high resolution seismic and sediment cores.  Fig.8 also shows the connection of the Ironton fault plane to the Lafitte and Lake Hermitage salt domes.  These are elements of the linked tectonic system described by Gagliano et al (2003).  There is likely to be a genetic relationship between the structural evolution of the fault and the salt domes.  It is probable that halokentic salt movement on the domes has affected the history of movement on the fault.  It is also possible that the dissolution of salt at the domes may have contributed to triggering fault slip events.  Halokenesis and salt dissolution should be considered as potential causes for future slip events on the Ironton fault.

A near-surface expression of the Ironton fault can be seen in a profile of cone penetrometer tests (CPT) taken along the proposed conveyance structure for the MBSD.  Interpreted images of a profile including the same tests were published in the Geotechnical Baseline Report for 30% Design (CPRA, 2014).  These interpreted images indicated the correlation of the top of the Pleistocene interval with an overlying deposit of “Near Shore Gulf” sand shown in orange in Fig. 9.  This sand layer is likely to represent a basal Holocene transgressive unit associated with the first submergence of the Pleistocene surface by rising sea level.

Picture9.png
Picture10.png

Frazier (1967) interpreted borings taken along the Barataria Waterway, and defined the top of deposits associated with the Bayou Families Delta at between 30 and 40 feet below the surface in this area.  The sandy, silty layer in this depth range on Fig. 9 is likely to be associated with this delta lobe.  These approximate ages may be used to generally approximate rates of fault slip on the Ironton fault, however much more accurate age determinations from sediment cores used in association with high resolution seismic profiles would be necessary to provide estimates accurate enough to be used in predictive modeling. 

Fig. 9 indicates 6 feet of throw on the Ironton fault at the top of the Pleistocene surface.  If the more detailed geological evaluation recommended here supports this interpretation, it may concluded that the 6 feet of throw is consistent with the pattern of throw versus depth seen on the cross section in Fig. 3.  It is likely that the magnitude of this throw is the result of the cumulative effect of multiple fault slip events that have occurred since the deposition of the basal transgressive sand.  Fig. 10 shows the fault plane contours of the Ironton fault crossing the CPT profile at the interpreted location of the fault on the profile.  Oil and gas wells in the area are identified by their Louisiana DNR serial number, and the depth of the Ironton fault is indicated in blue.

Characteristics of other faults in the area

Each of the nine other faults shown in Fig. 1 will be considered here to provide a collective characterization that may be used as a framework within which to develop the ability to forecast a probabilistic distribution for the frequency and magnitude of fault slip events on the Ironton fault.  Each of these faults exhibits some indication of recent fault movement.  Some reveal the potential impacts of faults on infrastructure, and some document differential rates of Holocene sediment accumulation across the fault.

Picture11.png

The Vacherie fault (Fig. 11) provides the highest quality data point for the magnitude of a recent individual fault slip event.  The event occurred between April 12th and 15th, 1943, and the magnitude of the vertical displacement was reported by Fisk (1944) to be 8 inches.  There was some reported local ground shaking, but there was no detection of seismic activity at the Loyola University seismometer only 50 miles away.  This appears to have been an aseismic fault slip event.  Fisk used a boring profile across the fault escarpment to measure 3 feet of vertical displacement across the fault at the top of the Pleistocene.  He interpreted this to indicate recurrent movement on the fault.  Fisk also noted “Of particular interest is the fact that faulting occurred at a place along the river where repeated crevassing took place.”

It is probable that the 3 feet of displacement at the top of the Pleistocene is representative of the cumulative total of all of the recurrent fault slip events throughout the Holocene.  If each recurring event produced a displacement of 8 inches, then there could have been three previous events over the 7,000 year span that is likely to be represented by the Holocene sediments.  Some of the displacement could have been taken up by compaction of the poorly consolidated sediments, which may allow for more events of the same magnitude, but it may also be true that larger magnitude events in the past would have been necessary to cause a crevasse of the river.  If the history of movement on the Vacherie fault could be accurately reconstructed, it would likely show a log-normal distribution in the magnitude of the fault slip events.  Subsurface mapping indicates that the fault has been active since at least the early Miocene.  Movement since the deposition of the first Holocene sediments would probably be represented by many small slip events that would have propagated a creeping movement on the fault, and a few large events, perhaps with magnitudes of a foot or more vertical displacement adequate to cause a crevasse of the river. 

The correlation of the fault and the site of the historical crevasses of the river is significant.  The 1943 fault slip event occurred during flood stage on the river.  It is likely that the dilation of near-surface aquifers near the river channel during flood stages may alter the stress fields within the near surface sedimentary layers, and provide a potential trigger mechanism for fault slip events.  The 6 feet of vertical displacement at the top of the Pleistocene seen on the CPT profile in Fig. 9 is twice the value seen on the Vacherie fault.  This may suggest that the Ironton fault has had more and higher magnitude fault slip events during the Holocene.  There is no evidence of historical crevasses of the river at the Ironton fault, but some effort should be made to evaluate the magnitude of an event that would be necessary to cause a crevasse of the river, and the probability that such an event could occur using predictive modeling based on a likely distribution of the magnitude of historical events.  This is particularly true given that differential sedimentary loading across the fault, as intended by the MBSD project, may provide a trigger mechanism for a fault slip event.

Picture12.png

The Highway 11 fault (Fig. 12) provides the most tangible evidence for the impact of a fault on infrastructure in the area.  The fault also provides valuable insights because it is the single place in south Louisiana where a fault is crossed by a high resolution seismic line capable of measuring near-surface fault throw on multiple horizons (Lopez et al, 1997).  Vertical displacement on the fault at the surface has also been measured by elevation surveys across the Highway 11 Bridge (Hopkins et al, 2016).

The seismic profile across the fault exhibits high-quality imaging of Pleistocene sedimentary layers to a depth of about 150 feet with a vertical resolution of less than a few feet.  This type of imaging capability is necessary to evaluate faults in the near-surface, and is recommended for a geological evaluation of the MBSD project area.  Lopez et al (1997) provided interpreted correlation of seismic horizons associated with sedimentary layers across the fault.  The fault exhibits increasing throw with depth and thickening of stratigraphic intervals on the hanging wall side up to the surface.  It is likely that this pattern extends to all faults delineated on Fig. 1, and it suggests continual episodic movement on the faults throughout the Quaternary.  If this seismic profile had been combined with high resolution dating of sediment cores, it would have been possible to reconstruct a detailed history of fault movement.  A reconstruction of fault movement from this type of data is likely to reveal averaged movement over intervals of time rather than individual episodic fault slip events, but it would be valuable for putting reasonable ranges on the magnitudes of individual events within a given interval. 

The elevation profiles from Hopkins et al (2016) indicate the magnitude and span of displacement across the fault, and valuable estimates of average rates of fault slip.  The lower of the two elevation profiles in Fig. 12 shows the entire span of the Highway 11 Bridge.  It is important to note that while there is a localized increase in offset immediately adjacent to the fault, the entire hanging wall side of the fault has a lower elevation than the footwall side.  This indicates that subsidence due to cumulative fault slip has affected an area of up to 100 square miles or more.  The average rate of subsidence associated with fault slip is between 1.3 and 3.0 mm year over the time span measured.  This is likely to be the result of cumulative slow slip movement, as there have been no reported episodic events on the fault over the past few decades.  Lopez et al (1997) did however document the occurrence of two small earthquakes associated with the South Point fault near the end of the bridge in 1987.

The St. Rose fault (Fig. 13) also appears to exhibit the impact of faults on infrastructure.  The surface trace of the fault is clearly delineated by a sharp tree line in the cypress swamp.  Trees on the hanging wall side of the fault have mostly died.  The plane of the fault can be mapped in the subsurface with well log correlation and 3-D seismic data interpretation.  One possible explanation for the death of the trees is the migration of saline fluids to the surface along the fault plane.  A similar configuration exists at the Montegut fault (Fig. 16) where Kuecher et al (2001) measured a distinct anomaly in total dissolved solids in soils adjacent to the fault.  Saline fluid migration has also been documented on faults in the Baton Rouge fault system. 

The clearly defined surface trace of the St. Rose fault very closely coincides with two apparent impacts on infrastructure.  An elevation survey on Highway 626 shows a vertical offset on the road bed.  The trace of the fault also coincides with a crack and vertical offset in a T-wall in the St. Rose drainage structure. 

Picture13.png

The Gentilly fault (Fig. 14) was identified in a boring profile by Fisk (1944).  The vertical offset at the top of the Pleistocene (blue on the profile) is nearly 20 feet.  As a result, the Holocene sediments, and in particular the organic clay and peat deposits near the surface are thicker on the hanging wall side of the fault than they are on the foot wall side.  A map of peat thickness across the area by Gould and Morgan (1962) shows that

Picture14.png

the thickest accumulation of peat in the area is on the hanging wall side of the fault.

Zou et al (2015) used US Coast and Geodetic Survey data from repeat surveys of elevation benchmarks in the area to delineate a high subsidence anomaly that coincides the thick accumulation of peat on the hanging wall side of the Gentilly fault.  This is likely to be a result of a feedback loop between fault slip and sediment accumulation.  Subsidence rates are logically higher across an area of thick peat accumulation because sedimentary layers composed primarily of poorly consolidated organic material are likely to exhibit higher rates of compaction resulting in subsidence expressed at the surface.  A consideration of why the peat layer is thick in the first place would lead to the recognition that peat accumulation is itself a response to subsidence.  The accumulation of peat has both a cause and effect relationship with subsidence in the area.  Marsh and swamp ecosystems on a delta plain maintain elevation through organic growth with some input of mineral sediment from river flooding or tidal flux.  As plants in the ecosystem die and are submerged new plants grow on top of them.  The rate of the accumulation of organic material is controlled by the rate of subsidence.  Areas with the high subsidence rates, such as the hanging walls of faults, tend to have thicker peat accumulations.

The total vertical displacement of 20 feet at the top of the Pleistocene suggests that the Gentilly fault has had higher magnitude and possibly more frequent individual slip events than either the Vacherie or Ironton faults.  Possible ranges for values of a distribution of magnitude and/or frequency of fault slip events can be considered.  If a distinct fault slip event occurred an average of once a century over the 7,000 year span representing the accumulated Holocene sediments, then there would have been 70 individual slip events and the average magnitude of each event would have been about 3.4 inches of vertical displacement.  If distinct slip events occurred once a decade, then there would have been 700 events with an average magnitude of about 1/3 inch.  Given the size distribution patterns of other natural phenomenon like earthquakes, floods and hurricanes, it is more likely that fault slip events have not occurred at evenly spaced intervals of time with equal magnitudes.  It is more likely that the historical distribution the magnitude of slip events on the Gentilly fault exhibits a log-normal pattern, and that high-end magnitudes with less than a 1% chance of occurrence may have values of 2 feet or greater, while low-end magnitudes may be fractions of an inch.  The cumulative effect of all of the individual fault slip events since the beginning of Holocene deposition accounts for the 20 feet of vertical offset at the top of the Pleistocene.

The Central Wetlands Unit, which overlies the thick peat accumulation on the hanging wall side of the Gentilly fault is recognized as a hot spot of wetlands loss.  It is likely that subsidence associated with fault slip and differential compaction of the peat layer across the fault has contributed significantly to the submergence of the wetlands.  The magnitude of the rate of subsidence is evidenced by the submergence of Old Paris Road, which as constructed in the 1930s, and is now 2 feet below the water’s surface.  This is an obvious impact on infrastructure.  It also appears likely that the surface trace of the Gentilly fault can be extrapolated to connect with a fault segment documented by Dokka (2011).  He used a LIDAR digital elevation model to reveal the escarpment of a fault crossing Gentilly Blvd. in eastern New Orleans.  Dokka also documented shear and extension fractures in the streets that coincided with the location of the fault.

Picture15.png

The Barataria fault (Fig 15.) is a critical component of the linked tectonic system that includes the Ironton fault.  The fault plane connects the Bay de Chene and Lake Washington salt domes at the western and eastern ends of the fault plane respectively.  It is probable that the domes and the fault have interacted throughout the Cenozoic Era, and that halokenesis and salt dissolution have played a causal role in fault slip events.  Expansion index evaluation across the fault indicates continual episodic movement of the fault throughout the Pleistocene, but there is not adequate data to determine the vertical offset at the top of the Pleistocene. A boring profile from Frazier (1967) crosses the fault near its western edge.  Frazier interpreted the sedimentary layers associated with the Bayou des Families (#7), Bayou Blue (#10) and early Mississippi (#13) deltas.  The vertical trace of the fault has been superimposed on the profile to highlight the offset of the sedimentary layers across the fault.  It is likely that each of the Holocene deltas defined by Frazier interacted with the components of the linked tectonic system, and that the architecture of each delta lobe reflects its interaction with the structural system.  As will be discussed in a subsequent section, there is good evidence that sedimentary loading associated with the Miocene deltas that were active across this area caused accelerated slip on certain faults.  It is probable that the deposition of deltaic sediments in the Holocene also contributed to fault slip events, and fault slip probably played an important role in the submergence of the delta lobes.  The same depth interval associated with the Bayou des Families delta on the boring profile in Fig. 15 can be seen in the CPT profile in Fig. 9, and it appears that this interval has been offset by the Ironton fault. 

Picture16.png

The Montegut fault (Fig 16.) is along a trend of faults that extend westward from the Barataria fault that Gagliano et al (2003) called the Golden Meadow fault trend.  The fault is one of several faults that were mapped in the subsurface with 3-D seismic data and extrapolated to surface traces that coincide with lineations in the marsh (Culpepper et al, 2019).  Gagliano et al (2003) used a boring profile to document thicker peat accumulations on the hanging wall side of the fault, and related the fault to a lineation in the marsh surface and an area of dead cypress trees.  Keucher et al (2001) identified a distinct anomaly in total dissolved solids in the soils adjacent to the fault.  It is likely that saline fluid migration along the fault plane is associated with the soil salinity anomaly, and may have contributed to the death of the cypress trees. 

The Bastian Bay fault (Fig. 17) parallels the Empire fault (Fig. 18), and both extend from the eastern flank of the Lake Washington salt dome.  The fault plane has been mapped in the subsurface with well logs and 3-D seismic data, and it can be extrapolated to a surface trace that coincides with a distinct lineation in the marsh.  A profile of well log correlations across the fault shows increasing fault throw with depth throughout the Pleistocene, and a vertical displacement at the top of the Pleistocene can be reasonably estimated at about 40 feet.  Martin (2006) also documented expansion indices across the fault in the same depth range that indicate continual fault movement.  The magnitude of the vertical offset at the top of the Pleistocene also suggests that it may be the result of larger and/or more frequent fault slip events that those previously discussed here.  The nature of one of these events may be revealed in a comparison of the marsh surface across the area in images from ProPublica (2014).  In the period between 1973 and 1975 a large area across the hanging wall of the Bastian Bay fault submerged and converted from mostly saline marsh to open water.  Gagliano et al (2003) estimated that the vertical displacement on the fault at the surface was 3.0 to 3.5 feet, and reported anecdotal evidence from oyster fishermen and camp owners.

Picture17.png

This pattern of rapid subsidence across a large area would be consistent with a fault slip event, but it may be possible that the “event” lasted for several weeks or months.  In this type of scenario the entire magnitude of the fault slip would not be contained in a single catastrophic slip event, but it may have distributed over time as a series of smaller slip events or an extended period in which the rate of a continuous slow-slip motion was accelerated.  Much more investigation is needed to understand the potential for this type of fault movement in the delta plain environment, but it should be considered in an evaluation of the Ironton fault.

The Empire fault (fig 18) was named by Gagliano et al (2003).  The same fault has also recently been referred to as the Adams Bay fault, but the original name is used here to avoid confusion.  The Empire fault is parallel to and just north of the Bastian Bay fault.  The Empire fault plane has been mapped in the subsurface with well logs and 3-D seismic data, and it can be extrapolated to one of the most distinct surface escarpments on the delta plain.  Expansion indices across the fault indicate continual movement throughout the Pleistocene (Martin, 2006)  Gagliano et al (2003) constructed an auger boring profile across the fault to illustrate a vertical offset across the fault of 3.0 to 3.5 feet at the surface.  Martin (2006) also quoted anecdotal evidence from Greg Linscombe of the Louisiana Department of Wildlife and Fisheries that the marsh break associated with the fault escarpment formed between 1976 and 1977.  The magnitude and pattern of this fault slip event are similar to that associated with the Bastian Bay fault, but they appear to be separated in time by 3 or 4 years.  The relationship between these two apparent fault slip events may provide insights into the nature of the connection between elements in the linked tectonic system.  It may be possible that the Bastian Bay fault slip event changed the stress fields in the near-surface in a way that affected and contributed to the cause of the Empire fault slip event.  Much more work is needed to understand these relationships, but the possibility for other components of the linked tectonic system to affect the Ironton fault should be considered in the evaluation.

Picture18.png
Picture19.png

The Magnolia fault (Fig. 19) appears to be linked to the Ironton fault through a connection with the Lake Hermitage salt dome.  The Magnolia faults extends out of the eastern flank of the dome toward the Potash salt dome, but it not clear if it connects directly to the Potash dome from available data coverage.  The plane of the Magnolia fault has been mapped with well logs and 3-D seismic data (Bullock et al, 2018) and the fault plane has been extrapolated to a surface trace which coincides with a fairly distinct lineation in the marsh.  Bullock et al (2018) used a sediment core profile to illustrate the near surface offset of the fault and increased thickness of sedimentary layers on the hanging wall side of the fault.  The upper layer of accumulated peat is thicker on the hanging wall of the fault, and higher average subsidence rates over about the last 1,000 years have also been interpreted by the incorporation of carbon dating of the organic sediments. 

The surface trace of the Magnolia fault crosses the Lake Hermitage Marsh Creation Project.  A reasonable peat thickness map can be constructed from the geotech cores from this project and cores taken by Bullock et al (2018).  Peat thickness is greater on the hanging wall side of the fault, and it is likely that this reflects a cause and effect relationship with fault-induced subsidence.  This configuration also provides valuable insights into the effect of sediment loading across the fault, as emplaced fill from marsh creation project extends across portions of the footwall and hanging wall sides of the fault.  Simoneaux et al (2016) reported that “the magnitude of lateral displacement of these soft organic soils after fill material placement (aka ‘mudwaving’) was grossly underestimated.”  It is not clear if lateral displacement of organic soils was greater on the hanging wall side of the fault, but that would be a reasonable assumption given the evidence of their thickness.  This phenomenon should be further investigated, and some effort should be made to evaluate the potential for the lateral displacement of thick deposits of highly organic soils in the vicinity of the MBSD Project the as a result of mineral sediment deposited by the MBSD Project.  This will be discussed in more detail in a subsequent section.

 

Subregional geological structure

To meet the objective of providing insights into the historical relationships between sediment loading, fault slip events, and fault-induced subsidence that may be incorporated into predictive models for land elevation and land area gain in the MBSD project area, this evaluation will review the sub-regional geological structure in the vicinity of the project area. It will also examine the interactions between the structural system and four major Miocene delta systems.  Fig. 20 shows an expanded scope of the linked tectonic system around the MBSD project area with a map of fault plane contours in yellow and salt dome contours in green.  The contours are not labeled, but the contour interval is 1,000’ feet, and each of the faults extends to the surface with a 0’-depth contour.  Figs. 21-23 show the extent of the mapped horizons at the Mid-Miocene, Pliocene and Pleistocene respectively.  The Pleistocene map as derived from subsurface depth contours on the maps in the Coastal Plain section of the USACE Engineering Geology Mapping website.  Faults that appear to reach the surface were integrated with the contours, and the contours were adjusted to reflect reasonable amounts of throw at the depth of the Pleistocene.  Four sub-regional profiles are indicated on each map.  These profiles are shown in block diagram form in Figs. 24-27. 

Picture20.png
Picture21.png
Picture22.png
Picture23.png
Picture24.png

Profile 1 crosses the Lafitte and Bay de Chene salt domes at the western ends of the Ironton and Barataria faults.  The profile is just west of the surface trace of the Ironton fault, and the fault can be seen on the profile extending up from an intersection with the Lafitte salt dome, but not reaching the Pliocene mapped horizon.  The salt domes exhibit a pattern that reflects a history in which they were probably squeezed up from an original Jurassic layer of salt under the weight of late Cretaceous and Paleogene sediments.  Their tilted alignment may also suggest some association with down to the north faults that are related to “toe structures” that are antithetic to the Baton Rouge fault system.

The Lafitte graben is expressed above the top of the salt dome.  “Graben” is derived from the German word for grave, and it refers to the down-dropped area between sets of opposing faults.  Movement on these faults is directly associated with halokenetic movement on the salt dome.  The fact that these faults appear to extend to the surface suggests some degree of geologically recent salt movement.  Profile 1 very closely coincides with the boring profile constructed by Frazier (1967).  The inset in the upper left corner of Fig. 24 shows the boring profile with the interior faults of the Lafitte graben superimposed.  The black layer at the top of the graben indicates accumulated peat within the graben.  This would further suggest recent fault movement and related subsidence at the Lafitte graben.

The profile also shows the dramatically expanded upper Miocene section on the hanging wall of the Barataria fault.  This is reflected in the expansion index graph for the fault in Fig. 31.

Picture25.png

Profile 2 crosses near the center of the Ironton and Barataria Bay fault planes.  It clearly shows the listric nature of the faults.  With the ability to seismically image below the 20,000’ depth of the profile, it would be possible to see the deeper extents of these fault planes.  It is likely that the faults would continue to flatten and to eventually merge into a nearly horizontal glide plane or decollment surface, as suggested by the dashed lines.  This probable deep connection between the faults is another aspect of the linked tectonic system that is not presented in the maps and profiles of this evaluation.  Dokka et al (2006) referred to this configuration as the “southeast Louisiana allochthon”, an extensional complex driven by gravity instabilities.  They used GPS data to assign rates of vertical and lateral motion associated with movement on this tectonic system.  (Also see Karegar et al (2015), Fig. 33)

Picture26.png

Profile 3 crosses the Magnolia fault and the Lake Washington salt dome.  Lake Washington exhibits a pattern that suggests a strong component of lateral movement through the middle Miocene followed by a period of vertical diapiric movement.  The top of the dome is at a depth of about 2,000’ and it intersects the Pliocene mapped horizon.  There appears to be a genetic relationship between the Lake Washington salt dome and the Magnolia fault.  Expansion on the fault appears to be greatest in the middle Miocene.  The thickness of the upper Miocene increases into the center of structural syncline that is probably related to salt withdrawal in the early stages of vertical diapiric movement on the dome.  The extent of this syncline can be seen in map view on Fig. 21.  It is probable that the downward movement of the syncline due to salt withdrawal would have contributed to movement on the Magnoila fault throughout the upper Miocene and Pliocene.

Throughout the Cenozoic Era the faults moved in response to sedimentary loading.  Some of the smaller faults on the profile that were active in the upper and middle Miocene ceased their movement, and became covered with enough sedimentary overburden to prevent future movement.  The major faults continued to move episodically throughout the span of geological time, which has allowed them to propagate to the surface.  It is likely that, if properly imaged, these faults would exhibit evidence of movement in the recent geologic past and in the historical past over the last few centuries.  It is more probable that these faults are still capable of active movement in the present, than they that they lost that capability sometime in the recent past.

Picture27.png

Profile 4 crosses the Delacroix Island fault, the Potash salt dome, the Empire fault and the Bastian Bay fault.  Interval thickening across the fault is greatest in the middle Miocene at Delacroix Island and upper Miocene at Bastian Bay.  The Bastian Bay fault also significantly expands the Pliocene interval.  As will be discussed in the next section, the changes in expansion index for any of these faults can be related to their relationships with the active deltas that prograded across southeast Louisiana in the Miocene.  There is a very clear relationship between sediment loading during active deltaic deposition and the magnitude of fault movement.  Insights from these relationships could be used to formulate a predictive framework to assess the impacts of sediment loading at the MBSD Project.

 

Fault-Delta Interactions

Curtis (1970) provided detailed delineations of the Miocene delta systems of southeast Louisiana and a generalized conception about their relationships with faults that were active on the continental margin at the time.  The general pattern over the period of several million years was a progressive advancement of the continental shelf edge to the south.  Delta systems tended to cluster into phases, and each phase had a relationship with a different set of faults.  Fig. 28 from Curtis (1970) begins with the advancement of deltas of the Mid-Miocene transgressive phase.  The block diagram on the right shows the engagement of the delta systems with a theoretical fault system in the center of the block.  Mid-Miocene regressive phase deltas fully engaged with this fault system, and deltaic sediments are thicker on the hanging wall side of the fault.  Fault activity is implicitly caused by sedimentary loading from the delta, and fault-induced subsidence creates accommodation space that allows for the accumulation of more sediment.  Deltas in the Late Miocene regressive phase prograded beyond the area of influence of the fault system, and engaged with another fault system further to the south and more capable of accepting the sedimentary load.  Through this process each fault system has an expansion index that reflects its engagement with a delta system of a certain age.

Picture28.png

Six faults in the vicinity of the MBSD Project are evaluated here to consider the relationships between sedimentary loading associated with Miocene deltas and fault movement. Each of these faults labeled in Fig. 29 is a major fault that is considered to be a part of a linked tectonic system with other faults and salt domes. Each fault has been active since at least the middle Miocene, and all appear to extend to the surface and show some evidence of activity in the recent past. Fig. 29 represents a generalized map of fault traces and salt domes in the subsurface. It is not intended to represent any particular subsurface stratigraphic horizon, but it is a representation of an approximate location of the fault trace at a depth of about 10,000 feet. The map is a compilation of interpretations from numerous sources. It is primarily intended to provide more detail than fault maps currently available in the technical literature, such as

Picture29.png

Wallace (1962).  The exact location of faults in the subsurface should be determined from horizon structure maps, such as those in Figs. 21 and 22.  This generalized representation is accurate enough to be used in comparison to the paleogeographic reconstructions in Figs. 30-32.  Expansion index graphs for each fault were constructed in exactly the same manner as that for the Ironton fault in Fig. 4.  On each graph in the succeeding figures age in million years before present is on the vertical axis and expansion index is on the horizontal axis.

Picture30.png

The Delta Farms, Ironton and Delacroix Island faults align with the boundary between the continental shelf and the upper continental slope in middle Miocene between 12 and 13 million years before present.  The time interval is highlighted in yellow on the stratigraphic chart in the upper right and on each of the expansion index graphs.  These shelf edge faults are an integral part of the paleogeographic reconstruction.  The outlines of the active delta systems are derived from Curtis (1970) by relating the transgressive-regressive sequences to global sea level cycles.  Eustatic curves for global sea levels were not established at the time of the publication of Curtis (1970), but the patterns of the transgressive-regressive cycles used at the time can be correlated to the eustatic curves with reasonable accuracy.  A regional interval isopach (thickness) map for the same stratigraphic interval reveals two discrete depositional centers in which deltaic deposits of that age are measurably thicker than in surrounding areas.  A generalized cross section across any of these faults at the time of deposition would be similar to the center block diagram in Fig. 28.

The expansion index graph for each fault reveals maximum expansion on the fault during this time interval.  The faults are most active at the time of depositional loading, and fault-induced subsidence provides the accommodation capacity that keeps the delta engaged with the tectonic system during the span of deposition.  These relationships underscore the observation that faults are as much a part of a delta system as the channels and sand bars.  Fault movement caused by depositional loading should be an expected response in an active delta environment.

Picture31.png

Fig. 31 shows the paleogeographic reconstruction for the end of the middle Miocene between 10 and 12 million years before present.  The shelf edge and the active delta complexes have prograded to the south and the active depositional centers have shifted to another set of faults that include the Barataria and Magnolia faults.  The expansion index graphs on these faults indicate maximum fault movement during this time interval.  The graphs in Fig. 30 exhibit expansion indices greater than 1 during this time period This indicates that those faults also interacted with delta systems of this time, but with significantly lower magnitudes of fault activity.

Picture32.png

The paleogeographic reconstruction for the upper (or late) Miocene in Fig. 32 shows the major depositional centers positioned along the northern boundary of the Terrebonne Trough.  The Trough is defined by the alinement of down to the south faults along its northern boundary (including the Barataria and Bastian Bay and Montegut faults), and down to the north faults and salt domes along its southern boundary.  The accommodation capacity of this structural system kept the major delta systems engaged throughout the remainder of the Miocene and into the Pliocene.  The magnitude of accumulation of relatively thick and dense sedimentary layers within this basin created a regional density anomaly that is obvious on maps of the Earth’s gravitational field.  The density differential across the area has been a contributing factor to fault movement and subsidence since the end of the Miocene. 

 

Framework assessments for fault slip and subsidence rate models

The objectives of this evaluation are to provide a framework for predictive models for the magnitude and frequency of fault slip events and the relationships between fault slip, subsidence, elevation and land gain associated with the MBSD Project.  This type of framework can now be considered in the context of the characteristics other faults presented here.  The Ironton fault is a part of a linked tectonic system which allows for comparison with other faults in the system to help derive some preliminary values for the magnitude and frequency of fault slip, the rates of subsidence associated with the fault, and the probable response of loading sediment across the fault.  The principal metric that allows for comparison among the faults is the throw of the fault at the top of the Pleistocene horizon.  This value is estimated to be 6 feet on the Ironton fault based on the interpretation of the Conveyance Structure CPT profile presented here in Fig. 9.  The value of throw at this horizon on other faults in the area established here is 3 feet on the Vacherie fault (Fig. 11), about 20 feet on the Gentilly fault (Fig. 14), and about 40 feet on the Bastian Bay fault (Fig. 17).  Fisk (1944) documented a recent slip event on the Vacherie fault with 8 inches of vertical displacement, and Gagliano provided evidence for a recent slip event on the Bastian Bay fault of 3.0 to 3.5 feet of vertical displacement.

The intention of the recommendations made in McLindon et al (2017) was to collect the data necessary to evaluate the potential for fault slip in the vicinity of the MBSD Project.  In the absence of the data necessary to fully develop a probabilistic model for future fault slip events, the values provided in this evaluation can be used to make some framework estimates.  The expansion index graph on the Ironton fault in Fig. 5 gives some insight into the historical distribution of fault slip events.  If the full historical distribution could be known, it would be very likely to exhibit a log-normal pattern for the magnitude of the events.  This would mean that the distribution would consist of very many small events and very few large events.  This type of size distribution is commonly associated with magnitude of earthquakes, which are the result of fault slip events in tectonic regions.  The 6 feet of vertical offset at the top of Pleistocene on the Ironton fault is the cumulative offset of all slip events that have occurred since the beginning of Holocene deposition in the area.  Frazier (1967) dated the basal transgressive deposit on the top of the Pleistocene in St. Bernard Parish at about 7,000 years before present.  The orange sandy layer on the CPT profile in Fig. 9 is likely to represent a similar basal transgressive unit, and it is reasonable to assume that it has a similar age.  This would mean that 6 feet of vertical offset has occurred over the last 7,000 years.  If fault slip events occurred an average of once every century over that time span, then there would have been 70 discrete events with an average magnitude of about 1 inch of offset.  If the magnitude of the slip events were log-normally distributed, as should be expected, then the distribution would reflect many small events of perhaps a fraction of an inch of offset and a few large events with magnitudes of possibly 1 foot of offset or more.  This type of distribution would be consistent with the implications of the relationships between slip events on the Vacherie fault and crevasses of the river at the site of the fault (Fisk, 1944).  Given that 8 inches of vertical offset in 1943 was not adequate to crevasse the river, it is reasonable to assume that a larger magnitude event would have been required, which would support a high-end value of about 1 foot vertical displacement at Ironton. 

These framework estimates indicate that if a probabilistic distribution for future slip events were developed, it would include a non-zero value for the probability of the occurrence of a slip event on the Ironton fault during the time span of the MBSD project, and a smaller, but still non-zero value for the occurrence of an event with a magnitude of 1 foot of vertical offset.  The acquisition of high resolution seismic data and high resolution dating of sediment cores across the fault should contribute to a significant improvement in the development of probabilistic models, and the recommendations of McLindon et al (2017) are reiterated here. 

Estimates for subsidence rate models for the vicinity of the MBSD Project can be derived from observations about the relationships between subsidence rates, the thickness of the Holocene interval, and the organic content of the Holocene sediments in the areas examined here.  The depth to the top of the Pleistocene (or the thickness of the Holocene) can be estimated at each of the faults examined here.  The organic content of the soils has been estimated by the thickness of accumulated peat at the Gentilly fault, the Montegut fault and the Magnolia fault.  A fairly detailed map of the thickness of accumulated peat and organic clays has been made with data from geotech cores taken for the MBSD Project and the Bayou Dupont Marsh Creation Project (Fig. 35).  Byrnes et al (2019a & 2019b) found a “compelling relationship between subsidence and the age, composition and thickness of Holocene deltaic deposits” in the Breton and Barataria hydrologic basins, which includes the MBSD Project area.  An estimate of total Holocene thickness could be derived from the top of Pleistocene map (Fig. 23), however the Holocene “topstratum” contains the more highly organic and compactable sediments of the Holocene and the topstratum isopach map (Kulp, 2000) is more informative in evaluating subsidence.

Picture33.png

Fig. 33 combines Holocene topstratum thickness contours from Kulp (2000) with measurements of horizontal and vertical (subsidence) velocities from Karegar et al (2015) and surface fault traces where blue are down to the south faults and red are down to the north faults.  Land motion velocities are measured at stations in the Continuously Operating Reference System using GPS technology.  The relationship between the values of vertical velocities measured at these stations and the thickness of the Holocene is obvious.  The three stations with subsidence values near 7 mm/yr lie along the axis of the Terrebonne Trough which coincides with the axis of thick Holocene sediments.  Karegar et al (2015) also found that the horizontal vectors of land motion “may reflect slow downslope movement on a series of listric normal faults due to gravitational sliding”.  Some consideration should be given to the potential impacts of horizontal land motion due to movement on the Ironton fault on the MBSD Project. 

Picture34.png

A patchwork compilation of interpretations by Gould & Morgan (1962), Kosters (1989) and in this evaluation using geotech cores from the MBSD Project, the Bayou Dupont Marsh Creation Project, and the Lake Hermitage Marsh Creation Project in Fig. 34 provides insight into the probable thickness distribution of highly organic soils across the rest of the area.  There is a general tendency for thicker peats to accumulate off the natural levee flanks of distributary channels and on the hanging wall of faults that appear to extend to the surface.  Keucher et al (2001) indicated that the consolidation of highly organic and clay rich-facies in recent delta deposits is a primary component of subsidence.  Zou et al (2015) measured a subsidence rate anomaly from geodetic data that coincided with the thick peats between the natural levees of the Mississippi River and the Bayou Sauvage distributary channel and on the hanging wall of the Gentilly fault.  Maximum values of subsidence in this area were over 15 mm/yr.  Dixon et al (2006) documented similar values in the same area using InSAR technology.  A closer examination of peat thickness in the vicinity of the MBSD Project in Fig. 35 shows thickness contours constructed from measurements of thickness of the upper layer of organic clay and peat in geotech cores.  The pattern of the contours conforms to the general concept of thick peats accumulating off the natural levee flanks of distributary channels (taken from maps on the USACE Engineering Geology Mapping website) and on the hanging wall of a fault.

Picture35.png

Thick peats have accumulated between the distributary channels Cheniere Traverse Bayou to the north Bayou Dupont in the middle and Bayou des Families to the southwest and the hanging walls of the Ironton fault and in the Lafitte graben.  There are no direct measurements of subsidence in these areas.  Predictive modeling of elevation and land gain for the MBSD Project should allow for subsidence rate values as high as 15 mm/yr  above the areas of thick accumulations of organic clay and peat based on the analogy with the Gentilly fault area.  Some qualitative insights into the patterns of subsidence in this area can be gained from the patterns of wetlands loss shown in Fig. 36.  Reconstructions of what the marsh surface looked like in 1956 and 2009 from ProPublica (2014) indicate areas of wetlands loss due to subsidence in relationship to the thick peat accumulations.  The compelling correlation suggests that compaction of the upper layer of organic clays and peats is the primary mechanism of subsidence.  It is probable that the episodic slip events on the faults has contributed to historical subsidence that allowed for the thicker accumulation of peat on the hanging wall of the fault, but the contribution of fault slip to the current subsidence rate is unknown.  Predictive modeling for the MBSD Project should provide for an accelerated rate of compaction of the upper organic soils in response to the deposition of mineral sediment caused by the MBSD Project, which is currently in the models described in the Draft EIS (indicated by the red outline on Fig. 35).  The models should also allow for lateral displacement of the soft organic soils (mudwaving) in response to

Picture36.png

sediment accumulation based on analogous configurations of a fault and thick peats at the Lake Hermitage Marsh Creation Project (Fig. 19).

The profile in Fig. 36 is constructed from geotech cores from the MBSD Project and the Bayou Dupont Marsh Creation Project it shows the thickening of the organic clay and peat across the Ironton fault.  The thick accumulation of highly organic soils on the hanging wall is very similar to the pattern at the Gentilly fault.

Picture37.png

Morton et al (2009) concluded that the extraction of oil and gas from the Lafitte Field has been the primary cause of subsidence over the last few decades in the vicinity of the MBSD Project.  They provided no discussion of a causal process by which fluid extraction at the field would have induced subsidence, nor did they provide any evidence that the conditions necessary for such a process existed at the field.  Although never stated in Morton et al (2009), the implicit mechanism for extraction-induced subsidence at the field was compaction of the produced reservoirs caused by the change in volume due to the extraction of fluids.  Chan and Zobak (2007) documented that the magnitude of subsidence at the surface due to the extraction of fluids in the subsurface is a function of pressure change in the reservoir at depth.  The reservoirs that they studied in Lapeyrouse Field in Terrebonne Parish had all experienced significant reductions in reservoir pressure due to the extraction of natural gas.  Their calculations of reservoir compaction were based on the reduction in reservoir pressure.  Their models showed that if reservoir compaction did occur, the expression of subsidence at the surface would be in the form of a bowl-shaped depression immediately above the reservoir.

The reservoirs at Lafitte Field are almost entirely oil reservoirs that are portions of very large saline aquifers.  The oil is more buoyant than the saline pore fluid, and so oil accumulations tend to be in reservoirs around anticlinal structures, such as the Lafitte salt dome.  The size of any individual reservoir is very small relative to the size of the entire aquifer.  As oil is extracted from these reservoirs, the natural expansion of the aquifer is allowed by a slight compressibility of the pore fluids due to dissolved gases.  The saline fluids naturally flow into the pore spaces from which oil has been extracted in a process called “water drive”.  There is no significant change in pore pressure, nor is there any meaningful change in the volume of the aquifer.  The submerged marshes in this area are well-removed from the produced reservoirs at Lafitte field.  It is considered to be highly improbable that oil and gas extraction at Lafitte Field contributed to subsidence in the vicinity of the MBSD project.  There is compelling evidence that subsidence in this vicinity is due primarily to the compaction of the upper layer of organic clays and peats.

Reiteration of Recommenations

McLindon et al (2017) recommended:

  1. A review the subsurface geology using oil and gas industry 3-D seismic data. This may be accomplished through a collaborative engagement with owners, licensees and interpreters of the 3-D seismic surveys in the area. Such a collaborative engagement may be facilitated with the assistance of the New Orleans Geological Society, the Louisiana Mid-Continent Oil and Gas Association, or the Louisiana Oil and Gas Association.

  2. The acquisition of high resolution seismic data in the immediate vicinity of the diversion structure. This should necessarily include land-based acquisition along both banks of the river and marine acquisition in the river channel, as indicated in Figure 6.

  3. The acquisition of sediment core profiles across potential faults. The arrangement of these core profiles should be of adequate density to allow for the interpretation of faults by the vertical offset and variations in thickness of the sedimentary layers. The evaluation of core profiles should include detailed stratigraphic analysis and age-dating of the sedimentary layers to allow for estimates of historical subsidence rates and rates of fault movement.

  4. The addition of subsidence measurement capabilities similar to those of the Myrtle Grove Superstation at several additional locations in the vicinity of the diversion. These stations should be positioned with advance knowledge of the location of faults in the area to allow for the direct measurement of variations in subsidence velocities across the faults.

  5. The integration of detailed variations in subsidence rate and estimates of fault slip rate into predictive subsurface geological models including models for the response to sediment loading associated with diversion operations.

 

These recommendations are reiterated here and should be enhanced to insure that the acquisition of high resolution seismic data is based on survey design that provides for quality imaging to at least 200 feet and vertical resolution of at most a few feet.  Sediment cores should be subjected optically stimulated luminescence and carbon dating techniques.

It is further recommended here that the Draft EIS should consider:

1. The potential impacts of an episode of fault slip on the MBSD Project infrastructure based on a predictive model for the magnitude and frequency of future episodes

2. The potential for accelerated rates of subsidence due to sedimentary loading associated with the MBSD Project,

3. The potential for an induced fault slip event due to sedimentary loading associated with the MBSD Project.

REFERENCES

 

Akintomide, A. O. and Dawers, N. H., 2019, Spatial and Temporal Variation of Fault Activity in the Terrebonne Salt Withdrawal Basin, Southeastern Louisiana: Response to Salt Evacuation and Sediment Loading, presentation, AAPG Annual Meeting, San Antonio, Texas.

Armstrong, C., Mohrig, D., Hess, T., George, T., Straub, K.M.,  2014, Influence of growth faults on coastal fluvial systems: Examples from the late Miocene to Recent Mississippi River Delta, Sedimentary Geology, v. 301, p. 120-132

Bullock, J. S., Kulp, M. A., McLindon, C. D., 2018, Evaluation of the Magnolia growth fault, Plaquemines Parish, southeastern Louisiana, poster session, GSA Annual Meeting, Indianapolis, Indiana.

Byrnes, M. R., Britsch, L. D., Burlinghoff, J. L., Johnson, R., Khalil, S., 2019a, Recent Subsidence rates for Barataria Basin, Louisiana, Geo-Marine Letters, 14 p. doi.org/10.1007/s00367-019-00573-3.

Byrnes, M. R., Britsch, L. D., Burlinghoff, J. L., Johnson, R., Khalil, S., 2019b. Determining Recent Subsidence Rates for Breton Sound and Eastern Pontchartrain Basins, Louisiana: Implications for Engineering and Design of Coastal Restoration Projects. Final Report prepared for Louisiana Coastal Protection and Restoration Authority. Contract 4400009020, Task 8, 58 p.

Chan, A.W. & Zoback, M.D., 2007. The role of hydrocarbon production on land subsidence and fault reactivation in the Louisiana coastal zone. Journal of Coastal Research, 23(3), 771–786

Coastal Protection and Restoration Authority (CPRA), 2014a, Mid-Barataria Sediment Diversion - Final Draft Executive Summary Report 30% Basis of Design, State of Louisiana.

Culpepper, D.B., McDade, E. C., Dawers, N. H., Kulp, M. A., Zhang, R., 2019, Synthesis of Fault Traces in SE Louisiana Relative to Infrastructure,  TranSET Project No. 17GTLSU12

Curtis, D. M., 1970, Miocene Deltaic Sedimentation, Louisiana Gulf Coast, S.E.M.P. Special Publication 15, Deltaic Sedimentation, Modern and Ancient, p. 293-308

Dixon, T.H., Amelung, F., Ferretti, A., Novali, F., Rocca, F., Dokka, R., Sella, G., Sang-Wan, K., Wdowinski, S., Whitman, D., 2006, Subsidence and flooding in New Orleans, Nature, v. 441, p. 587-588

Dokka, R. K., G. Sella, and T. H. Dixon, 2006, Tectonic control of subsidence and southward displacement of southeast Louisiana with respect to stable North America, Geophys. Res. Lett., 33, L23308, 5 p.

Dokka, R.K., 2011, The role of deep processes in late 20th century subsidence of New Orleans and coastal areas of southern Louisiana and Mississippi, Journal of Geophysical Research, v. 16, 25 pgs.

Fisk, H.N., 1944, Geological Investigation of the Alluvial Valley of the Lower Mississippi River, Mississippi River Commission.

Frank, J. P., 2017, Evidence of fault movement during the Holocene in Southern Louisiana: integrating 3-D seismic data with shallow high resolution seismic data, MS Thesis, University of New Orleans, 91 p.

Frazier, D.E., 1967, Recent deltaic deposits of the Mississippi River: their development and chronology, Trans. G.C.A.G.S., v. 17, p. 287‐315

Gagliano, S.M., Kemp III, E. B., Wicker, K. M., Wiltenmuth, K. S., 2003. Active Geological Faults and Land Change in Southeastern Louisiana. Prepared for U.S. Army Corps of Engineers, New Orleans District, Contract No. DACW 29-00-C-0034.

Gould, H. R. and Morgan, J.P., 1962 Coastal Swamps and Marshes, Geology of the Gulf Coast and Central Texas, p 287-341

Hopkins, M., Lopez, J., Songy, A. 2018, Subsidence rates from faulting determined by real-time kinematic (RTK) elevation surveys of bridges in Lake Pontchartrain,  presentation, State of the Coast Conference 2018, New Orleans, Louisiana.

Johnston, A., Zhang, R., Gottardi, R., Dawers, N. H., 2017, Investigating the relationships between tectonics and land loss near Golden Meadow, Louisiana by utilizing 3-D seismic and well log data, poster session, GSA Annual Meeting, Seattle, Washington.

Karegar, M.A., Dixon, T.H., Malservisi, R., 2015, A three-dimensional surface velocity field for the Mississippi River Delta: Implications for coastal restoration and flood potential, Geology, v. 43, p. 519-522

Kosters, E. C., 1989, Organic-Clastic Facies Relationships and Chronostratigraphy of the Barataria Interlobe Basin, Mississippi Delta Plain, Journal of Sedimentary Petrology, v. 59, no. 1, p. 98-113

Kuecher, G.J., Roberts, H.H., Thompson, M.D. and Matthews, L., 2001, Evidence for Active Growth Faulting in the Terrebonne Delta Plain, South Louisiana: Implications for Wetland Loss and the Vertical Migration of Petroleum., Environmental Geosciences, v. 8, p. 77‐94

Kulp M. 2000. Holocene stratigraphy, history, and subsidence of the Mississippi River delta region, north-central Gulf of Mexico. PhD thesis. Univ. Kentucky, Lexington. 283 pp.

Levesh, J. L,, Kulp, M. A., McLindon, C. D., 2019, Fault-slip history of the Delacroix Island fault system and its effect on Holocene salt marshes of the Mississippi River delta plain,  presentation, GSA Annual Meeting, Charleston, South Carolina

Lopez, J., Penland, S., Williams, J., 1997, Confirmation of Active Geologic Faults in Lake Pontchartrain in Southeast Louisiana, GCAGS Transactions, v. 47, p. 299-303

Martin, E., 2006, Fault Induced Subsidence Near Empire and Bastian Bay, Louisiana, MS Thesis, Tulane University, 163 p.

McLindon, C.D., 2017, History of fault slip and interaction with deltaic deposition from the middle Miocene to the Present – Barataria fault, coastal Louisiana, poster session, American Geophysical Union, annual meeting, New Orleans, Louisiana.

McLindon, C.D., Dawers, N.A., Culpepper, D., Kulp, M.A., and McDade, E., 2017, Comments to the U.S. Army Corps of Engineers, New Orleans Division in reference to the The Environmental Impact Statement for the Mid-Barataria Sediment Diversion.  11 pg.

Morton, R.A., Bernier, J.C., and Kelso, K.W., 2009, Recent subsidence and erosion at diverse wetland sites in the southeastern Mississippi delta plain: U.S. Geological Survey Open-File Report, 2009–1158, 39 p., plus app. (p. 41-221).

ProPublica, 2014, Losing Ground

Simoneaux, R., Richard, B., Cantu, K., and Roy K., 2016, The Impacts of In-Situ Soil Displacements on Settlement and Consolidation of Marsh Creation Fill – Lake Hermitage Marsh Creation, Abstract – State of the Coast Conference, New Orleans, LA

Wallace, W.E., 1962, Fault Map of South Louisiana (Scale 1:500,000), Gulf Coast Association of Geological Societies Transactions Vol. XII, p. 194.

Zou, L., Kent, J., Lam, N. S.-N., Cai, H., Qjang, Y., and Li, K., 2016, Evaluating Land Subsidence Rates and Their Implications for Land Loss in the Lower Mississippi River Basin, Water, v. 8, 15 p.

Louisiana's Oil and Gas Industry - The Missing Link In Coastal Sustainability

This article was originally published in the Quarterly journal of the Society of Independent Professional Earth Scientists

The sustainability of coastal zones is becoming a global priority.  Nearly half of the world’s population lives within 100 miles of the ocean, and half of them live within 50 miles.  Population centers tend to be in cities close to the ocean that are built on coastal plains or river deltas.  South Louisiana is often recognized as a model for sustainability in coastal zones.  What is generally not recognized is the critical role that the oil and gas industry must play in that model.  The concept of sustainability has evolved over the past few decades.  Many of the early ideas of trying to restore coastal ecosystems to a pristine state that preceded the industrial revolution have given way to a managed sustainability that balances ecology with economy.  Real sustainability allows for the utilization of the natural resources of the coast while working to protect and preserve the essential ecological components.  Fisheries, transportation, ports and oil and gas production can all coexist with a productive natural environment, if effectively managed.  The truth is that the viability of the economic component of a coastal system in many ways determines the viability of the ecological component.  The worst areas of pollution and coastal degradation tend to be in countries with the poorest economies.

Louisiana’s plan for sustainability is laid out in the document “Integrated Ecosystem Restoration & Hurricane Protection: Louisiana’s Comprehensive Master Plan for a Sustainable Coast”.   The 2021 Fiscal Budget for the Master Plan was recently published by the Louisiana Coastal Protection and Restoration Authority (CPRA).  It is based on $3.34 billion in revenue in the successive three-year period 2021-23.  81% of that revenue is derived from the oil and gas industry.  A significant portion of this revenue comes from settlements from the Deepwater Horizon disaster, but almost $400 million comes from the Gulf of Mexico Energy Security Act (GOMESA), which provides revenue from offshore leasing, rentals and royalties.  Another $100 million comes from the Coastal Protection and Restoration Trust Fund, which collects revenues from leases, rentals, royalties and severance taxes in onshore Louisiana.  The economic viability of Louisiana’s plan for coastal sustainability is directly tied to the economic viability of the oil and gas industry along the coast of Louisiana.

The real value of the oil and gas industry’s contribution may lie in the science of sustainability.  Regardless of the magnitude of its expenditure, the success of the Master Plan will ultimately be determined by the degree to which it is based on an accurate assessment of the natural processes in the wetlands.  Earth science lies at the foundation of understanding the natural processes of the coast, and the oil and gas industry can provide a wealth of earth science data and knowledge base.  The application of industry data to studying coastal Louisiana began in the first years of this century.  Dr. Woody Gagliano started his career at the LSU Coastal Studies Institute before he founded his own company, Coastal Environments, Inc.  Dr. Gagliano is a recipient of the Coalition to Restore Coastal Louisiana’s Lifetime Achievement Award.  In a 2003 publication that accompanied his report to the US Army Corps of Engineers he stated “A century of oil and gas activity has made the subsurface of Southeastern Louisiana, which lies in the Gulf Coast Salt Dome Basin (Gulf Basin), one of the most understood geological provinces on earth.  Thousands of wells have been drilled; hundreds of thousands of seismic lines have been run; and the geological literature is voluminous.  The success of the oil and gas industry in this and other similar sedimentary basins throughout the world attest to the validity of the process-response models of sediment loading, compaction, salt movement and fault adjustments developed by geologists and geophysicists working in the Gulf Basin.   A working premise of the present study was to use this solid foundation of research as a basis for understanding and predicting land form and environmental change.”  Gagliano had recognized that subsurface geological processes played a fundamental role in controlling surface processes including subsidence and wetlands loss, and that oil and gas industry knowledge base was essential to understanding the geology.

Figure1.png

Gagliano’s study pre-dated the use of 3-D seismic for coastal research.  It was not until 2014 that the first 3-D survey was used in a university research project on coastal geology.  Geologists at Tulane University and the University of Texas at Austin used a 530-square mile 3-D survey in Plaquemines Parish which was donated by Western Geophysical.  Their research produced the publication “Influence of growth faults on coastal fluvial systems: Examples from the Late Miocene to the Recent Mississippi River Delta”, which documented 28 faults mapped with the 3-D survey.  The study found that “Most of the seismically imaged faults appear to extend up to the modern land surface and some affect the modern delta morphology. … several of these faults correspond to abrupt shifts from emergent wetlands to fully submerged areas of open water on the delta surface.”  Their study found that the submergence of wetlands along the downthrown sides of faults that reached the surface was major factor in causing wetlands loss.  Since then, seven 3-D surveys and one 2-D survey have been made available for university research through direct donations and internships at oil and gas companies.  These data have been used in four MS theses at the University of New Orleans, three MS theses at the University of Louisiana at Lafayette and one PhD dissertation at Tulane University (Figure 1).  Each of these research projects mapped faults that were projected to extend to the surface and to affect surface coastal processes.  Publications from this research will be entering the scientific literature over the next few years.  The first summary of this research “Synthesis of Fault Traces in Southeast Louisiana Relative to Infrastructure” was published in 2019 by the Transportation Consortium of South-Central States (TranSET).  The surface fault traces are now available in GIS format on the Louisiana Department of Transportation and Development website.  The TranSET study stated that “Highly detailed 3D seismic reflection data from the energy sector was donated for use by local universities to map geologic features identified on seismic data to help predict areas of future land loss. It is hypothesized that the locations of geologic features as deep-seated faults and salt domes will prove to be closely related to many areas of coastal land-loss in south Louisiana. As a result of this collaboration, coastal restoration, flood-control and sustainability initiatives that are based upon this knowledge should be optimized to deal with future relative sea-level rise.”

Ironically, no meaningful consideration is given to the impact of faults in modern coastal sustainability planning.  This is the missing link in the sustainability models.  The paradigm under which the concepts about the causes of change in Louisiana’s coastal wetlands were developed did not include geology.  Nobody has described this better than Dr. Len Bahr, who ran Louisiana’s coastal restoration program from 1993 to 2008.  He wrote is in his blog post of May 2013  “Since its very inception in 1989 Louisiana’s coastal restoration program has been dominated by coastal wetland ecologists like me, folks who deal in relatively short-term surface processes, not the long term geophysical and riverine processes that underlie the delta. In other words, the planning expertise has been dominated by those who deal primarily with surface processes on the visible veneer of the delta, not the riverine hydrodynamics and sedimentary processes that created the delta and the underlying tectonic processes and shallow and deep subsidence to which the delta ultimately responds.”  Louisiana’s coastal sustainability planning has been very slow to integrate a consideration of the subsurface geological processes to which the delta ultimately responds.

It is the nature of science that new data is constantly changing our collective understanding of natural processes.  This is true for coastal Louisiana.  It turns out that many of the initial conceptions about wetlands loss were incorrect, but many of the popular narratives that developed early in the paradigm still persist.  The average person in coastal Louisiana would be likely to tell you that the wetlands are in a state of ecological crisis.  They would note that this is evidenced by the fact that the coast is losing land at the rate of a football field an hour.   They would also note that the principal causes of loss are the lack of a sediment supply caused by the levees and erosion that is caused by “saltwater intrusion”.  None of these contentions are scientifically accurate.  Correcting the misconceptions about the coast is essential to the long-term success of sustainability.  New data is beginning to challenge the popular narratives of the old paradigm, but in order to be fully integrated into coastal sustainability planning new data must be understood in the context of subsurface geology.

Figure2.png

CPRA and the USGS monitor a wide array of data from the coastal wetlands.  These data include rates of subsidence, sediment accretion, surface elevation change and land area change.  Data are measured at nearly 400 stations in the Coastwide Reference Monitoring System (CRMS).  Figure 2 shows a map of the CRMS network and data being collected at some typical sites.  Most CRMS stations measure rates of sediment accretion that are greater than subsidence rates, and in many cases greater than subsidence plus sea level rise rates.  CPRA has reported that 70% of the stations in the delta region measure “positive elevation trajectories”, meaning that the marsh surface is gaining elevation relative to a fixed rod at the station.  Land loss rates were high in the mid-1970s when the popular narratives were formulated, but rates have been continually decreasing.  In fact, a 2019 study from the University of Arkansas found that “Remarkably, during the recent period … coastal Louisiana gained a modest amount of land”.  CPRA determined in a 2019 study that subsidence rates used in the formulation of the 2017 Master Plan were over 1.5 times higher than the rates actually being measured by modern technology.  This is the scientifically factual assessment of coastal Louisiana.  The coastal marshes actually have a fairly robust sediment supply.  Sediment accretion rates are high relative to subsidence, and much of the wetlands have been able to maintain surface elevation (Figure 3).  The result has been over a decade of stability in total net wetlands area.

Figure3.png

Geological interpretation using oil and gas industry data offers the best explanation for the revelations of the new body of data.  The university research has produced a fairly extensive coverage of surface fault trace maps.  These have shown conclusively that historical wetlands loss over the past few decades has been concentrated along the downthrown sides of the faults (Figure 4).  Historical data from several tide gauges along the coast have been used to document a subsidence event in the 1970s that coincided with the land loss event (Figure 5).  The most probable explanation for the coincidence of the timing of these events and the concentration of loss along the faults is that there was an episode of fault slip that spanned about two decades.  This is what Dr. Woody Gagliano believed.  LSU Geology professor Dr. Roy Dokka also believed that fault slip was a significant contributor to subsidence along the coast.  He referred to the gradual progressive slip over time as being like “decades-long earthquakes”.  This type of slow-slip movement on faults has been documented in tectonic regions.  It will require more detailed study to fully document this phenomenon in the gravity-driven faulting of the delta region.

Figure4.png
Figure5.png

The rapid submergence of large areas of coastal marsh over these decades would have resulted in a substantial addition of mineral sediment to the coastal hydrologic system.  While high rates of land loss were being measured at the marsh surface, sediment was being redistributed in the coastal bays and estuaries.  Several studies have shown that significant sediment accretion happens during tropical storms.  Tidal flux carries sediment from the bays and estuaries onto the marsh platforms.  The high rates of sediment accretion and surface elevation change measured by CRMS can be best explained as a redistribution of sediment from marsh material that was lost to subsidence during the fault slip event.

The integration of interpretations from oil and gas industry seismic data with data from CRMS points to a positive outlook for the long-term sustainability of the Louisiana coast.  It appears that the natural processes of accretion are working to maintain the elevation of resilient marsh platforms in the face of subsidence and sea level rise.  Many of the edges of these marsh platforms are coincident with the surface traces of faults.  The most successful sustainability plans will work to supplement the natural processes of sediment accretion – perhaps by depositing dredged material from the Mississippi River in the bays adjacent to the marsh.  The most effective of these plans can only be developed with a more complete understanding of the subsurface geology that is controlling the surface process.  This can only be achieved by finding a way to integrate oil and gas industry data and knowledge base into the planning process.  That is the missing link in coastal sustainability.

Wells to Monuments Project Proposal - Phase 1

Wells to Monuments Program

Phase 1 Proposal

Purpose and Justification

Picture19.png

Introduction

The underlying concept for the Wells to Monuments Program was conceived by Dr. Elizabeth McDade of Chinn-McDade Associates, LLC.  Dr. McDade is also a co-author of the TranSET study referenced here.  This proposal is made by the group of collaborators listed below.  The intent of the program is to repurpose orphaned oil and gas wells to be used as geodetic monuments equipped with Global Navigation Satellite System (GNSS) equipment in a manner consistent with protocols of the U.S. National Geodetic Survey.  The most prevalent GNSS is the Global Positioning System (GPS) is a U.S.-owned utility that provides users with positioning, navigation, and timing (PNT) services.  It may also be possible that these monuments could be attached to corner reflectors design to reflect radar signals from an interferometric synthetic aperture radar (InSAR) satellite.  Geodetic monuments provide a survey base for the state and support mapping, boundary determination, property delineation, infrastructure development, resource evaluation surveys and scientific applications.  In addition, monuments equipped with GPS and/or InSAR reflectors allow for the measurement of vertical and horizontal velocities of the Earth’s surface to which the monument is anchored.  Orphaned oil and gas wells generally consist of steel casing that is anchored many thousands of feet below the surface.  This means that monuments created from repurposed wells will measure velocities that are not impacted by near-surface compaction (e.g., Keogh and Törnqvist, 2019).  The ability to compare vertical and horizontal velocities between new deep-anchored, as well as existing shallow-anchored GNSS and tide gauge monuments would provide new insights that both increase the spatial resolution of coastal monitoring, but that also enable separation of shallow and deep contributions to coastal subsidence (vertical velocity) and all forms of infrastructure design and planning, which may be impacted by both vertical and horizontal velocities.  Spatial patterns in the new data may also enable quantification of slow creep along coastal fault systems.

Picture1.png

Collaborators

Ahmed Abdalla

Assistant Professor Research, LSU Center for GeoInformatics

Reda M. Amer

Adjunct Professor, Tulane University School of Science & Engineering

Asst. Professor, GIS Remote Sensing Director, Lamar University



Mark Byrnes

Principal Coastal Scientist, Applied Coastal Research and Engineering, Inc.



J. Anthony Cavell

Surveyor, LSU Center for GeoInformatics



Cynthia J. Ebinger

Marshall-Heape Chair Professor, Tulane University School of Science & Engineering



Jeffrey Freymueller

Endowed Chair for Geology of the Solid Earth,  Michigan State University Department of Earth and



Karen Luttrell

Associate Professor, LSU Department of Geology & Geophysics



Chris McLindon

Geologist, McLindon Geosciences, LLC



Cliff Mugnier

Chief of Geodesy, LSU Center for GeoInformatics



Randy L. Osborne

Network Manager, LSU Center for GeoInformatics



Diana Carolina Hurtado Pulido

Ph.D. Candidate, Tulane University, School of Science & Engineering

 

The Wells to Monuments Program Proposal consists of four phases:

1.       Purpose and Justification

2.       Engineering Design and Cost Estimation

3.       Bureaucratic Procedure

4.       Implementation

This document is Phase 1 of the program proposal.  It is intended to outline the scientific purpose and justification for the program, and it will be followed by separate documents that define the subsequent phases.  Increasing the density of geodetic monuments in Louisiana is considered to have implicit value beyond the scientific purposes discussed here.  Louisiana currently has relatively sparse coverage of monuments in the National Geodetic Survey CORS (Continually Operating Reference Station) program.  The State of North Carolina, by comparison, has made a substantial investment in its network of geodetic monuments, and has much more robust coverage.  The benefits of the increased density of coverage in Louisiana that would result from the implementation of this proposed program extend to all phases of infrastructure design and planning including transportation, navigation, ports and coastal restoration.

Picture20.png

In addition to providing measurements of subsidence and enhancing positioning applications, other useful information can be extracted from the GPS signals. For the wells that are located over open water, reflectometry techniques can be applied to estimate the height of the GPS antenna above the water surface from the GPS receiver’s recorded signal to noise ratio, effectively turning the GPS site into a tide gauge (e.g., Larson et al., 2013; Larson et al., 2017; Larson, 2019). In addition, all of the sites may be usable to estimate precipitable water vapor (PWV) from the atmospheric delays on the GPS signals (Bevis et al., 1992; Bevis et al., 1994; Moore et al., 2015; Shuanggen et al., 2015). PWV observations from GPS (e.g., https://www.suominet.ucar.edu) have been incorporated in improved weather forecasts, and have been used to study severe weather events and intense rainfall events (e.g., Basivi et al., 2015; Sapucci et al., 2018).

It is the intention of this proposal that the data generated by the monuments created through its implementation should be open access.  This offers the opportunity to build on current in-state data collection that is managed by the LSU Center for GeoInformatics.

Proposed Wells

Six wells were chosen from the Louisiana Department of Natural Resources (DNR) Orphaned Well Program for consideration in the proposal.  The selection of these wells was based on three factors:  1) a distribution of locations that would effectively increase the density of geodetic monuments in the coastal zone, 2) locations that are likely to provide a range of velocity signals based on subsurface geology, and 3) the utility of the well in the program based on the conditions of the well and associated infrastructure at the surface.  The wells are shown here referenced by their American Petroleum Institute (API) well serial number, which is a unique identifier for each well, and to which all records of the well are referenced.  The well locations are shown in context to the existing CORS network.

According to the DNR program the term “orphaned” refers to the current status of a particular well site and indicates that the operator of record is no longer a viable responsible party.  Following a specific notification procedure, well sites are declared officially orphaned after being sent to the State Register to be published.  Wells are removed from the Orphaned Well Program through the Louisiana Oilfield Site Restoration Program, which was created in 1993. The focus of the program is to properly plug and abandon orphaned wells, and to restore the site to approximate pre-well site conditions by cutting well casing below the surface and removing all surface facilities and infrastructure associated with the well.  Revenue for the Oilfield Site Restoration Program is generated from a fee on oil and gas production in the state which is paid quarterly by Louisiana oil and gas operators.  The fund has collected on average about $4.5 million each year.  The program has plugged and abandoned and restored 2306 wells since 1993 at a cost of $64 million.  There are currently 2833 wells remaining in the program.  The intention of this proposal is to use funds from the Oilfield Site Restoration Program to properly plug and abandon each well below the surface but to leave the well head and surface infrastructure in place to be repurposed as a deeply-anchored geodetic monument.

Subsidence

The primary objective of this proposal is to provide new infrastructure for the direct measurement of subsidence velocities in the Louisiana coastal zone.  The National Academies Consensus Study Report (2018) recognized subsidence as an essential research gap in understanding the long-term evolution of the coastal system “The causes, rates, and patterns of subsidence along the Gulf Coast are not sufficiently well understood to allow for accurate prediction at the local to regional scale.”  The reasons for this insufficiency are both a paucity of subsidence measurement infrastructure and a poor comprehension of the relative contribution of the various causal processes.  According to the most recent study of subsidence published by the Louisiana Coastal Protection and Restoration Authority (Byrnes et al, 2019) “Understanding the causes and rates of subsidence is critical to successful planning and implementation for Louisiana Coastal Master Plan Projects”.  The study recognized that subsidence measured at the earth’s surface is the result of natural causes such as the consolidation of Holocene, Pleistocene, and Tertiary age sediments, fault-induced elevation changes due to basin tectonics, and down-warping of the underlying lithosphere due to sediment loading, and human-induced causes such as lowering the groundwater table, overburden associated with flood protection levees, and marsh settling due to altered hydrology.  Other studies have recognized Glacial Isostatic Adjustment (GIA) as a contributing cause of subsidence (Love et al, 2016).  Some of these processes are rooted in geological processes that are happening deep within the earth, such as down-warping of the lithosphere, and some are rooted in very shallow processes, such as near-surface consolidation.  The study found that previous attributions of deep-seated geologic processes as the primary mechanisms of subsidence measured at the surface by Dokka (2006,2011) assumed that all measured elevation changes reflected movement deeper than the base of piling or rods.  Byrnes et al found that unanchored and/or unsleeved benchmark rods or pilings in Holocene sediment are affected by the downward force exerted by the consolidating deltaic sediment.

Keogh and Törnqvist (2019) highlighted the importance of understanding the depth at which benchmark monuments are anchored.  Tide gauges, which have been are the primary source of data used to calculate multi-decadal- to century-scale rates of relative sea-level change, are each associated with an elevation benchmark.   Keogh and Törnqvist showed that depth at which these monuments are anchored varies substantially across the coast relative to the top of the Pleistocene surface.  Holocene sediments above this surface are highly compactable.

Picture21.png

Using deeply anchored wells as monuments would provide a true direct measurement of subsidence velocities that are not affected by the consolidation of Holocene deltaic sediments.  This proposal recommends that each deeply-anchored monument created from a repurposed well should be co-located with a shallow-anchored monument.  Each of these monuments could be equipped with both a GPS device and an InSAR corner reflector.  This configuration would allow for the most accurate assessment of subsidence signals.

Picture22.png

Horizontal velocities

These figures from Karegar et al (2015) show values of vertical and horizontal velocities measured by CORS stations across Louisiana.  The inset in the graph on the left shows the relationship between vertical velocities and the thickness of Holocene sediments, as measured by Kulp (2000).  Byrnes et al (2019) similarly found that the spatial variability in subsidence has a “compelling relationship” with the thickness of Holocene deltaic deposits.  Karegar et al state that the horizontal velocities in the graph on the right “may reflect slow downslope movement on a series of listric normal faults due to gravitational sliding but could also represent the horizontal component of differential compaction”.  Velocity data collected from deeply-anchored monuments could provide valuable insights into horzontal motion in coastal Louisiana.  This vector of land motion is rarely considered in infrastructure design and planning.

Picture2.png

Subsurface Geology

Every significant recent study of vertical and horizontal land motion in coastal Louisiana (Byrnes et al, 2019; Dokka et al (2006); Jones et al (2016); Karegar et al (2015); Zou et al (2015)) has recognized the potential contribution of fault slip.  However, there has been limited access to accurate maps of subsurface fault planes or surface fault traces.  The first accurate map of surface fault traces in coastal Louisiana was published by Armstrong et al (2014).  Researchers  at  Tulane  and UT-Austin used an oil and gas industry 3-D  seismic survey in Plaquemines Parish to project the surface traces of about 20 faults.  Since then, on-going work by research groups at Tulane, UNO and ULL have significantly expanded the coverage of surface fault trace interpretation using seismic data.  This research has been funded in part by grants from the U.S. Department of Transportation’s University Transportation Centers Program and through the Louisiana CPRA Center of Excellence Research Grants Program (the RESTORE Act) administered by The Water Institute of the Gulf.  Preliminary interpretations have been published by the Transportation Consortium of South Central States (Culpepper et al, 2019), and are available in GIS form on the LA Department of Transportation Website.  These published interpretations have been integrated here with unpublished interpretations by McLindon Geosciences that are the result of subsurface geological interpretation in coastal Louisiana utilizing seismic and well log data over the past 40 years.  The surface fault trace maps presented here follow the color scheme used by Armstrong et al and are considered to be the most comprehensive and accurate interpretations available.

Picture3.png

These maps illustrate the strong relationship between vertical velocities measured by existing monuments and the thickness of Holocene deltaic sediments from Kulp (2000) that was noted by Karegar et al (2015) and Byrnes et al (2019).  It is important to recognize, however, that the axis of maximum thickness of the Holocene coincides with that of the regional geological feature commonly called the Terrebonne Trough.  Maximum Holocene thickness lies between sets of conjugate faults that form the northern and southern boundaries of the trough.  Progressive fault slip throughout the Holocene could have contributed to the increased thickness of the Holocene along this axis.  Recent fault movement could be contributing to the subsidence velocities that may otherwise appear to be solely dependent on Holocene thickness.  Deeply-anchored monuments could provide data that could help to differentiate the relative contributions of these two factors.  Horizontal velocities of the CORS sites reported by Karegar et al (2015) appear to be consistent with expected fault movement based on these interpretations.  If the horizontal velocities were entirely due to a horizontal component of compaction, the station at Grand Isle should measure a more northerly vector of motion, toward the center of the compacting basin.

The CORS network and the proposed wells are shown in context to stations in the Coastwide Reference Monitoring System (CRMS), which is operated by CPRA and USGS.  CRMS has collected a wealth of valuable data over the past two decades.  A typical station measures vertical sediment accretion, surface elevation change, land area change and soil salinity.  Currently, the highly accurate surface elevation change data does not record measurements relative to a fixed datum like sea level.  Elevation change measurements are made relative to a fixed Rod Surface Elevation Table (RSET), which is probably subsiding at an unknown rate.  An increased density of geodetic monuments among the CRMS stations would improve the estimations of subsidence rate at each RSET and would thereby allow for a more accurate estimation of the surface elevation change relative to sea level.  Eventually InSAR corner reflectors may be attached to the RSET at CRMS stations and calibrated to nearby geodetic monuments.  A direct measurement of subsidence rates at CRMS stations would allow for detailed evaluation of the relationships between accretion, subsidence and land area change.  Many CRMS stations that are measuring high sediment accretion rates are also measuring net gains in land area.  It is likely that the rate of accretion is greater than the rate of relative sea level rise (subsidence plus eustatic sea level rise) in these areas.

The geological interpretations presented here indicate that the horizontal vectors of land motion in coastal Louisiana are primarily affected by movement associated with the network of faults and salt features including salt domes and salt welds.  Hedec and Jackson (2011) modified the genetic block diagram from Rowan (1999), which shows the relationships among these features.  Schuster (1995) mapped outlines of the ridges of salt bodies that parallel the coast and the directions of horizontal movement along the salt welds (evacuated allochthonous salt).  The implied vectors of horizontal motion are consistent with those being measured by the existing geodetic monuments.  The increased density of monuments provided by the proposed program would allow for a more detailed understanding of these relationships.

Picture17a.png

Geological overview of proposed well locations

Three of the proposed well locations are shown in relationship to the surface fault traces published by Culpepper et al (2019).  Subsurface depth contours of the associated salt domes are shown in green.  These locations offer a range of settings relative to the subsurface geology that are likely to produce a range of land motion velocity signals.  The higher vertical velocities should be expected at wells 170578814100 to the south, which is near the center axis of the Terrebonne Trough and well 170572144800 which is in the Golden Meadow graben.  Well 170572297200 is north of the bounding faults of the trough and adjacent to the Clovelly salt dome.  This location should measure lower vertical velocity than those to the south.  It may also measure some component of velocity that is influenced by the dome, which rises to less than 400 feet below the surface.  Research in this area continues at Tulane with the addition of a 3-D seismic survey to the south.  Future publications resulting from this research should significantly expand the scope and detail of geological interpretations in this area.  The addition of land motion data from these proposed locations would offer a significant enhancement to understanding the results of the ongoing geological research.  The ultimate goal is to better understand relationships between subsurface geological processes and surface morphological and ecological processes.

Picture4.png

These three proposed locations are shown in relationship to surface fault traces and salt domes mapped by McLindon Geosciences.  A portion of these fault traces are included in the publication by Armstrong et al (2014).  The eventual publication of ongoing research at UNO will provide more detailed interpretation of faults in this area.  Maximum land motion velocities values would be expected from well 170750164600 to the south, which is anchored on the hanging wall block of the Bastian Bay fault.  The lowest values would be expected from well 1707523036200 to the north, which is on the footwall block of the Ironton fault.  Well 170512090000 in the center is in a complex location involving the Barataria fault, the Lake Five fault and the Bay de Chene salt dome and would be expected to produce intermediate velocity values.  Each of these locations will be considered in the context of more detailed geological interpretations in the pdf version of the full proposal document.

Picture5.png

Well 170750164600 is of particular interest because of its proximity to Pete Hebert’s camp. Gagliano et al (2003) collected anecdotal evidence of subsidence from local oyster fishermen and camp owners.  Pete Hebert that told Gagliano about his camp built in the 1960s in Plaquemines Parish west of Buras.  The camp was located on the banks of Bayou Ferrand, but by the mid-1970s it was standing in open water.  Hebert reported that the land surface under the camp sank by 3 to 3 ½ feet over that time span.  Gagliano later referred to the Bastian Bay fault as the “Rosetta Stone” for deciphering relationships among fault movement, subsidence and land loss in southeast Louisiana.  Subsequent scientific investigations of the Bastian Bay fault using seismic data were conducted by Armstrong et al (2014) and Dawers and Martin (2006).  The episode of fault movement implied by the rapid sinking of Pete Hebert’s camp would coincide with the formation of the Fort St. Phillip crevasse during the 1973 flood.  It may be the case that the crevasse was caused by the fault movement.  Historical crevasses at Darrow, Vacherie and Davis Pond appear to be associated with faults.  It may also be the case that movement on the fault during a major flood event was not coincidental.  The dilation of shallow aquifers during high river stages is likely to change the stress field across the surface and may contribute to triggering fault slip.  Episodic fault movement would be likely to produce variable velocity signals from a monument in this location.

Picture6.png
Slide1.PNG
Slide2.PNG

Summary and Conclusions

The Wells to Monuments Program offers an opportunity to significantly increase the density of geodetic monuments in the Louisiana coastal zone by making use of existing infrastructure that would otherwise incur costs to be removed.  The array of proposed well locations offers several sites that appear to be well suited for the co-location of deep- and shallow-anchored monuments that could make use of both GPS and InSAR technologies.  Collecting velocity data from locations that were specifically chosen for the geologic setting is the most expedient way to determine the relationships between subsurface geological features and land motion at the surface – and ultimately between subsurface geology and surface coastal processes.

 

References

Akintomide, A. O. and Dawers, N. H., 2019, Spatial and Temporal Variation of Fault Activity in the Terrebonne Salt Withdrawal Basin, Southeastern Louisiana: Response to Salt Evacuation and Sediment Loading, presentation, AAPG Annual Meeting, San Antonio, Texas.

Armstrong, C., Mohrig, D., Hess, T., George, T., Straub, K.M.,  2014, Influence of growth faults on coastal fluvial systems: Examples from the late Miocene to Recent Mississippi River Delta, Sedimentary Geology, v. 301, p. 120-132

Basivi, R., Fabry, F., Braun, J. J., & Vanhove, T. M. (2015). Precipitable water from GPS over the continental United States: Diurnal cycle, intercomparisons with NARR, and link with convective initiation. Journal Of Climate, 28, 2584-2599. doi:10.1175/JCLI-D-14-00366.1

Bevis, M., S. Businger, T. A. Herring, C. Rocken, R. A. Anthes, and R. H. Ware, 1992: GPS meteorology: Remote sensing of atmospheric water vapor using the Global Positioning System. J. Geophys. Res., 97, 15 787–15 801, doi:10.1029/92JD01517.

Bevis, M., S. Businger , S. Chiswell, T. A. Herring, R. Anthes, C. Rocken, and R. H. Ware, 1994: GPS meteorology: Mapping zenith wet delays onto precipitable water. J. Appl. Meteor., 33, 379–386, doi:10.1175/1520-0450(1994)033,0379:GMMZWD.2.0.CO;2.

Bullock, J. S., Kulp, M. A., McLindon, C. D., 2018, Evaluation of the Magnolia growth fault, Plaquemines Parish, southeastern Louisiana, poster session, GSA Annual Meeting, Indianapolis, Indiana.

Byrnes, M. R., Britsch, L. D., Burlinghoff, J. L., Johnson, R., Khalil, S., 2019, Recent Subsidence rates for Barataria Basin, Louisiana, Geo-Marine Letters, 14 p. doi.org/10.1007/s00367-019-00573-3.

Byrnes, M. R., Britsch, L. D., Burlinghoff, J. L., Johnson, R., Khalil, S., 2019. Determining Recent Subsidence Rates for Breton Sound and Eastern Pontchartrain Basins, Louisiana: Implications for Engineering and Design of Coastal Restoration Projects. Final Report prepared for Louisiana Coastal Protection and Restoration Authority. Contract 4400009020, Task 8, 58 p.

Culpepper, D., E.C. McDade, N. Dawers, M. Kulp, R. Zhang, 2019, Synthesis of Fault Traces in SE Louisiana Relative to Infrastructure, Transportation Consortium of South-Central States, Project No. 17GTLSU 12, 54 p.

Dixon, T.H., Amelung, F., Ferretti, A., Novali, F., Rocca, F., Dokka, R., Sella, G., Sang-Wan, K., Wdowinski, S., Whitman, D., 2006, Subsidence and flooding in New Orleans, Nature, v. 441, p. 587-588

Dawers, N.H. and E. Martin. 2005. Fault-related changes in Louisiana coastal geometry, Louisiana Governor’s Applied Coastal Research and Development Program, GACRDP Technical Report Series 05-000, 21 p.

Dokka, R.K., 2006, Modern-day tectonic subsidence in coastal Louisiana, Geology, v. 34, p. 281-284.

Dokka, R. K., G. Sella, and T. H. Dixon, 2006, Tectonic control of subsidence and southward displacement of southeast Louisiana with respect to stable North America, Geophys. Res. Lett., 33, L23308, 5 p.

Dokka, R.K., 2011, The role of deep processes in late 20th century subsidence of New Orleans and coastal areas of southern Louisiana and Mississippi, Journal of Geophysical Research, v. 116, 25 p.

Frank, J. P., 2017, Evidence of fault movement during the Holocene in Southern Louisiana: integrating 3-D seismic data with shallow high resolution seismic data, MS Thesis, University of New Orleans, 91 p.

Gagliano, S.M., Kemp III, E. B., Wicker, K. M., Wiltenmuth, K. S., Sabate, R. W., 2003 Neo-Tectonic Framework of Southeast Louisiana and Applications to Coastal Restoration, Transactions G.C.A.G.S., v. 53, p. 262-272

Gagliano, S.M., Kemp III, E. B., Wicker, K. M., Wiltenmuth, K. S., 2003. Active Geological Faults and Land Change in Southeastern Louisiana. Prepared for U.S. Army Corps of Engineers, New Orleans District, Contract No. DACW 29-00-C-0034.

Gagliano, S.M., et.al., 2005, Effects of Earthquakes, Fault Movements, and Subsidence on the South Louisiana Landscape, The Louisiana Civil Engineer Journal of the Louisiana Section of The American Society of Civil Engineers Baton Rouge, Louisiana Volume 13, Number 2, pp. 5-7, 19-22

Hopkins, M., Lopez, J., Songy, A. 2018, Subsidence rates from faulting determined by real-time kinematic (RTK) elevation surveys of bridges in Lake Pontchartrain,  presentation, State of the Coast Conference 2018, New Orleans, Louisiana.

Hudec, M. R. and M.P.A. Jackson, 2011, The Salt Mine, AAPG Memoir 95, 305 p.

Ivins, E.R., Dokka, R.K., Blom, R.G., 2007, Post-glacial sediment load and subsidence in coastal Louisiana, Geophysical Research Letters, v. 34, L16303, 5 p.

Johnston, A., Zhang, R., Gottardi, R., Dawers, N. H., 2017, Investigating the relationships between tectonics and land loss near Golden Meadow, Louisiana by utilizing 3-D seismic and well log data, poster session, GSA Annual Meeting, Seattle, Washington.

Jones, C.E., An, K., Blom, R.G., Kent, J.D., Ivins, E.R., and Bekaert, D., 2016, Anthropogenic and geologic influences on subsidence in the vicinity of New Orleans, Louisiana, Journal of Geophysical Research: Solid Earth, v. 121, DOI: 10.1002/2015JB012636.

Karegar, M.A., Dixon, T.H., Malservisi, R., 2015, A three-dimensional surface velocity field for the Mississippi River Delta: Implications for coastal restoration and flood potential, Geology, v. 43, p. 519-522

Keogh, M. and T. Törnqvist, 2019, Measuring rates of present-day relative sea-level rise in low-elevation coastal zones:  a critical evaluation, Ocean Science, 15, 61-73.

Kulp M. 2000. Holocene stratigraphy, history, and subsidence of the Mississippi River delta region, north-central Gulf of Mexico. PhD thesis. Univ. Kentucky, Lexington. 283 pp.

Larson, K. M., R. D. Ray, F. G. Nievinski, and J. T. Freymueller, The Accidental Tide Gauge: A GPS Reflections Case Study from Kachemak Bay, Alaska, IEEE GRSL, Vol 10(5), 1200-1204, doi:10.1109/LGRS.2012.2236075, 2013.

Larson, K.M., R. Ray, and S.P. Williams, A ten-year comparison of water levels measured with a geodetic GPS receiver versus a conventional tide gauge, J. Atmos. Ocean Tech., Vol. 34(2), 295-307, doi:10.1175/JTECH-D-16-0101.1, 2017.

Larson, K.M. Unanticipated Uses of the Global Positioning System, Annual Rev. Earth and Planet. Sci., Vol. 47, 19-40, doi:10.1146/annurev-earth-053018-060203, 2019.

Levesh, J. L,, Kulp, M. A., McLindon, C. D., 2019, Fault-slip history of the Delacroix Island fault system and its effect on Holocene salt marshes of the Mississippi River delta plain,  presentation, GSA Annual Meeting, Charleston, South Carolina

Love, R., Milne, G.A., Tarasov, L., Engelhart, S.E., Hijma, M.P., Latychev, K., Horton, B.P. and Törnqvist, T.E., 2016. The contribution of glacial isostatic adjustment to projections of sea‐level change along the Atlantic and Gulf coasts of North America. Earth's Future, 4(10), pp.440-464.

McLindon, C.D., 2017, History of fault slip and interaction with deltaic deposition from the middle Miocene to the Present – Barataria fault, coastal Louisiana, poster session, American Geophysical Union, annual meeting, New Orleans, Louisiana.

Meckel, T. A. 2008, An attempt to reconcile subsidence rates determined from various techniques in southern Louisiana, Quaternary Science Review, v. 27, p. 1517–1522

Moore, A., Small, I., Gutman, S., Bock, Y., Dumas, J., Fang, P., Haase, J., Jackson, M.and Laber, J. (2015) National Weather Service Forecasters use GPS precipitable water vapor for enhanced situational awareness during the Southern California Summer Monsoon. Bulletin of the American Meteorological Society, 96, 1867–1877.

National Academies of Sciences, Engineering, and Medicine. 2018, Understanding the Long-term Evolution of the Coupled Natural-Human Coastal System. The Future of the U.S. Gulf Coast: Washington, DC: The National Academies Press. doi: 10.17226/25108

Rowan, M.G., M.P.A. Jackson and B.D. Trudgill, 1999, Salt-Related Fault Families and Fault Welds in the Northern Gulf of Mexico, AAPG Bulletin v. 83 no. 9, p. 1454-84

Schuster, D. C., 1995, Deformation of Allochthonous Salt and Evolution of Salt-Structural Systems, Eastern Louisiana Gulf Coast, in M.P.A. Jackson, D.G. Roberts, and S. Snelson, eds., Salt Tectonics: a global perspective AAPG Memoir 65, p. 177-198

Shinkle, K. D., and R. K. Dokka, 2004, Rates of vertical displacement at benchmarks in the Lower Mississippi Valley and the northern Gulf Coast, NOAA Tech. Rep. 50

Sapucci LF, Machado LAT, de Souza EM, Campos TB. Global Positioning System precipitable water vapour (GPS-PWV) jumps before intense rain events: A potential application to nowcasting. Meteorol Appl. 2019;26:49–63. https:// doi.org/10.1002/met.1735

Shangguan, M., Heise, S., Bender, M., Dick, G., Ramatschi, M. and Wickert, J. (2015) Validation of GPS atmospheric water vapor with WVR data in satellite tracking mode. Annales de Geophysique, 33, 55–61.

Zou, L., Kent, J., Lam, N. S.-N., Cai, H., Qjang, Y., and Li, K., 2016, Evaluating Land Subsidence Rates and Their Implications for Land Loss in the Lower Mississippi River B

 

The Bastian Bay fault

in
AdamsBayCamp2.jpg

The best inputs to a scientific investigation always come from the people in the field.  Woody Gagliano spent many years in the marshes of coastal Louisiana studying the patterns of change.  He collected scientific data, but also spoke to the local oyster fishermen and camp owners.  It was Pete Hebert that told Gagliano about his camp built in the 1960s in Plaquemines Parish west of Buras.  The camp was located on the banks of Bayou Ferrand, but by the mid-1970s it was standing in open water.  Hebert reported that the land surface under the camp sank by 3 to 3 ½ feet over that time span.  Gagliano realized that the submergence of Hebert’s camp was due to movement on the Bastian Bay fault, which he had found by constructing profiles of the shallow subsurface from auger borings.  Gagliano later referred to the Bastian Bay fault as the “Rosetta Stone” for deciphering relationships among fault movement, subsidence and land loss in southeast Louisiana.  Subsequent scientific investigation of the Bastian Bay fault by geologists at Tulane and the University of Texas at Austin confirmed the location of the fault using seismic data (Armstrong et al, 2013; Dawers and Martin, 2006)

Camp1.PNG
Camp2.PNG

The Bastian Bay fault is near the eastern end of the Terrebonne Trough.  The fault is one of a series of down-to-the-south faults (highlighted in blue) that bound the northern rim of the trough.  They appear to have a conjugate relationship with a set of northern-dipping faults (highlighted in red) and salt domes that bound the southern rim.  Miocene, Pliocene and Quaternary sedimentary layers are all thicker within the trough, and the weight differential created by these thick piles of sediment appears to have driven subsidence and fault movement for at least the past 15 million years. 

Regional map.png

Most of the major faults bounding the Terrebonne Trough appear to extend to the surface where recent fault movement has played a critical role in causing wetlands loss.  The figure below shows the subsurface depth contours of the Bastian Bay fault plane in yellow and their connection to the Lake Washington salt dome shown in green.  The zero-depth contour coincides with the surface trace of the Bastian Bay fault.  Similar fault plane contour maps for each of the other faults that reach the surface could be constructed using seismic data and well logs. 

Contours2.PNG

Like most of the other faults along the Terrebonne Trough, the Bastian Bay fault controlled the formation of geological structures at depth.  Most of these structures control the accumulation of oil and gas in the Miocene, Pliocene and Pleistocene sands.  Bastian Bay Field is an example of a “rollover fault structure” in which an anticlinal structure is created by the thickening of sedimentary layers into the fault at the time of deposition.   At Bastian Bay Field the producing gas reservoirs are found in sand layers that were deposited by ancestral deltas of the Mississippi River during the upper Miocene.  The subsurface structure map on the “O” Sand published by the New Orleans Geological Society shows the concentric contours of the anticlinal structure, which has a crest at depth of 12,400 feet below the surface.  The red area indicates the extent of gas accumulation at the top of the rollover structure.  The Bastian Bay fault runs along the northern edge of the anticline and forms a boundary of the gas reservoir.  A small synthetic fault running parallel to the Bastian Bay fault forms the southern boundary of this gas reservoir.  Wells drilled to a depth of 12,400 feet or greater within the red area encountered gas in the “O” Sand.

OSandMap.png

A cross section of well logs from the field shows the arrangement of the fault and the producing gas reservoirs.  The “throw” of the fault at any horizon is the difference in elevation across the fault.  At the upper blue Miocene horizon the throw is about 1,500 feet.  At the “X” Sand the throw is over 3,000 feet.  The interval thickness between each of the labeled sand layers is greater on the downthrown (or hanging wall) side of the fault than it is on the upthrown (or footwall) side.  This indicates continuous fault movement during the deposition of the deltaic sands, and it is why faults like this are referred to as “growth faults”.  

CrossSection.png

A simplified model for the history of the Bastian Bay fault can be constructed by isolating the fault in a depositional sequence in which all other faults and salt domes have been removed and sedimentary layers are considered to be flat and continuous away from the fault.  The reality is that the complex arrangement of structural features in the subsurface have all been interacting with each other throughout the deposition of the sedimentary layers, but this is very difficult to model.

BlockDiagram1.png

The Bastian Bay fault probably began as slide structure on the lower continental slope during the lower Miocene.   A study of the mechanics of faults by geologists W. Crans, G. Mandl and J. Harembore at Shell Research Laboratory found that when loose sediment is deposited on a gently inclined slope (about 3 degrees), its weight has a surface-parallel component that tends to pull it down the slope.  As long as this force is balanced by the reactive shear stress along the slope, the sediment will stay on the slope in stable equilibrium.  As sedimentation continues and the height of the sediment layer increases, both the “pulling” weight component and the reactive shear stress increase proportionally.  The equilibrium will remain stable as long as the slope-parallel shear stress remains below the limit determined by the shear strength of the sediment.

The shear strength of the sediment may be reduced when pore fluids become “geopressured” meaning that they support a portion of the overburden.  The greater the magnitude of geopressure, the lower the shear strength of the sediments.  At a critical point when the driving weight component exceeds the friction along the base of a sedimentary layer it may begin to slide along the basal slip plane.  In this model the initial lateral slip of the fault results in a negative vertical component of slip in the active plastic region at the head of the fault and a positive vertical component in the passive plastic region at the toe of the fault.  The accommodation capacity created at the head of the fault allows for a differential accumulation of sediment on the hanging wall of the fault, which will continue to drive fault movement in a feedback loop.

CransDiagram.png
BlockDiagram2.png

The influx of terrigenous sediments into the Terrebonne Trough during the middle Miocene triggered episodes of dramatic fault movement on all of the bounding faults including the Bastian Bay fault.  It is significant to note that these faults tended to act as an interconnected system of faults throughout this time.  It is probable that a triggering event on one fault could result movement that rippled throughout the system.  The simplified model on the Bastian Bay fault during this period shows that the original lateral slip plane of the fault begins to arc upward to a steeper angle.  This is the result of more competent normally pressured sediments. The curving arc of the fault plane defines the nature of a listric fault.  The model also shows that fluids migrate along the fault plane, including hydrocarbons.  It is almost certain that the gas accumulated in the Miocene reservoirs at Bastian Bay Field migrated up the Bastian Bay fault and into the sand layers at some time after their deposition. 

BlockDiagram3.png

Movement on the Bastian Bay fault continued to propagate throughout the Pliocene and the Quaternary.  The movement of fluids up the fault plane continued throughout this time period and into the present.  It is probable that saline fluids reach the surface along these major faults.  This may be responsible for some of the salinity anomalies that have been measured in the coastal wetlands (Keucher et al, 2001).

BlockDiagram4.png

A subsurface structure map on the top of the Pliocene stratigraphic interval shows the trace of the Bastian Bay fault.  The fault trace is constructed on this map by marking the intersection of the depth of the mapped horizon and the depth of the fault plane.  In the vicinity of the Bastian Bay Field the top of the Pliocene horizon is at a depth of about 2350 feet on the hanging wall and about 2250 feet on the footwall.  The two sides of the fault trace on the map are coincident with the equivalent 2350-foot and 2250-foot contours on the fault plane map.  This means that the throw of the fault in this area at the top of the Pliocene is about 100 feet.  Because the elevation of the surface in south Louisiana is effectively zero, a map on the depth of the top of the Pliocene is also a map of the thickness of the Quaternary (the Pleistocene and the Holocene).  The Quaternary is therefore 100 feet thicker on the downthrown side of the Bastian Bay fault than it is on the upthrown sides.  This indicates that the fault continued to be active in the Quaternary, which generally qualifies it as an “active fault”. 

Contours1.PNG

A detailed well log cross section of the Quaternary across the Bastian Bay fault shows that the fault remained continually active throughout the interval.  Individual stratigraphic intervals in the Pleistocene are vertically offset by the fault and they are generally thicker on the hanging wall side than they are on the footwall side.  This pattern extends into a tentative interpretation of the base of the Holocene interval from well log correlation.  It is apparent from this interpretation that there is stratigraphic evidence that the Bastian Bay fault was active during the Holocene. It is probable that late Holocene activity on the fault would have been associated with the deposition of sediments in the Bayou Robinson delta lobe between 950 and 300 years ago (Kulp et al, 2005).  Bayou Ferrand was originally formed as a distributary channel in this delta lobe.

QuaternaryProfile.png

Gagliano realized that the submergence of the marshes around Pete Hebert’s camp was a result of the latest episode of movement on the Bastian Bay fault.  It is the nature of all faults to move episodically in (geologically) short bursts of activity.  In the tectonic regions of California these episodes are generally very brief, and fault slip happens dramatically in a matter of seconds.  In a delta system where the forces are all driven by the weight of the deltaic sediments fault slip can happen much more slowly and will generally not create a detectable seismic signal.  The most recent episode of fault slip on the Bastian Bay fault appears to have happened over a period of about a decade between the early 1970s and the early 1980s.  Imagery of the area taken from the ProPublica publication “Losing Ground” is used in this short video sequence to show the submergence of the marsh on the hanging wall side of the Bastian Bay fault during that time period.  The interconnected nature of the system of faults bounding the Terrebonne Trough throughout their geological history strongly suggests that similar patterns of movement happened on most of these faults during the same time period.  This is likely to be the principal cause of the episode of wetlands loss that occurred between the 1970s and 1980s in southeastern Louisiana.

REFERENCES

Armstrong, C., Mohrig, D., Hess, T., George, T., Straub, K.M., 2014, Influence of growth faults on coastal fluvial systems: Examples from the late Miocene to Recent Mississippi River Delta, Sedimentary Geology, v. 301, p. 120-132

Crans, W., Mandl, G., Haremboure, J., 1980, On the Theory of Growth Faulting: A Geomechanical Delta Model Based on Gravity Sliding, Journal of Petroleum Geology, v. 2, p. 265-307

Dawers, N. H. and Martin, E., 2006, Fault Related Changes in Louisiana Coastal Geometry, Louisiana Governor’s Applied Coastal Research and Development Program, GACRDP Technical Report Series 05-000, 21 p.

Gagliano, S.M., Kemp III, E. B., Wicker, K. M., Wiltenmuth, K. S., Sabate, R. W., 2003 Neo-Tectonic Framework of Southeast Louisiana and Applications to Coastal Restoration, Transactions G.C.A.G.S., v. 53, p. 262-272

Gagliano, S.M., Kemp III, E. B., Wicker, K. M., Wiltenmuth, K. S., 2003. Active Geological Faults and Land Change in Southeastern Louisiana. Prepared for U.S. Army Corps of Engineers, New Orleans District, Contract No. DACW 29-00-C-0034

Kulp, M. A., Fitzgerald, D., Penland, S., 2005, Sand-Rich Lithosomes of the Holocene Mississippi River Delta Plain, SEPM Special Publication No. 83, SEPM (Society for Sedimentary Geology), ISBN 1-56576-113-8, p. 279–293.

Land area gains in the Louisiana wetlands

Several stations in the Coastwide Reference Monitoring System (CRMS) have measured net gains in land area in a 1-square kilometer area around the station over the past two decades. These same stations have also measured significant rates of sediment accretion and positive surface elevation change as well.

Four example stations are represented here. The location of each is shown with respect to the USGS Land Area Change Map (Scientific Investigations Map 3381, Couvillion et al, 2017). Each of these stations is in a belt the runs just north of a belt of land loss hot spots discussed in the last blog. It appears likely that the land area gains measured at these stations is due in part to a redistribution of sediment that was lost to submergence in the belt to the south.

0395_slide.png

JuxtaposeJS Embed

 
0400_slide.png

JuxtaposeJS Embed

 
0287_slide.png

JuxtaposeJS Embed

 
4355_slide.png

JuxtaposeJS Embed

The values of sediment accretion rates and surface elevation change rates recorded at these four CRMS stations are higher than both the rates of subsidence and sea level rise for coastal Louisiana. A majority of the nearly 400 CRMS stations across the coast have measured rates of accretion and surface elevation change greater than either subsidence or sea level rise. New wetlands area appear to be created in many of these naturally accreting areas. In other areas accretion is helping to mitigate the rate of wetlands loss.

CRMS was establish under the Coastal Wetlands Planning, Protection and Restoration Act (CWPPRA), which passed in 1990.  Significant data from CRMS has accumulated within the past two decades, and is now allowing for detailed evaluation of accretion and other important processes.  A typical CRMS marsh location includes elevated boardwalks in an H-shaped configuration.  Accretion data is measured by a set of procedures in which a white feldspar powder is spread on the surface of the marsh at designated locations around the boardwalks.  At regular intervals sediment samples are recovered by inserting a tipped slim copper tube into the soil column and filling it with liquid nitrogen.  This process freezes the sediment around the tube and allows for the extraction of an undisturbed core-like sample.  The surface created by the feldspar powder is obvious as a white layer within the core.  The thickness of sediment accreted at the site can be measured above the white layer.  Over time multiple feldspar markers at multiple locations around each site provide enough data to allow for accurate estimates of long and short term accretion rates through statistical analysis.  The figure below illustrates the essential elements of accretion measurement process and a typical graphical representation of the data from one site.

Slide1.PNG

The methodology for measuring surface elevation change was established by Cahoon et al, 2000.  This methodology has been employed by CRMS to record data over the past two decades.  As shown in figure below, a rod surface elevation table (R-SET) is established by driving the rod to refusal and surveying its elevation.  A set of fiberglass rods are suspended from an arm with is perpendicular to the rod.  This allows for the direct measurement of the elevation of the surface of the marsh relative to the R-SET.  Changes in surface elevation are not measurements of positive or negative velocities relative to a fixed datum such as sea level, but only to the surveyed elevation of the R-SET.  In order to relate changes in surface elevation measured at a CRMS site to subsidence the site would have to be surveyed on a repeated basis to determine changes in elevation relative to the fixed datum, or the site would have to be equipped with a GPS device.

An evaluation of CRMS surface elevation change data was presented by Leigh Anne Sharp and Camille Stagg at the 2018 State of the Coast Conference.  They found “… 332 [CRMS sites] have surface elevation data that can be assessed along with land change data to infer which processes are influencing recent land change … In the Deltaic Plain land loss, in general, is not associated with elevation loss, suggesting that erosion, not subsidence, may be responsible for continuing land loss there. Over the last decade, wetland surface elevation trajectories have been positive at 75% of sites, including those in coastal areas near the Mississippi River delta that have seen the most historic land loss. Land change data reveal continuing land loss at just 14% of CRMS sites and land gain at 10% of CRMS sites including those in the birdsfoot delta and in the upper Mermentau basin. Furthermore, at least eight sites that were classified as floating marsh at project inception are now attached, including sites downstream from diversions and near the Atchafalaya delta. In general, CRMS data reveal that the coast of Louisiana is dynamic, that much of it has been stable in recent years, and that restoration efforts in concert with natural processes have increased stability in some areas”.

Slide2.PNG

The belt of the wetlands in which gains in land area have been measured follows the trend of what Shea Penland called the “Teche Shoreline”. This boundary is a near perfect coincidence of ecological, geological and soil composition boundaries. Geologically, the Teche Shoreline coincides with a line of faults across the marsh surface. These faults appear to have played a significant role in causing the wetlands loss event described in the previous blog. Hot spots of wetlands loss from this event line up along the hanging wall or downthrown side of the faults. Ecologically, the Teche Shoreline coincides with a sharp transition from freshwater marsh to intermediate, brackish and saline marshes. The succession of marsh type over time appears to be associated with subsidence. The coincidence of the faults and the ecosystem boundary suggests that subsidence caused by fault slip is driving the succession of ecosystems as the marshes are progressively submerged in more saline waters. The Teche Shoreline also coincides with a distinct boundary in composition of the marsh soils. Marshes to the north of the boundary generally have organic contents greater than 90%, while those to the south have lower organic compositions. It appears likely that long term accretion to the north of the boundary has promoted the retention of organic content over many centuries. To the south of the boundary marshes have been repeatedly submerged and removed by erosion, then rebuilt by new deltas. This cyclical pattern of marsh creation and submergence south of the Teche Shoreline has resulted in a lower organic content relative to the more stable area to the north.

The CRMS data indicates that areas north of the Teche Shoreline are naturally building new wetlands area in the face of subsidence and sea level rise. It seems logical that the sediment source with which new land area is being created is the same marsh material that was lost to submergence and erosion in the last wetlands loss event. This sediment is naturally carried to the north over time by tidal flux. The rates of sediment accretion and surface elevation change associated with land area gains suggests that these marshes can continue to maintain elevation under future scenarios of subsidence and sea level rise for many decades to come.

Teche Shoreline.png