scholarly journals An 11-year record of mass balance of Brewster Glacier, New Zealand, determined using a geostatistical approach

2016 ◽  
Vol 63 (238) ◽  
pp. 199-217 ◽  
Author(s):  
NICOLAS J. CULLEN ◽  
BRIAN ANDERSON ◽  
PASCAL SIRGUEY ◽  
DOROTHEA STUMM ◽  
ANDREW MACKINTOSH ◽  
...  

ABSTRACTRecognising the scarcity of glacier mass-balance data in the Southern Hemisphere, a mass-balance measurement programme was started at Brewster Glacier in the Southern Alps of New Zealand in 2004. Evolution of the measurement regime over the 11 years of data recorded means there are differences in the spatial density of data obtained. To ensure the temporal integrity of the dataset a new geostatistical approach is developed to calculate mass balance. Spatial co-variance between elevation and snow depth allows a digital elevation model to be used in a co-kriging approach to develop a snow depth index (SDI). By capturing the observed spatial variability in snow depth, the SDI is a more reliable predictor than elevation and is used to adjust each year of measurements consistently despite variability in sampling spatial density. The SDI also resolves the spatial structure of summer balance better than elevation. Co-kriging is used again to spatially interpolate a derived mean summer balance index using SDI as a co-variate, which yields a spatial predictor for summer balance. The average glacier-wide surface winter, summer and annual balances over the period 2005–15 are 2484, −2586 and −102 mm w.e., respectively, with changes in summer balance explaining most of the variability in annual balance.

2019 ◽  
Vol 65 (251) ◽  
pp. 422-439 ◽  
Author(s):  
KUNPENG WU ◽  
SHIYIN LIU ◽  
ZONGLI JIANG ◽  
JUNLI XU ◽  
JUNFENG WEI

ABSTRACTTo obtain information on changes in glacier mass balance in the central Nyainqentanglha Range, a comprehensive study was carried out based on digital-elevation models derived from the 1968 topographic maps, the Shuttle Radar Topography Mission DEM (2000) and TerraSAR-X/TanDEM-X (2013). Glacier area changes between 1968 and 2016 were derived from topographic maps and Landsat OLI images. This showed the area contained 715 glaciers, with an area of 1713.42 ± 51.82 km2, in 2016. Ice cover has been shrinking by 0.68 ± 0.05% a−1 since 1968. The glacier area covered by debris accounted for 11.9% of the total and decreased in the SE–NW directions. Using digital elevation model differencing and differential synthetic aperture radar interferometry, a significant mass loss of 0.46 ± 0.10 m w.e. a−1 has been recorded since 1968; mass losses accelerated from 0.42 ± 0.20 m w.e. a−1 to 0.60 ± 0.20 m w.e. a−1 between 1968–2000 and 2000–2013, with thinning noticeably greater on the debris-covered ice than the clean ice. Surface-elevation changes can be influenced by ice cliffs, as well as debris cover and land- or lake-terminating glaciers. Changes showed spatial and temporal heterogeneity and a substantial correlation with climate warming and decreased precipitation.


2011 ◽  
Vol 65 ◽  
pp. 214-217
Author(s):  
Yao Ge Wang ◽  
Peng Yuan Wang

Interpolation is the core problem of Digital Elevation Model (DEM). The Coons DEM model is better than bilinear interpolation and moving surface fitting. It is constructed by grid boundary curve, the curve interpolates by some adjoining grid points. Its spatial pattern of error is random in global area, there is no significant global spatial autocorrelation, but it is an increasing trend along with the terrain average gradient increases.There is significant local spatial autocorrelation, the spatial pattern of error converges strongly in local areas.


2020 ◽  
Vol 12 (14) ◽  
pp. 2297 ◽  
Author(s):  
Mingyang Lv ◽  
Huadong Guo ◽  
Jin Yan ◽  
Kunpeng Wu ◽  
Guang Liu ◽  
...  

The Karakoram has had an overall slight positive glacier mass balance since the end of 20th century, which is anomalous given that most other regions in High Mountain Asia have had negative changes. A large number of advancing, retreating, and surging glaciers are heterogeneously mixed in the Karakoram increasing the difficulties and inaccuracies to identify glacier surges. We found two adjacent glaciers in the eastern Karakoram behaving differently from 1995 to 2019: one was surging and the other was advancing. In order to figure out the differences existing between them and the potential controls on surges in this region, we collected satellite images from Landsat series, ASTER, and Google Earth, along with two sets of digital elevation model. Utilizing visual interpretation, feature tracking of optical images, and differencing between digital elevation models, three major differences were observed: (1) the evolution profiles of the terminus positions occupied different change patterns; (2) the surging glacier experienced a dramatic fluctuation in the surface velocities during and after the event, while the advancing glacier flowed in a stable mode; and (3) surface elevation of the surging glacier decreased in the reservoir and increased in the receiving zone. However, the advancing glacier only had an obvious elevation increase over its terminus part. These differences can be regarded as standards for surge identification in mountain ranges. After combining the differences with regional meteorological conditions, we suggested that changes of thermal and hydrological conditions could play a role in the surge occurrence, in addition, geomorphological characteristics and increasing warming climate might also be part of it. This research strongly contributes to the literatures of glacial motion and glacier mass change in the eastern Karakoram through remote sensing.


2000 ◽  
Vol 31 ◽  
pp. 164-170 ◽  
Author(s):  
Michael P. Bishop ◽  
Jeffrey S. Kargel ◽  
Hugh H. Kieffer ◽  
David J. MacKinnon ◽  
Bruce H. Raup ◽  
...  

AbstractA large number of multispectral and stereo-image data are expected to become available as part of the Global Land Ice Measurements from Space project. We investigate digital elevation model extraction, anisotropic reflectance correction and selected glacier analysis tasks that must be developed to achieve full utility of these new data. Results indicate that glaciers in the Karakoram and Nanga Parbat Himalaya, northern Pakistan, exhibit unique spectral, spatial and geomorphometric patterns that can be exploited by various models and algorithms to produce accurate information regarding glacier extent, supraglacial features and glacier geomorphology The integration of spectral, spatial and geomorphometric features, coupled with approaches for advanced pattern recognition, can help geoscientists study glacier mass balance, glacier erosion, sediment-transfer efficiency and landscape evolution.


2013 ◽  
Vol 54 (63) ◽  
pp. 41-51 ◽  
Author(s):  
J.P. Conway ◽  
N.J. Cullen

AbstractThe turbulent sensible and latent heat fluxes are important components of the surface energy balance over glaciers in the Southern Alps of New Zealand, contributing over half the energy available for ablation during large melt events. To calculate these terms confidently in glacier mass-balance models it is essential to use appropriate parameterizations for surface roughness and atmospheric stability. Eddy covariance measurements at Brewster Glacier were obtained over an ice surface to help facilitate an assessment of the calculation of the turbulent heat fluxes. The roughness length for momentum was found to be 3.6 x 10−3m, while the roughness lengths for temperature and humidity were two orders of magnitude smaller, in agreement with surface renewal theory. A Monte Carlo approach was used to assess the uncertainty in turbulent heat fluxes calculated using the bulk aerodynamic method. It was found that input-data and roughness-length uncertainty could not explain underestimates of observed sensible heat fluxes during periods with low wind speed and large temperature gradients. During these periods a katabatic wind speed maximum alters the formulation of the turbulent exchange coefficient to that typically observed in a neutral atmosphere and this has implications for glacier mass-balance sensitivity.


2015 ◽  
Vol 11 (2) ◽  
pp. 603-636 ◽  
Author(s):  
C. Bravo ◽  
M. Rojas ◽  
B. M. Anderson ◽  
A. N. Mackintosh ◽  
E. Sagredo ◽  
...  

Abstract. Glacier behaviour during the mid-Holocene (MH, 6000 year BP) in the Southern Hemisphere provides observational data to constrain our understanding of the origin and propagation of palaeo-climatic signals. We examine the climatic forcing of glacier expansion in the MH by evaluating modelled glacier equilibrium line altitude (ELA) and climate conditions during the MH compared with pre-industrial time (PI, year 1750) in the mid latitudes of the Southern Hemisphere, specifically in Patagonia and the South Island of New Zealand. Climate conditions for the MH are obtained from PMIP2 models simulations, which in turn force a simple glacier mass balance model to simulate changes in equilibrium-line altitude during this period. Climate conditions during the MH show significantly (p ≤ 0.05) colder temperatures in summer, autumn and winter, and significantly (p ≤ 0.05) warmer temperatures in spring. These changes are a consequence of insolation differences between the two periods. Precipitation does not show significant changes, but exhibits a temporal pattern with less precipitation from August to September and more precipitation from October to April during the MH. In response to these climatic changes, glaciers in both analysed regions have an ELA that is 15–33 m lower than PI during the MH. The main causes of this difference are the colder temperature during the MH, reinforcing previous results that mid-latitude glaciers are more sensitive to temperature change compared to precipitation changes. Differences in temperature have a dual effect on mass balance. First, during summer and early autumn less energy is available for melting. Second in late autumn and winter, lower temperatures cause more precipitation to fall as snow rather than rain, resulting in more accumulation and higher surface albedo. For these reasons, we postulate that the modelled ELA changes, although small, may help to explain larger glacier extents observed in the mid Holocene in both South America and New Zealand.


2013 ◽  
Vol 7 (1) ◽  
pp. 55-101
Author(s):  
T. Goulden ◽  
C. Hopkinson ◽  
M. N. Demuth

Abstract. Glacial mass balance estimated through the geodetic method requires glacial surface coordinate observations from historical and contemporary sources. Contemporary observations and historical topographic maps are typically referenced to separate horizontal and vertical datums and observed with different sampling intervals. This research demonstrates the sensitivity of glacial change detection to the datum considerations and sampling schemes through case studies of Andrei, Bridge and Peyto glaciers in Western Canada. To simulate the procedure of observing the glacial surfaces, profile lines were sampled from Digital Elevation Model (DEMs) on contour intervals for historical data and horizontal intervals for contemporary data. Profile lines from the following scenarios were compared: (1) different horizontal and vertical sampling schemes; (2) the horizontal datum was correctly reconciled but the vertical datum was not; (3) the vertical datum was correctly reconciled but the horizontal datum was not; (4) both the horizontal and vertical datums were correctly reconciled; and (5) both the horizontal and vertical datums were incorrectly reconciled. Vertical errors of up to 6.9 m, 6.0 m and 5.0 m were observed due to sampling effects and vertical errors of 22.2 m, 9.9 m and 55.0 m were observed due to datum inconsistencies on Bridge, Andrei and Peyto glacier respectively. Horizontal datum inconsistencies manifested as erratic levels of growth or downwasting along the glacial surface profile and vertical datum errors manifested as a consistent vertical offset. Datum inconsistencies were identified to contribute errors of up to 257.2 × 106 m3 (or 87%) and 54.6 × 106 m3 (or 580%) of estimated volume change below and above the equilibrium line respectively on Peyto Glacier. The results of this study provide an estimate of typical errors due to sampling constraints or datum inconsistencies as well as guidance for identifying where these error sources have contaminated mass balance results.


Author(s):  
S. D. Jawak ◽  
A. J. Luis

Available digital elevation models (DEMs) of Antarctic region generated by using radar altimetry and the Antarctic digital database (ADD) indicate elevation variations of up to hundreds of meters, which necessitates the generation of local DEM and its validation by using ground reference. An enhanced digital elevation model (eDEM) of the Schirmacher oasis region, east Antarctica, is generated synergistically by using Cartosat-1 stereo pair-derived photogrammetric DEM (CartoDEM)-based point elevation dataset and multitemporal radarsat Antarctic mapping project version 2 (RAMPv2) DEM-based point elevation dataset. In this study, we analyzed suite of interpolation techniques for constructing a DEM from RAMPv2 and CartoDEM-based point elevation datasets, in order to determine the level of confidence with which the interpolation techniques can generate a better interpolated continuous surface, and eventually improves the elevation accuracy of DEM from synergistically fused RAMPv2 and CartoDEM point elevation datasets. RAMPv2 points and CartoDEM points were used as primary data for various interpolation techniques such as ordinary kriging (OK), simple kriging (SK), universal kriging (UK), disjunctive kriging (DK) techniques, inverse distance weighted (IDW), global polynomial (GP) with power 1 and 2, local polynomial (LP) and radial basis functions (RBF). Cokriging of 2 variables with second dataset was used for ordinary cokriging (OCoK), simple cokriging (SCoK), universal cokriging (UCoK) and disjunctive cokriging (DCoK). The IDW, GP, LP, RBF, and kriging methods were applied to one variable, while Cokriging experiments were employed on two variables. The experiment of dataset and its combination produced two types of point elevation map categorized as (1) one variable (RAMPv2 Point maps and CartoDEM Point maps) and (2) two variables (RAMPv2 Point maps + CartoDEM Point maps). Interpolated surfaces were evaluated with the help of differential global positioning system (DGPS) points collected from study area during the Indian Scientific Expedition to Antarctic (ISEA). Accuracy assessment of the RAMPv2 DEM, CartoDEM, and combined eDEM (RAMPv2 + CartoDEM) by using DGPS as ground reference data shows that eDEM achieves much better accuracy (average elevation error 8.44 m) than that of existing DEM constructed by using only CartoDEM (13.57 m) or RAMPv2 (41.44 m) alone. The newly constructed eDEM achieves a vertical accuracy of about 7 times better than RAMPv2 DEM and 1.5 times better than CartoDEM. After using accurate DGPS data for accuracy assessment, the approximation to the actual surface of the eDEM extracted here is much more accurate with least mean root mean square error (RMSE) of 9.22 m than that constructed by using only CartoDEM (RMSE = 14.15 m) point elevation data and RAMPv2 (RMSE = 69.48 m) point elevation data. Our results indicate that, the overall trend of accuracy for the interpolation methods for generating continuous elevation surface from CartoDEM + RAMPv2 point elevation data, based on RMSE, is as follows: GP1 > IDW > GP2 > OK > LP2 > DK > LP1 > RBF > SK > UK. In case of cokriging interpolation methods, OCoK yields more accurate eDEM with the least RMSE of 8.16 m, which can be utilized to generate a highly accurate DEM of the research area.. Based on this work, it is inferred that GP2 and OCok interpolation methods and synergistic use of RAMPv2 and CartoDEM-based point elevation datasets lead to a highly accurate DEM of the study region. This research experiment demonstrates the stability (w.r.t multi-temporal datasets), performance (w.r.t best interpolation technique) and consistency (w.r.t all the experimented interpolation techniques) of synergistically fused eDEM. On the basis of average elevation difference and RMSE mentioned in present research, the newly constructed eDEM may serve as a benchmark for future elevation models such as from the ICESAT-II mission to spatially monitor ice sheet elevation.


2019 ◽  
Vol 8 (4) ◽  
pp. 180 ◽  
Author(s):  
Jiang Wu ◽  
Junjie Fang ◽  
Jiangbo Tian

Digital elevation model (DEM) resolution is closely related to the degree of expression of real terrain, the extraction of terrain parameters, and the uncertainty of statistical models. Therefore, based on DEMs with various resolutions, this paper explores the representation and distinguishing ability of different roughness algorithms to measure terrain parameters. Fuyang, a district of Hangzhou City with various landform types, was selected as the research area. Slope, root mean squared height, vector deviation, and two-dimensional continuous wavelet transform were selected as four typical roughness algorithms. The resolutions used were 5, 10, 25, and 50 m DEM on the scale for plains, hills, and mountainous areas. The statistical criteria of effect size and entropy were used as indicators to evaluate and analyze the different roughness algorithms. The results show that in terms of these measures: (1) The expression ability of the SLOPE and root mean squared height (RMSH) algorithms is better than that of the vector deviation method, while the two-dimensional continuous wavelet method based on frequency analysis emphasizes the terrain information within a certain range. (2) The terrain distinguishing ability of the SLOPE and RMSH is not sensitive to the changes in resolution, with the other two algorithms varying with the changes in resolution.


2020 ◽  
Vol 66 (260) ◽  
pp. 927-937
Author(s):  
Mingyang Lv ◽  
Duncan J. Quincey ◽  
Huadong Guo ◽  
Owen King ◽  
Guang Liu ◽  
...  

AbstractGlaciers in the eastern Pamir have reportedly been gaining mass during recent decades, even though glaciers in most other regions in High Mountain Asia have been in recession. Questions still remain about whether the trend is strengthening or weakening, and how far the positive balances extend into the eastern Pamir. To address these gaps, we use three different digital elevation models to reconstruct glacier surface elevation changes over two periods (2000–09 and 2000–15/16). We characterize the eastern Pamir as a zone of transition from positive to negative mass balance with the boundary lying at the northern end of Kongur Tagh, and find that glaciers situated at higher elevations are those with the most positive balances. Most (67% of 55) glaciers displayed a net mass gain since the 21st century. This led to an increasing regional geodetic glacier mass balance from −0.06 ± 0.16 m w.e. a−1 in 2000–09 to 0.06 ± 0.04 m w.e. a−1 in 2000–15/16. Surge-type glaciers, which are prevalent in the eastern Pamir, showed fluctuations in mass balance on an individual scale during and after surges, but no statistical difference compared to non-surge-type glaciers when aggregated across the region.


Sign in / Sign up

Export Citation Format

Share Document