Data Mining & Professional Cartography

Mapping our changing climate

Basics

The focus of the exercise is to further explore the importance of accurate display techniques and communication of technical and complex data to a general audience. You will create a professional map or figure comparing temperature change maps for two different periods of time relative to a baseline period.

You will download synthesized global climate data from NASA  /Goddard Institute for Space Studies’ Surface Temperature Analysis site.

The main challenge of this exercise is to wrangle the data values, symbolize them accurately, and communicate their meaning clearly.

When preparing a map to present information, the map maker must make choices about what information to include or leave off the map in order to best support the map’s purpose.

This tutorial was originally conceived by Dr. Peter Howe at Utah State University.

Data

  • Climate data downloaded from NASA Goddard Institute for Space Studies
  • Continental boundaries downloaded from ArcGIS Online

NetCDF Files

You will be working with global temperature data in this lab.  Many organizations are making global temperature, precipitation, atmospheric pressure, and wind data available for scientists and the public to access, visualize, and analyze. 

The format of the data is NetCDF, short for Network Common Data Format.  This is a type of raster file that allows you to include multiple raster layers in the same file.  It is used in the fields of meteorology and climatology to display large sets of weather and climate data at different points in time.  NetCDF files are an open source file format that have metadata embedded in the file, which makes them easy to share and use on different computers with different software packages.

Download Data

  • Go to the NASA GISS website (Google NASA GISS).  
  • You are looking for Global Maps of Surface Temperatures
    • Navigate to Datasets > Surface Temperature (the upper left portion of their home page)
    • Navigate to Global Maps (link on the right hand side)

Global Map in infared from November 2019

 

You are looking for this form

NASA GISS Surface Temperature Anyalysis version 4 form

 

Set-up the Form

You will produce at least two datasets. Each will display a calculation of global temperature change between two periods of time.

Use the same base period for each calculation. For example, you might show the change in surface temperatures from the 1900s to the the 1980s and then in a second map how global temperatures have changed between the 1900s and the 2010s.


  • Map type = Anomalies
    • This calculates the difference in mean temp between the two time periods you will input
  • Mean Period = You decide, but you must include the mean period you choose on your map
    • This filters the part of the year over which the temperatures are averaged
    • You can choose a summer month, a winter month, a calendar year or meteorological year or season.
  • Time Interval = You decide, you will set up this form twice
  • Base Period = You decide, but keep it the same between your two runs of the model.
  • Smoothing Radius = Use the default
  • Map Projection = Use the default

Press Make Map

Video demonstration and discussion of the data mining options:

The map will show temperature differences between a baseline period and a more current period of time. 

You are mapping the difference in temperature (the change) not temperature.

GISS infared of temperature for the world from 1890-1919

Model Run June 2022

The number in the upper right (0.91 C in the example) is the calculated global mean temperature (increase in this case) between the mean periods; the difference (degrees C) between the global mean of the earlier time period and the global mean of the later time period. Globally, there has been a 0.91 C increase in mean annual temperature between the average temperature for the 30 year period 1890-1919 and the 30 year period 1992-2021.

 

The model also produces a graph. You can see that the largest increases in temperature occur in the northern latitudes. At around 60 degrees latitude in the southern hemisphere, there is a region showing very little change.

Graph of the model run from January 2021. Zonal Mean (in celsius) versus Latitude

Download Data for later use

  • Left click netCDF at the bottom of the page.
  • Save the file to a deliberate location.
  • The file will have a generic name: nmaps.

I've name this file to represent current preiod 1989-1990 against baseline period 1951-1980. File folder when downloading with "amaps.nc" in file name

fas fa-thumbtack
Pro Tip
All the maps come with the same name. Name the files as you download them so you don't lose track.

Adding NetCDF files to ArcGIS Pro

NetCDF files cannot be added to ArcGIS Pro using the Add Data button like you can with other data. The NetCDF file must first be converted to raster format.

Open the tool Make NetCDF Raster Layer


 

netcdf tool

Set up the tool as seen above with your file names and run.

fas fa-sticky-note
Note
The NetCDF Raster Layers are temporary files. Later I will show you how to save them so they are permanent in your files.

Another way to add NetCDF data is to use the “Add Multidimensional Raster Layer” option in the Add Data drop down.

Here's how:

Set up the Add Multidimensional Raster Layer window like above.

This is continuous raster data. Nearest Neighbor is not an appropriate interpolation method.

Change to Bilinear.

Repeat with your second netcdf file to add both datasets to the contents pane.

 

Download Country Boundaries from ArcGIS Online

File folder with "Relevance" sort feature selected

Notice you can sort the results by relevance (mine defaults to “date modified” which I don’t find very helpful).

Change the symbology so that the polygons aren’t filled and the line is unobtrusive.

 

Acquaint yourself with these new datasets

First – look at the range of values in the table of contents for each layer. What are the min and max values? What are the units of measure? What do these values mean? You really need to know this to make effective maps.

Zoom in until you can see individual cells. Use the Explore tool to sample the cell values.

Should positive values be displayed with warm colors or cool colors?

Open the Properties for the new raster layer > Source tab

"Explore" option in top ribbon selected with "Contents" of various colors shown

  • The Spatial Reference (coordinate system) is listed as WGS 84. This is a geographic coordinate system (GCS).  
  • It is a floating point raster. This means that the cell values can contain decimals. This also means it won’t have an attribute table. (Each row is a unique value. Floating point could have an infinite number of rows... The table would be too big.)
  • The default units are degrees (think lat/lon). The cell size is 2,2. This means each cell of the raster spans 2 degrees by 2 degrees.  

Evaluate the spatial distribution of the temperature changes across the map

Why? Because that’s what GIS experts do.

What is the map telling you about change in global surface temperatures over time (how are the values spatially distributed)? (Your map will look different depending on which periods of time you are comparing, of course.)

Infared of the world.

In my map, there is an area in the north that seems to be experiencing a high increase in average annual temperatures relative to the base period (1900-1929) but the rest of the map seems to suggest cooling across most latitudes with exceptions in Brazil and near the Antarctic Peninsula.

 

Does that jive with what we know about increasing global temperatures? The draft map created on the GISS website for these time periods showed consistent warming and a global average increase of 0.9 degrees Celsius:

GISS infared of temperature for the world from 1890-1919

fas fa-thumbtack
Pro Tip
The anomaly values are displayed with a default color ramp. While the warm and cool colors make sense for values that measure the amount of warming or cooling since the base period, you cannot assume that warm colors are being properly mapped to positive ('warming') values.

 

 

What is the range in values for your anomaly layer?

My experimental run produced the following:

Temperature Anomally values. Your values will be different

The color ramp goes from cool to warm colors (blue to red), but the range of values is most likely not centered evenly around 0 degrees. (Depending on your time periods, you may not have any negative values at all.)  Additionally, the color scheme appears to be dominated by cool colors. Notice how the yellow is well above the midpoint. Much of the ‘warming’ temps are displayed with cool colors (greens and blues).

That creates a misleading visualization of the change in temperature values as we tend to associate yellows oranges and reds with warmth and greens, blues, and purples as cool colors.

 

Wrangling the Symbology 

While there is a lot of room for uniqueness and individuality in cartography, there are some visualization rules that must be followed.

With this data, those rules include:

  • There must be a neutral color that represents directionality from a meaningful midpoint. In our case “0” temperature change between the two periods of time.
  • Intuitive colors must be used to represent warming and cooling.
  • The symbology on the two maps must be united. If “dark red” equals 2 degrees C increase in temperature on one map, the same shade of red must represent 2 degrees C warming on the second map. In other words, apply the same range of values to both color schemes.

You need to present these temperature change values in an intuitive and accurate manner.

Use a diverging color scheme with warm and cool colors ensuring the “0” value (no modeled temp change between time periods) is set to the neutral color.

color hue

Be aware of the "Neutral" color. If it is a 'warm' color, then you are imparting bias by displaying areas with no change in temperature as areas that apprear to have a slight warning trend over time. Both examples above have a yellow mid-color on which the 'no change' values will pivot. That will make your map look warmer than it. 

fas fa-sticky-note
Note
If your raster does not contain negative values (cooling) then a sequential color ramp is appropriate. However, if you are displaying two maps and one does contain negative values, you can still apply one color scheme to both maps.

Matching color schemes

I'm going to use these maps to walk through the symbology work:

 

With the crazy distribution of positive and negative values, this is most easily done by altering the symbology’s min and max value to be equal opposites, to ensure 0 is in the “middle” of the value distribution.

Color hue red to blue

Open the symbology of the new random raster.

Set up the parameters:

  • Primary Symbology = Stretch
  • Color Scheme – choose a divergent color ramp, one with a neutral mid-color.
  • Inspect the values to determine if your colors need to be inverted
    • Warm colors are positive value, negative values and cool colors
  • Label boxes: Round and include units (C ).
  • Change Stretch Type to Minimum Maximum

In the Statistics tab:

  • Change the Statistics option to Custom
  • Replace the Min and Max values with the rounded ‘greatest absolute’ value.

Desired range selector

In my example, 4.49893 was the greatest absolute value. I rounded the value to slightly above in order to include that value. Then entered the corresponding opposite (-4.499) in the minimum box.

Desired range scheme advanced options

Import symbology to each raster layer

Symbology import "Import from layer file" option

Note: you might need to toggle the statistics to Custom

Evaluate the resulting symbology

  • You should be able to detect a visible trend between the two time periods given what we know about global climate warming.
  • Use the Explorer tool to sample locations to determine that positive temperature change values are drawing with warm colors and negative with cool colors.

Save Raster Symbology

These rasters are temporary layers. You should make them permanent so you can work with them later without having to recreate them from the netcdf files. At the same time, we will save the symbology work you have done.

  • Export Rasters: Right click in contents pane > Data > Export Raster.
  • Save Symbology: Right click in contents pane > Sharing > Save as layer file.

Save the symbology of both new raster layers.

Now it’s all about making a professional presentation.

Map Layout

In ArcGIS Pro you will need a Map for each Raster Layer.

  • Insert > New Map

Add one of your saved anomaly raster layer files to each map.

Rename the maps to reflect the different time periods.

How to rename a map

Insert New Layout to create the digital page the maps will be set up on.

Insert Map Frame > Choose one of your maps > Draw a space on the layout

How to insert a new map frame

How to draw a map frame. A dotted line on a white page in a rectangular shape

Repeat to draw a space for the second map

Change the display coordinate system for the map frames to get a more natural look:

  • Table of Contents > right click on the MAP title (not the map frame or raster layer) > Properties > Coordinate system tab
  • Use a Projected World coordinate system.  I fancy the Natural Earth projections.

Drawing order of a map. Different editable frames shown

Legends

You can insert legends from the Insert tab.

Draw a space for it on the page like you did with the map frames.

Legend formatting takes a lot of work in ArcGIS. I’d recommend building a legend in some user friendly graphics software.

This is one legend I built in PowerPoint:

Legend of a temperature map

Why PowerPoint? Because it would take 2 months to do this in ArcGIS. You can also use a free graphics website called “Canva.” Canva and PowerPoint are user friendly and efficient for simple graphic manipulations (resizing, rotating, adding marks and text).

Color bar

How to insert a picture on ArcGIS Pro

Here’s how you might proceed.

  1. Take a screengrab of the color bar in the contents pane of ArcPro. 
    1. Or for better results, insert a legend and expand the color bar first to conserve the quality (I’ve got a video demonstrating this) and/or export the layout to something high resolution to add to PowerPoint or Canvas.
  2. Paste it into a PowerPoint slide and enlarge it and crop it so you have one clean image of the color bar.
    1. You might have a balanced color bar (dark red to dark blue instead of light blue like above). That’s fine. You decide how you want to present your values.
  3. Add text boxes to label the bar with your minimum and maximum temperature change values and zero.
  4. Add tick marks and units of measure for the values
  5. Enlarge the final product, screen-grab it or export/save the image as a png file. In ArcGIS Pro, used the insert picture option to put the new legend on your map.

There are many ways you can approach the legend for a map like this. 

Here are some examples:

Map showing difference of temperature over 1952-2005. Created by J. Roughgarden

Map showing temperature difference. Created by T. Prichard

Fully descriptive temperature map for the globe. Created by S. Wing