SUN-DISSERTATION-2019.Pdf (5.590Mb)
Total Page:16
File Type:pdf, Size:1020Kb
© Copyright by Zhiyue Sun, 2019 All Rights Reserved ESTIMATING GLACIER ICE THICKNESS AND MASS BALANCE USING MULTIPLE SPACEBORNE GEODETIC TECHNIQUES A Dissertation Presented to the Faculty of the Department of Civil and Environmental Engineering University of Houston In Partial Fulfillment of the Requirements for the Degree Doctor of Philosophy in Geosensing Systems Engineering and Sciences by Zhiyue Sun May 2019 ESTIMATING GLACIER ICE THICKNESS AND MASS BALANCE USING MULTIPLE SPACEBORNE GEODETIC TECHNIQUES __________________________ Zhiyue Sun Approved: ________________________________ Chair of the Committee Hyongki Lee, Ph.D. Associate Professor Civil & Environmental Engineering Committee Members: ________________________________ Ramesh L. Shrestha, Ph.D. Hugh Roy & Lillie Cranz Cullen Distinguished Professor Civil & Environmental Engineering ________________________________ Craig L. Glennie, Ph.D. Associate Professor Civil & Environmental Engineering ________________________________ Yushin. Ahn, Ph.D. Assistant Professor Geomatics Engineering California State University, Fresno, CA ________________________________ Alexander. Braun, Ph.D. Professor Geological Sciences and Geological Engineering Queen’s University, Kingston, Canada ________________________________ ________________________________ Suresh K. Khator, Ph.D. Craig L. Glennie, Ph.D. Associate Dean Geosensing Systems Engineering & Cullen College of Engineering Science Program Director ESTIMATING GLACIER ICE THICKNESS AND MASS BALANCE USING MULTIPLE SPACEBORNE GEODETIC TECHNIQUES An Abstract of a Dissertation Presented to the Faculty of the Department of Civil and Environmental Engineering University of Houston In Partial Fulfillment of the Requirements for the Degree Doctor of Philosophy in Geosensing Systems Engineering and Sciences by Zhiyue Sun May 2019 Abstract The melting of glaciers and ice caps makes a significant contribution to present- day sea level rise. Recently, multiple space-borne remote sensing techniques have been successfully used to obtain geodetic observations, which were used to estimate the contribution of global glacier melt to sea level rise. However, estimates from different approaches yield large discrepancies in certain glacierized areas such as the Eastern Nyainqen Tanglha. Moreover, other glacierized regions such as the Novaya Zemlya in the Russian Arctic have been under-studied. Therefore, to characterize and quantify accurate glacier mass balance estimates over Eastern Nyainqen Tanglha and Novaya Zemlya, an iterative velocity-based method is proposed to estimate glacier thickness. First, a new iterative method is presented, which estimates ice thickness using surface velocity and surface topography. The temperature-related rate factor for temperate glaciers is empirically obtained based on multiple in-situ measurements, while the rate factor for non-temperate glaciers follows the assumption made by previous studies. A validation was performed with 15 previous methods over 8 glaciers. Based on the comparative results, the proposed method in glacierized areas where direct observations are limited is promising. Second, multiple traditional spaceborne techniques for observing surface mass balance were tested in the Novaya Zemlya. The mass variation trend was obtained based on observations from Gravity Recovery and Climate Experiment (GRACE) and radar altimeter over 4 selected glaciers. The glacier outflow and influx were also determined using surface velocity and snowfall observations. Additionally, the contributions of outflow and influx to mass change for Novaya Zemlya were identified. vi Finally, the proposed thickness estimate method is applied to 4 largest glaciers in Nyainqen Tanglha. Surface velocity estimates were calculated from Advanced Land Observing Satellite (ALOS)/ALOS-2 Phased Array type L-band SAR (PALSAR) image pairs via speckle matching. Surface slope distribution was calculated from the Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) Digital Elevation Model (DEM) data. Consequently, the thickness distribution map from 2008 to 2016 was obtained. Furthermore, the mass balance and thickness changes were determined. These estimates have an agreement with previous ones from GRACE and ASTER DEM differencing approaches. vii TABLE OF CONTENTS ABSTRACT .............................................................................................................................................. VI TABLE OF CONTENTS ...................................................................................................................... VIII LIST OF FIGURES .................................................................................................................................. XI LIST OF TABLES ................................................................................................................................... XV 1 INTRODUCTION ................................................................................................................................... 1 1.1 BACKGROUNDS ................................................................................................................................ 2 1.1.1 Glacier Studies Using In-Situ Measurements ...................................................................... 2 1.1.2 Remote Sensing of Glaciers ................................................................................................. 3 1.2 MOTIVATION AND OBJECTIVES ........................................................................................................ 5 1.2.1 Motivation ............................................................................................................................ 5 1.2.2 Objectives and Contribution ................................................................................................ 7 1.3 DISSERTATION STRUCTURE .............................................................................................................. 8 2 THEORETICAL BACKGROUND ....................................................................................................... 9 2.1 SPECKLE MATCHING & FEATURE TRACKING ................................................................................... 9 2.2 PRINCIPLE OF RADAR ALTIMETRY ...................................................................................................12 2.3 GRACE DATA PROCESSING ............................................................................................................14 2.4 ERA-INTERIM REANALYSIS DATA PROCESSING .............................................................................19 3 NEW VELOCITY-BASED GLACIER THICKNESS ESTIMATION METHOD ..........................21 3.1 BACKGROUND .................................................................................................................................21 3.2 STUDY SITE & DATA .......................................................................................................................24 3.3 METHODOLOGY ...............................................................................................................................25 3.3.1 Equations for Thickness Estimates .....................................................................................25 3.3.2 Iterative Process to Determine 흓 .......................................................................................30 3.3.3 Optimal Englacial Temperature .........................................................................................32 3.4 RESULT ............................................................................................................................................39 3.4.1 Initial Inputs .......................................................................................................................39 viii 3.4.2 Comparison with Existing Methods ....................................................................................41 3.4.3 Selection of Contour Interval ..............................................................................................48 3.4.4 Error Analysis.....................................................................................................................50 3.5 CONCLUSION ...................................................................................................................................53 4 TRADITIONAL TECHNIQUE OF ESTIMATING ICE LOSS IN GLACIERS: CASE STUDY IN NOVAYA ZEMLYA. ................................................................................................................................55 4.1 BACKGROUND & DATA ...................................................................................................................55 4.2 TECHNIQUES & DATA ......................................................................................................................59 4.2.1 Ice Mass Changes from GRACE .........................................................................................59 4.2.2 Elevation Changes from Envisat RA-2 Radar Altimetry .....................................................60 4.2.3 Glacier Velocity from PALSAR Speckle-matching .............................................................62 4.2.4 Glacier Velocity from Landsat Feature Offset-tracking .....................................................63 4.2.5 Atmospheric and Oceanic Observations from ERA-Interim ...............................................64 4.3 RESULTS AND DISCUSSION ...............................................................................................................65