{"id":16434,"date":"2011-08-10T10:24:40","date_gmt":"2011-08-10T16:24:40","guid":{"rendered":"http:\/\/rankexploits.com\/musings\/?p=16434"},"modified":"2011-08-10T16:46:14","modified_gmt":"2011-08-10T22:46:14","slug":"nh-ice-update-vulnerable-ice-method","status":"publish","type":"post","link":"https:\/\/rankexploits.com\/musings\/2011\/nh-ice-update-vulnerable-ice-method\/","title":{"rendered":"NH Ice Update &#038; Vulnerable Ice Method"},"content":{"rendered":"<p>Time for the scheduled Wednesday NH ice update which includes the new method I concocted to predict the minimum.  Based on my full new method, my best estimate of the 7 day minimum ice extent is currently 4.64 million km^2 and the probability of breaking the breaking the 2007 sea ice minimum is approximately 9%.  Because the method is involves weighting results of two candidate models by their relative probabilities based on the corrected Aikake criteria, the explanation of how I estimate this requires multiple blog post. One post will explain the estimate based on the candidate model that appears most likely; other posts the less likely model and also explains how I weight the predictions from two models and estimate the uncertainty in the weighted. <\/p>\n<p>The new method was motivated mostly by my desire to try to have better and better short term models as we move toward the time of the NH sea ice minimum. More specifically, the choice of predictor a discussion with  who suggested <a href=\"http:\/\/neven1.typepad.com\/blog\/2011\/07\/sie-2011-update-14-two-steps-back-one-step.html?cid=6a0133f03a1e37970b01543440f4fa970c#comment-6a0133f03a1e37970b01543440f4fa970c\">CRANDALS<\/a> who suggested that current ice <em>area<\/em> would be a better predictor or the ultimate ice extent minimum than current ice extent.  <!--more--><\/p>\n<p>I admitted this might well be true. In my first cut for a short term fit, I&#8217;d limited myself to finding the best fit based on daily sea ice extent only and using only 1 predictor.  This results in fits like the ones discussed in my previous <a href=\"http:\/\/rankexploits.com\/musings\/2011\/wed-nh-ice-update\/\">NH Ice Extent Update<\/a>.  <\/p>\n<p>As I would like to improve that model by considering current ice area, I asked for a link to <i>daily<\/i> NH ice area data. Neven pointed me to the <a href=\"http:\/\/arctic.atmos.uiuc.edu\/cryosphere\/timeseries.anom\">Cryosphere Today<\/a> NH ice area data starting in 1970. Yesterday, I fiddled around to find the best predictor of the 7-day smooth JAXA minimum ice extent using limiting analysis to data since 1979.  For ice extent, I used a merge of JAXA and NASA-GSFC. For ice area I used the Cryosphere today.<\/p>\n<p>Once I had the data, I needed to think of candidate models to test.  Currently, loads of people are trying to predict the NH Ice Extent minimum, and each discussed statistical models with the predictors based on different notions of what <i>ought<\/i> to plausibly predict ice extent or ice melt rate. My initial goal is to consider candidate models that use <i>one<\/i> predictor to predict the ice extent.  That is: For now, I am avoiding multiple regressions.  <\/p>\n<p>With the available data, I computed the corrected Aikake criteria (AICc) for simple linear regressions for the following candidate models: <\/p>\n<ol>\n<li>Minimum 7-day Ice Extent predicted by Current 7-day Ice Extent. Motivation: for short term predictions, the best predictor will be the current value of the thing predicted. AICc=901.9<\/li>\n<li>Minimum 7-day Ice Extent predicted by Current 7-day Ice Area. Motivation: For slightly longer term predictions, the area gives a best estimate of the &#8220;less vulnerable&#8221; ise.  AICc=894.2<\/li>\n<li>Remaining &#8211; Extent Losses  predicted by the <i>difference<\/i> between the Current 7-day Ice Extent and the Current 7-day Ice Area.  Here I define &#8220;Remaining &#8211; Extent Losses&#8221; as the difference between the 7 day-average NH ice extent on the current day number counted from Jan 1 and the minimum for a particular year. I will call the difference between the ice extent and ice area the amount of &#8220;vulnerable ice&#8221;.   Note that using this model, the best estimate of the ice minimum is the difference between the current NH ice extent&#8211; which is known&#8211; and the predicted melt.   Motivation: The difference between extent and area describes the amount of &#8220;vulnerable ice&#8221;. AICc =892.0 <\/li>\n<\/ol>\n<p>Based on AICc, of the three candidate models, predicting remaining extent loss as a function of the amount of vulnerable ice is the most efficient. Note that I find CRANDALS is correct that the current area gives a better prediction of the NH Ice Extent minimum&#8211; though I suspect as we near the actual minimum this must change. In the limit that we reach the minimum in the ice extent, the current ice extent must give a very good prediction!  Still, currently: the best predictor is amount of vulnerable ice, the second best is area. <\/p>\n<p>The fit based on vulnerable ice is shown below:<br \/>\n<a href=\"http:\/\/rankexploits.com\/musings\/wp-content\/uploads\/2011\/08\/PredictNHIceLoss.png\"><img loading=\"lazy\" decoding=\"async\" src=\"http:\/\/rankexploits.com\/musings\/wp-content\/uploads\/2011\/08\/PredictNHIceLoss-500x500.png\" alt=\"\" title=\"PredictNHIceLoss\" width=\"500\" height=\"500\" class=\"aligncenter size-medium wp-image-16450\" srcset=\"https:\/\/rankexploits.com\/musings\/wp-content\/uploads\/2011\/08\/PredictNHIceLoss-500x500.png 500w, https:\/\/rankexploits.com\/musings\/wp-content\/uploads\/2011\/08\/PredictNHIceLoss-300x300.png 300w, https:\/\/rankexploits.com\/musings\/wp-content\/uploads\/2011\/08\/PredictNHIceLoss.png 1008w\" sizes=\"auto, (max-width: 500px) 100vw, 500px\" \/><\/a><br \/>\nThe open black circles are data for remaining losses as a function of the amount of vulnerable ice. Note that greater losses are correlated with larger amounts of vulnerable ice.  The solid blue line indicates the best fit estimate for the remaining extent loss based on the current amount of vulnerable ice; the dashed lines indicate the \u00c2\u00b195% scatter about the mean assuming the residuals are normally distributed (which they appear not to be).  The filled orange circle indicates best estimate of the ice loss rate based on the current amount of vulnerable ice (i.e. a remaining loss of 1.55 million square km); the open circles indicate the \u00c2\u00b195% confidence intervals based on scatter about the best fit curve and including the uncertainty in the trend and constant term for the model. <\/p>\n<p>The  (+) indicates extent loss, that if achieved will result in a tie with the 2010 ice extent. The (x) indicates the extent loss that will result in a tie with the 2007 minimum.  Both fall inside the \u00c2\u00b195% uncertainty intervals indicating that, according to this model, we may well see a new record low NH ice extent, or we may see the NH ice extent above that seen in 2010.  Based on this model, the probability of a new low is estimated at 7% while the probability of a  NH ice extent that exceeds 2010&#8217;s value is 28%. These are noted in black in the lower right hand corner of the figure. <\/p>\n<p>Some of you will note that other probabilities are listed in blue. The second set of probabilities takes into account that while the &#8220;vulnerable ice&#8221; model appears to explain the past data better than the other two models, it is not sufficiently better to discount that possibility that the other models are in reality better predictors. Taking predictions of the other models into account, I estimate the probability of a new minimum at 9%.  I will discuss the estimate based on the other two fits and how I weighted the two estimates in a future post planned within 2 days.<\/p>\n<p>If you have other methods of predicting the current ice, or know of other troves of daily data, point me to them. Right now I know lots of people are fascinated by the horserace, and I like to quantify what I think is going to happen.  Even if I guess wrong, I enjoy this more than just looking at the JAXA chart. <\/p>\n<p><b>Links to data sources:<\/b><\/p>\n<ol>\n<li><a href=\"http:\/\/www.ijis.iarc.uaf.edu\/seaice\/extent\/plot.csv\">JAXA<\/a> ice extent<\/li>\n<li><a href=\"http:\/\/polynya.gsfc.nasa.gov\/datasets\/nh_daily_observed_or_interpolated_sie_1972_2002.txt\">gsfc.nasa.gov<\/a><\/li>\n<li><a href=\"http:\/\/arctic.atmos.uiuc.edu\/cryosphere\/timeseries.anom\">Cryosphere Today<\/a>.<\/li>\n<\/ol>\n","protected":false},"excerpt":{"rendered":"<p>Time for the scheduled Wednesday NH ice update which includes the new method I concocted to predict the minimum. Based on my full new method, my best estimate of the 7 day minimum ice extent is currently 4.64 million km^2 and the probability of breaking the breaking the 2007 sea ice minimum is approximately 9%. &hellip; <a href=\"https:\/\/rankexploits.com\/musings\/2011\/nh-ice-update-vulnerable-ice-method\/\" class=\"more-link\">Continue reading <span class=\"screen-reader-text\">NH Ice Update &#038; Vulnerable Ice Method<\/span> <span class=\"meta-nav\">&rarr;<\/span><\/a><\/p>\n","protected":false},"author":2,"featured_media":0,"comment_status":"closed","ping_status":"closed","sticky":false,"template":"","format":"standard","meta":{"footnotes":""},"categories":[49],"tags":[334],"class_list":["post-16434","post","type-post","status-publish","format-standard","hentry","category-betting","tag-nh-ice"],"_links":{"self":[{"href":"https:\/\/rankexploits.com\/musings\/wp-json\/wp\/v2\/posts\/16434","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/rankexploits.com\/musings\/wp-json\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/rankexploits.com\/musings\/wp-json\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/rankexploits.com\/musings\/wp-json\/wp\/v2\/users\/2"}],"replies":[{"embeddable":true,"href":"https:\/\/rankexploits.com\/musings\/wp-json\/wp\/v2\/comments?post=16434"}],"version-history":[{"count":0,"href":"https:\/\/rankexploits.com\/musings\/wp-json\/wp\/v2\/posts\/16434\/revisions"}],"wp:attachment":[{"href":"https:\/\/rankexploits.com\/musings\/wp-json\/wp\/v2\/media?parent=16434"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/rankexploits.com\/musings\/wp-json\/wp\/v2\/categories?post=16434"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/rankexploits.com\/musings\/wp-json\/wp\/v2\/tags?post=16434"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}