Project

General

Profile

« Previous | Next » 

Revision 5600791d

Added by Adam M. Wilson almost 11 years ago

adding figures for initial LST temporal window evaluation 5

View differences:

climate/research/LST_evaluation.Rmd
89 89 
```{r,fig.height=3}
90 90 
histogram(dlst2$lst-dlst2$lst2m,col="grey",xlab="Anomolies (1 month - 2 month means)")
91 91 
```
92 
The 95th quantile of that is only `r round(quantile(abs(dlst2$lst-dlst2$lst2m),na.rm=T,0.95),1)`, so the differences are quite small. From this analysis, it appears that broadening the temporal window window will maintain relativly consistent estimate of LST (R^2 ranged from 0.88-0.9) even when using only data from the surrounding months.
92 
The 95th quantile of the absolute value is only `r round(quantile(abs(dlst2$lst-dlst2$lst2m),na.rm=T,0.95),1)`, so the differences are quite small. From this analysis, it appears that broadening the temporal window will maintain a relatively consistent estimate of LST (R^2 ranged from 0.88-0.9) even when using only data from the surrounding months.
93 93 

  		




  94
  94
  
    
Let's see how the seasonal cycle is represented by these proxies for a few randomly selected stations.  Here the red line is the observed TMax data, the heavy black line represents the mean LST in that month, and the green line is a three-month moving window, while the purple line is the 2-month window (not including the month of interest).


  		




  95
  95
  
    

  		












  

    
      climate/research/LST_evaluation.html
    
  





  181
  181
  
    
<body>


  		




  182
  182
  
    
<h1>LST Temporal Aggregation Evaluation</h1>


  		




  183
  183
  
    

  		




  184
  
  
    
<h3>Adam M. Wilson and Giuseppe Amatulli (Compiled on Mon Mar 31 17:07:34 2014  using code hash: 4e82a75)</h3>


  		




  
  184
  
    
<h3>Adam M. Wilson and Giuseppe Amatulli (Compiled on Mon Mar 31 17:18:07 2014  using code hash: fdf72e9)</h3>


  		




  185
  185
  
    

  		




  186
  186
  
    
<p>A short script to explore the implications of a 3-month moving window Land Surface Temperature Climatology.</p>


  		




  187
  187
  
    

  		




  ......




  259
  259
  
    
<h1>TMax~LST relationship seasonal variability</h1>


  		




  260
  260
  
    

  		




  261
  261
  
    
<p>First let&#39;s explore the variability of the TMax~LST relationship by month, grouped by tile.  In this plot grey lines indicate a Tmax~LST regression within each tile (stations may be present in multiple tiles). Variability among grey lines represents spatial variability in the fitted lapse rate and the heavy black line is overall mean relationship (by month).


  		




  262
  
  
    
<img src="http://i.imgur.com/dLVOWhY.png" alt="plot of chunk unnamed-chunk-4"/> </p>


  		




  
  262
  
    
<img src="http://i.imgur.com/fQUQNDv.png" alt="plot of chunk unnamed-chunk-4"/> </p>


  		




  263
  263
  
    

  		




  264
  264
  
    
<h1>Comparison of monthly means with 3-month moving window</h1>


  		




  265
  265
  
    

  		




  266
  266
  
    
<p>Here we compare the monthly means with a three month moving window (e.g. January mean LST with December-February mean LST).  Note that the relationship is very good (R<sup>2</sup> &gt;0.95) but slightly weaker in winter months, likely due primarily to seasonal minimums.  Heavy black line is 1:1, and thin line is the fitted regression.</p>


  		




  267
  267
  
    

  		




  268
  
  
    
<p><img src="http://i.imgur.com/FdyccKW.png" alt="plot of chunk unnamed-chunk-5"/> </p>


  		




  
  268
  
    
<p><img src="http://i.imgur.com/d2tZz4H.png" alt="plot of chunk unnamed-chunk-5"/> </p>


  		




  269
  269
  
    

  		




  270
  270
  
    
<h1>Comparison of monthly means with 2-month moving window that does not include the month</h1>


  		




  271
  271
  
    

  		




  272
  272
  
    
<p>Here we compare the monthly means with a two month moving window that does not include the month of interest (e.g. January mean LST with December and February mean LST, but not including the January LST).  Heavy black line is 1:1, and thin line is the fitted regression.</p>


  		




  273
  273
  
    

  		




  274
  
  
    
<p><img src="http://i.imgur.com/PGSvyuj.png" alt="plot of chunk unnamed-chunk-6"/> </p>


  		




  
  274
  
    
<p><img src="http://i.imgur.com/C79D3K2.png" alt="plot of chunk unnamed-chunk-6"/> </p>


  		




  275
  275
  
    

  		




  276
  276
  
    
<p>Now let&#39;s look at the differences between the 1-month and 2-month LST values.  These represent a measure of how wrong we would be if we only had data from the two surrounding months and not the month in question.  </p>


  		




  277
  277
  
    

  		




  278
  278
  
    
<pre><code class="r">histogram(dlst2$lst - dlst2$lst2m, col = &quot;grey&quot;, xlab = &quot;Anomolies (1 month - 2 month means)&quot;)


  		




  279
  279
  
    
</code></pre>


  		




  280
  280
  
    

  		




  281
  
  
    
<p><img src="http://i.imgur.com/60NdxK7.png" alt="plot of chunk unnamed-chunk-7"/> </p>


  		




  
  281
  
    
<p><img src="http://i.imgur.com/TsHkdVv.png" alt="plot of chunk unnamed-chunk-7"/> </p>


  		




  282
  282
  
    

  		




  283
  
  
    
<p>The 95th quantile of that is only 4.2, so the differences are quite small. From this analysis, it appears that broadening the temporal window window will maintain relativly consistent estimate of LST (R<sup>2</sup> ranged from 0.88-0.9) even when using only data from the surrounding months.</p>


  		




  
  283
  
    
<p>The 95th quantile of the absolute value is only 4.2, so the differences are quite small. From this analysis, it appears that broadening the temporal window will maintain a relatively consistent estimate of LST (R<sup>2</sup> ranged from 0.88-0.9) even when using only data from the surrounding months.</p>


  		




  284
  284
  
    

  		




  285
  285
  
    
<p>Let&#39;s see how the seasonal cycle is represented by these proxies for a few randomly selected stations.  Here the red line is the observed TMax data, the heavy black line represents the mean LST in that month, and the green line is a three-month moving window, while the purple line is the 2-month window (not including the month of interest).</p>


  		




  286
  286
  
    

  		




  ......




  289
  289
  
    
ss = sample(dlst2l$station, 10)


  		




  290
  290
  
    
</code></pre>


  		




  291
  291
  
    

  		




  292
  
  
    
<p><img src="http://i.imgur.com/1JNQ2vn.png" alt="plot of chunk unnamed-chunk-9"/> </p>


  		




  
  292
  
    
<p><img src="http://i.imgur.com/XvvNDr5.png" alt="plot of chunk unnamed-chunk-9"/> </p>


  		




  293
  293
  
    

  		




  294
  294
  
    
<h1>Remaining Questions</h1>


  		




  295
  295
  
    

  		












  

    
      climate/research/LST_evaluation.md
    
  





  4
  4
  
    

  		




  5
  5
  
    

  		




  6
  6
  
    

  		




  7
  
  
    
### Adam M. Wilson and Giuseppe Amatulli (Compiled on Mon Mar 31 17:07:34 2014  using code hash: 4e82a75)


  		




  
  7
  
    
### Adam M. Wilson and Giuseppe Amatulli (Compiled on Mon Mar 31 17:18:07 2014  using code hash: fdf72e9)


  		




  8
  8
  
    

  		




  9
  9
  
    
A short script to explore the implications of a 3-month moving window Land Surface Temperature Climatology.


  		




  10
  10
  
    

  		




  ......




  28
  28
  
    
# TMax~LST relationship seasonal variability 


  		




  29
  29
  
    

  		




  30
  30
  
    
First let's explore the variability of the TMax~LST relationship by month, grouped by tile.  In this plot grey lines indicate a Tmax~LST regression within each tile (stations may be present in multiple tiles). Variability among grey lines represents spatial variability in the fitted lapse rate and the heavy black line is overall mean relationship (by month).


  		




  31
  
  
    
![plot of chunk unnamed-chunk-4](http://i.imgur.com/dLVOWhY.png) 


  		




  
  31
  
    
![plot of chunk unnamed-chunk-4](http://i.imgur.com/fQUQNDv.png) 


  		




  32
  32
  
    

  		




  33
  33
  
    

  		




  34
  34
  
    

  		




  35
  35
  
    
# Comparison of monthly means with 3-month moving window


  		




  36
  36
  
    
Here we compare the monthly means with a three month moving window (e.g. January mean LST with December-February mean LST).  Note that the relationship is very good (R^2 >0.95) but slightly weaker in winter months, likely due primarily to seasonal minimums.  Heavy black line is 1:1, and thin line is the fitted regression.


  		




  37
  37
  
    

  		




  38
  
  
    
![plot of chunk unnamed-chunk-5](http://i.imgur.com/FdyccKW.png) 


  		




  
  38
  
    
![plot of chunk unnamed-chunk-5](http://i.imgur.com/d2tZz4H.png) 


  		




  39
  39
  
    

  		




  40
  40
  
    

  		




  41
  41
  
    
# Comparison of monthly means with 2-month moving window that does not include the month


  		




  42
  42
  
    
Here we compare the monthly means with a two month moving window that does not include the month of interest (e.g. January mean LST with December and February mean LST, but not including the January LST).  Heavy black line is 1:1, and thin line is the fitted regression.


  		




  43
  43
  
    

  		




  44
  
  
    
![plot of chunk unnamed-chunk-6](http://i.imgur.com/PGSvyuj.png) 


  		




  
  44
  
    
![plot of chunk unnamed-chunk-6](http://i.imgur.com/C79D3K2.png) 


  		




  45
  45
  
    

  		




  46
  46
  
    

  		




  47
  47
  
    
Now let's look at the differences between the 1-month and 2-month LST values.  These represent a measure of how wrong we would be if we only had data from the two surrounding months and not the month in question.  


  		




  ......




  51
  51
  
    
histogram(dlst2$lst - dlst2$lst2m, col = "grey", xlab = "Anomolies (1 month - 2 month means)")


  		




  52
  52
  
    
```


  		




  53
  53
  
    

  		




  54
  
  
    
![plot of chunk unnamed-chunk-7](http://i.imgur.com/60NdxK7.png) 


  		




  
  54
  
    
![plot of chunk unnamed-chunk-7](http://i.imgur.com/TsHkdVv.png) 


  		




  55
  55
  
    

  		




  56
  
  
    
The 95th quantile of that is only 4.2, so the differences are quite small. From this analysis, it appears that broadening the temporal window window will maintain relativly consistent estimate of LST (R^2 ranged from 0.88-0.9) even when using only data from the surrounding months.


  		




  
  56
  
    
The 95th quantile of the absolute value is only 4.2, so the differences are quite small. From this analysis, it appears that broadening the temporal window will maintain a relatively consistent estimate of LST (R^2 ranged from 0.88-0.9) even when using only data from the surrounding months.


  		




  57
  57
  
    

  		




  58
  58
  
    
Let's see how the seasonal cycle is represented by these proxies for a few randomly selected stations.  Here the red line is the observed TMax data, the heavy black line represents the mean LST in that month, and the green line is a three-month moving window, while the purple line is the 2-month window (not including the month of interest).


  		




  59
  59
  
    

  		




  ......




  65
  65
  
    
```


  		




  66
  66
  
    

  		




  67
  67
  
    

  		




  68
  
  
    
![plot of chunk unnamed-chunk-9](http://i.imgur.com/1JNQ2vn.png) 


  		




  
  68
  
    
![plot of chunk unnamed-chunk-9](http://i.imgur.com/XvvNDr5.png) 


  		




  69
  69
  
    

  		




  70
  70
  
    

  		




  71
  71
  
    

  		












Also available in:  Unified diff