Open Mind

Red River

November 1, 2007 · 23 Comments

I recently mentioned that the timing of snowmelt runoff affects drought conditions in the American southwest. This harkens back to early this year, when I had communications with Pat Neuman, a retired hydrologist working in the midwest. He communicated to me that his studies of snowmelt runoff timing in the Red River (and other midwestern) basin indicated that it was coming earlier and earlier in the year. I even posted on the topic.

I thought I’d show an example, using streamflow for the Red River as measured at Fargo, ND. This is not the Red River of movie fame, it’s more properly called the “Red River of the North,” and forms the boundary between North Dakota and Minnesota. For each year, the data are the day of the year for the beginning of snowmelt runoff measured at Fargo, ND.


fargo1.jpg

It looks like over the last few decades, snowmelt runoff has been beginning earlier. If we perform a linear regression, we can indeed confirm that the day of beginning snowmelt has been getting earlier, and the result is statistically significant.

The changes are a lot easier to see if we plot moving averages, showing the average for each ten-year interval spaced one year apart:

fargo2.jpg

From this it looks like the changes in snowmelt runoff timing have themselves changed. In fact it rather looks like the average timing was reasonably steady until the mid-1960s, but has been getting rapidly earlier since then. A good way to characterize this is to find the moment in time at which fitting two separate straight lines (connected at their endpoints) gives the overall best fit:

fargo3.jpg

For the first time interval (1910 to 1964), the trend rate is not statistically significant, so it gives evidence of no change. For the second time interval (1965 to 2003), the trend rate is strongly significant. It appears that up to 1964, snowmelt runoff timing for the Red River at Fargo, ND, had its ups and downs but remained constant over the long haul. Since 1965, snowmelt runoff into the Red River has been happening sooner. In fact, the slope of the line for the recent data indicates that snowmelt is coming earlier by about 0.62 days per year. From 1965 to 2003, its timing is about three and a half weeks earlier.

That’s not a trivial difference, it’s quite significant. It affects the supply of water to human society, and to native living things. In southern California, it affects both the severity and the length of the summer drought season. And it’s not just affecting the Red River of the North; such changes have been documented in New England, the northern midwest, and extensively in the American West (see here for some references).

One of the most worrisome aspects of global warming is that it will bring undesirable changes in the hydrological cycle. This doesn’t just mean changes in rainfall amounts and distribution; it also depends on how the timing of the seasons is itself being affected. Fortunately, some of these changes have been studied and can be partly planned for. Let’s hope we can limit the amount of global warming the planet experiences, so the changes don’t overwhelm us.

Categories: Global Warming · climate change

23 responses so far ↓

  • JesusChristHimself // November 1, 2007 at 7:00 pm | Reply

    It dumps into some lake in Canada.

    I grew up in northeast South Dakota. The winters there now little resemble the winters described in the Laura Ingalls Wilder books. The harsh winters I experienced as a boy in the late 50s and early 60s were a blast. We had fun.

    When I go home to see my parents now the winters are decidedly more mild.

  • Eli Rabett // November 1, 2007 at 9:20 pm | Reply

    That is essentially the hinged spline that GISS uses to correct for UHI

    [Response: Pretty much, except they find the one that gives the best match to the trends in nearby rural stations, while I'm just looking for the one that gives a best fit to a "hinged spline."]

  • Alexander Ac // November 1, 2007 at 9:46 pm | Reply

    Tamino,

    earlier snowmelt has definitely significant impact on water availability and length of wildfire season in (Western) US. Regarding your recent post on Calif. wildfires, there is one interesting paper on *history* of wildfires in Calif. (NewScientist):

    http://environment.newscientist.com/article/mg19626283.900;jsessionid=LHADOECDMNNN

    in that paper they claim, that just before the European came to America, area of wildfires was much larger, than could be observed for recent wildfires in Calif. or even US. I am sure, if any climatesceptic would find it, he would pick up on it. :-)

    Best,

  • Ohg Rea Tone // November 2, 2007 at 4:05 pm | Reply

    Hey Tamino,

    Just checking in and wandering around your blog. This is very well done and appreciated by many of us. Thanks,
    Ohg Rea Tone

  • tristram shandy // November 2, 2007 at 5:24 pm | Reply

    Tamino.

    I like your approach ( and thanks for correcting Rabbett. Hansen used to hinge in 1950 and has now moved to something more flexible) Would you endorse this approach with other time series? Also, when you say “best match” to the trends what error did you minimize? I’m curious because I have a similiar problem.

    [Response: For the start-of-snowmelt data, I identify the hinge point (t,x)=(time, temperature) so that best-fit straight lines intersecting the hinge point, one to the left the other to the right, minimizes the sum of the squares of the residuals for *both* segments, with residual defined as (data value) - (straight-line value).

    If the hinge time is given, it's a straightforward (linear) problem to find the hinge temperature value, as well as the slopes of the lines to the left and the right. Finding the best hinge time is a nonlinear problem, so I just test all values on a suitably fine grid and choose the best one.

    Another method I've used often is this: make an initial guess for the hinge time. Fit straight lines separately to the data to the left, and to the right, of that time. Then determine the intersection of those two straight lines; that defines a new estimate of the hinge time. Iterate. If the procedure converges (which it generally will do, or fail to do, rather quickly) you've got it. Convergence is perfect in the sense that an iteration gives exactly the same intersection as the previous iteration.]

  • JesusChristHimself // November 2, 2007 at 7:36 pm | Reply

    When I was kid the farmers of South Dakota banded together and defeated a huge Missouri River irrigation project. They are dedicated dryland farmers. It matters a great deal when the moisture occurs. It doesn’t take much of a drought to kill off most/all of the crop. Too much moisture can also present major problems.

    Note that North Dakota is a major wheat producer.

    “Snow is an important source of moisture, and a slow snowmelt is a very beneficial event. Winds that blow snow off the fields or into drifts and low areas are not desired. This is true for adjoining grazing lands too. …”

    http://geography.about.com/library/misc/ucwheat.htm

  • Hank Roberts // November 2, 2007 at 9:04 pm | Reply

    Alexander, the point of the New Scientist article is that fire suppression recently in California is what causes the extremely hot and damaging fires to occur. They compare Baja California where fires still are tolerated — and so happen more often over larger areas but at very low intensity because there’s little available fuel.

    That used to be the norm further north; I recall a 10-year average between natural fires in N. California coast range was estimated.

    I’ve worked on fire restoration and had one site burn naturally about 15 years after a devastating fire, and it did quite well, the fire went between the trees, burned the brush that had already been piled after clearing, didn’t get into the treetops. With luck it’ll happen again in another decade. IF that goes on for a century, the Ponderosa pine and oak canopy will be back shading out the brush, and after that it can burn on its own repeatedly. But it’ll take another century of hand-clearing and removal of fuels to let the trees get as big as they used to be before the area was first logged and that shady canopy lost.

    The trick is to keep the fast-buck loggers out, they want the biggest trees cut under the excuse of improving the site.

  • George // November 2, 2007 at 11:55 pm | Reply

    It’s interesting to compare the snowmelt time graph with theannual mean temperature data from Fargo.

  • EliRabett // November 3, 2007 at 2:30 am | Reply

    TS, Tamino agreed that the method was pretty much the same, but the criteria for picking the hinge point was different.

    “The urban adjustment in the current GISS analysis is a similar two-legged adjustment, but the date of the hinge point is no longer fixed at 1950, the maximum distance used for rural neighbors is 500 km provided that sufficient stations are available, and “small-town” (population 10,000 to 50,000) stations are also adjusted. The hinge date is now also chosen to minimize the difference between the adjusted urban record and the mean of its neighbors. In the United States (and nearby Canada and Mexico regions) the rural stations are now those that are “unlit” in satellite data, but in the rest of the world, rural stations are still defined to be places with a population less than 10,000. The added flexibility in the hinge point allows more realistic local adjustments, as the initiation of significant urban growth occurred at different times in different parts of the world.”

  • tristram shandy // November 3, 2007 at 2:38 pm | Reply

    Eli,

    Thanks, I have the GISSTemp source. I think
    it’s a 1000km radius for the UHI adjustment, unless it’s ben changed recently

  • Kelly O'Day // November 3, 2007 at 3:59 pm | Reply

    Tamino:

    Great post. I looked at P Neuman’s original work and noticed that his moving average trend chart showed that the Red River trend was quite different than the St Louis or St Croix rivers.

    This prompted me to put together an Excel workbook of Neuman’s beginning snowmelt days for the 3 basins and compare them with a CuSum/ change point analysis.

    Here’s a link to my post.

    link

    Any thougths on why St Louis and St Croix river beginning day of snowmelt trends are so different than the Red River basin trend?

    Is it possible for Patrick Neuman to update his snowmelt data with 2004 – 2007 results to see if the trends are still headed in the same directions?

    PS – I particularly like your posts that let me get my hands on the actual data and see if I can reproduce your results. Any of your readers who want to play with the snowmelt data are welcome to download my workbook at the link above.

  • NeuvoLiberal // November 4, 2007 at 2:58 pm | Reply

    Hi Tamino,

    just came across something rather startling about October temperature means for the Arctic region. This brief post claims that the temperatures this October have been up to 15°C warmer than the averages for 1971-2000 for the region. Is this verifiably true? I am looking around for other news and data sources for this, but wanted to pass this along to you and the readers here. Any pointers would be much appreciated. Thank you. -NL

  • Hank Roberts // November 4, 2007 at 4:39 pm | Reply

    Neuvo, that’s a picture, can’t tell much without a cite, can you read the text to give a link to its source?

    This may help:
    http://nsidc.org/arcticmet/patterns/climate_change.html

  • NeuvoLiberal // November 4, 2007 at 6:28 pm | Reply

    Hi Hank Roberts: I don’t know Dutch either and so I used babelfish to get the following. The text doesn’t give clues as to the source:

    Quote:
    That is almost warmer then here!!

    The largest deviation was true it average over October -2° was above the marine water instead of normal -19 à -20° above a ijslaag (with sneeuwdekje).

    Greet,
    Alwin

    I happened upon this from blogger blog post. I’ll try to post a comment and/or PM and see if I can track down the original source. Thanks.

  • MrPete // November 4, 2007 at 7:21 pm | Reply

    I’m curious how they determine the date of the beginning of snowmelt runoff? Doesn’t seem like a simple thing to measure.

  • MrPete // November 4, 2007 at 7:30 pm | Reply

    Very interesting stuff. From here, they calculate SMR (SnowMelt Runoff) date as the date when 50% of volume for April-July (AMJJ) has taken place.
    The cited paper concludes:
    “The shift in SMR is related to a regional
    step increase in AMJJ temperatures during the mid-1980s. A small portion of the variability in the timing of SMR and AMJJ temperatures is attributable to PDO and ENSO; however, these climate indices are not sufficient to explain the changes in SMR timing. In addition, it appears that the observed change in the timing of SMR in the western United States is a regional response to natural climate variability and may not be related to global trends in temperature. It is unknown whether the shift to earlier SMR timing will continue or if there will be a change to later SMR timing.”

  • Gareth // November 4, 2007 at 9:03 pm | Reply

    Hank & NL: That anomaly looks to be over the ice-free part of the Arctic Ocean – N of the Bering Strait. Given that there wass no ice there during October, it’s hardly surprising that the difference with respect to climatology is that large – the baseline period would have had ice all over that region. It won’t get back to “normal” until there’s complete ice cover – end of November?

  • george // November 5, 2007 at 12:13 am | Reply

    The trends for the snowmelt date for the Red River do seem to be in keeping with the temperature trends for Fargo over the same time period (and the hinge point seems to be right around the same year)

    But spring “runoff” for some rivers is sometimes controlled by dams, which might explain differences in the trends for different rivers in the same region of the country.

    I don’t know whether that is the source of the discrepancy referred to above: “the Red River trend was quite different than the St Louis or St Croix rivers”.

    But it is at least possible that the snowmelt runoff date in question could also depend on man rather than nature alone (ie, temperature). One also must be careful to consider tributaries because the dam(s) could actually be a on a tributary and not the river in question.

  • NeuvoLiberal // November 5, 2007 at 1:19 pm | Reply

    Tamino and folks: I was able to get in touch with the person that posted those October temperature contour maps for the Arctic, and he was kind enough to send me the link where one can produce the data and plots. The datasets are from NCEP/NCAR. Here is the link (tiny URL: http://tinyurl.com/36l9dl ).

    I produced the following Oct’07 map for the Arctic circle (Lat range: 66N-90N) from that link: seen here. But, to get October means for the 1971-2000 period, it seems that I may need to learn how to manipulate netcdf files. Will report if and once I manage to get do that.

    Gareth, thanks for your thought on this. Would be interesting to see how the means for the upcoming ice-covered months would compare with the corresponding monthly means for the base period of ‘71-2k.

    Thanks. -NL

  • Tristram Shandy // November 8, 2007 at 3:30 am | Reply

    It might be instructive and just plain fun to look at a Hurst component.

    This might be fun to correlate with Temperatures

    http://nwis.waterdata.usgs.gov/nwis/peak/?site_no=05054000

  • JesusChristHimself // November 8, 2007 at 11:48 am | Reply

    “Any thougths on why St Louis and St Croix river beginning day of snowmelt trends are so different than the Red River basin trend? …”

    When I was growing up that area was perceived to be colder than the Dakotas.

  • George // November 8, 2007 at 4:37 pm | Reply

    “Any thoughts on why St Louis and St Croix river beginning day of snowmelt trends are so different than the Red River basin trend? …”

    It may be because the temperature at the headwaters of the St Louis and St Croix river’s are influenced by nearby Lake Superior.

    Lake superior has been warming in recent decades , but lake surface temp and air temp in the region do not stay in lock step. The relationship between the air temp and the surface water temp is a complex one – and very nonlinear (see below)

    But the temperature of the lake surface water nonetheless does have an effect on the air temp in the surrounding regions and it is at least plausible that this could affect the snowmelt date in the region.

    There is an interesting effect related to lake warming that shows a similar trend to that of snowmelt date over the past few decades .

    The temp of lake superior is heavily influenced in early summer by what is know as “Spring Overturn.”

    Up until that time (beginning of July) lake surface temp ramps up fairly slowly. After overturn it warms very quickly after spring overturn.
    spring overturn has been occurring about 1/2 day earlier each year over the past few decades, just slightly less than the 0.62 days/yr that Tamino found for the Red River above. )

    It should be noted that Spring overturn occurs much later than the beginning of snowmelt, so the overturn effect is not itself causing the snowmelt.

  • Kelly O'Day // December 4, 2007 at 3:14 pm | Reply

    Tamino:

    I found your discussion on segmented regression very interesting. So interesting, in fact, that I developed an Excel workbook to provide segmented – piecewise regression. It includes a manual approach similar to the method you described to Tristram Shandy above .

    I also developed a minimization solution using Excel’s Solver.

    The workbook, with the Red River data, is available here

    Kelly

    [Response: I'm impressed. And this is not the first time you've shown considerable ingenuity using available tools to accomplish analytical goals.]

Leave a Comment