Abstract Lake surface conditions are critical for representing lake‐atmosphere interactions in numerical weather prediction. The Community Land Model's 1‐D lake component (CLM‐lake) is part of NOAA's High‐Resolution Rapid Refresh (HRRR) 3‐km weather/earth‐system model, which assumes that virtually all the two thousand lakes represented in CONUS have distinct (for each lake) but spatially uniform depth. To test the sensitivity of CLM‐lake to bathymetry, we ran CLM‐lake as a stand‐alone model for all of 2019 with two bathymetry data sets for 23 selected lakes: the first had default (uniform within each lake) bathymetry while the second used a new, spatially varying bathymetry. We validated simulated lake surface temperature (LST) with both remote and in situ observations to evaluate the skill of both runs and also intercompared modeled ice cover and evaporation. Though model skill varied considerably from lake to lake, using the new bathymetry resulted in marginal improvement over the default. The more important finding is the influence bathymetry has on modeled LST (i.e., differences between model simulations) where lake‐wide LST deviated as much as 10°C between simulations and individual grid cells experienced even greater departures. This demonstrates the sensitivity of surface conditions in atmospheric models to lake bathymetry. The new bathymetry also improved lake depths over the (often too deep) previous value assumed for unknown‐depth lakes. These results have significant implications for numerical weather prediction, especially in regions near large lakes where lake surface conditions often influence the state of the atmosphere via thermal regulation and lake effect precipitation.
more »
« less
Landsat-derived bathymetry of lakes on the Arctic Coastal Plain of northern Alaska
Abstract. The Pleistocene sand sea on the Arctic Coastal Plain (ACP) ofnorthern Alaska is underlain by an ancient sand dune field, a geologicalfeature that affects regional lake characteristics. Many of these lakes,which cover approximately 20 % of the Pleistocene sand sea, are relativelydeep (up to 25 m). In addition to the natural importance of ACP sand sealakes for water storage, energy balance, and ecological habitat, the needfor winter water for industrial development and exploration activities makeslakes in this region a valuable resource. However, ACP sand sea lakes havereceived little prior study. Here, we collect in situ bathymetric data totest 12 model variants for predicting sand sea lake depth based on analysisof Landsat-8 Operational Land Imager (OLI) images. Lake depth gradients weremeasured at 17 lakes in midsummer 2017 using a Humminbird 798ci HD SI Comboautomatic sonar system. The field-measured data points were compared tored–green–blue (RGB) bands of a Landsat-8 OLI image acquired on 8 August2016 to select and calibrate the most accurate spectral-depth model for eachstudy lake and map bathymetry. Exponential functions using a simple bandratio (with bands selected based on lake turbidity and bed substrate)yielded the most successful model variants. For each lake, the most accuratemodel explained 81.8 % of the variation in depth, on average. Modeled lakebathymetries were integrated with remotely sensed lake surface area toquantify lake water storage volumes, which ranged from 1.056×10-3 to 57.416×10-3 km3. Due to variations in depthmaxima, substrate, and turbidity between lakes, a regional model iscurrently infeasible, rendering necessary the acquisition of additional insitu data with which to develop a regional model solution. Estimating lakewater volumes using remote sensing will facilitate better management ofexpanding development activities and serve as a baseline by which toevaluate future responses to ongoing and rapid climate change in the Arctic.All sonar depth data and modeled lake bathymetry rasters can be freelyaccessed at https://doi.org/10.18739/A2SN01440 (Simpson and Arp, 2018) andhttps://doi.org/10.18739/A2HT2GC6G (Simpson, 2019), respectively.
more »
« less
- Award ID(s):
- 1806213
- PAR ID:
- 10218061
- Date Published:
- Journal Name:
- Earth System Science Data
- Volume:
- 13
- Issue:
- 3
- ISSN:
- 1866-3516
- Page Range / eLocation ID:
- 1135 to 1150
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Melt and supraglacial lakes are precursors to ice shelf collapse and subsequent accelerated ice sheet mass loss. We used data from the Landsat 8 and Sentinel-2 satellites to develop a threshold-based method for detection of lakes found on the Antarctic ice shelves, calculate their depths and thus their volumes. To achieve this, we focus on four key areas: the Amery, Roi Baudouin, Nivlisen, and Riiser-Larsen ice shelves, which are all characterized by extensive surface meltwater features. To validate our products, we compare our results against those obtained by an independent method based on a supervised classification scheme (e.g., Random Forest algorithm). Additional verification is provided by manual inspection of results for nearly 1000 Landsat 8 and Sentinel-2 images. Our dual-sensor approach will enable constructing high-resolution time series of lake volumes. Therefore, to ensure interoperability between the two datasets, we evaluate depths from contemporaneous Landsat 8 and Sentinel-2 image pairs. Our assessments point to a high degree of correspondence, producing an average R2 value of 0.85, no bias, and an average RMSE of 0.2 m. We demonstrate our method’s ability to characterize lake evolution by presenting first evidence of drainage events outside of the Antarctic Peninsula on the Amery Ice shelf. The methods presented here pave the way to upscaling throughout the Landsat 8 and Sentinel-2 observational record across Antarctica to produce a first-ever continental dataset of supraglacial lake volumes. Such a dataset will improve our understanding of the influence of surface hydrology on ice shelf stability, and thus, future projections of Antarctica’s contribution to sea level rise.more » « less
-
Abstract Depth regulates many attributes of aquatic ecosystems, but relatively few lakes are measured, and existing datasets are biased toward large lakes. To address this, we used a large dataset of maximum (Zmax;n = 16,831) and mean (Zmean;n = 5,881) depth observations to create new depth models, focusing on lakes < 1,000 ha. We then used the models to characterize patterns in lake basin shape and volume. We included terrain metrics, water temperature and reflectance, polygon attributes, and other predictors in a random forest model. Our final models generally outperformed existing models (Zmax; root mean square error [RMSE] = 8.0 m andZmean; RMSE = 3.0 m). Our models show that lake depth followed a Pareto distribution, with 2.8 orders of magnitude fewer lakes for an order of magnitude increase in depth. In addition, despite orders of magnitude variation in surface area, most size classes had a modal maximum depth of ~ 5 m. Concave (bowl‐shaped) lake basins represented 79% of all lakes, but lakes were more convex (funnel‐shaped) as surface area increased. Across the conterminous United States, 9.8% of all lake water was within the top meter of the water column, and 48% in the top 10 m. Excluding the Laurentian Great Lakes, we estimate the total volume in the conterminous United States is 1,057–1,294 km3, depending on whetherZmaxorZmeanwas modeled. Lake volume also exhibited substantial geographic variation, with high volumes in the upper Midwest, Northeast, and Florida and low volumes in the southwestern United States.more » « less
-
Nearly 25% of all lakes on earth are located at high latitudes. These lakes are formed by a combination of thermokarst, glacial, and geological processes. Evidence suggests that the origin of periglacial lake formation may be an important factor controlling the likelihood of lakes to drain. However, geospatial data regarding the spatial distribution of these dominant Arctic and subarctic lakes are limited or do not exist. Here, we use lake-specific morphological properties using the Arctic Digital Elevation Model (DEM) and Landsat imagery to develop a Thermokarst lake Settlement Index (TSI), which was used in combination with available geospatial datasets of glacier history and yedoma permafrost extent to classify Arctic and subarctic lakes into Thermokarst (non-yedoma), Yedoma, Glacial, and Maar lakes, respectively. This lake origin dataset was used to evaluate the influence of lake origin on drainage between 1985 and 2019 in northern Alaska. The lake origin map and lake drainage datasets were synthesized using five-year seamless Landsat ETM+ and OLI image composites. Nearly 35,000 lakes and their properties were characterized from Landsat mosaics using an object-based image analysis. Results indicate that the pattern of lake drainage varied by lake origin, and the proportion of lakes that completely drained (i.e., >60% area loss) between 1985 and 2019 in Thermokarst (non-yedoma), Yedoma, Glacial, and Maar lakes were 12.1, 9.5, 8.7, and 0.0%, respectively. The lakes most vulnerable to draining were small thermokarst (non-yedoma) lakes (12.7%) and large yedoma lakes (12.5%), while the most resilient were large and medium-sized glacial lakes (4.9 and 4.1%) and Maar lakes (0.0%). This analysis provides a simple remote sensing approach to estimate the spatial distribution of dominant lake origins across variable physiography and surficial geology, useful for discriminating between vulnerable versus resilient Arctic and subarctic lakes that are likely to change in warmer and wetter climates.more » « less
-
Abstract As ice sheets load Earth's surface, they produce ice‐marginal depressions which, when filled with meltwater, become proglacial lakes. We include self‐consistently evolving proglacial lakes in a glacial isostatic adjustment (GIA) model and apply it to the Laurentide ice sheet over the last glacial cycle. We find that the locations of modeled lakes and the timing of their disappearance is consistent with the geological record. Lake loads can deflect topography by >10 m, and volumes collectively approach 30–45 cm global mean sea‐level equivalent. GIA increases deglaciation‐phase lake volume up to five‐fold and average along‐ice‐margin depth ≤90 m compared to glaciation‐phase ice volume analogs—differences driven by changes in the position and size of the peripheral bulge. Since ice‐marginal lake depth affects grounding‐line outflow, GIA‐modulated proglacial lake depths could affect ice‐sheet mass loss. Indeed, we find that Laurentide ice‐margin retreat rate sometimes correlates with proglacial lake presence, indicating that proglacial lakes aid glacial collapse.more » « less
An official website of the United States government

