Advanced Silicon Avalanche Photodiodes on NASA's Global Ecosystem Dynamics Investigation (GEDI) Mission

Silicon Avalanche Photodiodes (APDs) are used in NASAs Global Ecosystem Dynamics Investigation (GEDI) which was launched in December 2018 and is currently measuring the Earths vegetation vertical structure from the International Space Station. The APDs were specially made for space lidar with a much lower hole-to-electron ionization coefficient ratio (k-factor ~0.008) than that of commercially available silicon APDs in order to reduce the APD excess noise from the randomness of the avalanche gain. A silicon heater resistor was used under the APD chip to heat the device up to 70C and improve its quantum efficiency at 1064 nm laser wavelength while maintaining a low dark current such that the overall signal to noise ratio is improved. Special APD protection circuits were used to raise the overload damage threshold to prevent device damage from strong laser return by specular surfaces, such as still water bodies, and space radiation events. The APD and a hybrid transimpedance amplifier circuit were hermetically sealed in a package with a sufficiently low leak rate to ensure multi-year operation lifetime in space. The detector assemblies underwent a series of pre-launch tests per NASA Goddard Environmental Verification Standard for space qualification. They have performed exactly as expected with GEDI in orbit. A detailed description of the GEDI detector design, signal and noise model, and test results are presented in this paper

ICESat GLAS Altimetry Measurements: Received Signal Dynamic Range and Saturation Correction

NASAs Ice, Cloud, and land Elevation Satellite (ICESat), which operated between 2003 and 2009, made the first satellite-based global lidar measurement of Earths ice sheet elevations, sea-ice thickness and vegetation canopy structure. The primary instrument on ICESat was the Geoscience Laser Altimeter System (GLAS), which measured the distance from the spacecraft to Earths surface via the roundtrip travel time of individual laser pulses. GLAS utilized pulsed lasers and a direct detection receiver consisting of a silicon avalanche photodiode (SiAPD) and a waveform digitizer. Early in the mission, the peak power of the received signal from snow and ice surfaces was found to span a wider dynamic range than planned, often exceeding the linear dynamic range of the GLAS 1064-nm detector assembly. The resulting saturation of the receiver distorted the recorded signal and resulted in range biases as large as 50 cm for ice and snow-covered surfaces. We developed a correction for this saturation range bias based on laboratory tests using a spare flight detector, and refined the correction by comparing GLAS elevation estimates to those derived from Global Positioning System (GPS) surveys over the calibration site at the salar de Uyuni, Bolivia. Applying the saturation correction largely eliminated the range bias due to receiver saturation for affected ICESat measurements over Uyuni and significantly reduced the discrepancies at orbit crossovers located on flat regions of the Antarctic ice sheet

Space Flight LiDARs, Navigation & Science Instrument Implementations: Lasers, Optoelectronics, Integrated Photonics, Fiber Optic Subsystems and Components

For the past 25 years, the National Aeronautics and Space Administration (NASA) Goddard Space Flight Center's Photonics Group in the Engineering Directorate has been substantially contributing to the flight design, development, production, testing and integration of many science and navigational instruments. The Moon to Mars initiative will rely heavily upon utilizing commercial technologies for instrumentation with aggressive schedule deadlines. The group has an extensive background in screening, qualifying, development and integration of commercial components for spaceflight applications. By remaining adaptable and employing a rigorous approach to component and instrument development, they have forged and fostered relationships with industry partners. They have been willing to communicate lessons learned in packaging, part construction, materials selection, testing, and other facets of the design and production process critical to implementation for high-reliability systems. As a result, this successful collaboration with industry vendors and component suppliers has enabled a history of mission success from the Moon to Mars (and beyond) while balancing cost, schedule, and risk postures. In cases where no commercial components exist, the group works closely with other teams at Goddard Space Flight Center and other NASA field centers to fabricate and produce flight hardware for science, remote sensing, and navigation applications. Summarized here is the last ten years of instrumentation development lessons learned and data collected from the subsystems down to the optoelectronic component level

The Ice, Cloud, and Land Elevation Satellite-2 (ICESat-2): Science Requirements, Concept, and Implementation

The Ice, Cloud, and land Elevation Satellite (ICESat) mission used laser altimetry measurements to determine changes in elevations of glaciers and ice sheets, as well as sea ice thickness distribution. These measurements have provided important information on the response of the cryosphere (Earths frozen surfaces) to changes in atmosphere and ocean condition. ICESat operated from 2003-2009 and provided repeat altimetry measurements not only to the cryosphere scientific community but also to the ocean, terrestrial and atmospheric scientific communities. The conclusive assessment of significant ongoing rapid changes in the Earths ice cover, in part supported by ICESat observations, has strengthened the need for sustained, high accuracy, repeat observations similar to what was provided by the ICESat mission. Following recommendations from the National Research Council for an ICESat follow-on mission, the ICESat-2 mission is now under development for planned launch in 2018. The primary scientific aims of the ICESat-2 mission are to continue measurements of sea ice freeboard and ice sheet elevation to determine their changes at scales from outlet glaciers to the entire ice sheet, and from 10s of meters to the entire polar oceans for sea ice freeboard. ICESat carried a single beam profiling laser altimeter that produced approximately 70 m diameter footprints on the surface of the Earth at approximately 150 m along-track intervals. In contrast, ICESat-2 will operate with three pairs of beams, each pair separated by about 3 km across-track with a pair spacing of 90 m. Each of the beams will have a nominal 17 m diameter footprint with an along-track sampling interval of 0.7 m. The differences in the ICESat-2 measurement concept are a result of overcoming some limitations associated with the approach used in the ICESat mission. The beam pair configuration of ICESat-2 allows for the determination of local cross-track slope, a significant factor in measuring elevation change for the outlet glaciers surrounding the Greenland and Antarctica coasts. The multiple beam pairs also provide improved spatial coverage. The dense spatial sampling eliminates along-track measurement gaps, and the small footprint diameter is especially useful for sea surface height measurements in the often narrow leads needed for sea ice freeboard and ice thickness retrievals. The ICESat-2 instrumentation concept uses a low energy 532 nm (green) laser in conjunction with single-photon sensitive detectors to measure range. Combining ICESat-2 data with altimetry data collected since the start of the ICESat mission in 2003, such as Operation IceBridge and ESAs CryoSat-2, will yield a 15+ year record of changes in ice sheet elevation and sea ice thickness. ICESat-2 will also provide information of mountain glacier and ice cap elevations changes, land and vegetation heights, inland water elevations, sea surface heights, and cloud layering and optical thickness

W42 - a scalable spatial database system for holding Digital Elevation Models

The design of a scalable system for holding spatial data in general and digital elevation models (DEMs) in specific has to account for the characteristics of data from various application fields. The data can be heterogeneous in coverage, as well as in resolution, information content and quality. A database aiming at the representation of world-wide DEMs has to consider these differences in the design of the system with respect to the structure and the algorithms. The database system W42, which is presented in the work at hand, is a scalable spatial database system capable of holding, extracting, mosaicking, and fusing spatial data represented in raster- as well as in vector-format. Design aspects for this task can be specified as holding spatial data in unique data structures and providing unique access functions to the data. These are subject of this work as well as first experiences gained from the implementation of part of the extensions made for the TanDEM-X mission

Elevation and Deformation Extraction from TomoSAR

3D SAR tomography (TomoSAR) and 4D SAR differential tomography (Diff-TomoSAR) exploit multi-baseline SAR data stacks to provide an essential innovation of SAR Interferometry for many applications, sensing complex scenes with multiple scatterers mapped into the same SAR pixel cell. However, these are still influenced by DEM uncertainty, temporal decorrelation, orbital, tropospheric and ionospheric phase distortion and height blurring. In this thesis, these techniques are explored. As part of this exploration, the systematic procedures for DEM generation, DEM quality assessment, DEM quality improvement and DEM applications are first studied. Besides, this thesis focuses on the whole cycle of systematic methods for 3D & 4D TomoSAR imaging for height and deformation retrieval, from the problem formation phase, through the development of methods to testing on real SAR data. After DEM generation introduction from spaceborne bistatic InSAR (TanDEM-X) and airborne photogrammetry (Bluesky), a new DEM co-registration method with line feature validation (river network line, ridgeline, valley line, crater boundary feature and so on) is developed and demonstrated to assist the study of a wide area DEM data quality. This DEM co-registration method aligns two DEMs irrespective of the linear distortion model, which improves the quality of DEM vertical comparison accuracy significantly and is suitable and helpful for DEM quality assessment. A systematic TomoSAR algorithm and method have been established, tested, analysed and demonstrated for various applications (urban buildings, bridges, dams) to achieve better 3D & 4D tomographic SAR imaging results. These include applying Cosmo-Skymed X band single-polarisation data over the Zipingpu dam, Dujiangyan, Sichuan, China, to map topography; and using ALOS L band data in the San Francisco Bay region to map urban building and bridge. A new ionospheric correction method based on the tile method employing IGS TEC data, a split-spectrum and an ionospheric model via least squares are developed to correct ionospheric distortion to improve the accuracy of 3D & 4D tomographic SAR imaging. Meanwhile, a pixel by pixel orbit baseline estimation method is developed to address the research gaps of baseline estimation for 3D & 4D spaceborne SAR tomography imaging. Moreover, a SAR tomography imaging algorithm and a differential tomography four-dimensional SAR imaging algorithm based on compressive sensing, SAR interferometry phase (InSAR) calibration reference to DEM with DEM error correction, a new phase error calibration and compensation algorithm, based on PS, SVD, PGA, weighted least squares and minimum entropy, are developed to obtain accurate 3D & 4D tomographic SAR imaging results. The new baseline estimation method and consequent TomoSAR processing results showed that an accurate baseline estimation is essential to build up the TomoSAR model. After baseline estimation, phase calibration experiments (via FFT and Capon method) indicate that a phase calibration step is indispensable for TomoSAR imaging, which eventually influences the inversion results. A super-resolution reconstruction CS based study demonstrates X band data with the CS method does not fit for forest reconstruction but works for reconstruction of large civil engineering structures such as dams and urban buildings. Meanwhile, the L band data with FFT, Capon and the CS method are shown to work for the reconstruction of large manmade structures (such as bridges) and urban buildings

WORLDDEM – A NOVEL GLOBAL FOUNDATION LAYER

Airbus Defence and Space's WorldDEM™ provides a global Digital Elevation Model of unprecedented quality, accuracy, and coverage. The product will feature a vertical accuracy of 2m (relative) and better than 6m (absolute) in a 12m x 12m raster. The accuracy will surpass that of any global satellite-based elevation model available. WorldDEM is a game-changing disruptive technology and will define a new standard in global elevation models. The German radar satellites TerraSAR-X and TanDEM-X form a high-precision radar interferometer in space and acquire the data basis for the WorldDEM. This mission is performed jointly with the German Aerospace Center (DLR). Airbus DS refines the Digital Surface Model (e.g. editing of acquisition, processing artefacts and water surfaces) or generates a Digital Terrain Model. Three product levels are offered: WorldDEMcore (output of the processing, no editing is applied), WorldDEM™ (guarantees a void-free terrain description and hydrological consistency) and WorldDEM DTM (represents bare Earth elevation). Precise elevation data is the initial foundation of any accurate geospatial product, particularly when the integration of multi-source imagery and data is performed based upon it. Fused data provides for improved reliability, increased confidence and reduced ambiguity. This paper will present the current status of product development activities including methodologies and tool to generate these, like terrain and water bodies editing and DTM generation. In addition, the studies on verification & validation of the WorldDEM products will be presented

White paper – On the use of LiDAR data at AmeriFlux sites

Our aim is to inform the AmeriFlux community on existing and upcoming LiDAR technologies (atmospheric Doppler or Raman LiDAR often deployed at flux sites are not considered here), how it is currently used at flux sites, and how we believe it could, in the future, further contribute to the AmeriFlux vision. Heterogeneity in vegetation and ground properties at various spatial scales is omnipresent at flux sites, and 3D mapping of canopy, understory, and ground surface can help move the science forward

Post-drought decline of the Amazon carbon sink

Amazon forests have experienced frequent and severe droughts in the past two decades. However, little is known about the large-scale legacy of droughts on carbon stocks and dynamics of forests. Using systematic sampling of forest structure measured by LiDAR waveforms from 2003 to 2008, here we show a significant loss of carbon over the entire Amazon basin at a rate of 0.3 ± 0.2 (95% CI) PgC yr−1 after the 2005 mega-drought, which continued persistently over the next 3 years (2005–2008). The changes in forest structure, captured by average LiDAR forest height and converted to above ground biomass carbon density, show an average loss of 2.35 ± 1.80 MgC ha−1 a year after (2006) in the epicenter of the drought. With more frequent droughts expected in future, forests of Amazon may lose their role as a robust sink of carbon, leading to a significant positive climate feedback and exacerbating warming trends.The research was partially supported by NASA Terrestrial Ecology grant at the Jet Propulsion Laboratory, California Institute of Technology and partial funding to the UCLA Institute of Environment and Sustainability from previous National Aeronautics and Space Administration and National Science Foundation grants. The authors thank NSIDC, BYU, USGS, and NASA Land Processes Distributed Active Archive Center (LP DAAC) for making their data available. (NASA Terrestrial Ecology grant at the Jet Propulsion Laboratory, California Institute of Technology)Published versio

Biomass estimation from simulated GEDI, ICESat-2 and NISAR across environmental gradients in Sonoma County, California

Estimates of the magnitude and distribution of aboveground carbon in Earth's forests remain uncertain, yet knowledge of forest carbon content at a global scale is critical for forest management in support of climate mitigation. In light of this knowledge gap, several upcoming spaceborne missions aim to map forest aboveground biomass, and many new biomass products are expected from these datasets. As these new missions host different technologies, each with relative strengths and weaknesses for biomass retrieval, as well as different spatial resolutions, consistently comparing or combining biomass estimates from these new datasets will be challenging. This paper presents a demonstration of an inter-comparison of biomass estimates from simulations of three NASA missions (GEDI, ICESat-2 and NISAR) over Sonoma county in California, USA. We use a high resolution, locally calibrated airborne lidar map as our reference dataset, and emphasize the importance of considering uncertainties in both reference maps and spaceborne estimates when conducting biomass product validation. GEDI and ICESat-2 were simulated from airborne lidar point clouds, while UAVSAR's L-band backscatter was used as a proxy for NISAR. To estimate biomass for the lidar missions we used GEDI's footprint-level biomass algorithms, and also adapted these for application to ICESat-2. For UAVSAR, we developed a locally trained biomass model, calibrated against the ALS reference map. Each mission simulation was evaluated in comparison to the local reference map at its native product resolution (25 m, 100 m transect, and 1 ha) yielding RMSEs of 57%, 75%, and 89% for GEDI, NISAR, and ICESat-2 respectively. RMSE values increased for GEDI's power beam during simulated daytime conditions (64%), coverage beam during nighttime conditions (72%), and coverage beam daytime conditions (87%). We also test the application of GEDI's biomass modeling framework for estimation of biomass from ICESat-2, and find that ICESat-2 yields reasonable biomass estimates, particularly in relatively short, open canopies. Results suggest that while all three missions will produce datasets useful for biomass mapping, tall, dense canopies such as those found in Sonoma County present the greatest challenges for all three missions, while steep slopes also prove challenging for single-date SAR-based biomass retrievals. Our methods provide guidance for the inter-comparison and validation of spaceborne biomass estimates through the use of airborne lidar reference maps, and could be repeated with on-orbit estimates in any area with high quality field plot and ALS data. These methods allow for regional interpretations and filtering of multi-mission biomass estimates toward improved wall-to-wall biomass maps through data fusion.</p