Search NASA⌕ Search

SEARCH · Search NASA

Results for “Surface hydrology”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
NASA NTRS2013

Connecting Satellite Observations with Water Cycle Variables Through Land Data Assimilation: Examples Using the NASA GEOS-5 LDAS

A land data assimilation system (LDAS) can merge satellite observations (or retrievals) of land surface hydrological conditions, including soil moisture, snow, and terrestrial water storage (TWS), into a numerical model of land surface processes. In theory, the output from such a system is superior to estimates based on the observations or the model alone, thereby enhancing our ability to understand, monitor, and predict key elements of the terrestrial water cycle. In practice, however, satellite observations do not correspond directly to the water cycle variables of interest. The present paper addresses various aspects of this seeming mismatch using examples drawn from recent research with the ensemble-based NASA GEOS-5 LDAS. These aspects include (1) the assimilation of coarse-scale observations into higher-resolution land surface models, (2) the partitioning of satellite observations (such as TWS retrievals) into their constituent water cycle components, (3) the forward modeling of microwave brightness temperatures over land for radiance-based soil moisture and snow assimilation, and (4) the selection of the most relevant types of observations for the analysis of a specific water cycle variable that is not observed (such as root zone soil moisture). The solution to these challenges involves the careful construction of an observation operator that maps from the land surface model variables of interest to the space of the assimilated observations.

Land data assimilation↗
NASA NTRS1999

Coupled Land Atmosphere Predictability

We have designed and executed a set of predictability experiments, designed around the driest and wettest June soil moisture anomalies from a CCM3 simulation forced by observed SST for the period from 1958 through 1998. Each set contains an ensemble of five runs, all begun on June I radiation date. One set of these experiments helps to assess the extent to which the wet or dry conditions depend solely on the initial state of the atmosphere. The other set of experiments helps to assess the extent to which the wet or dry conditions depend sole on the initial state of the land surface. Preliminary analysis of these experiments suggests that the initial atmospheric state is more important than the initial state of the surface soil moisture in predicting the occurrence of wet or dry periods. These results suggest that when the atmosphere is inclined to generate dry surface conditions (through reduced moisture availability and increased evaporation) it matters little what initial levels of soil water are at; the soil will rapidly dry out. The ensemble forcing the ensembles with dry soil conditions, but utilizing 'normal' atmospheric conditions show little if any indication of the sharp reduction in soil moisture experienced in the control, indicating that the 'normal' atmospheric state is more important than the initial state of the soil moisture in predicting the occurrence of drought. Our work to date has documented the response of surface hydrologic variability, particularly over North America, to atmospheric forcing. While we have seen some suggestion that pre-existing surface anomalies can affect atmospheric circulation over this region, by far the strongest signal is the atmospheric anomalies leading those of soil moisture and surface energy balance changes. This is not an unexpected result since wintertime NA precipitation links to remote atmospheric forcing by ENSO are known to be statistically significant. While warm season links to remote forcing are more tenuous and not well explored, the present results encourage us to examine in more detail the SST forcing from the tropical and North Pacific.

Roads, John↗
NASA NTRS1997

Soil Moisture: The Hydrologic Interface Between Surface and Ground Waters

A hypothesis is presented that many hydrologic processes display a unique signature that is detectable with microwave remote sensing. These signatures are in the form of the spatial and temporal distributions of surface soil moisture. The specific hydrologic processes that may be detected include groundwater recharge and discharge zones, storm runoff contributing areas, regions of potential and less than potential evapotranspiration (ET), and information about the hydrologic properties of soils. In basin and hillslope hydrology, soil moisture is the interface between surface and ground waters.

Engman, Edwin T.↗
NASA NTRS2017

The SMAP Level-4 ECO Project: Linking the Terrestrial Water and Carbon Cycles

The SMAP (Soil Moisture Active Passive) Level-4 projects aims to develop a fully coupled hydrology-vegetation data assimilation algorithm to generate improved estimates of modeled hydrological fields and carbon fluxes. This includes using the new NASA Catchment-CN (Catchment-Carbon-Nitrogen) model, which combines the Catchment land surface hydrology model with dynamic vegetation components from the Community Land Model version 4 (CLM4). As such, Catchment-CN allows a more realistic, fully coupled feedback between the land hydrology and the biosphere. The L4 ECO project further aims to inform the model through the assimilation of Soil Moisture Active Passive (SMAP) brightness temperature observations as well as observations of Moderate Resolution Imaging Spectroradiometer (MODIS) fraction of absorbed photosynthetically active radiation (FPAR). Preliminary results show that the assimilation of SMAP observations leads to consistent improvements in the model soil moisture skill. An evaluation of the Catchment-CN modeled vegetation characteristics showed that a calibration of the model's vegetation parameters is required before an assimilation of MODIS FPAR observations is feasible.

MODIS↗
NASA NTRS1998

A Technique for Assimilating GOES-Derived Land Surface Products into Regional Models to Improve the Representation of Land Surface Forcing

As the parameterizations of surface energy budgets in regional models have become more complete physically, models have the potential to be much more realistic in simulations of coupling between surface radiation, hydrology, and surface energy transfer. Realizing the importance of properly specifying the surface energy budget, many institutions are using land-surface models to represent the lower boundary forcing associated with biophysical processes and soil hydrology. However, the added degrees of freedom due to inclusion of such land-surface schemes require the specification of additional parameters within the model system such as vegetative resistances, green vegetation fraction, leaf area index, soil physical and hydraulic characteristics, stream flow, runoff, and the vertical distribution of soil moisture. A technique has been developed for assimilating GOES-IR skin temperature tendencies into the surface energy budget equation of a mesoscale model so that the simulated rate of temperature change closely agrees with the satellite observations. A critical assumption of the technique is that the availability of moisture (either from the soil or vegetation) is the least known term in the model's surface energy budget. Therefore, the simulated latent heat flux, which is a function of surface moisture availability, is adjusted based upon differences between the modeled and satellite-observed skin temperature tendencies. An advantage of this technique is that satellite temperature tendencies are assimilated in an energetically consistent manner that avoids energy imbalances and surface stability problems that arise from direct assimilation of surface shelter temperatures. The fact that the rate of change of the satellite skin temperature is used rather than the absolute temperature means that sensor calibration is not as critical. An advantage of this technique for short-range forecasts is that it does not require a complex land-surface formulation within the atmospheric model. As a result, the need to specify poorly known soil and vegetative characteristics is eliminated. The GOES assimilation technique has been incorporated into the PSU/NCAR MM5. Results will be presented to demonstrate the ability of the assimilation scheme to improve short-term simulations of near-surface air temperature and mixing ratio during the warm season for several selected cases during June of 1995 and June 1997 which exhibit a variety of atmospheric and land-surface conditions. In addition, the simulations produced by the assimilation technique will be compared with those produced by the MM5/BATS coupled model initialized with soil water content derived from a soil hydrology model.

Lapenta, William M.↗
NASA NTRS1993

Coupled land surface/hydrologic/atmospheric models

The topics covered include the following: prototype land cover characteristics data base for the conterminous United States; surface evapotranspiration effects on cumulus convection and implications for mesoscale models; the use of complex treatment of surface hydrology and thermodynamics within a mesoscale model and some related issues; initialization of soil-water content for regional-scale atmospheric prediction models; impact of surface properties on dryline and MCS evolution; a numerical simulation of heavy precipitation over the complex topography of California; representing mesoscale fluxes induced by landscape discontinuities in global climate models; emphasizing the role of subgrid-scale heterogeneity in surface-air interaction; and problems with modeling and measuring biosphere-atmosphere exchanges of energy, water, and carbon on large scales.

Pielke, Roger↗
NASA NTRS1999

Real-Time Assimilation of Goes-Derived Products into A Mesoscale Model and It's Impact on Short-Term (06-36hr) Forecasts from 17 October 1998 through the Present

As the parameterizations of surface energy budgets in regional models have become more complete physically, models have the potential to be much more realistic in simulations of coupling between surface radiation, hydrology, and surface energy transfer. Realizing the importance of properly specifying the surface energy budget, many institutions are using land-surface models to represent the lower boundary forcing associated with biophysical processes and soil hydrology. However, the added degrees of freedom due to inclusion of such land-surface schemes require the specification of additional parameters within the model system such as vegetative resistances, green vegetation fraction, leaf area index, soil physical and hydraulic characteristics, stream flow, runoff, and the vertical distribution of soil moisture. Spatial heterogeneity of these parameters makes correct specification problematic since measurements are not routinely available. A technique has been developed for assimilating GOES-IR skin temperature tendencies, solar insolation, and surface albedo into the surface energy budget equation of a mesoscale model so that the simulated rate of temperature change closely agrees with the satellite observations. The technique has been successfully employed in a number of mesoscale models in case-study mode. We have taken the next step and developed a study to determine if assimilating these types of data into mesoscale models in real-time can improve short-term (648h) forecasts of temperature, relative humidity, and QPF on a daily basis over relatively large regions. Therefore, an operational modeling/assimilation system has been developed at the GHCC during the past summer that allows us to produce simulations out to 48 hours in a timely manor. The PSU/NCAR MM5 is used in a nested configuration with a 25 km grid covering the southeastern third of the US. The model has been on-line since 1 July 1998 and forecast products are posted on our web site. The satellite algorithms that generate data to be assimilated came on-line 17 October 1998. Quantitative assessment of the forecast quality is performed via traditional verification statistics. In addition, invaluable qualitative information is obtained through close collaboration with several NWSFO's who are using the MM5 products in real-time on a daily basis. The assimilation technique has been applied in an off-line mode since 17 October. Results based on bulk statistical verification of surface meteorology over the entire Southeastern US show that assimilating the GOES-derived land surface tendencies and solar radiation results in a significant reduction of the shelter air temperature and RH bias on a daily basis. In fact, the assimilation technique has produced improved temperature and RH forecasts for 97% of the 100 simulations performed to date. Work is currently underway to determine the sensitivity of the assimilation procedure to the availability of satellite data, length of assimilation period, model initialization, and synoptic-scale meteorological conditions. In addition, results from a detailed energy budget analysis using the Early Eta, our operational MM5, and the assimilation runs will help us to better understand the satellite assimilation the land-surface energy budge. Research during the spring-summer of 1999 will focus on the impact of the assimilation technique during the warm season where it is hypothesized that it can have a positive impact on QPF during conditions of weak synoptic-scale forcing.

Lapenta, William M.↗
NASA NTRS1976

An aid to the development of Botswana's resources: Section on hydrology

The author has identified the following significant results. It is proved that FCC's can be used for a simple estimate of the total evaportranspiring area of the Okavango Delta, sufficiently accurate for preliminary inputs for the development of mathematical model of the surface hydrology of the delta. The color coded matrix has shown as interesting inverse correlation with an array on the same grid prepared by ecologists from air photography study, for percent liable to flood.

Dickson, W. L.↗
NASA NTRS2003

From the South Pole to the Northern Plains: The Argyre Planitia Story

Parker (1985, 1994) first described evidence for catastrophic flooding from a large lake or sea within Argyre Planitia through the Uzboi-Holden- Ladon-Margaritifer Valles system during the Noachian. The channel connection to Argyre had been recognized during the mid-1970s, based primarily on Russian orbiter images taken at that time. The most critical reviews of these inferences related to the relative timing of the plains materials, sinuous ridges, and debris aprons in southern Argyre, and the connection, via Uzboi Vallis, of ponding within Argyre to flooding through the Chryse Trough. The prevailing "competing" hypothesis for formation of materials within Argyre is that they are a result of south circumpolar glacial processes, with glacial scour and stagnation producing the pitting and sinuous ridges (eskers) on the basin floor rather than lacustrine erosion and deposition followed much later by a process akin to rock glacier formation of the debris aprons in a colder Amazonian climate. Argyre was part of a larger surface hydrological system that also included two large valley networks draining the Margaritifer Sinus region northwest of Argyre. The morphometry of these systems suggest a combination of precipitation and groundwater sapping, with surface runoff for their formation.

T J Parker↗
NASA NTRSDate not supplied

HydroSAR: A Cloud-based SAR Data Analysis Service to Monitor Hydrological Disasters and their Impact on Population and Agriculture

Weather-related hazards are ubiquitous around the world including: 1) hurricane storm surges, 2) rapid snowmelt and heavy rainfall, 3) severe weather leading to flash floods, and 4) seasonal freeze and thaw of rivers that may lead to ice jams. Each of these hazards affects human settlements and has the potential to impact agricultural productivity. In each setting, end-users in disaster management need access to data processing tools helpful in mapping past and current disasters. Analysis of past events supports risk mitigation by understanding what has already occurred and how to alleviate those impacts in the future. Having capabilities to generate the same products in a response setting means that lessons learned from risk analysis will carry forward to event response. Synthetic aperture radar (SAR) data are particularly useful for these activities due to their all-weather 24/7 monitoring capabilities. In this effort we present HydroSAR, a cloud-based SAR data analysis service for the mapping of meteorological and hydrological disasters as well as their impact on population and agriculture. As part of this project we have developed a series of SAR-based value added products for the monitoring of surface hydrology (image time series, change detection, flood extent, flood depth) and the assessment of impacts on population (flood depth) and agriculture (active agriculture, inundated agriculture, flood duration). We also developed a cloud-based platform for generating these products over affected areas and are working with end-users to integrate derived product into decision-making workflows The paper will briefly introduce the SAR-based products that were developed for this effort. We describe the cloud-based production pipeline that was built to automatically generate these products in near-real time over extended regions. The integration of SAR-based information into hazard preparation and response activities is described for a number of recent disasters including the 2019 forest fires in Alaska, 2019 flooding in the U.S. Midwest, the 2020 U.S. severe weather easter outbreak, 2020 tropical storm Christobal, 2020 cyclone Amphan, 2020, Alaska Spring breakup flooding, and the 2020 flood season in Eastern India, Bangladesh, and Nepal.

Franz Josef Meyer↗
NASA NTRS2005

The Use of Mascons to Resolve Time-Variable Gravity from GRACE

We have analyzed GRACE Level 1B data to resolve time-variable gravity using a local mascon approach. The spherical harmonic solutions released to date resolve the signal from surface hydrology over land areas at spatial scales of 750 to 1000 km over one month intervals [Wahr et al., 2004; Tapley et al., 2004]. In our local approach, we solve explicitly for the mass of water in surface blocks using only the KBRR data collected as GRACE overflies the region of interest. The local representation of gravity minimizes leakage of errors from other areas due to aliasing or mismodelling. In this paper, we report on the analysis of GRACE data from January 2003 through August 2004 over three regions: the Amazon, the Indian subcontinent, and the continental United States. We solve for mass change at 10-day intervals using 4 deg x 4 deg blocks. We give an overview of our latest results, and we present the results of error analyses, and comparisons to both hydrology models and in-situ data.

Lemoine, Frank G.↗
NASA NTRS2014

"Efficiency Space" - A Framework for Evaluating Joint Evaporation and Runoff Behavior

At the land surface, higher soil moisture levels generally lead to both increased evaporation for a given amount of incoming radiation (increased evaporation efficiency) and increased runoff for a given amount of precipitation (increased runoff efficiency). Evaporation efficiency and runoff efficiency can thus be said to vary with each other, motivating the development of a unique hydroclimatic analysis framework. Using a simple water balance model fitted, in different experiments, with a wide variety of functional forms for evaporation and runoff efficiency, we transform net radiation and precipitation fields into fields of streamflow that can be directly evaluated against observations. The optimal combination of the functional forms the combination that produces the most skillful stream-flow simulations provides an indication for how evaporation and runoff efficiencies vary with each other in nature, a relationship that can be said to define the overall character of land surface hydrological processes, at least to first order. The inferred optimal relationship is represented herein as a curve in efficiency space and should be valuable for the evaluation and development of GCM-based land surface models, which by this measure are often found to be suboptimal.

GCM↗
NASA NTRS2019

Chapter 14: 100 Years of Progress in Hydrology

The focus of this chapter is progress in hydrology for the last 100 years. During this period, we have seen a marked transition from practical engineering hydrology to fundamental developments in hydrologic science, including contributions to Earth system science. The first three sections in this chapter review advances in theory, observations, and hydrologic prediction. Building on this foundation, the growth of global hydrology, land-atmosphere interactions and coupling, ecohydrology, and water management are discussed, as well as a brief summary of emerging challenges and future directions. Although the review attempts to be comprehensive, the chapter offers greater coverage on surface hydrology and hydrometeorology for readers of this American Meteorological Society (AMS) Monograph.

Peters-Lidard, Christa D.↗
NASA NTRS1994

Structural analysis of airborne flux estimates over a region

Aircraft-based observations of turbulence fields of velocity, moisture, and temperature are used to study coherent turbulent structures that dominate turbulent transfer of moisture and heat above three different eco-systems. Flux traces are defragmented, to reconstruct the presumed full size (along the sampled transect) of these structures, and flux traces are simplified by elimination of those that contribute negligibly to the flux estimate. Structures are analyzed in terms of size, spatial distribution, and contribution to the flux, in the four 'quadrant' modes of eddy-covariance transfer (excess up/down and deficit up/down). The effect of nonlinear detrending of moisture and temperature data on this 'structural analysis,' over surfaces with heterogeneous surface wetness, is also examined. Results over grassland, wetland, and moist and dry agricultural land, show that nonlinear detrending may provide a more physically realistic description of structures. Significant differences are observed between structure size and associated relative flux contribution, between moist and dry areas, with smaller structures playing a more important role over the moist areas. Structure size generally increases with height, as spatial reorganization from smaller structures into larger ones takes place. This coincides with a gradual loss of surface 'signature' (position and clustering of plumes above localized source areas). The data are expected to provide a basis for an eventual statistical description of boundary-layer transfer events , and help to interpret the link between boundary-layer transfer and hydrological surface conditions.

Caramori, Paulo↗
NASA NTRS2015

Technical Report Series on Global Modeling and Data Assimilation, Volume 43: Initial Evaluation of the Climate - MERRA-2

The years since the introduction of MERRA have seen numerous advances in the GEOS-5 Data Assimilation System as well as a substantial decrease in the number of observations that can be assimilated into the MERRA system. To allow continued data processing into the future, and to take advantage of several important innovations that could improve system performance, a decision was made to produce MERRA-2, an updated retrospective analysis of the full modern satellite era. One of the many advances in MERRA-2 is a constraint on the global dry mass balance; this allows the global changes in water by the analysis increment to be near zero, thereby minimizing abrupt global interannual variations due to changes in the observing system. In addition, MERRA-2 includes the assimilation of interactive aerosols into the system, a feature of the Earth system absent from previous reanalyses. Also, in an effort to improve land surface hydrology, observations-corrected precipitation forcing is used instead of model-generated precipitation. Overall, MERRA-2 takes advantage of numerous updates to the global modeling and data assimilation system. In this document, we summarize an initial evaluation of the climate in MERRA-2, from the surface to the stratosphere and from the tropics to the poles. Strengths and weaknesses of the MERRA-2 climate are accordingly emphasized.

MERRA-2↗
NASA NTRS2012

Projections of Declining Surface-Water Availability for the Southwestern United States

16 of the CMIP5 models had all the data needed for this work for at least one simulation that was continuous from 1950 to 2040. Details of the models analyzed here are provided in Table S1. The model data analyzed here are available at http://strega.ldeo.columbia.edu:81/expert/home/.naomi/.AR5/.v2/.historical:rcp85/.mmm16/ a. Assessing the climatology of the models Despite increases in horizontal resolution of many models compared to their CMIP3 counterparts none of these models can adequately resolve the topography of the south west United States, such as the Sierra Nevada and Rocky Mountains and the associated orographic precipitation. This requires that caution be used when interpreting the results presented here. To assess the ability of the models to simulate the current hydroclimate, in Figure S1 we show the observed (from the Global Precipitation Climatology Centre gridded rain gauge data, (1)) monthly climatology of precipitation and the same for all the models and the multimodel mean for the California-Nevada, Colorado headwaters and Texas regions. The GPCC data uses rain gauges only and interpolates to regular grids of which we used the 1◦ by 1◦ one. Details of the data set can be found in (2). While the models apparently overestimate precipitation in California and Nevada the seasonal cycle with wet winters and dry summers is very well represented. It is also possible that the rain gauge observations are biased low by inadequately sampling the higher mountain regions. How ever the models might also be expected to underestimate orographic precipitation due to inadequate horizontal resolution. The 25 models are also too wet in the Colorado headwaters region but correctly represent the quite even distribution though the year. The bimodal distribution of precipitation in Texas, with peaks in May and September, and the absolute amounts, are well modeled but with the September peak too weak. The positive precipitation bias translates into a positive runoff bias for the Colorado headwaters as also shown in Figure S1. Here the observed runoff values are taken from simulations of the Variable Infiltration Capacity (VIC) land surface-hydrology model (3) forced by observed meteorology (5) that were conducted as part of the North American Land Data Assimilation System project phase 2 ( (NLDAS-2), http://www.emc.ncep.noaa.gov/mmb/nldas/. Runoff for California-Nevada is better simulated but there is a positive bias over Texas despite no strong precipitation bias. To check whether regional climate models better simulate P and runoff in these regions we analyzed the historical simulation with the Regional Climate Model version 3 driven by the National Centers for Environmental Prediction-Department of Energy Reanalysis 2 available from the North American Regional Climate Change Assessment Program (http://www.narccap.ucar.edu). This model configuration retained these biases in P and runoff although they were reduced in amplitude. Given these varying biases we plot P and P − E changes in actual values but apply the simplest bias correction possible to the runoff and soil moisture values and show the modeled changes in terms of percentages of the 20th Century model climatologies. A thorough assessment of the simulation of North American climate in CMIP5 models is conducted in Sheffield at al. (North American Climate in CMIP5 Experiments. Part I: Evaluation of 20th Century Continental and Regional Climatology, manuscript submit ted to J. Climate, available at http://www.climate.noaa.gov/index.jsp?pg=./cpo pa/ mapp/cmip5 publications.html). Sheffield et al. analyze the climatology of precipitation, surface air temperature, low level winds, moisture fluxes, runoff etc. and conclude that the main features of the hydrological cycle, including characteristics of the atmospheric moisture balance and its seasonality, are captured in the CMP5 models subject to biases in total precipitation amounts. We chose to use all available models instead of selecting some and rejecting others based on an assessment of model realism. This is in accord with the suggestions of Mote et al. for CMIP3 (4) but future work needs to revisit this matter for the case of the CMIP5 ensemble.

Seager, Richard↗
NASA NTRSDate not supplied

On the Prediction of Surface Melt Over the Greenland Ice Sheet at NWP and S2S Timescales

Greenland surface melt “events” – encompassing large swaths of the ice sheet over approximately 4-10 day periods – have occurred with increasing frequency in the 21st Century. Surface processes are significant to the ice sheet’s overall contribution to sea level. In addition, the sudden onset of widespread surface melt produces hazardous conditions for researchers and for communities adjacent to the ice sheet. The mechanisms producing melt events, their associated large-scale forcing, and timescales of interest are highly variable and poorly understood. Here we examine the predictability of melt events in a numerical weather prediction system (NWP) and a subseasonal-to-seasonal system (S2S). The NASA Global Modeling and Assimilation Office forward processing (GMAO-FP) is an NWP system that incorporates an evolving, global hybrid 4D ensemble–variational (4DEnVar) data assimilation system and a state-of-the-art model that is run twice-daily with forecasts out to five and ten days. The GMAO S2S system is a hindcast-anomaly, coupled-model prediction system initialized with an atmospheric analysis, and in situ and satellite ocean data assimilation. Both systems represent ice sheet surface hydrology and snowpack processes. We examine the performance of the systems over the period 2018-2023, a span that includes approximately nine significant events, against analyses and satellite-derived melt extent. In addition to extent, we examine related conditions and their role in the melt forecast, including the onset of both omega and Rex atmospheric blocking patterns, sea surface temperatures, sea ice extent, water vapor transport, precipitation, and surface preconditioning. We further examine the tendency for false-alarm events in each system and assess reasons for the forecast performance. North Atlantic atmospheric blocking patterns leading to surface melt are well forecast by the NWP system out to around eight days, while the S2S system generally underestimates the magnitude of such events at monthly lead times. These systems suggest that while short term forecasts of ice sheet melt and melt events are relatively robust, S2S predictions of surface conditions are highly dependent on the ability of a coupled system to generate the short-term atmospheric conditions associated with surface melt events.

Richard Cullather↗
NASA NTRS2021

Development of Ice-Shelf Estuaries Promotes Fractures and Calving

As the global climate warms, increased surface meltwater production on ice shelves may trigger ice-shelf collapse and enhance global sea-level rise. The formation of surface rivers could help prevent ice-shelf collapse if they can efficiently evacuate meltwater. Here we present observations of the evolution of a surface river into an ice-shelf estuary atop the Petermann Ice Shelf in northwest Greenland and identify a second estuary at the nearby Ryder Ice Shelf. This surface-hydrology process can foster fracturing and enhance calving. At the Petermann estuary, sea ice was observed converging at the river mouth upstream, indicating a flow reversal. Seawater persists in the estuary after the surrounding icescape is frozen. Along the base of Petermann estuary, linear fractures were initiated at the calving front and propagated upstream along the channel. Similar fractures along estuary channels shaped past large rectilinear calving events at the Petermann and Ryder ice shelves. Increased surface melting in a warming world will enhance fluvial incision, promoting estuary development and longitudinal fracturing orthogonal to ice-shelf fronts, and increase rectilinear calving. Estuaries could develop in Antarctica within the next half-century, resulting in increased calving and accelerating both ice loss and global sea-level rise.

Cryospheric science↗