Lab 5: Analyzing Glacier Change with Satellite Imagery

Crystal Huscroft

In this lab you will be using your understanding of glacial processes to make observations and measurements documenting and predicting the consequences of climate change for a Canadian alpine glacier. You will learn how to analyze the characteristics of glaciers and glacial landforms from a variety of types and sources of satellite imagery.

Learning Objectives

After completion of this lab, you will be able to

  • Make observations, identify, and appreciate the consequences of global warming for alpine glaciers.
  • Analyze the characteristics of alpine glaciers in order to make predictions regarding the future survival of a glacier.
  • Summarize and collate Geographic Information System (GIS) point, line, and area information from different platforms into a single GIS project.
  • Acquire and interpret remotely sensed imagery.
  • Communicate effectively regarding the impacts of global warming for British Columbian glaciers.

Pre-readings

Introduction to Glacial Processes and Landforms

Understanding of the following key terms is required for this laboratory activity:

  • Ablation zone
  • Accumulation zone
  • Annual ice horizon
  • Cirque glacier
  • Crevasses
  • Equilibrium line
  • Firn
  • Glacial ice
  • Glacier toe
  • Moraine
  • Snowline
  • Trimline
  • Valley glacier

If you are unfamiliar with any of these terms, look them up before continuing.

How to Determine Whether a Glacier Will Survive Current Climate

You will be asked to predict whether or not the glacier that you choose to study will survive within the current climate.  Glaciers form and persist in areas where snow accumulation is consistently equal to or greater than snow and ice melt or sublimation (ablation).  Glaciers that are in danger of disappearing will show evidence of insufficient snow accumulation to persist through the melt season.   As described by Dr. Mauri Pelto on Glacier Survival website (please read) and Forecasting temperate alpine glacier survival from accumulation zone observations [PDF] (Pelto, 2010), the likelihood of a glacier’s survival can be assessed using the following visual clues in satellite imagery:

  1. Newly exposed rock outcrops in the accumulation zone. These rock outcrops may appear lighter than surrounding rock due to a lack of lichen growth on their surface.
  2. Lowering of the surface and a reduction in size of the upper margins of the glacier. This might also be indicated by light coloured rock margins around the edges of the glacier due to exposure of unvegetated (mostly lichen) rock and debris.
  3. Discontinuous snow cover in the accumulation zone. This may appear as patches of bare ice surrounded by snow.
  4. Less than a third snow cover over the glacier late in the melt season.

Press here to connect directly to a a video that captures a glaciologist, Dr. Mauri Pelto, applying the concepts above and performing the same type of analysis that you are asked to perform in this lab.

Video 10.1

Video 10.1. Analyzing Sacagawea Glacier, Wind River Range Wyoming Survival. 

Introduction to Sentinel 2 Imagery

In this lab, you will be analyzing false colour images captured by the Sentinel 2 satellite program. The Sentinel 2 program consists of twin satellites launched by the European Space Agency that are able to capture imagery of Earth at a resolution of 10-60 m. It consistently captures images between 56° S to 84° N every 5 days. Unlike our eyes, the satellite system can detect energy in the visible, near-infrared, and shortwave infrared portions of the electromagnetic spectrum.

Although the satellites measure 13 spectral bands of reflected and emitted energy, false colour images combining three specific bands of energy are excellent for studying glaciers because they are able to detect and differentiate ice and snow (Bands 12, 8A, 4). In these false colour images, snow appears light teal blue, ice appears dark teal blue, and water appears a very dark shade of blue.

Calculating Rates of Glacier Advance or Retreat

A rate is the amount of change in some quantity over a given period of time. In this lab we are interested in rates of advance or retreat of glaciers. This rate of movement is the rate of change (Δ) in distance per change in time. Glaciers move slowly, and so in this context the rate is commonly expressed in metres per year (m/y). The rate may be expressed as (Equation 10.1):

Equation 10.1

    \[ Rate = \frac{\Delta Distance}{\Delta Time} \]

Glaciers generally advance each winter and retreat each summer. However, these annual advances and retreats may not be equal and so, over time, the glacier toe may, on average, advance or retreat. The average rate of advance or retreat is calculated by taking the total distance that the glacier toe has moved from the defined starting point/time and dividing it by the number of years that have elapsed.

For example, if the glacier has retreated 70 m in 35 years, using Equation 10.1,

    \[ Rate = \frac{\Delta Distance}{\Delta Time} \]

    \[ Rate = \frac{70m}{35y} = 2m/y. \]

For a more detailed explanation of calculating rates, please refer to SERC Carleton’s How do I calculate rates? Calculating changes through time in the geosciences.

Lab Exercises

In the exercises in this lab, you will choose a Canadian glacier and create a report consisting of satellite imagery and commentary regarding the likelihood of your glacier’s survival in the current climate. You will be accessing topographic and satellite information regarding your glacier from the following platforms so that you can discuss factors that have affected your glacier’s behaviour in terms of retreat and the likelihood that the glacier will survive current and future warming:

  1. Google Earth (Web)
  2. Google Earth Timelapse
  3. Sentinel Playground

You will be guided to complete a lab report, which you will save as a PDF.

EX2: Analyzing Glacier Change and Processes With Satellite Imagery

Step 1: Capture Google Earth Engine Timelapse imagery (Figures EX2.1 & EX2.2).

  1. Visit Google Earth Engine Timelapse and navigate to a large scale view of your glacier. It may load slowly, so give it a few moments. Zoom in as far as possible with your glacier in the centre of the view. Press the pause button on the lower left of the screen, then choose the earliest view (e.g., 1984) where the position of the toe of your glacier is visible. The earliest view should be before 1990.
  2. Create a screen capture of this view and save it with the file name <5_Lastname YYYY imagery> With YYYY indicating the year. Add the image to your document as Figure EX2.1.
  3. Write a caption for Figure EX2.1 indicating what the image is of, the date of the imagery and include a link to this view by using the Share Current View button and pasting it as the image source.
  4. Keeping the exact same view, repeat the previous three steps to produce Figure EX2.2 using the most recent imagery and save the screen capture with the file name <6_year_Lastname_Glacier_extent>. Figure 10.5 provides an example. Write a caption for Figure EX2.2 describing any changes you observe.
Figure 10.5. Sample Figure EX2.2 in a student’s assignment. This image depicts the glacier in 2018. You will want to choose the most recent image of your glacier.  As compared to the 1984 image (Figure EX2.1), we can see that the glacier toe has receded, the toe narrowed, a rock island has appeared in the upper regions of the glacier. Also, the upper area has started to divide more distinctly into three smaller cirques. This could suggest the glacier is thinning in the accumulation zone. Source: Google Earth Engine Timelapse – glacier Sampled from Google Earth by C. Huscroft. Used in accordance with Google Earth terms and conditions.

Step 2: Estimate the approximate location of the glacier toe in the earliest and most recent years (creating a point feature and Figure EX2.3).

  1. (Note this step may not work on some slower computers.  Please see the next step if not successful).  Open Google Earth (Web) and ensure that the settings are set to display locations in decimal degrees.  Then return to Google Earth Engine Timelapse and click the Timelapse In Earth menu at the top of the page.  If your computer processor is fast and you have enable GL acceleration of graphics on your browser, you will be able to hover your cursor over the toe of the glacier and record the location of the toe in the earliest and most recent years.
  2. If the previous step did not work, using landmarks visible in the earliest and most recent good imagery from Google Earth Engine Timelapse (this will most likely  be 1984), add a point feature to your Google Earth project and move it to your best estimate of the location of the toe of the glacier in the earliest and most recent years. Use as many landmarks as possible to situate your point. Make sure your view is in 2D and it helps if you view both images at the same scale.
  3. For comparison, measure the length of the glacier you estimate in 1984 and create a screen capture of your measurement while the value of the measurement is still visible. Save this screen capture of Google Earth (Web) as <7_Lastname_1984_approximate_glacier_toe> (example provided as Figure 10.6).
Figure 10.6. Sample Figure EX2.3 in a student’s assignment. By locating the approximate position of the toe in 1984 Google Earth Engine Timelapse imagery, we can see that the glacier was 2.8 km long in 1984. Source: Sampled from Google Earth by C. Huscroft. Used in accordance with Google Earth terms and conditions.

Step 3: Analyze the accumulation zone in the most recent satellite imagery (Figure EX2.4).

  1. Go to the Sentinel Playground for viewing recent satellite imagery and navigate to a large scale view of your glacier.
  2. View your chosen glacier by typing in the coordinates of your glacier or name of your glacier in the search pane (magnifying glass icon) and zoom in until the glacier fills a lot of your map area.
  3. In order to find the date and associated satellite displaying the best Sentinel 2  image capturing the most recent size and pattern of the true accumulation zone, use the calendar at the top of the page.  Inspect the summer and early fall images in order to find the year’s true accumulation zone.   This date coincides with the time of year when the size of the preserved snowpack from the previous winter will be at a minimum. Therefore, set the date range for your search to mid or late summer. For glaciers in low latitudes (southern and central BC or Alberta), inspect from early August to late September. For glaciers in higher latitudes (Yukon, Nunavut and northern BC or Alberta) inspect from early July to Late August.  A link to a video on moodle provides instructions on how to do this.  In order to refresh the displayed image, switch the satellite between the Sentinel-2 L1C and Sentinel 2 L2A units.
  4. Once you have found the best image, use the menu on the left, and set the menu to Custom and drag the following Band onto the output fields R:12  G:8a  B:4 .  This band combination accentuates snow surfaces from ice as well as water.
  5. Press the Generate button.  Copy the URL of this image (button on bottom) and download the image.
  6. Create a screen capture of the glacier and save it in a file named <8_Lastname_daymonthyear_Sentinel_2_image>.
  7. Add this file to your document as Figure EX2.4 (refer to Figure 10.7 for an example) and indicate the image source by adding the URL of your current view.
  8. Write a caption in which you describe and discuss the following questions:

    a. What is the likelihood of the glacier disappearing in the current climate regime? Please refer to the visual evidence described in the lab introduction (see top of page)  that supports your prediction.  You may also want to read and Forecasting temperate alpine glacier survival from accumulation zone observations [PDF] (Pelto, 2010),

    b. What observations are you basing your hypothesis on? This involves describing the visual evidence you used.

False colour image of unnamed glacier.

Figure 10.7. Sample Figure EX2.4 in a student’s assignment. Image source: https://www.sentinel-hub.com/explore/sentinelplayground/. Source: Source: Sampled from Sentinel Hub EO Browser by C. Huscroft. Materials are shared under CC BY 4.0 license

Step 4: Calculate the rate of change (Figure EX2.5).

  1. Using the coordinates recorded using GoogleEarth Timelapse, measure the between your earliest and most recent toe positions and calculate the average rate of retreat or advance of the glacier. Take a screen capture of the distance measurement and save the screenshot as <9_Lastname_glacier_position_measurement>.  Add this image to your document as Figure EX2.5. An example is provided as Figure 10.8.
  2. Write a caption that includes the total measured retreat (or advance) distance over the period from your earliest photo until the most recent satellite imagery, and the average rate of retreat.
Figure 10.8. Sample Figure EX2.5 in a student’s assignment. Source: Sampled from Google Earth by C. Huscroft. Used in accordance with Google Earth terms and conditions.

Step 5: Get curious and do a small scientific analysis of your own choosing.

  1. By exploring the tools and images available to you on the websites (Google Timelapse or Sentinel Playground)  you have used so far in this lab, do an extra analysis of your own creation and add extra screen capture(s) as figure(s) that support your analysis as Figure EX2.6 (and additional if needed). Try to keep it as simple as possible. Be sure to state the question you are trying to answer in your analysis as well as your method. Possible analyses could involve:
    • Doing an extra measurement or observation on your glacier.
    • Accessing imagery from a different time of your glacier or an entirely different subject (flood? forest fire, increase in clear-cutting) and doing comparisons.
  • Be sure to include coordinates and links and always credit the image sources.
  1. Review your text for grammatical and spelling mistakes.
  2. Export your document to a PDF. Your PDF file needs to be a reasonable size (generally less than 5MB, but ask your lab instructor for details). In order to do this, you will need to reduce (compress) the file size of all of your images to between 150-100 ppi. Be sure to check that you have not compressed the images so much that you are not able to see them clearly on your monitor. For PC users, the compression can be done with MSWord by clicking on the image, choosing Picture Tools – Compress.

Submit as directed by your instructor.

Worksheets

Return to Step 3, EX1

Sample Student Assignments Remember that Lab 10 will be just images 5 to 9.

References

GISGeography. (Last updated 2021, February 27). 5 Free Historical Imagery Viewers to Leap Back in the Past. https://gisgeography.com/free-historical-imagery-viewers/

Pelto, M. S. (2010). Forecasting temperate alpine glacier survival from accumulation zone observations. The Cryosphere 4, 67–75. https://doi.org/10.5194/tc-4-67-2010

Media Attributions

  • Sentinel-2 L2A image on 2022-10-17