Tuesday, August 25, 2015

The Evolution of Models - 12

Fig. 1 Jakobshavn Glacier, Greenland
I. Introduction

In the previous post of this series I presented some variations in software modelling which "length of doubling" produces in sea level change (SLC) projections.

"Doubling" is a description that scientists use to describe the observance of ice loss on ice sheets in the sense that if the sheet loses "x amount" of ice in one year, then a number of years later that "x amount" of ice loss per year doubles (x * 2, 2x), then that is considered a "doubling."

If it takes ten years for x to double, that is a "ten year doubling," and if it takes seven years for x to double, that is a "seven year doubling."

II. Recent Scientific Papers

One recent paper by 17 scientists, including James Hansen, indicates that with a ten year doubling, sea level rise (SLR) could be in the vicinity of ten feet by 2050 (A Paper From Hansen et al. Is Now Open For Discussion).

A more recent paper indicates that calving and melting of glaciers around the world is doubling faster than probably at any other time in recorded history:
"Observations show that glaciers around the world are in retreat and losing mass. Internationally coordinated for over a century, glacier monitoring activities provide an unprecedented dataset of glacier observations from ground, air and space. Glacier studies generally select specific parts of these datasets to obtain optimal assessments of the mass-balance data relating to the impact that glaciers exercise on global sea-level fluctuations or on regional runoff. In this study we provide an overview and analysis of the main observational datasets compiled by the World Glacier Monitoring Service (WGMS). The dataset on glacier front variations (~42 000 since 1600) delivers clear evidence that centennial glacier retreat is a global phenomenon. Intermittent readvance periods at regional and decadal scale are normally restricted to a subsample of glaciers and have not come close to achieving the maximum positions of the Little Ice Age (or Holocene). Glaciological and geodetic observations (~5200 since 1850) show that the rates of early 21st-century mass loss are without precedent on a global scale, at least for the time period observed and probably also for recorded history, as indicated also in reconstructions from written and illustrated documents. This strong imbalance implies that glaciers in many regions will very likely suffer further ice loss, even if climate remains stable."
(Historically Unprecedented Global Glacier Decline, emphasis added). Later in this post (Section V) I will try to put that into the context of "doubling."

III. You Were A Naughty Scientist

Generally, when scientists publish the type of information mentioned in Section II above, they are chastened for being "alarmist" in the unspoken Jim Crow law sort of way (Blind Willie McTell News).

The psychological denial of things civilization is out of control about, such as global warming induced climate change, has now wormed its way into the institutions of science.

Scientists are treated like teenagers with an allowance that is taken away if the scientists publish things that cannot be spoken:
"I suspect the existence of what I call the `John Mercer effect'. Mercer (1978) suggested that global warming from burning of fossil fuels could lead to disastrous disintegration of the West Antarctic ice sheet, with a sea level rise of several meters worldwide. This was during the era when global warming was beginning to get attention from the United States Department of Energy and other science agencies. I noticed that scientists who disputed Mercer, suggesting that his paper was alarmist, were treated as being more authoritative.

It was not obvious who was right on the science, but it seemed to me, and I believe to most scientists, that the scientists preaching caution and downplaying the dangers of climate change fared better in receipt of research funding. Drawing attention to the dangers of global warming may or may not have helped increase funding for relevant scientific areas, but it surely did not help individuals like Mercer who stuck their heads out. I could vouch for that from my own experience. After I published a paper (Hansen et al 1981) that described likely climate effects of fossil fuel use, the Department of Energy reversed a decision to fund our research, specifically highlighting and criticizing aspects of that paper at a workshop in Coolfont, West Virginia and in publication (MacCracken 1983).

I believe there is a pressure on scientists to be conservative. Papers are accepted for publication more readily if they do not push too far and are larded with caveats. Caveats are essential to science, being born in skepticism, which is essential to the process of investigation and verification. But there is a question of degree. A tendency for `gradualism' as new evidence comes to light may be ill-suited for communication, when an issue with a short time fuse is concerned."
(Scientific Reticence & SLR, emphasis added). This is how civilizations end up damaging themselves out of existence, or as Toynbee described it "committing suicide" (Civilization Is Now On Suicide Watch).

IV. Believing Your Eyes

There is no doubt that the ice sheets of Greenland (Fig. 1) and Antarctica are, like glaciers elsewhere, "in a hurry" to melt and calve:
"Satellite images show that the world's fastest moving glacier lost a piece of ice measuring nearly 5 square miles over two days. Scientists say it is one of the most significant calving events on record.

Radar images of the Jakobshavn glacier in Greenland, taken from the European Space Agency's (ESA) Sentinel-1A and Sentinel-2A weather satellites between July 27 and Aug. 19, show the massive glacier advancing westward before its front rapidly retreated as a huge chunk fell off it front, sometime between Aug. 14 and Aug. 16.

The ESA estimates that the ice is about 1,400 meters (4,600 feet) thick, meaning that the piece lost has a volume of about 17.5 cubic km, which is enough to cover the island of Manhattan in a layer of ice almost 1,000 feet deep."
(CBS News, emphasis added). So, what does this have to do with SLC modelling software that uses a doubling scheme?

V. Believing SLC Software Models

A. The Code

I wrote a simple SLR software model that does an "oscillating doubling" of 10,7,5,3,2 year sequences, as follows.

The code below is in bold, comments are in regular text between "/** */" C++ comment markers.

The SLR program is listed between the following horizontal lines:

#include <iostream>
#include <iomanip>
#include <fstream>

using namespace std;

/** URL for "80.32" is Table 1 @: http://pubs.usgs.gov/fs/fs2-00/ */
const double maxMeanSLR = 80.32;

/** csv comma delimited header (percent is % of maxMeanSLR) */
const char *pColumnNames = "Year,SLR,Rate,Percent";

/** reoeating doubling cycle */
const int doublingYrs[5] = {10,7,5,3,2};

/** current max est. SLR per year @ East Coast (0.866141732 in.) */
const double incPerYear = 0.022;

/** number of years to project (36 for yr. 2050; 86 for yr 2100) */
const int yearsToProject = 86;

/*************************************
This program uses a 10,7,5,3,2 year
repeating "doubling sequence."

This is designed to vary the ice
sheet melt/calving as different
zones are traversed and cycles of
global warming events wax and wane.

The rate of melt changes each
10 then 7 then 5 then 3 and then
2 years, and then repeats itself.
------------------------------------
all sea level values are in meters
*************************************/
int main()
{
    /** current sea level rise since 1750 @East Coast of U.S. */
    double meanSLR = 0.4572;

    /** set initial SLR per year @East Coast (0.022 m, 0.0722 ft, 0.866 in.) */
    double incThisYear = incPerYear;

    int doublingElement = 0, doublingCount = 1;

    /** print results to a .csv text file */
    freopen("slr-lite.csv", "w", stdout);

    /** csv header: write column names and commas at top of page */
    cout << pColumnNames << endl;

    for (int year = 0; year < yearsToProject; year++)
    {
        if (meanSLR < maxMeanSLR)
        {
            /** simulate ice sheet melt rate of increase */
            meanSLR += incThisYear;

            /** store year to year values in csv format */
            cout << year + 2015 << ","
                 << setprecision(5)
                 << meanSLR << ","
                 << incThisYear << ","
                 << (meanSLR / maxMeanSLR) * 100
                 << endl;

            /** check doubling settings */
            if (doublingCount == doublingYrs[doublingElement])
            {
                /** variable "incThisYear" records values in a 10,7,5,3,2 year cycle */
                incThisYear += incPerYear;

                /** test for end of 10,7,5,3,2 cycle */
                if (doublingYrs[doublingElement] == 2)
                {
                    /** reset to beginning, repeat cycle */
                    doublingElement = 0;

                    /** reset to original doubling amount */
                    incThisYear = incPerYear;
                }
                else /** go to next doubling value */

                {
                    doublingElement++;
                }
                /** reset the cycle counter */
                doublingCount = 1;
            }
            else
            {
                doublingCount++;
            }
        }
        else /** all ice has been depleted, so don't record data */
        {   
            cout << year + 2015
            << " - maximum SLR reached ("
            << maxMeanSLR
            << ")"
            << endl;
        }
    }

    return 0;
}


B. The Resulting Data and Graphs

When set to an ending in the year 2100, it prints the following csv ("comma separated values") file, which I used to do two graphs.

The data are on the left side, the two graphs (Fig. 2, Fig. 3) are on the right side.

Years = 2015-2100.
SLR = sea level rise in meters.
Rate is the value associated with the 10,7,5,3,2 doubling.
Percent is the percent of the total ice on Earth (equating to the 80.32 meters of potential global mean SLR) that has melted or calved by that year.

Year,SLR,Rate,Percent
2015,0.4792,0.022,0.59661
2016,0.5012,0.022,0.624
2017,0.5232,0.022,0.65139
2018,0.5452,0.022,0.67878
2019,0.5672,0.022,0.70618
2020,0.5892,0.022,0.73357
Fig. 2 2015-2050
2021,0.6112,0.022,0.76096
2022,0.6332,0.022,0.78835
2023,0.6552,0.022,0.81574
2024,0.6772,0.022,0.84313
2025,0.7212,0.044,0.89791
2026,0.7652,0.044,0.95269
2027,0.8092,0.044,1.0075
2028,0.8532,0.044,1.0623
2029,0.8972,0.044,1.117
2030,0.9412,0.044,1.1718
2031,0.9852,0.044,1.2266
2032,1.0512,0.066,1.3088
2033,1.1172,0.066,1.3909
2034,1.1832,0.066,1.4731
2035,1.2492,0.066,1.5553
2036,1.3152,0.066,1.6375
2037,1.4032,0.088,1.747
2038,1.4912,0.088,1.8566
2039,1.5792,0.088,1.9661
2040,1.6892,0.11,2.1031
2041,1.7992,0.11,2.24
2042,1.8212,0.022,2.2674
2043,1.8432,0.022,2.2948
2044,1.8652,0.022,2.3222
2045,1.8872,0.022,2.3496
2046,1.9092,0.022,2.377
2047,1.9312,0.022,2.4044
2048,1.9532,0.022,2.4318
2049,1.9752,0.022,2.4592
2050,1.9972,0.022,2.4866
2051,2.0192,0.022,2.5139
2052,2.0632,0.044,2.5687
2053,2.1072,0.044,2.6235
2054,2.1512,0.044,2.6783
2055,2.1952,0.044,2.7331
2056,2.2392,0.044,2.7878
2057,2.2832,0.044,2.8426
2058,2.3272,0.044,2.8974
2059,2.3932,0.066,2.9796
2060,2.4592,0.066,3.0618
2061,2.5252,0.066,3.1439
2062,2.5912,0.066,3.2261
2063,2.6572,0.066,3.3083
2064,2.7452,0.088,3.4178
2065,2.8332,0.088,3.5274
2066,2.9212,0.088,3.637
2067,3.0312,0.11,3.7739
2068,3.1412,0.11,3.9109
2069,3.1632,0.022,3.9382
2070,3.1852,0.022,3.9656
2071,3.2072,0.022,3.993
2072,3.2292,0.022,4.0204
2073,3.2512,0.022,4.0478
2074,3.2732,0.022,4.0752
2075,3.2952,0.022,4.1026
Fig. 3 2015-2100
2076,3.3172,0.022,4.13
2077,3.3392,0.022,4.1574
2078,3.3612,0.022,4.1848
2079,3.4052,0.044,4.2395
2080,3.4492,0.044,4.2943
2081,3.4932,0.044,4.3491
2082,3.5372,0.044,4.4039
2083,3.5812,0.044,4.4587
2084,3.6252,0.044,4.5134
2085,3.6692,0.044,4.5682
2086,3.7352,0.066,4.6504
2087,3.8012,0.066,4.7326
2088,3.8672,0.066,4.8147
2089,3.9332,0.066,4.8969
2090,3.9992,0.066,4.9791
2091,4.0872,0.088,5.0886
2092,4.1752,0.088,5.1982
2093,4.2632,0.088,5.3078
2094,4.3732,0.11,5.4447
2095,4.4832,0.11,5.5817
2096,4.5052,0.022,5.6091
2097,4.5272,0.022,5.6365
2098,4.5492,0.022,5.6638
2099,4.5712,0.022,5.6912
2100,4.5932,0.022,5.7186

C. The Analysis

The csv file printout shows the "0.22" Rate repeating for 10 yrs., then the "0.44" rate following for 7 yrs., the "0.66" rate following that for 5 yrs., and so on.

Then, when the first sequence (10,7,5,3,2) completes in the yr. 2041, the sequence begins anew to repeat itself.

One effect this technique has is that in 2050 the SLR is 1.9972 meters, or 6.6 ft. instead of the 10 feet using only the 10 yr. doubling (as shown here).

The point is that projection software is probably better if it has oscillating variation, because it is not realistic to expect the same rate for each year as the old style models tend to do ("linear" projections).

The graphs reflect the sequence end, and a new beginning, in a sawtooth type of look, reflecting what the zone effect can and should look like (e.g. The Evolution of Models - 5).

That is, the graphs reflect a sawtooth look which is reflective of uneven melt rates from year to year, with an upward trend nevertheless.

VI. Conclusion

The ice sheets continue to melt and calve, with an acceleration built into that melting and calving (because of the increase in fossil fuel use and green house gas emissions).

Software models and scientific papers need to take that into consideration.

Public officials also need to take that into consideration (You Are Here - 5).

In fact, everyone else needs to take that reality into consideration.

The next post in this series is here, the previous post in this series is here.

01:30 Warming is concentrated in the Arctic ... this pattern has been evident for decades
02:45 These same models that began development in the late 60's always had an Arctic amplified warming because the physics of that reflectivity is simple.
02:55 Removal of reflective cover absorb more sunlight. The primitive models got that and its of course validated by observation.
04:32 Within our lifetimes we will bear witness to this story unfolding rather rapidly.
04:40 The red line shows observations of sea ice area. It depicts how sea ice is retreating at twice the rate [doubling] that our best model projections show.
05:05 The models still don't get the true fidelity of climate, the rapidity of sea ice loss ...
05:30 ... retreating at a rate four times [doubling * 2] what model projections currently get ... the models are if anything underestimating ...
06:28 Greenland is absorbing an additional amount of solar energy that equals twice the energy use of the U.S. each year
07:53 Just in the past ten years science reveals more sensitive response than was previously encoded in models.


Researcher sweating it out at an Economic Summit:



3 comments:

  1. It ain't like makin' model air-planes.

    ReplyDelete
  2. Honing in on reality - that it's happening at all is the scare of course (because we're powerless to stop it) - but this gives us a good idea of where things are going (you could run the models with various doubling values at different times to get a range that the measured reality of sea level at any time might occur between). Nice work Dredd.

    Tom

    ReplyDelete
    Replies
    1. Yep.

      Copy the code, put it into your compiler, and generate many hypothetical scenarios.

      Delete