Assessing visual disturbance conditions on the Custer National Forest

Assessment of soil disturbance on the Custer National Forest was conducted during two summers to determine if the U.S. Forest Service Forest Soil Disturbance Monitoring Protocol (FSDMP) was able to distinguish post-harvest soil conditions in a chronological sequence of sites harvested using different ground-based logging systems. Results from the first year of sampling suggested that the FSDMP point sampling method may not be sensitive enough to measure post-harvest disturbance in stands with low levels of disturbance. Therefore, a revised random transect method was used during the second sampling season to determine the actual extent of soil disturbance in these cutting units. Using combined data collected from both summers I detected statistically significant differences (p < 0.05) in fine fraction bulk density measurements between FSDMP disturbance classes across all sites. Disturbance class 3 (most severe) had the highest reported bulk density, which suggest that the FSDMP visual class estimates are defined adequately allowing for correlations to be made between visual disturbance and actual soil physical characteristics. Forest site productivity can be defined by its ability to retain carbon and convert it to above- and belowground biomass. However, forest management activities that alter basic site characteristics have the potential to alter productivity. Soil compaction is one critical management impact that is important to understand; compaction has been shown to impede the root growth potential of plants, reduce water infiltration rates increasing erosion potential, and alter plant available water and nutrients, depending on soil texture. A new method to assess ground cover, erosion, and other soil disturbances was recently published by the U.S. Forest Service, as the Forest Soil Disturbance Protocol (FSDMP). The FSDMP allows soil scientists to visually assign a disturbance class estimate (0 – none, 1, 2, 3 – severe) from field measures of consistently defined soil disturbance indicators (erosion, fire, rutting, compaction, and platy/massive/puddled structure) in small circular (15 cm) plots to compare soil quality properties pre- and post- harvest condition. Using this protocol we were able to determine that ground-based timber harvesting activities occurring on the Custer National Forest are not reaching the 15% maximum threshold for detrimental soil disturbance outlined by the Region 1 Soil Quality Standards.

Patch selection by wintering white-tailed deer : consequences for eastern hemlock regeneration at local and landscape scales

Tsuga canadensis (eastern hemlock) is a highly shade-tolerant, late-successional, and long-lived conifer species found throughout eastern North America. It is most often found in pure or nearly pure stands, because highly acidic and nutrient poor forest floor conditions are thought to favor T. canadensis regeneration while simultaneously limiting the establishment of some hardwood species with greater nutrient requirements. Once a common species, T. canadensis is currently experiencing widescale declines across its range. The hemlock woolly adelgid (Adelges tsugae) is decimating the population across its eastern distribution. Across the Upper Great Lakes region, where the adelgid is currently being held at bay by cold winter temperatures, T. canadensis has been experiencing failures in regeneration attributed, in part, to herbivory by white-tailed deer (Odocoileus virginianus). Deer utilize T. canadensis stands as winter habitat in areas of high snow depth. Tsuga canadensis, once a major component of these forests, currently exists at just a fraction of its pre-settlement abundance due to historic logging and contemporary forest management practices, and what remains is found in small remnant patches surrounded by second- and third-growth deciduous forests. The deer population across the region, however, is likely double that of pre-European settlement times. In this dissertation I explore the relationship between white-tailed deer use of T. canadensis as winter habitat and the effect this use is having on regeneration and forest succession. For this research I quantified stand composition and structure and abiotic variables of elevation and snow depth in 39 randomly selected T. canadensis stands from across the western Upper Peninsula of Michigan. I also quantified composition and the configuration of the landscapes surrounding these stands. I measured relative deer use of T. canadensis stands as pellet group piles deposited in each stand during each of three consecutive winters, 2005-06, 2006-07, and 2007-08. The results of this research suggest that deer use of T. canadensis stands as winter habitat is influenced primarily by snow depth, elevation, and the composition and configuration of the greater landscapes surrounding these stands. Specifically, stands with more heterogeneous landscapes surrounding them (i.e., a patchy mosaic of conifer, deciduous, and open cover) had higher relative deer use than stands surrounded by homogenous deciduous forest cover. Additionally, the intensity of use and the number of stands used was greater in years with higher average snow depth. Tsuga canadensis regeneration in these stands was negatively associated with deer use and Acer saccharum (sugar maple) basal area. Of the 39 stands, 17 and 22 stands had no T. canadensis regeneration in small and large sapling categories, respectively. Acer saccharum was the most common understory tree species, and the importance of A. saccharum in the understory (stems < 10 cm dbh) of the stands was positively associated with overstory A. saccharum dominance. Tsuga canadensis establishment was associated with high-decay coarse woody debris and moss, and deciduous leaf litter inputs in these stands may be limiting access to these important microsites. Furthermore, A. saccharum is more tolerant to the effects of deer herbivory than T. canadensis, giving A. saccharum a competitive advantage in stands being utilized as winter habitat by deer. My research suggests that limited microsite availability, in conjunction with deer herbivory, may be leading to an erosion in T. canadensis patch stability and an altered successional trajectory toward one of A. saccharum dominance, an alternately stable climax species.

Dynamics and kinematics of eruptive activity at Fuego volcano, Guatemala 2005-2009

Volcanoes are the surficial expressions of complex pathways that vent magma and gasses generated deep in the Earth. Geophysical data record at least the partial history of magma and gas movement in the conduit and venting to the atmosphere. This work focuses on developing a more comprehensive understanding of explosive degassing at Fuego volcano, Guatemala through observations and analysis of geophysical data collected in 2005 – 2009. A pattern of eruptive activity was observed during 2005 – 2007 and quantified with seismic and infrasound, satellite thermal and gas measurements, and lava flow lengths. Eruptive styles are related to variable magma flux and accumulation of gas. Explosive degassing was recorded on broadband seismic and infrasound sensors in 2008 and 2009. Explosion energy partitioning between the ground and the atmosphere shows an increase in acoustic energy from 2008 to 2009, indicating a shift toward increased gas pressure in the conduit. Very-long-period (VLP) seismic signals are associated with the strongest explosions recorded in 2009 and waveform modeling in the 10 – 30 s band produces a best-fit source location 300 m west and 300 m below the summit crater. The calculated moment tensor indicates a volumetric source, which is modeled as a dike feeding a SW-dipping (35°) sill. The sill is the dominant component and its projection to the surface nearly intersects the summit crater. The deformation history of the sill is interpreted as: 1) an initial inflation due to pressurization, followed by 2) a rapid deflation as overpressure is explosively release, and finally 3) a reinflation as fresh magma flows into the sill and degasses. Tilt signals are derived from the horizontal components of the seismometer and show repetitive inflation deflation cycles with a 20 minute period coincident with strong explosions. These cycles represent the pressurization of the shallow conduit and explosive venting of overpressure that develops beneath a partially crystallized plug of magma. The energy released during the strong explosions has allowed for imaging of Fuego’s shallow conduit, which appears to have migrated west of the summit crater. In summary, Fuego is becoming more gas charged and its summit centered vent is shifting to the west - serious hazard consequences are likely.

Automation of shear-wave splitting parameter determination of local earthquakes at Yellowstone : application as indicator of crustal stress and temporal variation

Shear-wave splitting can be a useful technique for determining crustal stress fields in volcanic settings and temporal variations associated with activity. Splitting parameters were determined for a subset of local earthquakes recorded from 2000-2010 at Yellowstone. Analysis was automated using an unsupervised cluster analysis technique to determine optimum splitting parameters from 270 analysis windows for each event. Six stations clearly exhibit preferential fast polarization values sub-orthogonal to the direction of minimum horizontal compression. Yellowstone deformation results in a local crustal stress field differing from the regional field dominated by NE-SW extension, and fast directions reflect this difference rotating around the caldera maintaining perpendicularity to the rim. One station exhibits temporal variations concordant with identified periods of caldera subsidence and uplift. From splitting measurements, we calculated a crustal anisotropy of ~17-23% and crack density ~0.12-0.17 possibly resulting from stress-aligned fluid filled microcracks in the upper crust and an active hydrothermal system.

Body sensor network for in-home personal healthcare

A body sensor network solution for personal healthcare under an indoor environment is developed. The system is capable of logging the physiological signals of human beings, tracking the orientations of human body, and monitoring the environmental attributes, which covers all necessary information for the personal healthcare in an indoor environment. The major three chapters of this dissertation contain three subsystems in this work, each corresponding to one subsystem: BioLogger, PAMS and CosNet. Each chapter covers the background and motivation of the subsystem, the related theory, the hardware/software design, and the evaluation of the prototype’s performance.

Geophysical and geological analysis of the Collingwood Member of the Trenton Formation

The Collingwood Member is a mid to late Ordovician self-sourced reservoir deposited across the northern Michigan Basin and parts of Ontario, Canada. Although it had been previously studied in Canada, there has been relatively little data available from the Michigan subsurface. Recent commercial interest in the Collingwood has resulted in the drilling and production of several wells in the state of Michigan. An analysis of core samples, measured laboratory data, and petrophysical logs has yielded both a quantitative and qualitative understanding of the formation in the Michigan Basin. The Collingwood is a low permeability and low porosity carbonate package that is very high in organic content. It is composed primarily of a uniformly fine grained carbonate matrix with lesser amounts of kerogen, silica, and clays. The kerogen content of the Collingwood is finely dispersed in the clay and carbonate mineral phases. Geochemical and production data show that both oil and gas phases are present based on regional thermal maturity. The deposit is richest in the north-central part of the basin with thickest deposition and highest organic content. The Collingwood is a fairly thin deposit and vertical fractures may very easily extend into the surrounding formations. Completion and treatment techniques should be designed around these parameters to enhance production.

Management impacts on Carex assiniboinensis and Cardamine concatenata on the Ottawa National Forest

Monitoring of herbaceous plants on the Ottawa National Forest (ONF) is used to understand the impact of forest management on understory composition and site conditions. In their planning, national forests are required to take into account management impacts on diversity and ecosystem health. The effect of management on understory species is dependent on various factors, including the intensity of disturbance and the biology of the plant. In the first study in this report, a population of Carex assiniboinensis, a Michigan state threatened species, was monitored for seven seasons including before logging commenced, in order to determine the sedge’s response to a single-tree selection harvest. Analyses provided insights for management of C. assiniboinensis at the stand level over the short-term. In the second study in this report, the use of the cutleaf toothwort (Cardamine concatenata) as a Management Indicator Species on the ONF was reviewed. Data were analyzed to determine the suitability of using C. concatenata to monitor impacts of forest management on site conditions. The various factors that affect understory species population dynamics illuminated the challenges of using indicator species to monitor site conditions. Insights from the study provide a greater understanding of management impacts on understory species across the Ottawa National Forest.

Hydrogeochemical controls on uranium in aquifers of the Jacobsville Sandstone

The occurrence of elevated uranium (U) in sandstone aquifers was investigated in the Upper Peninsula of Michigan, focusing on aquifers of the Jacobsville Sandstone. The hydrogeochemical controls on groundwater U concentrations were characterized using a combination of water sampling and spectral gamma-ray logging of sandstone cliffs and residential water wells. 235U/238U isotope ratios were consistent with naturally occurring U. Approximately 25% of the 270 wells tested had U concentrations above the U.S. Environmental Protection Agency Maximum Contaminant Level (MCL) of 30 μg/L, with elevated U generally occurring in localized clusters. Water wells were logged to determine whether groundwater U anomalies could be explained by the heterogeneous distribution of U in the sandstone. Not all wells with relative U enrichment in the sandstone produced water with U above the MCL, indicating that the effect of U enrichment in the sandstone may be modified by other hydrogeochemical factors. Well water had high redox, indicating U is in its highly soluble (VI) valence. Equilibrium modeling indicated that aqueous U is complexed with carbonates. In general, wells with elevated U concentrations had low 235U/238U activity ratios. However, in some areas U concentrations and 235U/238U activity ratios were simultaneously high, possibly indicating differences in rock-water interactions. Limited groundwater age dating suggested that residence time may also help explain variations in well water U concentrations. Low levels of U enrichment (4 to 30 ppm) in the Jacobsville sandstone may make wells in its oxidized aquifers at risk for U concentrations above the MCL. On average, U concentrations were highest in heavy mineral and clay layers and rip up conglomerates. Uranium concentrations above 4 ppm also occurred in siltstones, sandstones and conglomerates. Uranium enrichment was likely controlled by deposition processes, sorption to clays, and groundwater flow, which was controlled by permeability variations in the sandstone. Low levels of U enrichment were found at deltaic, lacustrine and alluvial fan deposits and were not isolated to specific depositional environments.

The influence of concurrent disturbances on plant community dynamics in northern hemlock-hardwood forests

Throughout the Upper Great Lakes region, alterations to historic disturbance regimes have influenced plant community dynamics in hemlock-hardwood forests. Several important mesic forest species, eastern hemlock (Tsuga canadensis), yellow birch (Betula alleghaniensis), eastern white pine (Pinus strobus), and Canada yew (Taxus canadensis), are in decline due to exploitive logging practices used at the turn of the 20th century and the wave of intense fires that followed. Continued regeneration and recruitment failure is attributed to contemporary forest management practices and overbrowsing by white-tailed deer (Odocoileus virginianus). Therefore, I examined the influence of two concurrent disturbances, overstory removal and herbivory, on plant community dynamics in two hemlock-hardwood forests. I measured the post-disturbance regeneration response (herbaceous and woody species) inside and outside of deer exclosures in 20 artificial canopy gaps (50 – 450 m2) and monitored survival and growth for hundreds of planted seedlings. The results of this research show that interacting disturbances can play a large role in shaping plant community composition and structure in hemlock-hardwood forests. White-tailed deer herbivory homogenized the post-disturbance plant communities across the experimental gradient of gap areas, essentially making species compositions in small gaps “look like” those in large gaps. Deer browsing also influenced probability of survival for planted Canada yew cuttings; all else being equal an individual was nearly seven times more likely to survive if protected from herbivory (P < 0.001). In contrast, the ability of sugar maple (Acer saccharum) to persist under high levels of herbivory and respond rapidly to overstory release appears to be related to the presence of stem layering(i.e., portions of below-ground prostrate stem). Layering occurred in 52% of excavated saplings (n = 100) and was significantly associated with increased post-disturbance height growth. Understory light was also important to planted seedling establishment and height growth. Higher levels of direct under-canopy light negatively impacted survival for shade-tolerant hemlock and Canada yew, while an increase in diffuse light was linked to a higher probability of survival for yellow birch and height growth for hemlock and Canada yew. Increases in white pine height growth were also significantly associated with a decrease in canopy cover.

Triggering of volcanic activity by large earthquakes

Statistical analyses of temporal relationships between large earthquakes and volcanic eruptions suggest seismic waves may trigger eruptions even over great (>1000 km) distances, although the causative mechanism is not well constrained. In this study the relationship between large earthquakes and subtle changes in volcanic activity was investigated in order to gain greater insight into the relationship between dynamic stresses propagated by surface waves and volcanic response. Daily measurements from the Ozone Monitoring Instrument (OMI), onboard the Aura satellite, provide constraints on volcanic sulfur-dioxide (SO2) emission rates as a measure of subtle changes in activity. Time series of SO2 emission rates were produced from OMI data for thirteen persistently active volcanoes from 1 October 2004 to 30 September 2010. In order to quantify the affect of earthquakes at teleseismic distances, we modeled surface-wave amplitudes from the source mechanisms of moment magnitude (Mw) ≥7 earthquakes, and calculated the Peak Dynamic Stress (PDS). We assessed the influence of earthquakes on volcanic activity in two ways: 1) by identifying increases in the SO2 time series data and looking for causative earthquakes and 2) by examining the average emission rate before and after each earthquake. In the first, the SO2 time series for each volcano was used to calculate a baseline threshold for comparison with post-earthquake emission. Next, we generated a catalog of responses based on sustained SO2 emission increases above this baseline. Delay times between each SO2 response and each prior earthquake were analyzed using both the actual earthquake catalog, and a randomly generated catalog of earthquakes. This process was repeated for each volcano. Despite varying multiple parameters, this analysis did not demonstrate a clear relationship between earthquake-generated PDS and SO2 emission. However, the second analysis, which was based on the occurrence of large earthquakes indicated a response at most volcanoes. Using the PDS calculations as a filtering criterion for the earthquake catalog, the SO2 mass for each volcano was analyzed in 28-day windows centered on the earthquake origin time. If the average SO2 mass after the earthquake was greater than an arbitrary percentage of pre-earthquake mass, we identified the volcano as having a response to the event. This window analysis provided insight on what type of volcanic activity is more susceptible to triggering by dynamic stress. The volcanoes with very open systems included in this study, Ambrym, Gaua, Villarrica, Erta Ale and, Turrialba, showed a clear response to dynamic stress while the volcanoes with more closed systems, Merapi, Semeru, Fuego, Pacaya, and Bagana, showed no response.

Determining the role of environmental factors and disturbance in the distribution of reed canary grass within wetlands

Reed canary grass (Phalaris arundinacea L.) is an invasive species originally from Europe that has now expanded to a large range within the United States. Reed canary grass possesses a number of traits that allow it to thrive in a wide range of environmental factors, including high rates of sedimentation, bouts of flooding, and high levels of nutrient inputs. Therefore, the goals of our study were to determine if 1) certain types of wetland were more susceptible to Reed canary grass invasion, and 2) disturbances facilitated Reed canary grass invasion. This study was conducted within the Keweenaw Bay Indian Community reservation in the Upper Peninsula of Michigan, in Baraga County. We selected 28 wetlands for analysis. At each wetland, we identified and sampled distinct vegetative communities and their corresponding environmental attributes, which included water table depth, pH, conductivity, calcium and magnesium concentrations, and percent organic matter. Disturbances at each site were catalogued and their severity estimated with the aid of aerial photos. A GIS dataset containing information about the location of Reed canary grass within the study wetlands, the surrounding roads and the level of roadside Reed canary grass invasion was also developed. In all, 287 plant species were identified and classified into 16 communities, which were then further grouped into three broad groupings of wetlands: nonforested graminoid, Sphagnum peatlands, and forested wetlands. The two most common disturbances identified were roads and off-road recreation trails, both occurring at 23 of the 28 sites. Logging activity surrounding the wetlands was the next most common disturbance and was found at 18 of the sites. Occurrence of Reed canary grass was most common in the non-forested graminoid communities. Reed canary grass was very infrequent in forested wetlands, and almost never occurred in the Sphagnum peatlands. Disturbance intensity was the most significant environmental factor in explaining Reed canary grass occurrence within wetlands. Statistically significant relationships were identified at distances of 1000 m, 500 m, and 250 m from studied wetlands, between the level of road development and the severity of Reed canary grass invasion along roadsides. Further analysis revealed a significant relationship between roadside Reed canary grass populations and the level of road development (e.g. paved, graded, and ungraded).

Paleomagnetic Study of the Portage Lake Volcanics exposed in the Quincy Mine

A detailed paleomagnetic and rock-magnetic investigation was conducted on thirty six basaltic flows of the ~1095 Ma Portage Lake Volcanics. The flows were sampled along the East Adit of the Quincy Mine (Hancock, MI). Thirty two flows yielded well-defined primary magnetization directions carried by magnetite. A secondary magnetization component carried by hematite was also found in twenty nine flows. After correction for serial correlation between the flows, nineteen independent mean directions were calculated. The corresponding paleomagnetic pole is located at 25.5 °N, 182.1 °W (A95 = 3.5°). The new pole overlaps with the pole from the ~1087 Ma Lake Shore Traps suggesting a standstill of the North American plate during that time period. The low angular dispersion of virtual geomagnetic poles (S = 7.9°) suggests that the flows were erupted within a short time period, or that the strength of geomagnetic secular variation was lower than that of the recent field.

Seismic observations of regional gas blowdown at Teal South

A re-examination of seismic time-lapse data from the Teal South field provides support for a previously proposed model of regional pressure decline and the associated liberation of gas from nearby reservoirs due to the production from the only reservoir among them that is under production. The use of a specific attribute, instantaneous amplitude, and a series of time slices, however, provides increased detail in understanding fluid migration into or out of the reservoirs, and the path taken by pressure changes across faults. The regional decrease of pressure due to production in one reservoir has dramatic effects in nearby untapped reservoirs, one of which appears to exhibit evidence for the escape, and possible re-trapping nearby, of hydrocarbons from a spill point. The influx of water into the producing reservoir is also evidenced by a decrease in amplitude at one end of the oil-water contact.

DIFFERENTIATING PATH AND SOURCE PROCESSES IN COMPLEX GEOLOGIC STRUCTURE THROUGH PASSIVE SEISMIC STUDIES: BERING GLACIER, ALASKA AND VILLARRICA VOLCANO, CHILE

Measuring shallow seismic sources provides a way to reveal processes that cannot be directly observed, but the correct interpretation and value of these signals depend on the ability to distinguish source from propagation effects. Furthermore, seismic signals produced by a resonating source can look almost identical to those produced by impulsive sources, but modified along the path. Distinguishing these two phenomena can be accomplished by examining the wavefield with small aperture arrays or by recording seismicity near to the source when possible. We examine source and path effects in two different environments: Bering Glacier, Alaska and Villarrica Volcano, Chile. Using three 3-element seismic arrays near the terminus of the Bering Glacier, we have identified and located both terminus calving and iceberg breakup events. We show that automated array analysis provided a robust way to locate icequake events using P waves. This analysis also showed that arrivals within the long-period codas were incoherent within the small aperture arrays, demonstrating that these codas previously attributed to crack resonance were in fact a result of a complicated path rather than a source effect. At Villarrica Volcano, seismometers deployed from near the vent to ~10 km revealed that a several cycle long-period source signal recorded at the vent appeared elongated in the far-field. We used data collected from the stations nearest to the vent to invert for the repetitive seismic source, and found it corresponded to a shallow force within the lava lake oriented N75°E and dipping 7° from horizontal. We also used this repetitive signal to search the data for additional seismic and infrasonic properties which included calculating seismic-acoustic delay times, volcano acoustic-seismic ratios and energies, event frequency, and real-time seismic amplitude measurements. These calculations revealed lava lake level and activity fluctuations consistent with lava lake level changes inferred from the persistent infrasonic tremor.

TIME LAPSE SEISMIC OBSERVATIONS AND EFFECTS OF RESERVOIR COMPRESSIBILITY AT TEAL SOUTH OIL FIELD

One of the original ocean-bottom time-lapse seismic studies was performed at the Teal South oil field in the Gulf of Mexico during the late 1990’s. This work reexamines some aspects of previous work using modern analysis techniques to provide improved quantitative interpretations. Using three-dimensional volume visualization of legacy data and the two phases of post-production time-lapse data, I provide additional insight into the fluid migration pathways and the pressure communication between different reservoirs, separated by faults. This work supports a conclusion from previous studies that production from one reservoir caused regional pressure decline that in turn resulted in liberation of gas from multiple surrounding unproduced reservoirs. I also provide an explanation for unusual time-lapse changes in amplitude-versus-offset (AVO) data related to the compaction of the producing reservoir which, in turn, changed an isotropic medium to an anisotropic medium. In the first part of this work, I examine regional changes in seismic response due to the production of oil and gas from one reservoir. The previous studies primarily used two post-production ocean-bottom surveys (Phase I and Phase II), and not the legacy streamer data, due to the unavailability of legacy prestack data and very different acquisition parameters. In order to incorporate the legacy data in the present study, all three poststack data sets were cross-equalized and examined using instantaneous amplitude and energy volumes. This approach appears quite effective and helps to suppress changes unrelated to production while emphasizing those large-amplitude changes that are related to production in this noisy (by current standards) suite of data. I examine the multiple data sets first by using the instantaneous amplitude and energy attributes, and then also examine specific apparent time-lapse changes through direct comparisons of seismic traces. In so doing, I identify time-delays that, when corrected for, indicate water encroachment at the base of the producing reservoir. I also identify specific sites of leakage from various unproduced reservoirs, the result of regional pressure blowdown as explained in previous studies; those earlier studies, however, were unable to identify direct evidence of fluid movement. Of particular interest is the identification of one site where oil apparently leaked from one reservoir into a “new” reservoir that did not originally contain oil, but was ideally suited as a trap for fluids leaking from the neighboring spill-point. With continued pressure drop, oil in the new reservoir increased as more oil entered into the reservoir and expanded, liberating gas from solution. Because of the limited volume available for oil and gas in that temporary trap, oil and gas also escaped from it into the surrounding formation. I also note that some of the reservoirs demonstrate time-lapse changes only in the “gas cap” and not in the oil zone, even though gas must be coming out of solution everywhere in the reservoir. This is explained by interplay between pore-fluid modulus reduction by gas saturation decrease and dry-frame modulus increase by frame stiffening. In the second part of this work, I examine various rock-physics models in an attempt to quantitatively account for frame-stiffening that results from reduced pore-fluid pressure in the producing reservoir, searching for a model that would predict the unusual AVO features observed in the time-lapse prestack and stacked data at Teal South. While several rock-physics models are successful at predicting the time-lapse response for initial production, most fail to match the observations for continued production between Phase I and Phase II. Because the reservoir was initially overpressured and unconsolidated, reservoir compaction was likely significant, and is probably accomplished largely by uniaxial strain in the vertical direction; this implies that an anisotropic model may be required. Using Walton’s model for anisotropic unconsolidated sand, I successfully model the time-lapse changes for all phases of production. This observation may be of interest for application to other unconsolidated overpressured reservoirs under production.