Search NASA⌕ Search

SEARCH · Search NASA

Results for “Reservoir Modeling”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
NASA NTRS2018

Soil Moisture Active Passive (SMAP) Project Assessment Report for Version 4 of the L4_SM Data Product

This report provides an assessment of Version 4 of the SMAP Level 4 Surface and Root Zone Soil Moisture (L4_SM) product, released on 14 June 2018. The assessment includes comparisons of L4_SM soil moisture and temperature estimates with in situ measurements from core validation sites and sparse networks. The assessment further includes a global evaluation of the internal diagnostics from the ensemble-based data assimilation system that is used to generate the L4_SM product, including observation-minus-forecast (O-F) brightness temperature residuals and soil moisture analysis increments.Together, the core validation site comparisons and the statistics of the assimilation diagnostics areconsidered primary validation methodologies for the L4_SM product. Comparisons against in situ measurements from regional-scale sparse networks are considered a secondary validation methodology because such in situ measurements are subject to upscaling errors from the point-scale to the grid-cell scale of the data product.The Version 4 L4_SM product benefits from an improved land surface modeling system and from retrospective surface meteorological forcing data that are as consistent as possible with the present-day datain terms of their climatology. Specifically, the model changes include revised parameters and parameterizations for (i) the surface energy balance, (ii) recharge from below of the model's surface excess reservoir, and (iii) the snow depletion curve. Updated ancillary inputs include improved datasets for landcover, topography, and vegetation height. The Version 4 algorithm further includes a revised approach to precipitation corrections that improves the precipitation climatology in Africa and the high-latitudes. Moreover, for system calibration the model is forced retrospectively with MERRA-2 reanalysis data, which are more consistent with the near-real time GEOS forward processing (FP) data used during the SMAP period than the retrospective GEOS data that were available for previous L4_SM versions. An analysis of the time-average surface and root zone soil moisture shows that the global pattern ofarid and humid regions is captured by the Version 4 L4_SM estimates. Owing to the changes in the landsurface modeling system, surface soil moisture is typically drier by several volumetric percent in Version 4 compared to Version 3, whereas root zone soil moisture is wetter in Version 4 in some regions and drierin others. Because of these climatological differences, the Version 3 and Version 4 products should not be combined into a single dataset for use in applications.Results from the core validation site comparisons indicate that Version 4 of the L4_SM data product meets the self-imposed L4_SM accuracy requirement, which is formulated in terms of the RMSE after removal of the long-term mean difference (ubRMSE). The overall ubRMSE of the 3-hourly L4_SM dataat the 9 km scale is 0.039 m3 m-3 for surface soil moisture and 0.029 m3 m-3 for root zone soil moisture,below the 0.04 m3 m-3 requirement. The L4_SM estimates are an improvement over estimates from a model-only Nature Run version 7.2 (NRv7.2), which demonstrates the beneficial impact of the SMAP brightness temperature data. Overall, L4_SM surface and root zone soil moisture estimates are more skillful than NRv7.2 estimates, with statistically significant improvements at the 5% level for surface soil moisture R and anomaly R values. Results from comparisons of the L4_SM product to i

Reichle, Rolf H.↗
NASA NTRS2012

The Volatile Fraction of Comets as Quantified at Infrared Wavelengths - An Emerging Taxonomy and Implications for Natal Heritage

It is relatively easy to identify the reservoir from which a given comet was ejected. But dynamical models demonstrate that the main cometary reservoirs (Kuiper Belt, Oort Cloud) each contain icy bodies that formed in a range of environments in the protoplanetary disk, and the Oort Cloud may even contain bodies that formed in disks of sibling stars in the Sun s birth cluster. The cometary nucleus contains clues to the formative region(s) of its individual components. The composition of ices and rocky grains reflect a range of processes experienced by material while on the journey from the natal interstellar cloud core to the cometary nucleus. For that reason, emphasis is placed on classifying comets according to their native ices and dust (rather than orbital dynamics). Mumma & Charnley [1] reviewed the current status of taxonomies for comets and relation to their natal heritage.

Mumma, M. J.↗
NASA NTRS2024

Satellite-Based Tracking of Reservoir Operations for Flood Management During the 2018 Extreme Weather Event in Kerala, India

Uncoordinated management of hydropower dams during extreme and unexpected precipitation events in mountainous terrain can have disastrous consequences due to the competing nature of flood control and hydropower generation. Numerous cases of flooding events that have been exacerbated due to insufficient storage conditions in hydropower dams have been reported worldwide. There is a need for a scalable and publicly accessible monitoring framework that is capable of providing reliable, near-real time, and transparent reservoir operations data. A fully satellite-based framework is the most viable solution to build such capability. The Reservoir Assessment Tool (RAT 3.0), which utilizes high frequency remote sensing-based surface area and reservoir storage estimation alongside hydrological modelled inflow was applied here for the 2018 Kerala floods in India as a globally representative case for a mountainous river basin with high precipitation and hydropower dams. Application of satellite-based RAT 3.0 in monitoring the state of 19 reservoirs in Kerala during the flood event showed very promising results. In general, RAT 3.0, using satellite remote sensing, was found to be able to capture the temporal trend of the reservoir storage and pinpoint the sudden shift in filling or release decisions made by the dam operator. Inflow modelling in such regions was found to require careful calibration with identification of reservoirs that are heavily regulated being a critical aspect. The perennial high cloud cover in such regions necessitate and highlights the central role played by microwave and radar-based satellite sensors, such as the Surface Water and Ocean Topography (SWOT) mission, in tracking reservoir state. An operational version of RAT 3.0 for stakeholder agencies tailored for hydropower dams operating in high precipitation and mountainous environments is a real-world outcome of this study.

Sarath Suresh↗
NASA NTRS1977

Fluid-line math model

Simplified mathematical model simulates large hydraulic systems on either analog or digital computers. Models of pumps, servoactuators, reservoirs, accumulators, and valves are connected generating systems containing six hundred elements.

Kandelman, A.↗
NASA NTRS1981

Calibration of a turbidity meter for making estimates of total suspended solids concentrations and beam attenuation coefficients in field experiments

Management of water resources such as a reservoir requires using analytical models which describe such parameters as the suspended sediment field. To select or develop an appropriate model requires making many measurements to describe the distribution of this parameter in the water column. One potential method for making those measurements expeditiously is to measure light transmission or turbidity and relate that parameter to total suspended solids concentrations. An instrument which may be used for this purpose was calibrated by generating curves of transmission measurements plotted against measured values of total suspended solids concentrations and beam attenuation coefficients. Results of these experiments indicate that field measurements made with this instrument using curves generated in this study should correlate with total suspended solids concentrations and beam attenuation coefficients in the water column within 20 percent.

Usry, J. W.↗
NASA NTRS1987

Volatile reservoirs below the surface of the Elysium region of Mars: Geomorphic evidence

The Elysium volcanic province contains a variety of geomorphic evidence for the existence of large volatile reservoirs of subsurface volatiles. Study of these landforms yields insight into the distribution and size of these reservoirs and how they interact with the surface environment and will ultimately place constraints on the geometry, constitution, origin, time of formation, and temporal evolution of these important components of the Martian crust. Three principal types of landforms appear to be related to subsurface volatile reservoirs in the Elysium region of Mars: small outflow channels; large lahars; and vast expanses of knobby terranes around the margins of the Elysium dome. The evidence provided by these landforms is internally consistent with the presence of a large relatively shallow volatile reservoir in the Elysium region. If the geologic features described are reliable indicators of subsurface volatiles, they imply that: volatile reservoirs lie relatively close to the surface and underlie millions of sq km in this region; there is no apparent latitudinal variation in the depth or thickness of the volatile reservoirs; the precursors of the knobby terranes are or were important volatile reservoirs; volatiles may be lost in a variety of ways from these reservoirs; and volatiles were incorporated in an easily eroded surficial deposit in the middle history of Mars. The ultimate origin of water in this reservoir is uncertain. A model to explain the preferential entrapment of volatiles into the region's surface materials may be required.

Christiansen, Eric H.↗
NASA NTRS2015

Majorite-Garnet Partitioning of the Highly Siderophile Elements: New Results and Application to Mars

HSE and Os isotopes are used to constrain processes such as accretion, mantle evolution, crustal recycling, and core-mantle mixing, and to constrain the timing and depth of differentiation of Mars. Although showed that the HSE contents of the martian mantle could have been established by metal-silicate equilibrium in early Mars, the role of a cooling magma ocean and associated crystallization in further fractionating the HSEs is unclear. Garnet is thought to have played an important role in controlling trace element concentrations in the martian mantle reservoirs. However, testing these models, including Os isotopes, has been hindered by a dearth of partitioning data for the HSE in deep mantle phases - majorite, wadsleyite, ringwoodite, akimotoite - that may be present in the martian mantle. We examine the partitioning behavior of HSEs between majorite garnet (gt), olivine (oliv), and silicate liquid (melt).

Danielson, L. R.↗
NASA NTRS1989

A simple, mass balance model of carbon flow in a controlled ecological life support system

Internal cycling of chemical elements is a fundamental aspect of a Controlled Ecological Life Support System (CELSS). Mathematical models are useful tools for evaluating fluxes and reservoirs of elements associated with potential CELSS configurations. A simple mass balance model of carbon flow in CELSS was developed based on data from the CELSS Breadboard project at Kennedy Space Center. All carbon reservoirs and fluxes were calculated based on steady state conditions and modelled using linear, donor-controlled transfer coefficients. The linear expression of photosynthetic flux was replaced with Michaelis-Menten kinetics based on dynamical analysis of the model which found that the latter produced more adequate model output. Sensitivity analysis of the model indicated that accurate determination of the maximum rate of gross primary production is critical to the development of an accurate model of carbon flow. Atmospheric carbon dioxide was particularly sensitive to changes in photosynthetic rate. The small reservoir of CO2 relative to large CO2 fluxes increases the potential for volatility in CO2 concentration. Feedback control mechanisms regulating CO2 concentration will probably be necessary in a CELSS to reduce this system instability.

Garland, Jay L.↗
NASA NTRS2018

Improving the SMAP Level-4 Soil Moisture Product

The NASA Soil Moisture Active Passive (SMAP) mission generates, among other data sets, the Level 4 Soil Moisture (L4 SM) product. The L4 SM algorithm is based on the assimilation of SMAP radiometer brightness temperature observations into the NASA Catchment land surface model using a spatially distributed ensemble Kalman filter. The L4 SM data are published with a mean latency of approx. 2.5 days from the time of observation and provide global, three-hourly, 9 km resolution estimates of surface and root-zone soil moisture and related land surface states and fluxes. In 2018, the product was upgraded from Version 3 to Version 4. Underlying the new version is a revised modeling system that includes improved input parameter datasets for land cover, topography, and vegetation height that are based on recent, high quality space-borne remote sensing observations. Land cover inputs were updated to the GlobCover2009 product, which is based on satellite observations from the Medium Resolution Imaging Spectrometer. Topographic statistics now rely on observations from the Shuttle Radar Topography Mission. Finally, vegetation height inputs are derived from space-borne Lidar measurements. Additionally, SMAP Level-2 soil moisture retrievals and in situ soil moisture measurements were used to calibrate a particular Catchment model parameter that governs the recharge of soil moisture from the models root-zone excess reservoir into the surface excess reservoir. Specifically, the replenishment of soil moisture near the surface from below under non-equilibrium conditions was substantially reduced, which brings the models surface soil moisture more in line with the SMAP Level 2 and in situ soil moisture. Finally, the calibration of the assimilated SMAP brightness temperatures changed substantially from Version 3 to Version 4. Considerable effort went into the version upgrade, creating an expectation that the new version is improved over the old version. Indeed, some aspects of the new version are clearly better. However, other aspects are not, and on balance, the overall improvement is modest at best. In this presentation we summarize the skill of the new and old versions vs. independent in situ measurements and in terms of data assimilation diagnostics, including, for example, the statistics of the (soil moisture) analysis increments and the observation minus forecast (brightness temperatures) residuals. We share our experience with trying to improve to the L4 SM product and the lessons learned from the effort.

Reichle, Rolf↗
NASA NTRS2019

Quantifying Water Storage Change and Land Subsidence Induced by Reservoir Impoundment Using GRACE, Landsat, and GPS Data

The construction of hydropower dams is a common strategy to support a country's increasing need for electricityand river water management for industry and agriculture. Although the hydrological and geophysical impacts ofwater relocation are usually assessed prior to impoundment, their accuracy is generally limited due to the lack ofin situ observations, especially in a remote area. This study presents a workflow to quantify the terrestrial waterstorage change (TWS) and land subsidence induced by a reservoir's water impoundment using multiple satelliteobservations (GRACE, Landsat), land surface models (CABLE, GLDAS, NCEP, ECMWF), and GPS data. The studysite is the Bakun Dam, located in Sarawak, Malaysia, which is the largest hydropower dam in Southeast Asia.Commencing operation in late 2010, the dam induced a change of water mass and lake surface area that wasclearly observed by GRACE and Landsat observations, respectively. During the 17-month impounding period(from August 2010 to December 2011), GRACE observed a dramatic increase of approximately 200mmequivalent water height, while Landsat detected an increased lake extent of around 600 km2. In this paper, aforward model is developed to determine the increased water surface level corresponding to GRACE observations,estimated to be about 120 m. In contrast to GRACE, the TWS derived from land surface models cannotcapture the increased TWS, due to the lack of reservoir routing algorithms in the models. In addition, the landsubsidence was calculated using the disk load model constructed based on the GRACE-derived lake level andLandsat-derived lake extent; the result is validated with the GPS data from BIN1 station, located at the westerncoast of Borneo. The commencement stage of the Bakun Dam induces the large-scale land subsidence, whichcauses the GPS-BIN1 station to subside by ~9 mm, and move toward the Bakun Lake by ~4 mm. Computation ofthe surface displacements directly from GRACE spherical harmonic coefficient data fails to capture the subsidencefeature, mainly due to the truncation error. Overall, this study demonstrates that evaluating GRACE inconjunction with Landsat, LSMs, and GPS data allows the exploitation of the gravity signal at a much smallerspatial scale than its intrinsic resolution. Benefiting from global coverage, the newly developed satellite-basedalgorithm is a valuable tool for assessing the impacts of reservoir operation on hydrological and geophysicalchanges from local to regional scales.

Tangdamrongsub, Natthachet↗
NASA NTRS2026

Synchronized Eruptions on Io: Possible Evidence of Interconnected Subsurface Magma Reservoirs

On 27 December 2024, Juno's JIRAM (Jovian InfraRed Auroral Mapper) instrument observed an unprecedented volcanic event in Io's southern hemisphere, covering a vast region of ∼65,000 km2, near 73°S, 140°E. Within the imaged region, only one hot spot was previously known (Pfd454). This feature was earlier estimated to cover an area of 300 km2 with a total power output of 34 GW. JIRAM results show that the region produces a power output of 140–260 TW, over 1,000 times higher than earlier estimates and likely exceeding the brightest eruption ever recorded on Io, that of Surt in 2001 (∼80 TW). Three adjacent hot spots also exhibited dramatic power increases: P139, PV18, and an unnamed feature south of the main one that surged to ∼1 TW, placing all of them among the top 10 most powerful hot spots observed on Io. A temperature analysis of the features supports the simultaneous onset of these brightenings and suggests a single eruptive event propagating beneath the surface across hundreds of kilometers; this is the first time this has been observed on Io. This in turn would imply a connection among the hotspots' magma reservoirs, while other nearby hotspots that have been known to be active in the recent past, such as Kurdalagon Patera, appear unaffected. The simultaneity supports models of massive, interconnected magma reservoirs. The topology of these regional magma systems may resemble that of a large-scale sponge, in which the massive reservoirs are the pores, interconnected through a largely solid outer shell.

A Mura↗
NASA NTRS2002

Re-Os Isotopic Constraints on the Chemical Evolution and Differentiation of the Martian Mantle

The (187)Re-187Os isotopic systematics of SNC meteorites, thought to be from Mars, provide valuable information regarding the chemical processes that affected the Martian mantle, particularly with regard to the relative abundances of highly siderophile elements (HSE). Previously published data (Birck and Allegre 1994, Brandon et al. 2000), and new data obtained since these studies, indicate that the HSE and Os isotopic composition of the Martian mantle was primarily set in its earliest differentiation history. If so, then these meteorites provide key constraints on the processes that lead to variation in HSE observed in not only Mars, but also Earth, the Moon and other rocky bodies in the Solar System. Processes that likely have an effect on the HSE budgets of terrestrial mantles include core formation, magma ocean crystallization, development of juvenile crust, and the addition of a late veneer. Each of these processes will result in different HSE variation and the isotopic composition of mantle materials and mantle derived lavas. Two observations on the SNC data to present provide a framework for which to test the importance of each of these processes. First, the concentrations of Re and Os in SNC meteorites indicate that they are derived from a mantle that has similar concentrations to the Earth's mantle. Such an observation is consistent with a model where a chondritic late veneer replenished the Earth and Martian mantles subsequent to core formation on each planet. Alternative models to explain this observation do exist, but will require additional data to test the limitations of each. Second, Re-Os isotopic results from Brandon et al. (2000) and new data presented here, show that initial yos correlates with variations in the short-lived systems of (182)Hf- (182)W and (142)Sm-142Nd in the SNC meteorites (epsilon(sub W) and epsilon(sub 142Nd)). These systematics require an isolation of mantle reservoirs during the earliest differentiation history of Mars, and subsequent inefficient mixing between these reservoirs. These data show that models for the origin of isotopic variation for SNC meteorites require at least two long-lived mantle reservoirs, and possibly three. The range in the projected present day gamma(sub Os) of these reservoirs is from -5.4+/-2.6, to +4+/-1. The isotopic systematics of these reservoirs may be linked to development of cumulate crystal piles in a Martian magma ocean and variable amounts of late stage intercumulus melt. In this model, fractional crystallization of olivine and possibly other phases with slightly subchondritic Re/Os, from a solidifying magma ocean, resulted in a lower Re/Os ratio early cumulates, and a resultant low gamma(sub Os). Later cumulates or evolved melts crystallized with higher Re/Os ratios to produce the mantle reservoir(s) with consequent higher gamma(sub Os). Crystallization of the Martian magma ocean followed earliest core formation, as indicated by the correlation of epsilon(sub W) with epsilon(sub 142Nd) and initial gamma(sub Os).

Brandon, Alan D.↗
NASA NTRS1980

Trophic classification of Tennessee Valley area reservoirs derived from LANDSAT multispectral scanner data

LANDSAT MSS data from four different dates were extracted from computer tapes using a semiautomated digital data handling and analysis system. Reservoirs were extracted from the surrounding land matrix by using a Band 7 density level slice of 3; and descriptive statistics to include mean, variance, and ratio between bands for each of the four bands were calculated. Significant correlations ( 0.80) were identified between the MSS statistics and many trophic indicators from ground truth water quality data collected at 35 reservoirs in the greater Tennessee Valley region. Regression models were developed which gave significant estimates of each reservoir's trophic state as defined by its trophic state index and explained in all four LANDSAT frames at least 85 percent of the variability in the data. To illustrate the spatial variations within reservoirs as well as the relative variations between reservoirs, a table look up elliptical classification was used in conjunction with each reservoir's trophic state index to classify each reservoir on a pixel by pixel basis and produce color coded thematic representations.

Meinert, D. L.↗
NASA NTRS2022

SMAP Soil Moisture Assimilation To Enhance Streamflow Estimates Across South Asia

Streamflow estimation across areas facing water scarcity is important considering the evolving climatic conditions. In this study, the impact of assimilating Soil Moisture Active Passive (SMAP) soil moisture retrievals on the modeled streamflow across three main river basins in South Asia is explored. Model estimated runoff was hydraulically routed using the HYMAP routing model within the NASA Land Information System (LIS) framework. Streamflow estimates via soil moisture assimilation highlighted the potential improvements in streamflow across low-flow river tributaries in irrigated regions. Limited relative change was noted along high-flow river tributaries. Modeled streamflow tended to underestimate the flow magnitude at upstream stations and overestimated the streamflow at downstream stations within the Indus basin due to missing physics related to reservoir operation. Improvement of the modeling systems representativeness of the ground conditions by the inclusion of water management information would potentially improve the current streamflow modeling capability across South Asia.

Jawairia Ahmad↗
NASA NTRS1994

Evaluating GCM land surface hydrology parameterizations by computing river discharges using a runoff routing model: Application to the Mississippi basin

To relate general circulation model (GCM) hydrologic output to readily available river hydrographic data, a runoff routing scheme that routes gridded runoffs through regional- or continental-scale river drainage basins is developed. By following the basin overland flow paths, the routing model generates river discharge hydrographs that can be compared to observed river discharges, thus allowing an analysis of the GCM representation of monthly, seasonal, and annual water balances over large regions. The runoff routing model consists of two linear reservoirs, a surface reservoir and a groundwater reservoir, which store and transport water. The water transport mechanisms operating within these two reservoirs are differentiated by their time scales; the groundwater reservoir transports water much more slowly than the surface reservior. The groundwater reservior feeds the corresponding surface store, and the surface stores are connected via the river network. The routing model is implemented over the Global Energy and Water Cycle Experiment (GEWEX) Continental-Scale International Project Mississippi River basin on a rectangular grid of 2 deg X 2.5 deg. Two land surface hydrology parameterizations provide the gridded runoff data required to run the runoff routing scheme: the variable infiltration capacity model, and the soil moisture component of the simple biosphere model. These parameterizations are driven with 4 deg X 5 deg gridded climatological potential evapotranspiration and 1979 First Global Atmospheric Research Program (GARP) Global Experiment precipitation. These investigations have quantified the importance of physically realistic soil moisture holding capacities, evaporation parameters, and runoff mechanisms in land surface hydrology formulations.

Liston, G. E.↗
NASA NTRS2022

Reservoir Assessment Tool 2.0: Stakeholder driven improvements to satellite remote sensing based reservoir monitoring

In light of the rapidly increasing regulation of rivers due to planned and constructed reservoirs, monitoring reservoir operations has become very crucial. The Reservoir Assessment Tool (RAT) framework was developed to monitor reservoir operations globally, using hydrological modeling and satellite observations. With feedback from stakeholders, improvements in the RAT framework are demonstrated in this study using the Mekong River Basin as an example. A novel multi-sensor reservoir area mapping technique was developed using complementary strengths of optical and SAR sensors at a 1-5 day temporal resolution, allowing the quantification of sub-weekly reservoir operations. Additionally, the skill of radar altimeters in the RAT framework was tested using the Jason-3 altimeter. Using in-situ data from three Thai reservoirs in the Mekong Basin, consistent improvements were observed as compared to the original RAT framework. A functionality for visualizing forecasted outflow was also added using historically inferred reservoir operations or stakeholder-driven target storage expectations.

Reservoirs↗
NASA NTRS1994

Io's heat flow from infrared radiometry: 1983-1993

We report the following results from a decade of infrared radiometry of Io: (1) The average global heat flow is more than approx. 2.5 W/sq.m, (2) large warm (less than or equal to 200 K) volcanic regions dominate the global heat flow, (3) smal high-temperature (greater than or = 300 K) 'hotspots' contribute little to the average heat flow, (4) thermal anomalies on the leading hemisphere contribute about half of the heat flow, (5) a substantial amount of heat is radiated during Io's night, (6) high-temperature (greater than or = 600 K) 'outbursts' occurred during approx. 4% of the nights we observed, (7) 'Loki' is the brightest, persistent, infrared emission feature, and (8) some excess emission is always present at the longitude of Loki, but its intensity and other characteristics change between apparitions. Observations of Io at M(4.8 micrometer), 8.7 micrometer, N(10 micrometer), and Q(20 micrometer) with the Infrared Telescope Facility presented here were collected during nine apparitions between 1983 and 1993. These measurements provide full longitudinal coveraged as well as an eclipse observation and the detection of two outbursts. Reflected sunlight, passive thermal emission, and radiation from thermal anomalies all contribute to the observed flux densities. We find that a new thermophysical model is required to match all the data. Two key elements of this model are (1) a 'thermal reservoir' unit which lowers daytime temperatures, and (2) the 'thermal pedestal effect' which shifts to shorter wavelengths the spectral emission due to the reradiation of solar energy absorbed by the thermal anomalies. The thermal anomalies are modeled with a total of 10 source components at five locations. Io's heat flow is the sum of the power from these components.

Veeder, Glenn J.↗
NASA NTRS2019

Quantifying Water Storage Change and Land Subsidence Induced by Reservoir Impoundment Using GRACE, Landsat, and GPS Data

The construction of hydropower dams is a common strategy to support a country's increasing need for electricity and river water management for industry and agriculture. Although the hydrological and geophysical impacts of water relocation are usually assessed prior to impoundment, their accuracy is generally limited due to the lack of in situ observations, especially in a remote area. This study presents a workflow to quantify the terrestrial water storage change (TWS) and land subsidence induced by a reservoir's water impoundment using multiple satellite observations (GRACE, Landsat), land surface models (CABLE, GLDAS, NCEP, ECMWF), and GPS data. The study site is the Bakun Dam, located in Sarawak, Malaysia, which is the largest hydropower dam in Southeast Asia. Commencing operation in late 2010, the dam induced a change of water mass and lake surface area that was clearly observed by GRACE and Landsat observations, respectively. During the 17-month impounding period (from August 2010 to December 2011), GRACE observed a dramatic increase of approximately 200 mm equivalent water height, while Landsat detected an increased lake extent of around 600 km2. In this paper, a forward model is developed to determine the increased water surface level corresponding to GRACE observations, estimated to be about 120 m. In contrast to GRACE, the TWS derived from land surface models cannot capture the increased TWS, due to the lack of reservoir routing algorithms in the models. In addition, the land subsidence was calculated using the disk load model constructed based on the GRACE-derived lake level and Landsat-derived lake extent; the result is validated with the GPS data from BIN1 station, located at the western coast of Borneo. The commencement stage of the Bakun Dam induces the large-scale land subsidence, which causes the GPS-BIN1 station to subside by ~9 mm, and move toward the Bakun Lake by ~4 mm. Computation of the surface displacements directly from GRACE spherical harmonic coefficient data fails to capture the subsidence feature, mainly due to the truncation error. Overall, this study demonstrates that evaluating GRACE in conjunction with Landsat, LSMs, and GPS data allows the exploitation of the gravity signal at a much smaller spatial scale than its intrinsic resolution. Benefiting from global coverage, the newly developed satellite-based algorithm is a valuable tool for assessing the impacts of reservoir operation on hydrological and geophysical changes from local to regional scales.

Natthachet Tangdamrongsub↗