<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=utf-8">
</head>
<body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class="">
<div class=""><br class="">
</div>
Hi Leo,
<div class=""><br class="">
</div>
<div class="">Thanks for the update.  Were you surprised that melting decreased when you turned off LW downscaling?  This suggests that some of the lower-elevation glacier ECs rely on the downscaled LW to melt.  Would you like to go ahead and try the ice-column-only
 downscaling?</div>
<div class=""><br class="">
</div>
<div class="">I haven’t been able to think of a reason why the LW effect is much bigger in coupled runs than in offline runs.</div>
<div class=""><br class="">
</div>
<div class="">A couple of questions:  First, could you remind me of the definitions of QSNOMELT, QSNOMELT_ICE and QICE_MELT?  Is QSNOMELT_ICE the refreezing?  Second, what is our target for QICE_MELT?  It seems like we’re still way below where we want to be,
 even in B3 with the adjusted ramp, the lowered ice albedo, and the SCF bug fix.  Is there any significant increase in melting over time?</div>
<div class=""><br class="">
</div>
<div class="">Given that the surface forcing seems reasonable overall, I can’t help wondering if we’re missing something important.  What was the last simulation with a reasonable amount of surface melting?  Do you think we’re getting too much refreezing, maybe
 as a result of missing physics?</div>
<div class=""><br class="">
</div>
<div class="">Thanks,</div>
<div class=""><br class="">
</div>
<div class="">Bill L.</div>
<div class=""><br class="">
</div>
<div class=""><br class="">
<div>
<blockquote type="cite" class="">
<div class="">On Jan 24, 2017, at 2:13 PM, Kampenhout, L. van (Leo) <<a href="mailto:L.vanKampenhout@uu.nl" class="">L.vanKampenhout@uu.nl</a>> wrote:</div>
<br class="Apple-interchange-newline">
<div class="">
<div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class="">
Hi Bill and others,
<div class=""><br class="">
</div>
<div class="">I’ve been doing some more experiments with LW downscaling:</div>
<div class=""><a href="https://docs.google.com/document/d/1l6SU0PPba8eOuB8fiI3LYVDdyDMrwZkQf2lyh7Ztaqw/edit#" class="">https://docs.google.com/document/d/1l6SU0PPba8eOuB8fiI3LYVDdyDMrwZkQf2lyh7Ztaqw/edit#</a></div>
<div class=""><br class="">
</div>
<div class="">As I mentioned in a previous email, I found that it has little effect on the energy balance in offline runs (only about 5 Gt/yr). </div>
<div class="">However in coupled mode, there is a big response when I turn off LW downscaling : it is halved (cf. run B3 and B4). </div>
<div class=""><br class="">
</div>
<div class="">I can’t really explain this at the moment, it could even be caused by natural variability , but I suspect a real difference between coupled and offline runs. </div>
<div class="">
<div class="">I’ll keep working on this a bit more, any suggestions welcome! </div>
</div>
<div class=""><br class="">
</div>
<div class="">Leo</div>
<div class=""><br class="">
</div>
<div class=""><br class="">
</div>
<div class=""><br class="">
<div class="">
<blockquote type="cite" class="">
<div class="">On 23 Jan 2017, at 16:48, Lipscomb, William Henry <<a href="mailto:lipscomb@lanl.gov" class="">lipscomb@lanl.gov</a>> wrote:</div>
<br class="Apple-interchange-newline">
<div class="">
<div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class="">
<div class="">Hi Leo,</div>
<div class=""><br class="">
</div>
<div class="">Your finding about the downscaled LW is very interesting.  When we added the LW downscaling option, I was imagining giving greater LW (and greater melt) to lower-elevation ECs, relative to high-elevation ECs.  I think this suggestion originally
 came from Jeff Ridley at the Hadley Centre.  But I can see how the downscaling would redistribute LW energy from glaciated to bare-land columns, which could be a bad thing.</div>
<div class=""><br class="">
</div>
<div class="">I’d suggest that we first try turning off the LW radiation downscaling entirely, and see what happens.  Then we can come back to the idea of an ice-column-only downscaling.  I would expect that an ice-column-only downscaling would maximize melt
 in low-lying glacier ECs, though as you say the physical basis for this (without a corresponding treatment of SW) is a little shaky.</div>
<div class=""><br class="">
</div>
<div class="">Thanks,</div>
<div class=""><br class="">
</div>
<div class="">Bill L.</div>
<br class="">
<div class="">
<blockquote type="cite" class="">
<div class="">On Jan 23, 2017, at 8:13 AM, Bill Sacks <<a href="mailto:sacks@ucar.edu" class="">sacks@ucar.edu</a>> wrote:</div>
<br class="Apple-interchange-newline">
<div class="">
<div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class="">
Hi Leo,
<div class=""><br class="">
</div>
<div class="">Okay, thanks for the explanation of your thoughts. I'll let Bill L chime in here if he has anything more to add, but my understanding is that the current parameterization implicitly assumes that cloud fraction, cloud height and most other properties
 such as cloud liquid water content are the same for different elevation classes / bare land. In that case, downwelling longwave radiation is a function of cloud temperature. The current adjustment is based on differences in surface temperature, which won't
 be exactly the same as cloud temperature, but seems like a reasonable approximation given the limited information we have.</div>
<div class=""><br class="">
</div>
<div class="">Do you see problems with this?</div>
<div class=""><br class="">
</div>
<div class="">If you feel this is problematic, then I'd probably argue for just turning off longwave radiation downscaling entirely, rather than applying it only over ice-covered regions. This can be done via the namelist parameter, glcmec_downscale_longwave.</div>
<div class=""><br class="">
</div>
<div class="">Thanks,</div>
<div class="">Bill S</div>
<div class=""><br class="">
<div class="">
<blockquote type="cite" class="">
<div class="">On Jan 23, 2017, at 6:33 AM, Kampenhout, L. van (Leo) <<a href="mailto:L.vanKampenhout@uu.nl" class="">L.vanKampenhout@uu.nl</a>> wrote:</div>
<br class="Apple-interchange-newline">
<div class="">
<div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class="">
Hi Bill,
<div class=""><br class="">
</div>
<div class="">
<blockquote type="cite" class="">
<div class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;">
<div class="">I'm adding liwg-core to the cc list here, to include others in the discussion.<br class="">
</div>
</div>
</blockquote>
</div>
<div class="">
<div class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;">
<div class=""><br class="">
</div>
</div>
</div>
<div class="">that is probably a good idea, thanks. </div>
<div class=""><br class="">
</div>
<div class="">
<blockquote type="cite" class="">
<div class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;">
<div class="">
<div class="">Within the CISM domain, the elevation of the bare land portion of the grid cell comes from CISM. Jeremy and I discussed this extensively, and this felt like the most consistent way to handle downscaling. Are you suggesting continuing to do other
 aspects of the downscaling over bare land columns, but not the longwave downscaling? What would the physical justification be for not doing the longwave downscaling over the bare land columns?</div>
</div>
</div>
</blockquote>
</div>
<div class="">
<div class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;">
<div class="">
<div class=""><br class="">
</div>
<div class="">Downwelling longwave radiation is controlled primarily by the existence and thickness of clouds. So by downscaling longwave radiation you’re downscaling a cloud property. However, clouds also modulate shortwave radiation and this is not accounted
 for, which makes the current approach inconsistent. I believe this inconsistency currently leads to a net loss of radiation over the ice sheet, with a considerable magnitude (~ 0.5 W/m2 year round, equivalent to about 80 Gt/year of melt).  </div>
<div class=""><br class="">
</div>
<div class="">Now that I think about it, an alternative approach could be to introduce a lapse rate for SW transmittance, based on .e.g. Figure 6 of <a href="http://www.the-cryosphere.net/10/2379/2016/tc-10-2379-2016.pdf" class="">this paper</a>. This could
 compensate the LW radiation loss by an increase in SW radiation. Unsure how to relate the SW transmittance to SW downwelling though. </div>
<div class=""><br class="">
</div>
<div class="">Leo</div>
<div class=""><br class="">
</div>
</div>
</div>
</div>
<div class=""><br class="">
</div>
<div class=""><br class="">
</div>
<div class="">
<div class="">
<blockquote type="cite" class="">
<div class="">On 23 Jan 2017, at 13:43, Bill Sacks <<a href="mailto:sacks@ucar.edu" class="">sacks@ucar.edu</a>> wrote:</div>
<br class="Apple-interchange-newline">
<div class="">
<div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class="">
Hi Leo,
<div class=""><br class="">
</div>
<div class="">I'm adding liwg-core to the cc list here, to include others in the discussion.<br class="">
<div class=""><br class="">
</div>
<div class="">Within the CISM domain, the elevation of the bare land portion of the grid cell comes from CISM. Jeremy and I discussed this extensively, and this felt like the most consistent way to handle downscaling. Are you suggesting continuing to do other
 aspects of the downscaling over bare land columns, but not the longwave downscaling? What would the physical justification be for not doing the longwave downscaling over the bare land columns?</div>
<div class=""><br class="">
</div>
<div class="">Bill S</div>
<div class=""><br class="">
</div>
<div class="">
<div class="">
<blockquote type="cite" class="">
<div class="">On Jan 23, 2017, at 5:30 AM, Kampenhout, L. van (Leo) <<a href="mailto:L.vanKampenhout@uu.nl" class="">L.vanKampenhout@uu.nl</a>> wrote:</div>
<br class="Apple-interchange-newline">
<div class="">
<div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class="">
Hi Bill,
<div class=""><br class="">
</div>
<div class="">I was wondering about the topographic height of the bare land portion of the grid cell. Is this set by CISM or is is equal to the grid cell mean height? Either way, I think in general this height will be lower than that of the glaciated portion,
 so more radiation goes the the bare land (low elev) than to the glacier part (high elev). To my current understanding, this is a design decision with big implications ( -5 W/m2 over the entire ice sheet) so we may need to reconsider this. My proposal would
 be to turn off downscaling for the bare land columns, just for the glaciated columns. Has this ever been tested? If not, I would like to test the effect of this. </div>
<div class=""><br class="">
</div>
<div class="">Thanks,</div>
<div class="">Leo</div>
<div class=""><br class="">
</div>
<div class=""><br class="">
</div>
<div class=""><br class="">
<div class="">
<div class="">
<blockquote type="cite" class="">
<div class="">On 21 Jan 2017, at 14:41, Bill Sacks <<a href="mailto:sacks@ucar.edu" class="">sacks@ucar.edu</a>> wrote:</div>
<br class="Apple-interchange-newline">
<div class="">
<div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class="">
Hi Leo,
<div class=""><br class="">
</div>
<div class="">If I understand the problem right: You need to look at grid cell means for the longwave radiation downscaling to conserve: the adjustments are made over the bare land portion of the grid cell, too, so those need to be included in your averages
 to get conservation.</div>
<div class=""><br class="">
</div>
<div class="">Bill S</div>
<div class=""><br class="">
<div class="">
<blockquote type="cite" class="">
<div class="">On Jan 20, 2017, at 6:12 PM, Kampenhout, L. van (Leo) <<a href="mailto:L.vanKampenhout@uu.nl" class="">L.vanKampenhout@uu.nl</a>> wrote:</div>
<br class="Apple-interchange-newline">
<div class="">
<div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class="">
Hi Bill,
<div class=""><br class="">
</div>
<div class="">
<div class="">
<blockquote type="cite" class="">
<div class="">
<div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class="">
<div class=""><br class="">
</div>
<div class="">To make sure I understand:  It’s currently the case that in glacier columns with fractional snow cover, any shortwave radiation that ought to be going into the ice (or “soil”) is simply thrown away?  And the fix, at least for now, is to apply
 all this radiation to the snow cover until the snow is gone?</div>
</div>
</div>
</blockquote>
<div class=""><br class="">
</div>
<div class="">That is correct. Following this reasoning, the imbalance is greatest when the fractional snow cover is smallest, right before bare ice appears. </div>
<br class="">
<blockquote type="cite" class="">
<div class="">
<div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class="">
<div class=""></div>
<div class="">Looking at your pdf plot, it seems there is a mean ice-sheet-wide imbalance of ~0.1 W/m2 during summer and fall.  Based on a back-of-the-envelope calculation, that translates to a melting deficit of ~5 Gt/yr, which is less than I was hoping.  But
 obviously the imbalance is much bigger locally, which could lead to important feedbacks.  Also, with it being late on a Friday, I could have made a mistake somewhere on my envelope :-)</div>
<div class=""><br class="">
</div>
</div>
</div>
</blockquote>
<div class=""><br class="">
</div>
<div class="">I get the same value as you do, with the same caveat. </div>
<div class=""><br class="">
</div>
<div class="">There was a second finding in my exercise which I did not mention in my email because I’m unsure what to think of this yet. That is, the longwave downscaling also causes a SEB residual but only when integrating the _ICE energy flux variables,
 not with the grid-cell average variables. According to the code documentation, this cannot be true because there is a normalisation in place that forces the downscaled longwave to exactly equal the gridcell value. However, when I turn off LW downscaling the
 residual disappears. </div>
<div class=""><br class="">
</div>
<div class="">That said, the residual is about 0.5 W/m2 <i class="">throughout the year</i> with your Stefan-Boltzmann downscaling approach, and over 0.6 W/m2 with a lapse rate approach. Using the same envelope as before, this translates to about 80-100 Gt/yr
 of additional melt that we’re missing out on, excluding feedbacks. I don’t have time to look into this over the weekend, but I’m happy to do some more analysis next week.</div>
<div class=""><br class="">
</div>
<div class="">Cheers,</div>
<div class="">Leo</div>
<div class=""><br class="">
</div>
<br class="">
<blockquote type="cite" class="">
<div class="">
<div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class="">
<div class=""></div>
<div class="">Thanks again,</div>
<div class=""><br class="">
</div>
<div class="">Bill L.</div>
<div class=""><br class="">
</div>
<div class=""><br class="">
</div>
<div class="">
<blockquote type="cite" class="">
<div class="">On Jan 20, 2017, at 3:21 PM, David Lawrence <<a href="mailto:dlawren@ucar.edu" class="">dlawren@ucar.edu</a>> wrote:</div>
<br class="Apple-interchange-newline">
<div class="">
<div dir="ltr" class="">
<div class="">
<div class="">
<div class="">Thanks Leo,<br class="">
<br class="">
</div>
I read through your presentation and I think your logic is sound.  But, as Bill suggests, Sean is really the one who should review this because he knows this code really well.  I think he is in a meeting this afternoon, but I hope he can look at this on Monday
 at the latest.<br class="">
<br class="">
</div>
Cheers,<br class="">
<br class="">
</div>
Dave<br class="">
</div>
<div class="gmail_extra"><br class="">
<div class="gmail_quote">On Fri, Jan 20, 2017 at 2:55 PM, Bill Sacks <span dir="ltr" class="">
<<a href="mailto:sacks@ucar.edu" target="_blank" class="">sacks@ucar.edu</a>></span> wrote:<br class="">
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div style="word-wrap:break-word" class="">
<div class="">Hi Leo,</div>
<div class=""><br class="">
</div>
<div class="">Thank you very much for your careful analysis and detailed write-up of this problem. I have only skimmed through this and will let others – particularly Sean – comment in more detail.</div>
<div class=""><br class="">
</div>
<div class="">I'm not sure if this is related, but I have vague recollections of a discussion with Sean a few years ago about how the fractional snow parameterization should be treated for glacier_mec columns. I recall our feeling that at least some parts of
 this parameterization should not apply to glacier_mec columns, since they already have some treatment of subgrid snow variability through the multiple elevation class scheme. I can't recall any more details right now, though maybe Sean remembers more?</div>
<div class=""><br class="">
</div>
<div class="">Bill S</div>
<div class=""><br class="">
<div class="">
<blockquote type="cite" class="">
<div class="">
<div class="h5">
<div class="">On Jan 20, 2017, at 2:51 PM, Kampenhout, L. van (Leo) <<a href="mailto:L.vanKampenhout@uu.nl" target="_blank" class="">L.vanKampenhout@uu.nl</a>> wrote:</div>
<br class="m_825767671887645258Apple-interchange-newline">
</div>
</div>
<div class="">
<div class="">
<div class="h5">
<div class="">
<div style="word-wrap:break-word" class="">
<div class="">Forgot to attach the source code that contains the fix</div>
<div class=""><br class="">
</div>
</div>
<div style="word-wrap:break-word" class=""><br class="">
<div class="">
<blockquote type="cite" class="">
<div class="">On 20 Jan 2017, at 22:32, Kampenhout, L. van (Leo) <<a href="mailto:L.vanKampenhout@uu.nl" target="_blank" class="">L.vanKampenhout@uu.nl</a>> wrote:</div>
<br class="m_825767671887645258x_Apple-interchange-newline">
<div class="">
<div style="word-wrap:break-word" class="">
<div class=""><br class="">
</div>
<div class="">Hi Dave, Sean and Bill,</div>
<div class="">CC others,</div>
<div class=""><br class="">
</div>
<div class="">While doing the surface energy analysis for Greenland, I discovered that the surface energy budget was not fully closed. This holds both at the grid-cell integrated level as well as the glacier level (for this, I use the  variables ending in ‘_ICE’). </div>
<div class=""><br class="">
</div>
<div class="">See attachment <b class="">seb.pdf, </b>panel SEB residual<b class="">.
</b>Note that the imbalance occurs in summer, during the melt season. </div>
<div class=""><br class="">
</div>
<div class="">I reproduced this problem in a single column run, the results are attached. I believe the problem lies with the implementation of fractional snow and how glaciers are handled, namely by setting frac_sno_eff = 1.0. I hope you agree with the analysis.
 I also suggested a workaround which appears to work without problems in another ongoing BG run (2+ years). We are seeing slightly higher melt rates now over Greenland, which is good. How substantial the changes are is left to be seen.</div>
<div class=""><br class="">
</div>
<div class="">In the long-term a better solution may be desired. </div>
<div class=""><br class="">
</div>
<div class="">Cheers,</div>
<div class="">Leo</div>
<div class=""><br class="">
</div>
<div class=""><br class="">
</div>
<div class=""><br class="">
</div>
<div class=""></div>
</div>
<span id="m_825767671887645258x_cid:C5CE0560-E425-401C-A296-9CCE4F96CF90" class=""><swrad_bug_glaciers.pptx></span>
<div style="word-wrap:break-word" class="">
<div class=""></div>
<div class=""><br class="">
</div>
<div class=""></div>
</div>
<span id="m_825767671887645258x_cid:01B0A788-0195-4313-85D1-47A5630904A4" class=""><seb.pdf></span>
<div style="word-wrap:break-word" class="">
<div class=""></div>
</div>
</div>
</blockquote>
</div>
<br class="">
</div>
</div>
</div>
</div>
<span id="m_825767671887645258cid:E79F6E50257DD6459C0E9EDDD314726D@uu.nl" class=""><SurfaceRadiationMod.F90></span></div>
</blockquote>
</div>
<br class="">
</div>
</div>
</blockquote>
</div>
<br class="">
</div>
</div>
</blockquote>
</div>
<br class="">
<div class="">
<div style="letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px; word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class="">
<div class="">---</div>
<div class="">William Lipscomb</div>
<div class="">Los Alamos National Laboratory</div>
<div class="">Group T-3, MS B216</div>
<div class="">Los Alamos, NM 87545</div>
<div class="">505-667-0395</div>
<div class=""><br class="">
</div>
</div>
<br class="Apple-interchange-newline">
<br class="Apple-interchange-newline">
</div>
<br class="">
</div>
</div>
</blockquote>
</div>
<br class="">
</div>
</div>
</div>
</blockquote>
</div>
<br class="">
</div>
</div>
</div>
</blockquote>
</div>
<br class="">
</div>
</div>
</div>
</div>
</blockquote>
</div>
<br class="">
</div>
</div>
</div>
</div>
</blockquote>
</div>
<br class="">
</div>
</div>
</div>
</blockquote>
</div>
<br class="">
</div>
</div>
</div>
</blockquote>
</div>
<br class="">
<div class="">
<div style="letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px; word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class="">
<div class="">---</div>
<div class="">William Lipscomb</div>
<div class="">Los Alamos National Laboratory</div>
<div class="">Group T-3, MS B216</div>
<div class="">Los Alamos, NM 87545</div>
<div class="">505-667-0395</div>
<div class=""><br class="">
</div>
</div>
<br class="Apple-interchange-newline">
<br class="Apple-interchange-newline">
</div>
<br class="">
</div>
</div>
</blockquote>
</div>
<br class="">
</div>
</div>
</div>
</blockquote>
</div>
<br class="">
<div class="">
<div style="color: rgb(0, 0, 0); letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px; word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class="">
<div class="">---</div>
<div class="">William Lipscomb</div>
<div class="">Los Alamos National Laboratory</div>
<div class="">Group T-3, MS B216</div>
<div class="">Los Alamos, NM 87545</div>
<div class="">505-667-0395</div>
<div class=""><br class="">
</div>
</div>
<br class="Apple-interchange-newline">
<br class="Apple-interchange-newline">
</div>
<br class="">
</div>
</body>
</html>