A proposal for reservoir volume calculation in rainwater harvesting systems

This paper introduces a proposal for reservoir volume calculation in rainwater harvesting systems. The proposed method can be used for reservoir volume design in rainwater harvesting systems and is based on three important variables. These variables are water demand, system efficiency and repayment time. Several simulations were carried out in different scenarios considering typical values of both catchment area (for low-income and medium-income households) and water demand, with fixed water and tank costs. Results showed that the integrated analysis of demand, efficiency and repayment time may assist designers to determine a more adequate reservoir volume.

A bivariate return period based on copulas for hydrologic dam design: accounting for reservoir routing in risk estimation

A multivariate analysis on flood variables is needed to design some hydraulic structures like dams, as the complexity of the routing process in a reservoir requires a representation of the full hydrograph. In this work, a bivariate copula model was used to obtain the bivariate joint distribution of flood peak and volume, in order to know the probability of occurrence of a given inflow hydrograph. However, the risk of dam overtopping is given by the maximum water elevation reached during the routing process, which depends on the hydrograph variables, the reservoir volume and the spillway crest length. Consequently, an additional bivariate return period, the so-called routed return period, was defined in terms of risk of dam overtopping based on this maximum water elevation obtained after routing the inflow hydrographs. The theoretical return periods, which give the probability of occurrence of a hydrograph prior to accounting for the reservoir routing, were compared with the routed return period, as in both cases hydrographs with the same probability will draw a curve in the peak-volume space. The procedure was applied to the case study of the Santillana reservoir in Spain. Different reservoir volumes and spillway lengths were considered to investigate the influence of the dam and reservoir characteristics on the results. The methodology improves the estimation of the Design Flood Hydrograph and can be applied to assess the risk of dam overtopping

Avaliação da dinâmica temporal da evaporação no reservatório de Itumbiara, GO, utilizando dados obtidos por sensoriamento remoto

The object of this work was to study the dynamics of evaporation in the Itumbiara reservoir, located in Central Brazil, using MODIS-derived water surface temperature (product MOD11A1) and meteorological data acquired over the water surface. The evaporation rates were derived from latent heat flux, estimated through a mass transfer model. The estimates were carried out for the period between 1/1/2010 and 31/12/2010. The results showed that evaporation rate tends to increase from January to September and then decrease from September to December. The evaporation rate reached values near 20 mm day-1 in Itumbiara reservoir during the dry season in 2010. The mean evaporation rate for the wet season was 3.66 mm day-1 and 8.25 mm day-1 for the dry season. The total water volume evaporated from Itumbiara reservoir during 2010 was estimated at about 1.7 billion m³ (2,300 mm) which represents 10% of total reservoir volume. The results suggest that advection is the main transport mechanism which drives the evaporation in Itumbiara. The convective processes contribute secondarily to evaporation in Itumbiara reservoir.

Desarrollo de metodologías para mejorar la estimación de los hidrogramas de diseño para el cálculo de los órganos de desagüe de las presas

En la presente Tesis se ha llevado a cabo el contraste y desarrollo de metodologías que permitan mejorar el cálculo de las avenidas de proyecto y extrema empleadas en el cálculo de la seguridad hidrológica de las presas. En primer lugar se ha abordado el tema del cálculo de las leyes de frecuencia de caudales máximos y su extrapolación a altos periodos de retorno. Esta cuestión es de gran relevancia, ya que la adopción de estándares de seguridad hidrológica para las presas cada vez más exigentes, implica la utilización de periodos de retorno de diseño muy elevados cuya estimación conlleva una gran incertidumbre. Es importante, en consecuencia incorporar al cálculo de los caudales de diseño todas la técnicas disponibles para reducir dicha incertidumbre. Asimismo, es importante hacer una buena selección del modelo estadístico (función de distribución y procedimiento de ajuste) de tal forma que se garantice tanto su capacidad para describir el comportamiento de la muestra, como para predecir de manera robusta los cuantiles de alto periodo de retorno. De esta forma, se han realizado estudios a escala nacional con el objetivo de determinar el esquema de regionalización que ofrece mejores resultados para las características hidrológicas de las cuencas españolas, respecto a los caudales máximos anuales, teniendo en cuenta el numero de datos disponibles. La metodología utilizada parte de la identificación de regiones homogéneas, cuyos límites se han determinado teniendo en cuenta las características fisiográficas y climáticas de las cuencas, y la variabilidad de sus estadísticos, comprobando posteriormente su homogeneidad. A continuación, se ha seleccionado el modelo estadístico de caudales máximos anuales con un mejor comportamiento en las distintas zonas de la España peninsular, tanto para describir los datos de la muestra como para extrapolar a los periodos de retorno más altos. El proceso de selección se ha basado, entre otras cosas, en la generación sintética de series de datos mediante simulaciones de Monte Carlo, y el análisis estadístico del conjunto de resultados obtenido a partir del ajuste de funciones de distribución a estas series bajo distintas hipótesis. Posteriormente, se ha abordado el tema de la relación caudal-volumen y la definición de los hidrogramas de diseño en base a la misma, cuestión que puede ser de gran importancia en el caso de presas con grandes volúmenes de embalse. Sin embargo, los procedimientos de cálculo hidrológico aplicados habitualmente no tienen en cuenta la dependencia estadística entre ambas variables. En esta Tesis se ha desarrollado un procedimiento para caracterizar dicha dependencia estadística de una manera sencilla y robusta, representando la función de distribución conjunta del caudal punta y el volumen en base a la función de distribución marginal del caudal punta y la función de distribución condicionada del volumen respecto al caudal. Esta última se determina mediante una función de distribución log-normal, aplicando un procedimiento de ajuste regional. Se propone su aplicación práctica a través de un procedimiento de cálculo probabilístico basado en la generación estocástica de un número elevado de hidrogramas. La aplicación a la seguridad hidrológica de las presas de este procedimiento requiere interpretar correctamente el concepto de periodo de retorno aplicado a variables hidrológicas bivariadas. Para ello, se realiza una propuesta de interpretación de dicho concepto. El periodo de retorno se entiende como el inverso de la probabilidad de superar un determinado nivel de embalse. Al relacionar este periodo de retorno con las variables hidrológicas, el hidrograma de diseño de la presa deja de ser un único hidrograma para convertirse en una familia de hidrogramas que generan un mismo nivel máximo en el embalse, representados mediante una curva en el plano caudal volumen. Esta familia de hidrogramas de diseño depende de la propia presa a diseñar, variando las curvas caudal-volumen en función, por ejemplo, del volumen de embalse o la longitud del aliviadero. El procedimiento propuesto se ilustra mediante su aplicación a dos casos de estudio. Finalmente, se ha abordado el tema del cálculo de las avenidas estacionales, cuestión fundamental a la hora de establecer la explotación de la presa, y que puede serlo también para estudiar la seguridad hidrológica de presas existentes. Sin embargo, el cálculo de estas avenidas es complejo y no está del todo claro hoy en día, y los procedimientos de cálculo habitualmente utilizados pueden presentar ciertos problemas. El cálculo en base al método estadístico de series parciales, o de máximos sobre un umbral, puede ser una alternativa válida que permite resolver esos problemas en aquellos casos en que la generación de las avenidas en las distintas estaciones se deba a un mismo tipo de evento. Se ha realizado un estudio con objeto de verificar si es adecuada en España la hipótesis de homogeneidad estadística de los datos de caudal de avenida correspondientes a distintas estaciones del año. Asimismo, se han analizado los periodos estacionales para los que es más apropiado realizar el estudio, cuestión de gran relevancia para garantizar que los resultados sean correctos, y se ha desarrollado un procedimiento sencillo para determinar el umbral de selección de los datos de tal manera que se garantice su independencia, una de las principales dificultades en la aplicación práctica de la técnica de las series parciales. Por otra parte, la aplicación practica de las leyes de frecuencia estacionales requiere interpretar correctamente el concepto de periodo de retorno para el caso estacional. Se propone un criterio para determinar los periodos de retorno estacionales de forma coherente con el periodo de retorno anual y con una distribución adecuada de la probabilidad entre las distintas estaciones. Por último, se expone un procedimiento para el cálculo de los caudales estacionales, ilustrándolo mediante su aplicación a un caso de estudio. The compare and develop of a methodology in order to improve the extreme flow estimation for dam hydrologic security has been developed. First, the work has been focused on the adjustment of maximum peak flows distribution functions from which to extrapolate values for high return periods. This has become a major issue as the adoption of stricter standards on dam hydrologic security involves estimation of high design return periods which entails great uncertainty. Accordingly, it is important to incorporate all available techniques for the estimation of design peak flows in order to reduce this uncertainty. Selection of the statistical model (distribution function and adjustment method) is also important since its ability to describe the sample and to make solid predictions for high return periods quantiles must be guaranteed. In order to provide practical application of previous methodologies, studies have been developed on a national scale with the aim of determining a regionalization scheme which features best results in terms of annual maximum peak flows for hydrologic characteristics of Spanish basins taking into account the length of available data. Applied methodology starts with the delimitation of regions taking into account basin’s physiographic and climatic characteristics and the variability of their statistical properties, and continues with their homogeneity testing. Then, a statistical model for maximum annual peak flows is selected with the best behaviour for the different regions in peninsular Spain in terms of describing sample data and making solid predictions for high return periods. This selection has been based, among others, on synthetic data series generation using Monte Carlo simulations and statistical analysis of results from distribution functions adjustment following different hypothesis. Secondly, the work has been focused on the analysis of the relationship between peak flow and volume and how to define design flood hydrographs based on this relationship which can be highly important for large volume reservoirs. However, commonly used hydrologic procedures do not take statistical dependence between these variables into account. A simple and sound method for statistical dependence characterization has been developed by the representation of a joint distribution function of maximum peak flow and volume which is based on marginal distribution function of peak flow and conditional distribution function of volume for a given peak flow. The last one is determined by a regional adjustment procedure of a log-normal distribution function. Practical application is proposed by a probabilistic estimation procedure based on stochastic generation of a large number of hydrographs. The use of this procedure for dam hydrologic security requires a proper interpretation of the return period concept applied to bivariate hydrologic data. A standard is proposed in which it is understood as the inverse of the probability of exceeding a determined reservoir level. When relating return period and hydrological variables the only design flood hydrograph changes into a family of hydrographs which generate the same maximum reservoir level and that are represented by a curve in the peak flow-volume two-dimensional space. This family of design flood hydrographs depends on the dam characteristics as for example reservoir volume or spillway length. Two study cases illustrate the application of the developed methodology. Finally, the work has been focused on the calculation of seasonal floods which are essential when determining the reservoir operation and which can be also fundamental in terms of analysing the hydrologic security of existing reservoirs. However, seasonal flood calculation is complex and nowadays it is not totally clear. Calculation procedures commonly used may present certain problems. Statistical partial duration series, or peaks over threshold method, can be an alternative approach for their calculation that allow to solve problems encountered when the same type of event is responsible of floods in different seasons. A study has been developed to verify the hypothesis of statistical homogeneity of peak flows for different seasons in Spain. Appropriate seasonal periods have been analyzed which is highly relevant to guarantee correct results. In addition, a simple procedure has been defined to determine data selection threshold on a way that ensures its independency which is one of the main difficulties in practical application of partial series. Moreover, practical application of seasonal frequency laws requires a correct interpretation of the concept of seasonal return period. A standard is proposed in order to determine seasonal return periods coherently with the annual return period and with an adequate seasonal probability distribution. Finally a methodology is proposed to calculate seasonal peak flows. A study case illustrates the application of the proposed methodology.

Summary of the Snowmastodon Project Special Volume: A high-elevation, multi-proxy biotic and environmental record of MIS 6-4 from the Ziegler Reservoir fossil site, Snowmass Village, Colorado, USA

In North America, terrestrial records of biodiversity and climate change that span Marine Oxygen Isotope Stage (MIS) 5 are rare. Where found, they provide insight into how the coupling of the ocean-atmosphere system is manifested in biotic and environmental records and how the biosphere responds to climate change. In 2010-2011, construction at Ziegler Reservoir near Snowmass Village, Colorado (USA) revealed a nearly continuous, lacustrine/wetland sedimentary sequence that preserved evidence of past plant communities between similar to 140 and 55 lea, including all of MIS 5. At an elevation of 2705 m, the Ziegler Reservoir fossil site also contained thousands of well-preserved bones of late Pleistocene megafauna, including mastodons, mammoths, ground sloths, horses, camels, deer, bison, black bear, coyotes, and bighorn sheep. In addition, the site contained more than 26,000 bones from at least 30 species of small animals including salamanders, otters, muskrats, minks, rabbits, beavers, frogs, lizards, snakes, fish, and birds. The combination of macro- and micro-vertebrates, invertebrates, terrestrial and aquatic plant macrofossils, a detailed pollen record, and a robust, directly dated stratigraphic framework shows that high-elevation ecosystems in the Rocky Mountains of Colorado are climatically sensitive and varied dramatically throughout MIS 5 

The Cougar Point tuff: Implications for thermochemical zonation and longevity of high-temperature, large-volume silicic magmas of the Miocene Yellowstone hotspot

The 12.7-10.5 Ma Cougar Point Tuff in southern Idaho, USA, consists of 10 large-volume (>10²-10³ km³ each), high-temperature (800-1000 °C), rhyolitic ash-flow tuffs erupted from the Bruneau-Jarbidge volcanic center of the Yellowstone hotspot. These tuffs provide evidence for compositional and thermal zonation in pre-eruptive rhyolite magma, and suggest the presence of a long-lived reservoir that was tapped by numerous large explosive eruptions. Pyroxene compositions exhibit discrete compositional modes with respect to Fe and Mg that define a linear spectrum punctuated by conspicuous gaps. Airfall glass compositions also cluster into modes, and the presence of multiple modes indicates tapping of different magma volumes during early phases of eruption. Equilibrium assemblages of pigeonite and augite are used to reconstruct compositional and thermal gradients in the pre-eruptive reservoir. The recurrence of identical compositional modes and of mineral pairs equilibrated at high temperatures in successive eruptive units is consistent with the persistence of their respective liquids in the magma reservoir. Recurrence intervals of identical modes range from 0.3 to 0.9 Myr and suggest possible magma residence times of similar duration. Eruption ages, magma temperatures, Nd isotopes, and pyroxene and glass compositions are consistent with a long-lived, dynamically evolving magma reservoir that was chemically and thermally zoned and composed of multiple discrete magma volumes.

A modelling approach to the development of an active management strategy for the Queen Elizabeth II reservoir

This article outlines the outcome of work that set out to provide one of the specified integral contributions to the overarching objectives of the EU- sponsored LIFE98 project described in this volume. Among others, these included a requirement to marry automatic monitoring and dynamic modelling approaches in the interests of securing better management of water quality in lakes and reservoirs. The particular task given to us was to devise the elements of an active management strategy for the Queen Elizabeth II Reservoir. This is one of the larger reservoirs supplying the population of the London area: after purification and disinfection, its water goes directly to the distribution network and to the consumers. The quality of the water in the reservoir is of primary concern, for the greater is the content of biogenic materials, including phytoplankton, then the more prolonged is the purification and the more expensive is the treatment. Whatever good that phytoplankton may do by way of oxygenation and oxidative purification, it is eventually relegated to an impurity that has to be removed from the final product. Indeed, it has been estimated that the cost of removing algae and microorganisms from water represents about one quarter of its price at the tap. In chemically fertile waters, such as those typifying the resources of the Thames Valley, there is thus a powerful and ongoing incentive to be able to minimise plankton growth in storage reservoirs. Indeed, the Thames Water company and its predecessor undertakings, have a long and impressive history of confronting and quantifying the fundamentals of phytoplankton growth in their reservoirs and of developing strategies for operation and design to combat them. The work to be described here follows in this tradition. However, the use of the model PROTECH-D to investigate present phytoplankton growth patterns in the Queen Elizabeth II Reservoir questioned the interpretation of some of the recent observations. On the other hand, it has reinforced the theories underpinning the original design of this and those Thames-Valley storage reservoirs constructed subsequently. The authors recount these experiences as an example of how simulation models can hone the theoretical base and its application to the practical problems of supplying water of good quality at economic cost, before the engineering is initiated.

Simulations of phytoplankton dynamics in El Gergal Reservoir, Southern Spain (PROTECH)

The Mediterranean region is characterised by a variable climate with most of the rain falling during the winter and frequent summer droughts. Such warm, dry periods are ideal for the growth of large algal blooms that often consist of potentially toxic Cyanobacteria. This makes the management of water for human use particularly challenging in such a climate and it is important to understand how such blooms can be avoided or at least be reduced in size. PROTECH (Phytoplankton RespOnses To Environmental CHange) is a model that simulates the dynamics of different species of phytoplankton populations in lakes and reservoirs. Its distinct advantage over similar models is its ability to simulate the relative composition of the algal flora, allowing both quantitative and qualitative conclusions to be drawn e.g. whether Cyanobacteria could be a potential problem. PROTECH has been applied primarily to lakes and reservoirs in northern Europe. Recently, however, the model has been applied to water bodies in lower latitudes, including Australia to a water supply reservoir in the south of Spain, El Gergal. El Gergal is the last in a chain of reservoirs that supply water to the city of Seville. It was brought into service in April 1979 and has a maximum storage volume of 35 000 000 m3. This article summarises the application of PROTECH in order to simulate the following problems: • the effect of a large influx of Ceratium biomass into El Gergal from another reservoir • the effect of using alternative water sources instead of the Guadalquivir River (used occasionally to raise water levels in El Gergal) • the effect of installing tertiary sewage treatment on the Cala River • the effect of simulated drought conditions on phytoplankton in the reservoir.

Hydrological studies of Matatilla Reservoir, Uttar Pradesh

Matatilla Reservoir, located in semi arid region, (Lat. 25 degree 15'N and Long. 78 degree 23'E) has an area (at FRL) of 13,893 ha, volume and shore development 0.663 and 1.65, shoreline 73.6 km. Volume and shore development indicate that greater part of the reservoir is shallow, which is a favourable point for fish productivity. Temperature and dissolved oxygen gradually decreased with the increase in depth. Carbon dioxide was absent from the surface but invariably present in the bottom (3.6 ppm) pH remained alkaline (7.2-8.4 ppm) throughout the year. Alkalinity, chloride, calcium, magnesium, hardness and priductivity was maximum in pre-monsoon and minimum in monsoon except for calcium and manganesium in post-monsoon. Phosphate, nitrogen and ammonical nitrogen were found in traces. These variations may be due to influx and outflow of water and use of reservoir water for multipurpose activities.

The morphometrics and the proximate composition of the edible fresh water mussel Lamellidens lamellatus (Lea) from Bathalagoda tank, a man made reservoir in Sri Lanka

Morphometric characteristics of Lamellidens lamellatus were studied using a random sample of 138 individuals collected from Bathalagoda reservoir in Kurunegala district. The largest number of individuals in the collection was between 40-50 mm in length. They had a body weight between 60-80 gm. The most abundant bio-chemical component in the adductor and in the foot was protein. Carbohydrate and lipid quantities were almost equal while ash showed the least value. The total protein content ranged from 48.8% to 73.4%; the carbohydrate content from 12.64% to 23.8%; and lipid content from 12.0% to 28.1% of their dry weight. It was found that there were significant relationships between the length and the body weight; and between length and volume.

Studies on morphometry and hydrology of Gandhisagar Reservoir with special reference to its fisheries

Gandhisagar, the second largest reservoir of India is located in Mandsaur district of Madhya Pradesh at latitude 24°44'N and longitude 75°33'E at an altitude of 403.56m MSL in orientation from NE to SE. It has an extensive water spread area of 66000 ha at full reservoir level with a maximum and mean depth of 49.52 and 11.73 m respectively. The maximum length and width of the reservoir are 112 and 16km having a total shore line of 442km. Details of catchment area, bathymetry, standard hydrological data giving water level relation of the basin to water spread area, volume and fish production and the bottom topographical details of 11 experimental fishing stations and 6 fish landing centres are discussed.

How Porosity and Permeability Vary Spatially With Grain Size, Sorting, Cement Volume, and Mineral Dissolution In Fluvial Triassic Sandstones: The Value of Geostatistics and Local Regression

Although it is well known that sandstone porosity and permeability are controlled by a range of parameters such as grain size and sorting, amount, type, and location of diagenetic cements, extent and type of compaction, and the generation of intergranular and intragranular secondary porosity, it is less constrained how these controlling parameters link up in rock volumes (within and between beds) and how they spatially interact to determine porosity and permeability. To address these unknowns, this study examined Triassic fluvial sandstone outcrops from the UK using field logging, probe permeametry of 200 points, and sampling at 100 points on a gridded rock surface. These field observations were supplemented by laser particle-size analysis, thin-section point-count analysis of primary and diagenetic mineralogy, quantitiative XRD mineral analysis, and SEM/EDAX analysis of all 100 samples. These data were analyzed using global regression, variography, kriging, conditional simulation, and geographically weighted regression to examine the spatial relationships between porosity and permeability and their potential controls. The results of bivariate analysis (global regression) of the entire outcrop dataset indicate only a weak correlation between both permeability porosity and their diagenetic and depositional controls and provide very limited information on the role of primary textural structures such as grain size and sorting. Subdividing the dataset further by bedding unit revealed details of more local controls on porosity and permeability. An alternative geostatistical approach combined with a local modelling technique (geographically weighted regression; GWR) subsequently was used to examine the spatial variability of porosity and permeability and their controls. The use of GWR does not require prior knowledge of divisions between bedding units, but the results from GWR broadly concur with results of regression analysis by bedding unit and provide much greater clarity of how porosity and permeability and their controls vary laterally and vertically. The close relationship between depositional lithofacies in each bed, diagenesis, and permeability, porosity demonstrates that each influences the other, and in turn how understanding of reservoir properties is enhanced by integration of paleoenvironmental reconstruction, stratigraphy, mineralogy, and geostatistics.

Diet and feeding of fish from Grande River, located below the Volta Grande Reservoir, MG-SP

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Load models support sustainable: planning for Brazil s reservoir aquaculture

Cage aquaculture in hydroelectric reservoirs has great potential for expansion in Brazil, but there are concerns of negative environmental impacts. The environmental sustainability of cage culture depends on hamornization between farming practices and the hydrological peculiarities of the site. Mass balance modeling can estimate the amounts of nutrients that can be loaded without triggering eutrophication and resulting maximum allowable production volume. Careful climate zoning can also assist proper siting.

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.