Search NASA⌕ Search

SEARCH · Search NASA

Results for “lidar error”

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 325 records · Page 18

Merging IceSAT GLAS and Terra MODIS Data in Order to Derive Forest Type Specific Tree Heights in the Central Siberian Boreal Forest

Mapping of boreal forest's type, biomass, and other structural parameters are critical for understanding of the boreal forest's significance in the carbon cycle, its response to and impact on global climate change. We believe the nature of the forest structure information available from MISR and GLAS can be used to help identify forest type, age class, and estimate above ground biomass levels beyond that now possible with MODIS alone. The ground measurements will be used to develop relationships between remote sensing observables and forest characteristics and provide new information for understanding forest changes with respect to environmental change. Lidar is a laser altimeter that determines the distance from the instrument to the physical surface by measuring the time elapsed between the pulse emission and the reflected return. Other studies have shown that the returned signal may identify multiple returns originating from trees, building and other objects and permits the calculation of their height. Studies using field data have shown that lidar data can provide estimates of structural parameters such as biomass, stand volume and leaf area index and allows remarkable differentiation between primary and secondary forest. NASA's IceSAT Geoscience Laser Altimeter System (GLAS) was launched in January 2003 and collected data during February and September of that year. This study used data acquired over our study sites in central Siberia to examine the GLAS signal as a source of forest height and other structural characteristics. The purpose of our Siberia project is to improve forest cover maps and produce above-ground biomass maps of the boreal forest in Northern Eurasia from MODIS by incorporating structural information inherent in the Terra MISR and ICESAT Geoscience Laser Altimeter System (GLAS) instruments. A number of forest cover classifications exist for the boreal forest. We believe the limiting factor in these products is the lack of structural information, particularly in the vertical dimension. The emphasis of this project is to improve upon satellite maps of boreal forest structure parameters (i.e. height and biomass) using temporal, multi-angle, and vertical profile information of GLAS data. The existing and near future lidar data is useful for demonstrating these techniques and pursuing current estimates. Future lidar missions may be several years in the future, so we will work other new data sets that may aide in biomass estimates such as ALOS PALSAR We will continue this work to produce an accurate map of current above ground forest phytomass/carbon storage possible for the study area. We plan to develop, test, and integrate remote sensing methods for extracting forest canopy structure measures. We are compiling our field measurements and will compare them with the remote sensing methods where possible. We also be able to produce a realistic error bound on the remotely sensed carbon estimates.

Ranson, K. Jon↗

Shape from spectra

We introduce a new unified atmospheric–topographic correction approach that estimates surface geometry directly from the radiance measurement. Surface topography influences the at-sensor radiance measurement, making precise topography modeling critical in applications like vegetation or snow studies in mountainous terrain. Currently, elevation maps are used to derive topographic variables such as the slope and sky-view factor. This process is error-prone since static global digital elevation models do not generally achieve the accuracy required, and even minor mismatches in spatial resolution can introduce significant artifacts in downstream processing. Here we demonstrate that it is possible to estimate topographic parameters directly from spectral data, ensuring perfect physical consistency, temporal coincidence, and spatial alignment. We present experiments estimating topographic slope in two scenes in Southern California, with data from NASA’s Next Generation Airborne Visible/Near Infrared Imaging Spectrometer (AVIRIS-NG). We compared our radiance-based estimates against high-resolution lidar datasets. Our initial validation result showed a correlation of R 2 = 0.864 (n = 160) over the homogeneous surface of Beckman Auditorium’s cone-shaped roof on the Caltech campus in Pasadena, California. We then validate the model over a larger study site near Santa Clarita, California, finding R 2 = 0.923 (n = 40,000) in a 350 x 350m area. The accuracy of our model estimates, combined with its systematic advantages over the alternative, show the potential of the approach for use in both airborne campaigns and orbital missions.

Nimrod Carmon↗

Potential Basin-Scale Estimates of Arctic Snow Depth with Sea Ice Freeboards from CryoSat-2 and ICESat-2: An Exploratory Analysis

The potential of deriving snow depth estimates using differences in freeboard heights from CryoSat-2 (CS-2) and ICESat-2 (IS-2) is examined. In our analysis, we use lidar freeboard from the Airborne Topographic Mapper (ATM) on Operation IceBridge (OIB) as proxy of IS-2 total (snow plus ice) freeboard. Snow depths are estimates from the OIB snow radar. Differences in height between the total (ATM) and ice (CS-2) freeboards are related to snow depth by the refractive index of the snow layer, which is dependent on snow density. For two years (2014 and 2015), regression of the ATM and CS-2 freeboard differences against OIB snow depth gives correlations of approximately 0.80, estimated eta (sub s) of 1.21, and standard errors of 8 centimeters. The resulting refractive index, eta (sub s), can be compared to that expected of the Arctic snow cover in early spring (1.25 plus or minus 0.05). The expected biases and variability in the regression analysis are discussed. Results suggest that snow depth can be estimated from the freeboard differences. The benefits of adjusting the orbit of CS-2 for providing more optimized overlaps between IS-2 and CS-2 are considered.

Sea ice↗

Space-based Doppler lidar sampling strategies: Algorithm development and simulated observation experiments

Lidar Atmospheric Wind Sounder (LAWS) Simulation Models (LSM) were developed to evaluate the potential impact of global wind observations on the basic understanding of the Earth's atmosphere and on the predictive skills of current forecast models (GCM and regional scale). Fully integrated top to bottom LAWS Simulation Models for global and regional scale simulations were developed. The algorithm development incorporated the effects of aerosols, water vapor, clouds, terrain, and atmospheric turbulence into the models. Other additions include a new satellite orbiter, signal processor, line of sight uncertainty model, new Multi-Paired Algorithm and wind error analysis code. An atmospheric wind field library containing control fields, meteorological fields, phenomena fields, and new European Center for Medium Range Weather Forecasting (ECMWF) data was also added. The LSM was used to address some key LAWS issues and trades such as accuracy and interpretation of LAWS information, data density, signal strength, cloud obscuration, and temporal data resolution.

Emmitt, G. D.↗

Tolerancing the LITE optical system

This paper describes the optical system used in the Lidar In-Space Technology Experiment (LITE) and presents the results of a study designed to generate a tolerance budget for the LITE, using the combination of manual and Monte Carlo tolerancing techniques utilizing the SYNOPSYS optical analysis program. The tolerance budget derived for LITE, including the contributions of both the fabrication and the alignment of the optical-path elements, are presented. It is shown that the nature of the design allows for typical optical shop tolerances of about 5 mils of despace and decenter, 0.1 deg of tilt, and about 5 fringes of surface figure error.

Wilson, Mark E.↗

Global Monitoring of Clouds and Aerosols Using a Network of Micro-Pulse Lidar Systems

Long-term global radiation programs, such as AERONET and BSRN, have shown success in monitoring column averaged cloud and aerosol optical properties. Little attention has been focused on global measurements of vertically resolved optical properties. Lidar systems are the preferred instrument for such measurements. However, global usage of lidar systems has not been achieved because of limits imposed by older systems that were large, expensive, and logistically difficult to use in the field. Small, eye-safe, and autonomous lidar systems are now currently available and overcome problems associated with older systems. The first such lidar to be developed is the Micro-pulse lidar System (MPL). The MPL has proven to be useful in the field because it can be automated, runs continuously (day and night), is eye-safe, can easily be transported and set up, and has a small field-of-view which removes multiple scattering concerns. We have developed successful protocols to operate and calibrate MPL systems. We have also developed a data analysis algorithm that produces data products such as cloud and aerosol layer heights, optical depths, extinction profiles, and the extinction-backscatter ratio. The algorithm minimizes the use of a priori assumptions and also produces error bars for all data products. Here we present an overview of our MPL protocols and data analysis techniques. We also discuss the ongoing construction of a global MPL network in conjunction with the AERONET program. Finally, we present some early results from the MPL network.

Welton, Ellsworth J.↗

MODIS 3 Km Aerosol Product: Applications over Land in an Urban/suburban Region

MODerate resolution Imaging Spectroradiometer (MODIS) instruments aboard the Terra and Aqua satellites have provided a rich dataset of aerosol information at a 10 km spatial scale. Although originally intended for climate applications, the air quality community quickly became interested in using the MODIS aerosol data. However, 10 km resolution is not sufficient to resolve local scale aerosol features. With this in mind, MODIS Collection 6 is including a global aerosol product with a 3 km resolution. Here, we evaluate the 3 km product over the Baltimore/Washington D.C., USA, corridor during the summer of 2011, by comparing with spatially dense data collected as part of the DISCOVER-AQ campaign these data were measured by the NASA Langley Research Center airborne High Spectral Resolution Lidar (HSRL) and a network of 44 sun photometers (SP) spaced approximately 10 km apart. The HSRL instrument shows that AOD can vary by up to 0.2 within a single 10 km MODIS pixel, meaning that higher resolution satellite retrievals may help to characterize aerosol spatial distributions in this region. Different techniques for validating a high-resolution aerosol product against SP measurements are considered. Although the 10 km product is more statistically reliable than the 3 km product, the 3 km product still performs acceptably, with more than two-thirds of MODIS/SP collocations falling within the expected error envelope with high correlation (R > 0.90). The 3 km product can better resolve aerosol gradients and retrieve closer to clouds and shorelines than the 10 km product, but tends to show more significant noise especially in urban areas. This urban degradation is quantified using ancillary land cover data. Overall, we show that the MODIS 3 km product adds new information to the existing set of satellite derived aerosol products and validates well over the region, but due to noise and problems in urban areas, should be treated with some degree of caution.

MODIS↗

Calibration and Validation of Landsat Tree Cover in the Taiga-Tundra Ecotone

Monitoring current forest characteristics in the taiga-tundra ecotone (TTE) at multiple scales is critical for understanding its vulnerability to structural changes. A 30 m spatial resolution Landsat-based tree canopy cover map has been calibrated and validated in the TTE with reference tree cover data from airborne LiDAR and high resolution spaceborne images across the full range of boreal forest tree cover. This domain-specific calibration model used estimates of forest height to determine reference forest cover that best matched Landsat estimates. The model removed the systematic under-estimation of tree canopy cover greater than 80% and indicated that Landsat estimates of tree canopy cover more closely matched canopies at least 2 m in height rather than 5 m. The validation improved estimates of uncertainty in tree canopy cover in discontinuous TTE forests for three temporal epochs (2000, 2005, and 2010) by reducing systematic errors, leading to increases in tree canopy cover uncertainty. Average pixel-level uncertainties in tree canopy cover were 29.0%, 27.1% and 31.1% for the 2000, 2005 and 2010 epochs, respectively. Maps from these calibrated data improve the uncertainty associated with Landsat tree canopy cover estimates in the discontinuous forests of the circumpolar TTE.

Montesano, Paul Mannix↗

Lidar observation of ozone over Tsukuba (36 deg N, 140 deg E)

An ozone lidar system was installed at the National Institute for Environmental Studies (NIES) in Tsukuba, Japan in March 1988 and has been measuring vertical profiles of ozone (15 - 45 km) since September 1988. The lidar system consists of a XeCl (308 nm) excimer laser, its deuterium Raman shifter (339 nm), a XeF excimer laser (351 nm), a 2 m telescope, a receiving system and a data processing system. The precision of the derived ozone concentration is about 10 percent of an altitude of 40 km for a 4 hr observation. Temperature profiles (30 - 80 km) are also obtained from the Rayleigh scattering signals at 351 nm. Approximate 50 ozone measurements are carried out in a year and variations of vertical profiles of ozone such as seasonal variations and shorter-term variations are observed. Systematic errors due to aerosols had been negligible until the arrival of the stratospheric aerosols injected by the eruption of Mt. Pinatubo. Effects of the volcanic aerosols on ozone measurements depend on the differences between wavelengths used as the on- and off-resonance.

Nakane, H.↗

Surface Irradiances Consistent With CERES-Derived Top-of-Atmosphere Shortwave and Longwave Irradiances

The estimate of surface irradiance on a global scale is possible through radiative transfer calculations using satellite-retrieved surface, cloud, and aerosol properties as input. Computed top-of-atmosphere (TOA) irradiances, however, do not necessarily agree with observation-based values, for example, from the Clouds and the Earth's Radiant Energy System (CERES). This paper presents amethod to determine surface irradiances using observational constraints of TOA irradiance from CERES. A Lagrange multiplier procedure is used to objectively adjust inputs based on their uncertainties such that the computed TOA irradiance is consistent with CERES-derived irradiance to within the uncertainty. These input adjustments are then used to determine surface irradiance adjustments. Observations by the Atmospheric Infrared Sounder (AIRS), Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observations (CALIPSO), CloudSat, andModerate Resolution Imaging Spectroradiometer (MODIS) that are a part of the NASA A-Train constellation provide the uncertainty estimates. A comparison with surface observations from a number of sites shows that the bias [root-mean-square (RMS) difference] between computed and observed monthlymean irradiances calculated with 10 years of data is 4.7 (13.3) W/sq m for downward shortwave and 22.5 (7.1) W/sq m for downward longwave irradiances over ocean and 21.7 (7.8) W m22 for downward shortwave and 21.0 (7.6) W/sq m for downward longwave irradiances over land. The bias andRMS error for the downward longwave and shortwave irradiances over ocean are decreased from those without constraint. Similarly, the bias and RMS error for downward longwave over land improves, although the constraint does not improve downward shortwave over land. This study demonstrates how synergetic use of multiple instruments (CERES,MODIS, CALIPSO, CloudSat, AIRS, and geostationary satellites) improves the accuracy of surface irradiance computations.

Kato, Seiji↗

Asteroid approach covariance analysis for the Clementine mission

The Clementine mission is designed to test Strategic Defense Initiative Organization (SDIO) technology, the Brilliant Pebbles and Brilliant Eyes sensors, by mapping the moon surface and flying by the asteroid Geographos. The capability of two of the instruments available on board the spacecraft, the lidar (laser radar) and the UV/Visible camera is used in the covariance analysis to obtain the spacecraft delivery uncertainties at the asteroid. These uncertainties are due primarily to asteroid ephemeris uncertainties. On board optical navigation reduces the uncertainty in the knowledge of the spacecraft position in the direction perpendicular to the incoming asymptote to a one-sigma value of under 1 km, at the closest approach distance of 100 km. The uncertainty in the knowledge of the encounter time is about 0.1 seconds for a flyby velocity of 10.85 km/s. The magnitude of these uncertainties is due largely to Center Finding Errors (CFE). These systematic errors represent the accuracy expected in locating the center of the asteroid in the optical navigation images, in the absence of a topographic model for the asteroid. The direction of the incoming asymptote cannot be estimated accurately until minutes before the asteroid flyby, and correcting for it would require autonomous navigation. Orbit determination errors dominate over maneuver execution errors, and the final delivery accuracy attained is basically the orbit determination uncertainty before the final maneuver.

Ionasescu, Rodica↗

Mission Simulation of Space Lidar Measurements for Seasonal and Regional CO2 Variations

Results of mission simulation studies are presented for a laser-based atmospheric [82 sounder. The simulations are based on real-time carbon cycle process modeling and data analysis. The mission concept corresponds to the Active Sensing of [82 over Nights, Days, and Seasons (ASCENDS) recommended by the US National Academy of Sciences Decadal Survey of Earth Science and Applications from Space. One prerequisite for meaningful quantitative sensor evaluation is realistic CO2 process modeling across a wide range of scales, i.e., does the model have representative spatial and temporal gradients? Examples of model comparison with data will be shown. Another requirement is a relatively complete description of the atmospheric and surface state, which we have obtained from meteorological data assimilation and satellite measurements from MODIS and [ALIPS0. We use radiative transfer model calculations, an instrument model with representative errors ' and a simple retrieval approach to complete the cycle from "nature" run to "pseudo-data" CO2, Several mission and instrument configuration options are examined/ and the sensitivity to key design variables is shown. We use the simulation framework to demonstrate that within reasonable technological assumptions for the system performance, relatively high measurement precision can be obtained, but errors depend strongly on environmental conditions as well as instrument specifications. Examples are also shown of how the resulting pseudo - measurements might be used to address key carbon cycle science questions.

Kawa, Stephan↗

GEDI Launches a New Era of Biomass Inference from Space

Accurate estimation of aboveground forest biomass stocks is required to assess the impacts of land use changes such as deforestation and subsequent regrowth on concentrations of atmospheric CO 2 . The Global Ecosystem Dynamics Investigation (GEDI) is a lidar mission launched by NASA to the International Space Station in 2018. GEDI was specifically designed to retrieve vegetation structure within a novel, theoretical sampling design that explicitly quantifies biomass and its uncertainty across a variety of spatial scales. In this paper we provide the estimates of pan-tropical and temperate biomass derived from two years of GEDI observations. We present estimates of mean biomass densities at 1 km resolution, as well as estimates aggregated to the national level for every country GEDI observes, and at the sub-national level for the United States. For all estimates we provide the standard error of the mean biomass. These data serve as a baseline for current biomass stocks and their future changes, and the mission’s integrated use of formal statistical inference points the way towards the possibility of a new generation of powerful monitoring tools from space.

Ralph Dubayah↗

Measurements and Modeling of Aerosol Absorption and Single Scattering Albedo at Ambient Relative Hum

Uncertainties in the aerosol single scattering albedo have been identified to be an important source of errors in current large-scale model estimates of the direct aerosol radiative forcing of climate. A number of investigators have obtained estimates of the single scattering albedo from a variety of remote sensing and in situ measurements during aerosol field experiments. During the Tropospheric Aerosol Radiative Forcing Observational Experiment (TARFOX, 1996) for example, estimates of the aerosol single scattering albedo were obtained (1) as a best-fit parameter in comparing radiative flux changes measured by airborne pyranometer to those computed from independently measured aerosol properties; (2) from estimates of the aerosol complex index of refraction derived using a combination of airborne sunphotometer, lidar backscatter and in situ size distribution measurements; and (3) from airborne measurements of aerosol scattering and absorption using nephelometers and absorption photometers. In this paper, we briefly compare the results of the latter two methods for two TARFOX case studies, since those techniques provide height-resolved information about the aerosol single scattering albedo. Estimates of the aerosol single scattering albedo from nephelometer and absorption photometer measurements require knowledge of the scattering and absorption humidification (i.e., the increase in these properties in response to an increase in ambient relative humidity), since both measurements are usually carried out at a relative humidity different from the ambient atmosphere. In principle, the scattering humidification factor can be measured, but there is currently no technique widely available to measure the absorption of an aerosol sample as a function of relative humidity. Frequently, for lack of better knowledge, the absorption humidification is assumed to be unity (meaning that there is no change in aerosol absorption due to an increase in ambient relative humidity). This assumption then enters the estimate of the single scattering albedo at ambient relative humidity. To investigate the validity of this assumption we have carried out modeling studies of the absorption humidification factor, assuming that the aerosols contain an insoluble soot core and a coating which determines its hygroscopic growth behavior. The aerosol optical properties are then computed on the basis of the shell/core particle morphology using a Mie-code for concentric shells. From basic physical principles, it is conceivable that aerosol absorption increases when an atmospheric aerosol particle collects a non-absorbing shell, since the soot core is then exposed to an increased (focused) electric field strength. Indeed, our preliminary modeling studies show that the absorption of an atmospheric aerosol particle composed of a soot core and an aqueous sulfuric acid shell may increase by a factor of 50% due to a change in ambient relative humidity from 30 to 95%. We will show how this increased absorption is a function of the initial particle size and soot mass fraction.

Redemann, J.↗

Performance of the Lidar Design and Data Algorithms for the GLAS Global Cloud and Aerosol Measurements

The Geoscience Laser Altimeter System (GLAS) launched in early 2003 is the first polar orbiting satellite lidar. The instrument design includes high performance observations of the distribution and optical scattering cross sections of atmospheric clouds and aerosol. The backscatter lidar operates at two wavelengths, 532 and 1064 nm. For the atmospheric cloud and aerosol measurements, the 532 nm channel was designed for ultra high efficiency with solid state photon counting detectors and etalon filtering. Data processing algorithms were developed to calibrate and normalize the signals and produce global scale data products of the height distribution of cloud and aerosol layers and their optical depths and particulate scattering cross sections up to the limit of optical attenuation. The paper will concentrate on the effectiveness and limitations of the lidar channel design and data product algorithms. Both atmospheric receiver channels meet and exceed their design goals. Geiger Mode Avalanche Photodiode modules are used for the 532 nm signal. The operational experience is that some signal artifacts and non-linearity require correction in data processing. As with all photon counting detectors, a pulse-pile-up calibration is an important aspect of the measurement. Additional signal corrections were found to be necessary relating to correction of a saturation signal-run-on effect and also for daytime data, a small range dependent variation in the responsivity. It was possible to correct for these signal errors in data processing and achieve the requirement to accurately profile aerosol and cloud cross section down to 10-7 llm-sr. The analysis procedure employs a precise calibration against molecular scattering in the mid-stratosphere. The 1064 nm channel detection employs a high-speed analog APD for surface and atmospheric measurements where the detection sensitivity is limited by detector noise and is over an order of magnitude less than at 532 nm. A unique feature of the GLAS is a full acquisition of the surface return pulse, which has important application to the atmospheric transmission retrieval.

Spinhirne, James D.↗

MIR Measurements of Atmospheric Water Vapor Profiles

Three subjects related to atmospheric water vapor profiling using the 183.3 GHz absorption line are discussed in this paper. First, data acquired by an airborne millimeter-wave imaging radiometer (MIR) over ocean surface in the western Pacific are used to estimate three-dimensional (3-D) distribution of atmospheric water vapor. The instrument's radiometric measurements with mixed vertical and horizontal polarizations require modifications to the retrieval algorithm used in the past. It is demonstrated that, after the modifications, the new algorithm can provide adequate retrieval of water vapor profiles, even though the measured data are of mixed polarizations. Next, the retrieved profiles, in terms of water vapor mixing ratio rho (g/kg), are compared with those measured in near concurrence by dropsondes from a research aircraft in the western Pacific and by a ground-based Raman lidar at Wallops Island, Virginia. The ratio of the standard deviation to the mean rho is found to be 0.12 at 0.25 km altitude and gradually degraded to 0.67 at the highest altitude of the retrieval of 10.25 km. Finally, the effect of the "initial guess" relative humidity profile on the final retrieved product is analyzed with respect to the condition for the convergent retrieval. It is found that the effect is minimal if the initial profile is not unrealistically different from the true one. If the initial profile is very different from the true one, the final retrieved product could be subject to a substantial error. Tightening of the convergent condition in the retrieval helped reduce magnitude of the error, but not remove it totally. It is concluded that an initial profile based on climatology is likely to provide most reliable retrieval results.

Wang, J. R.↗

Remote Sensing of Cloud Top Heights Using the Research Scanning Polarimeter

Clouds cover roughly two thirds of the globe and act as an important regulator of Earth's radiation budget. Of these, multilayered clouds occur about half of the time and are predominantly two-layered. Changes in cloud top height (CTH) have been predicted by models to have a globally averaged positive feedback, however observational changes in CTH have shown uncertain results. Additional CTH observations are necessary to better and quantify the effect. Improved CTH observations will also allow for improved sub-grid parameterizations in large-scale models and accurate CTH information is important when studying variations in freezing point and cloud microphysics. NASA's airborne Research Scanning Polarimeter (RSP) is able to measure cloud top height using a novel multi-angular contrast approach. RSP scans along the aircraft track and obtains measurements at 152 viewing angles at any aircraft location. The approach presented here aggregates measurements from multiple scans to a single location at cloud altitude using a correlation function designed to identify the location-distinct features in each scan. During NASAs SEAC4RS air campaign, the RSP was mounted on the ER-2 aircraft along with the Cloud Physics Lidar (CPL), which made simultaneous measurements of CTH. The RSPs unique method of determining CTH is presented. The capabilities of using single and combinations of channels within the approach are investigated. A detailed comparison of RSP retrieved CTHs with those of CPL reveal the accuracy of the approach. Results indicate a strong ability for the RSP to accurately identify cloud heights. Interestingly, the analysis reveals an ability for the approach to identify multiple cloud layers in a single scene and estimate the CTH of each layer. Capabilities and limitations of identifying single and multiple cloud layers heights are explored. Special focus is given to sources of error in the method including optically thin clouds, physically thick clouds, multi-layered clouds as well as cloud phase. When determining multi-layered CTHs, limits on the upper clouds opacity are assessed.

Cloud height indicators↗

Deriving High-Resolution Reservoir Bathymetry from ICESat-2 Prototype Photon-Counting Lidar and Landsat Imagery

Knowledge of reservoir bathymetry is essential for many studies on terrestrial hydrological and biogeochemical processes. However, there are currently no cost-effective approaches to derive reservoir bathymetry at the global scale. This study explores the potential of generating high-resolution global bathymetry using elevation data collected by the 532-nm Advanced Topographic Laser Altimeter System (ATLAS) onboard the Ice, Cloud, and Land Elevation Satellite (ICESat-2). The novel algorithm was developed and tested using the ICESat-2 airborne prototype, the Multiple Altimeter Beam Experimental Lidar (MABEL), with Landsat-based water classifications (from 1982 to 2017). MABEL photon elevations were paired with Landsat water occurrence percentiles to establish the elevation-area (E-A) relationship, which in turn was applied to the percentile image to obtain partial bathymetry over the historic dynamic range of reservoir area. The bathymetry for the central area was projected to achieve the full bathymetry. The bathymetry image was then embedded onto the digital elevation model (DEM). Results were validated over Lake Mead against survey data. Results over four transects show coefficient of determination (R²) values from 0.82 to 0.99 and root-mean-square error (RMSE) values from 1.18 to 2.36 m. In addition, the E-A and elevation-storage (E-S) curves have RMSEs of 1.56 m and 0.08 km³, respectively. Over the entire dynamic reservoir area, the derived bathymetry agrees very well with independent survey data, except for within the highest and lowest percentile bands. With abundant overpassing tracks and high spatial resolution, the newly launched ICESat-2 should enable the derivation of bathymetry over an unprecedented number of reservoirs.

Yao Li↗