Evaluating Land Subsidence Rates and Their Implications for Land Loss in the Lower Mississippi River Basin

High subsidence rates, along with eustatic sea-level change, sediment accumulation and shoreline erosion have led to widespread land loss and the deterioration of ecosystem health around the Lower Mississippi River Basin (LMRB). A proper evaluation of the spatial pattern of subsidence rates in the LMRB is the key to understanding the mechanisms of the submergence, estimating its potential impacts on land loss and the long-term sustainability of the region. Based on the subsidence rate data derived from benchmark surveys from 1922 to 1995, this paper constructed a subsidence rate surface for the region through the empirical Bayesian kriging (EBK) interpolation method. The results show that the subsidence rates in the region ranged from 1.7 to 29 mm/year, with an average rate of 9.4 mm/year. Subsidence rates increased from north to south as the outcome of both regional geophysical conditions and anthropogenic activities. Four areas of high subsidence rates were found, and they are located in Orleans, Jefferson, Terrebonne and Plaquemines parishes. A projection of future landscape loss using the interpolated subsidence rates reveals that areas below zero elevation in the LMRB will increase from 3.86% in 2004 to 19.79% in 2030 and 30.88% in 2050. This translates to a growing increase of areas that are vulnerable to land loss from 44.3 km2/year to 240.7 km2/year from 2011 to 2050. Under the same scenario, Lafourche, Plaquemines and Terrebonne parishes will experience serious loss of wetlands, whereas Orleans and Jefferson parishes will lose significant developed land, and Lafourche parish will endure severe loss of agriculture land.


Introduction
The Lower Mississippi River Basin (LMRB) located in southeastern coastal Louisiana is a major producer of crude oil and natural gas in the U.S., containing a large portion (40%-45%) of the nation's coastal wetlands and acting as a buffer zone for in-land residents from hurricanes and storms [1].Since 1930, however, this area has lost more than 4921 km 2 (~1900 mi 2 ) of land, which accounts for 80% of the total coastal wetland loss in the U.S. [2].The land loss problem causes severe damage to local fishery industries, deteriorates wetland ecosystem balance and increases the risk of coastal hazards to both coastal residents and energy infrastructures, which may cause thousands of fatalities and billions of economic loss.The widespread land loss around coastal Louisiana stemmed from the combination of land subsidence, eustatic sea-level change and shoreline erosion [3].However, high subsidence rates have been considered a principal cause contributing to the ongoing extensive wetland loss in coastal Louisiana [4].
Subsidence is defined as the downward shifting of land surface relative to a reference datum.Due to its impacts on the economy, livelihood and culture, a number of researchers have attempted to evaluate the subsidence rates in coastal Louisiana using different approaches, including analyses based on discrete tide gauge records [5] or benchmarks surveys [6], numerical modeling [7] and RADARSAT satellite images [8].Similarly, different factors driving the subsidence process have been investigated quantitatively, including consolidation of Holocene sediments, sediment loadings, sea-level rise, movements on growth faults and hydrocarbon production [9].
Despite these prior studies, few investigations examining vertical change in coastal Louisiana have focused on creating an accurate subsidence rate surface to reveal the spatial patterns and potential impacts on the region's long-term sustainability.The objectives of this investigation are three-fold.First, we will create a subsidence rate surface in the study area using the empirical Bayesian kriging (EBK) interpolation method.Second, future elevation will be projected based on the calculated subsidence rates.Third, possible land loss impacts caused by subsidence induced elevation change will be estimated.This study is significant in that it provides scenarios of future environmental changes caused by land subsidence in the deltaic region.The knowledge gained can be used to support decision makers to formulate plans for mitigation, adaptation and restoration in advance.
The remainder of this paper is organized as follows.Section 2 describes the study area and previous related studies.Section 3 introduces the major data sources and methods to create the subsidence rate layer.Section 4 presents the results, and the land loss implications are discussed on Section 5. Section 6 presents the conclusions of this study.

Study Area
The Lower Mississippi River Basin (LMRB) is located in southeastern Louisiana, the United States (Figure 1), with elevation ranging from ´7.3 to 370 meters and a mean of 10.7 meters.The study region consists of 26 parishes, covering the cities of New Orleans and Baton Rouge, 63 watersheds and three basins-Lake Pontchartrain, Terrebonne and Barataria basins [10].Its total area is around 48,046 km 2 .The Mississippi River flows into the LMRB from West Feliciana Parish in the northwest to the Plaquemines Parish in the southeast and eventually discharges to the Gulf of Mexico.
The study region is extremely vulnerable to land subsidence, which is the result of complex interactions among natural, social and economic conditions.Frequently exposed to hurricanes and storms, southern Louisiana is particularly vulnerable to coastal inundation hazards.A study on business recovery in New Orleans after Hurricane Katrina confirmed that flood depth is significantly negatively correlated with business reopening probabilities [11].It is projected that Katrina-like storm surges will hit every other year if the climate warms 2 ˝C (~3.6 ˝F) [12].Ongoing sinking will exacerbate these hazards and further delay recovery processes.The consequences of subsidence extend beyond flooding.Research has closely linked subsidence with land loss observed across the coastal zone.Reed and Cahoon (1993) suggested that a slight downwards shifting in elevation can lead to frequent flooding, which would erode vegetation and further accelerate the loss of wetlands in these areas [13].Aerial monitoring showed that the statewide wetland loss rates in Louisiana were 36 km 2 /year, 100 km 2 /year and 65 km 2 /year during 1930-1950, 1960-1980 and 1980-1990, respectively [14].Unchecked erosion brings the Gulf of Mexico closer to human habitats, raising the public's concerns about the future sustainability of subsiding communities.Wetland losses also threaten energy infrastructure.As a major supplier of crude oil with many gas wells and pipelines beneath the low elevation area, the relative sea level rise associated with subsidence and global eustasy will expose facilities to seawater, which can immerse coastal soils, corrode the energy facilities, as well as increase the risk of damage to energy infrastructures.Several studies have attempted to estimate subsidence rates in coastal Louisiana.Based on the stationary drilling data, Roberts et al. (1994) found that the average subsidence rates increased an order of magnitude (0.43-3.96 mm/year) from west to east and doubled (2.17-4.29 mm/year) from north to south around the Mississippi River Delta [15].Surveys by Penland and others [5,[16][17][18]] concluded that subsidence rates for the Mississippi Deltaic plain were 9-13 mm/year.Using the 2710 benchmarks during 1920 and 1995, Shinkle and Dokka (2004) found that the subsidence rates were substantially higher than rates reported in previous studies with a maximum rate of 51.92 mm/year, and the rates increased in many areas during the latter half of the 20th century [6].Further, Kent and Dokka (2013) utilized ordinary kriging spatial interpolation to estimate subsidence rates in coastal Louisiana, resulting in rates with a range of 3-17 mm/year [19].
To thoroughly capture the mechanism of downward shifting phenomenon in the LMRB, factors from environmental process and human activities should be considered and quantitatively measured.Four factors have been frequently mentioned in the literature, including geological factors, such as consolidation of Holocene sediments, lithospheric flexure response to sediment loadings, tectonic movements on growth faults and anthropogenic factors, mostly human-induced fluid extraction and hydrocarbon production [20][21][22][23][24]. Yuill et al. (2009) found that the range of rates caused by the four factors were 1.0-8.0,1.0-5.0,0.1-20 and 0-3 mm/year [9].Chan (2007) and Mallman and Zoback (2007) concluded that the combination of the first three environmental processes explained only about a 3-mm/year subsidence rate, whereas vertical change caused by hydrocarbon production induced fault reactivation, and reservoir compaction was primary responsible for the high subsidence rates, as the change of hydrocarbon production volumes coincide with the changes of land loss rates during 1920 and 1995 [1,25].
These previous studies on subsidence rates in this region were mostly based on an analysis of discrete points.However, the influence of land submergence is statewide and area based.An accurate evaluation of the spatial pattern and trend of subsidence rates can give new insights into the land loss problem.This study aims to build a bridge between current field measurements and regional influences of subsidence by interpolating the sampled point values into a continuous surface and to comprehend its impacts at a broader scale.Since the coastal zone is adjacent to open water and has suffered significant land loss during the past decade, the study area is divided into the north and the south regions according to the boundary along the northern shoreline of Lake Pontchartrain (Figure 1), to compare the subsidence rates and their potential impacts caused by subsidence between inland and coastal zones [26].This imaginary boundary roughly divides areas that have experienced significant population growth (the north) and areas that have suffered long-term population decline (the south).

Data Sources
Technical Report #50 (TR50) provides the most comprehensive geodetic study of vertical change ever conducted for the Louisiana Gulf coast [6].The report, published by the National Geodetic Survey (NGS), assessed vertical change for 2710 discrete reference benchmarks observed during 96 first-order geodetic leveling surveys conducted between 1920 and 1995.Relative differences in height were measured for observations that coincided with two or more surveys.To account for decades of uncorrected subsidence, benchmark heights were validated by using vertical change estimates derived from relative sea-level rates measured at tide gauges associated with the leveling surveys [6].Because the heights measured at each benchmark were minimally constrained (i.e., tied only to one known benchmark), the relative displacements between survey epochs were free from errors and distortions that typically propagate within constrained network adjustments [6].The subsidence rate for each benchmark location was determined by the difference between the two most recent measurements divided by the number of gap years and then subtracting the global sea-level rise rate: In Equation ( 1), R subsidence and R sea´level-rise represent the yearly subsidence and sea-level-rise rates, respectively.Y begin and Y end stand for the beginning and ending years of the survey period, and H begin and H end are the corresponding elevation readings of the benchmark at the beginning and the ending years.Shinkle and Dokka adopted a constant sea-level rise rate of 1.25 mm/year.Of the 1178 benchmarks located in the LMRB (Figure 1), 385 sites are located in the north, and their rates ranged from ´0.42 mm/year (meaning no subsidence, but rather, minor uplifting) to 49.28 mm/year, with an average rate of 7.72 mm/year.The other 777 locations are in the south, and their subsidence rates ranged from 0.0 to 51.94 mm/year with a mean of 10.95 mm/year.
Other data involved in this study include elevation and land use/land cover products.Elevation data were obtained from the National Elevation Dataset (NED) of the U.S. Geological Survey (USGS), updated in 2004.Land use and land cover data were acquired from the National Land Cover Database (NLCD).Before analysis, the original 15 land cover types were re-aggregated to the 7 categories corresponding to Anderson's first-level LULC classification [27], and they include water, developed land, barren, forest, grass, agricultural land and wetlands.The resolutions for both datasets were 30 ˆ30 meters.

Empirical Bayesian Kriging
To create a subsidence rate surface accurately in the study area based on the sample points, it is necessary to choose a suitable spatial interpolation method to estimate values at unknown locations from observed values.Kriging is a popular geo-statistical interpolation technique that returns the best linear unbiased prediction of the intermediate data and avoids cluster effects, which is especially important for this investigation, since our data are quite spatially clustered.Kriging assumes that the spatial variation of an attribute consists of a trend m(Z) and a spatially-correlated component [28]: where Z 0 is the estimated value at one unknown location, m is the number of points in a neighborhood and W i and Z i are the weights and observed values of the i-th surrounding observation, respectively.Different treatments of the three components-the trend, the residual variogram and the error term-led to the development of various kriging methods [29].Here, we will briefly explain the principle of ordinary kriging assumed without a trend component as an example.To measure the spatially-correlated component, ordinary kriging uses isotropic semivariance, which is computed as: where γ phq is computed as half of the average squared difference between one location and another at a distance of h, with their values Z(0) and Z(h).The relationship between lag distance and corresponding semivariance is described as a model, such as spherical, exponential and Gaussian, with parameters nugget, sill and range [30].The weights of the surrounding observations are determined through the semivariance matrix between data points: where W is the vector of weights, V is the semivariance matrix between surrounding known data pairs and D is the semivariance between unknown and known points.Lagrange multiplier λ is added to guarantee the minimum estimation error subject to the constraint that the sum of weights is equal to 1.An advantage of this interpolation routine is that the standard error for each estimate can be calculated and mapped according to: where σ is the standard error and γ ph i0 q is the semivariance between the i-th known points and the point to be estimated.A kriging error map can be used to indicate where future sampling should be conducted to reduce the interpolation error [31].However, classical kriging assumes that data are generated from a Gaussian distribution with the correlation structure defined by the estimated semivariogram, which is difficult to establish in practice.
Bayesian statistics provides a statistically-robust approach for modeling the uncertainties with respect to the unknown distribution and parameters of the input samples [32].Empirical Bayesian kriging (EBK) combines Bayes' theorem and kriging interpolation and accounts for the error in estimating the true semivariogram through iterative simulations [33].Figure 2 explains the procedures of EBK interpolation used in this investigation.First, the input sample data are divided into I subsets with less than or equal to N samples.For each subset, observations are transformed to a Gaussian distribution, and a semivariogram model is derived.Based on the semivariogram, new data are simulated and back transformed at known points.This process is iterated for n times, and each time, the new data produces a new semivariogram.Using Bayes' rule, the weights for each semivariogram are computed: where θ i is the i-th set of semivariogram parameters nugget, sill and range.W pθ i |Zq is the weight for the i-th semivariogram; f pZ|θ i q evaluates the likelihood the observed data can be generated from the semivariogram; and P pθ i q stands for the probability of the i-th set of parameters θ i among the simulated semivariogram spectrum.The weighted sum of simulated semivariograms creates a "true" semivariogram model.For each location, the prediction is generated using the "new" semivariogram distribution in the point's neighborhood.Compared to classical kriging, EBK is more appropriate for interpolating non-stationary data for large areas and requires less sample data [32].

Results
In this investigation, the number N of data points in each subset was set to 100.For each subset of data, an empirical transformation was applied to satisfy the assumption of normality.The overlap factor x representing the number of subsets in which an observation will participate was set to 1.0 (meaning each observation is only being considered once), and the number of simulated semivariograms n was 100.The neighborhood search radius R was defined as 86 kilometers, which was determined by Moran's I autocorrelation analysis as the distance at which Moran's I is the highest [34].Since the output pixel size is 30 ˆ30 meters, which is very time consuming for a study area with nearly 50 million pixels, we used the parallel computing option available in ArcGIS to accelerate the processing.The parameter was set to 100%, meaning that all of the cores in a computer were used in parallel computing.It took 490 h to complete the whole process on an eight-core CPU workstation.
Figure 3 and Table 1 shows the interpolated subsidence rates.There are active areas with high subsidence rates above 20 mm/year.They are located at Terrytown in Jefferson parish (#1), southern Orleans parish along the Intracoastal Waterway (#2), the wetlands between Morgan City and Houma in Terrebonne parish (#3) and the deltaic zones around the Triumph community in Plaquemines parish (#4).The standard errors of the estimated interpolated rates were less than 2.5 mm/year in 53.12% and less than 5 mm/year in 91.84% of the research area (Figure 4).Higher estimation errors are generally located in the southeastern part of coastal Louisiana and selected areas in New Orleans.
We compared our estimation to previous investigations at the same location for further evaluation (Table 2).We also conducted an interpolation using ordinary kriging (OK) for comparison.All subsidence rates in the table are expressed as millimeters per year.In the New Orleans metropolitan area, interpolated rates from EBK are similar to the observations from Dixon et al.  [37][38][39], have used the same TR50 benchmark survey data from NGS and obtained similar subsidence estimations at various parts of the study region.Consistency between ranges in different sets of results demonstrated that our interpolation is reliable and accurate.Furthermore, a comparison between EBK and OK shows that interpolated subsidence rates from EBK are closer to previous studies than the values derived by OK.In general, subsidence rates are lower in the north inland area and higher in the south coastal zone (Figure 3 and Table 1).This north-south discrepancy is hypothesized as the outcome of both regional geophysical conditions and anthropogenic activities.The south zone contains numerous submerged vertical faults.For instance, Thibodaux fault crosses the city of New Orleans, and Theriot and Golden Meadow faults pass through St. Mary, Terrebonne, Lafourche and Plaquemines parishes [41,42].The structure of the Mississippi Delta is largely based on deposited salt structures derived from an underlying autochthonous Jurassic salt [43].The upward intrusion of salt into fault zones, known as "salt diapir", may trigger proximally located faults and create new radial fault zones, thus contributing to land sinking.Another significant factor is sediment compaction, which is strongly affected by human activities.Hydrocarbon production has had a significant presence in the South for more than a century.Withdrawal of subsurface petroleum, natural gas and significant quantities of groundwater triggers the loss of subsurface pore pressure, which in turn accelerates natural consolidation processes and increases the degree of land subsidence [7,21].Theoretically, both the tectonic and anthropogenic processes are localized, with the former being a longer term process.Previous investigations in coastal Louisiana suggested that long-term fault slip led to a 0.1-20 mm/year subsidence rate, and human activities aggravated the subsidence process by contributing up to a 23 mm/year subsidence rate in some areas [9,23,44].
All four high subsidence areas are located in the south zone, including two locations in greater New Orleans.The urban growth of New Orleans has led to extensive networks of pumps and canals designed to move water away from residential areas.These canals and pumps were mostly positioned on the natural levee of the Mississippi River.As the city grew, the pumps and canals were expanded to move water into Lake Pontchartrain.Seasonal high water levels and the threat from storm surge triggered levee construction along the river and related industrial navigation canals, disconnecting the soils from the natural hydrography and triggering sediment compaction and consolidation.Meanwhile, the wet and dry process caused by precipitation and pumping activities extracts water from the ground, making soil particles collapse onto each other and inducing more subsidence.Furthermore, tapping into subsurface aquifers for industrial purposes could activate fault slip and shallow surface deformations, particularly in New Orleans East [45].
Another two hotspots are located in the coastal wetlands in Terrebonne and Plaquemines parishes.The Terrebonne Basin located in Terrebonne parish consists of thick uncompacted sediment areas and is subject to high compaction, which contributes to high subsidence [46].Although human-related hydrocarbon production could induce surficial changes and subsidence, whether this factor has contributed to high subsidence in the Terrebonne Basin remains controversial.Coleman and Roberts (1989) and Boesch et al. (1994) [3,47] concluded that the production impacts on wetland subsidence are minimal, while Morton et al. (2003) [48] speculated that the high density of oil/gas wells and related energy infrastructures may be associated with higher levels of land subsidence.On the other hand, human activities were considered to have a huge impact on the wetlands subsidence in Plaquemines.Levees along the Mississippi River removed the natural, renewal sediment flows that are needed to keep the Mississippi River Delta from sinking.Plaquemines parish is the site of several oil refineries and a base for assistance to offshore oilrigs.Oil and gas extraction has also hastened land subsidence protected by the levees in Plaquemines [49].

Discussion
Elevation change is a complex process caused by subsidence, as well as changes in sedimentation and human activities.Decreasing elevation does not appear to affect land loss for areas with high elevation or long distances from open water.For coastal area, current elevation strongly determines the impact of subsidence on land loss in the LMRB, especially for the southern parishes adjacent to the Gulf of Mexico.The bubble graph in Figure 5 illustrates the averages of subsidence rates and elevations, together with the population for each parish to explore the relationships of the three elements at the parish scale.The circle size represents the total population within each parish in 2012.For ease of visualization, the north inland parishes of Washington, West Feliciana, East Feliciana and St. Helena are not shown because of their high mean elevation (40.96-61.78meters) and relatively low subsidence rates (4.1-6.6 mm/year).Parishes in the south zone are closer to the Gulf of Mexico and have lower elevation and higher subsidence rates.Subsidence-induced elevation change may cause severe land loss and hazard damages in the future, since some of the parishes in the south have large populations such as Orleans and Jefferson.Using the interpolated rates, we projected the future elevations in 2030 and 2050 assuming that elevation is under the influence of subsidence only, as shown in Figure 6.Areas below sea level are highlighted as red areas and defined as vulnerable regions, which have high potential for land loss.In 2004, only 3.86% of the whole study region, including 0.69% in the north and 7.17% in the south, is below sea level.The most vulnerable regions are located in Orleans parish.If subsidence continues at the current rate, the proportions of vulnerable regions in the north, the south and the whole LMRB will increase to 2.32%, 37.95% and 19.79% in 2030 and to 5.63%, 57.32% and 30.88% in 2050, respectively.In addition to Orleans parish, five more parishes, including Terrebonne, St. Bernard, Plaquemine, Jefferson and Lafourche, will have more than 40% vulnerable regions in their parishes in 2030.By 2050, three more parishes, St. Charles, St. John the Baptist and St. Mary parishes, will join the list of parishes that have more than 40% of the parish area below sea level.most immediate consequence of subsidence-induced elevation change is land loss.Land loss risks will be greater in the regions closest to open water.By 2050, the yearly growth of vulnerable lands will increase from 44.3 km 2 /year in 2011 to 240.7 km 2 /year in 2050 (Figure 7a).Further, land cover types were overlapped with potential land loss regions.The results indicate that by 2050, Lafourche, Plaquemines and Terrebonne parishes will suffer severe loss of wetlands; Orleans and Jefferson parishes will endure serious loss of developed land; and heavy agricultural land loss will take place in Lafourche Parish (Figure 7b).In reality, subsidence does not necessarily result in land loss.There are multiple, complex natural-human interactions that impact loss and alter local elevation, including sedimentation, hydrocarbon production and construction of levees and canals.Therefore, our estimation may exaggerate the real land loss during the evolution of coastal Louisiana, since the above scenario is based on a single-parameter (subsidence rate) input model.However, land subsidence will still impact population livelihood in profound ways.To protect communities from the consequences of subsidence, such as frequent flooding, governments will have to invest significantly in flood protection levees and pumping stations.Furthermore, Louisiana is a major agricultural state, producing large amounts of rice (ranked third nationally), corn, sugarcane, soybeans and others.Elevation shifting changes the ground's moisture and salinity, which will alter soil availability and decrease yields.Moreover, Louisiana coastal wetlands serve as buffering zones for storms and floods, habitats for aquatic creatures and the base for ecosystem food webs [50].The abundant aquaculture resources of fishes, oysters, shrimps and crabs contribute to 25% of all seafood in the U.S. and attract numerous tourists from out of state [51].Even a slight elevation change in wetlands will jeopardize the ecological balance of the wetlands environment, which not only raises the risk of destructive marine forces for human beings, but also threatens local fisheries and tourism.
Despite its contribution towards understanding coastal vertical change, TR50 has been criticized for over-estimating regional subsidence rates [52,53].Because leveling surveys primarily occurred along transportation corridors, critics contend that the findings in TR50 fail to represent vertical change within the intermediate coastal prairie and wetlands.However, the comparison between EBK, ordinary kriging and the estimates from previous studies (Table 2) reveals comparable findings.The city of New Orleans, which is adjacent to transportation corridors, shows the highest agreements.EBK estimated rates are slightly higher (1-5 mm/year) in most cases of coastal wetlands, except for the Mississippi Deltaic estuary in Plaquemines parish.In general, the subsidence rate layer obtained in this investigation is a reliable reference for modeling scenarios of future environmental changes and identifying regions vulnerable to land loss.

Conclusions
This study created a surface representing the general subsidence pattern and evaluated its potential impacts on land loss in the Lower Mississippi River Basin (LMRB).Using the 1178 benchmarks' elevation data for 1922-1995, subsidence rates at 30 m ˆ30 m were calculated under the consideration of global sea level rise using the ordinary kriging (OK) and empirical Bayesian kriging (EBK) interpolation methods.The results show that the subsidence rates derived by OK and EBK for the entire region ranged from 2.6 to 21.5 and from 1.7 to 29 mm/year, with an average rate of 9 and 9.4 mm/year, respectively.Compared to OK, subsidence rates obtained by EBK have higher agreements with previously-published findings.The subsidence rates increased from north (7.5 ˘5 mm/year) to south (11 ˘10 mm/year).By creating surface contours, four areas of high subsidence rates were observed, including two regions in New Orleans and wetlands in Terrebonne and Plaquemines parishes.The former two could be a result of human activities, such as frequent pumping of water and withdrawal of underground water, whereas the latter could be associated with local soil conditions, levee constructions and heavy hydrocarbon production.By combining the subsidence rates with elevation data in 2004, we projected the elevations in 2030 and 2050 and found that areas below sea level will increase from 3.86% in 2004 to 19.79% in 2030 and 30.88% in 2050.By tabulating the land loss projections with land use and land cover data, we found that heavy wetland loss would occur in Lafourche, Plaquemines and Terrebonne parishes.Orleans and Jefferson will suffer severe loss of developed land, and Lafourche will endure high agricultural land loss.
The significance of this investigation is that it produced a subsidence rate surface for the study region using a rigorous geostatistical interpolation procedure.It then provides scenarios of land subsidence that could threaten the region's sustainability.With sea level rising, the problem could be even more serious.The findings from this study will help decision makers formulate better plans for the region.By identifying the developed areas that have high potential of land loss in the following decades, planners could build more levees to protect this area or avoid the development in those areas.New development could be located further inland while allowing access to the coast for economic activities.Information on potential developed land loss will also be beneficial to residents so that they could plan ahead for protecting their properties and lives, such as lifting their houses or moving to higher ground.In future work, Global Positioning System (GPS)-measured subsidence rate and elevation data will be introduced to validate and calibrate our current estimation to provide a more accurate subsidence rate layer.In addition, frequent hurricanes and flooding in coastal Louisiana will distribute sediment to vulnerable regions and compensate the elevation decrease caused by subsidence.Constructed levees protect low elevation areas (e.g., New Orleans) from both river flooding and land loss.Hence, factors such as sedimentation rates, hydrocarbon production and levee distribution will be considered in modeling elevation change and estimating land loss to project a more realistic future landscape.

Figure 3 .
Figure 3. Interpolated subsidence-rate map in the Lower Mississippi River Basin.

Figure 4 .
Figure 4. Standard errors of interpolated subsidence rates.

Figure 5 .
Figure 5. Parish-level average elevation, subsidence rates and population in the LMRB in 2012 (Washington, West Feliciana, East Feliciana and St. Helena parishes are not included).

Figure 6 .
Figure 6.Projected elevation map with areas below sea level.

Figure 7 .
Figure 7. (a) Potential land loss rate from 2011 to 2050; (b) potential land loss types and areas by parish in the south region.

Table 1 .
Statistics of interpolated subsidence rates.

Table 2 .
Subsidence rates (in mm/year) from the literature and our results.