Urban-Meteorology-Reading / WRF-SUEWS

WRF-SUEWS coupling project
https://wrf-suews.readthedocs.org
MIT License
5 stars 2 forks source link

Urban haze: changing SUEWS to capture urban haze #54

Closed hamidrezaomidvar closed 4 years ago

sunt05 commented 5 years ago

better to be resolved at the SUEWS side.

suegrimmond commented 5 years ago

Yes – I have some suggestions

Best wishes Sue

Prof Sue Grimmond Dept. of Meteorology, University of Reading, Reading, RG6 6BB T: 44 118 378 6248 – messages get emailed to me O:Met Building (#58 on map) rm:1U14 E: c.s.grimmond@reading.ac.ukmailto:c.s.grimmond@reading.ac.uk W: http://micromet.reading.ac.uk/

From: Ting Sun notifications@github.com Sent: 02 August 2019 11:56 To: Urban-Meteorology-Reading/WRF-SUEWS WRF-SUEWS@noreply.github.com Cc: Subscribed subscribed@noreply.github.com Subject: Re: [Urban-Meteorology-Reading/WRF-SUEWS] Urban haze: changing SUEWS to capture urban haze (#54)

better to be resolved at the SUEWS side.

— You are receiving this because you are subscribed to this thread. Reply to this email directly, view it on GitHubhttps://github.com/Urban-Meteorology-Reading/WRF-SUEWS/issues/54?email_source=notifications&email_token=AE2KZXXVY4G4Z7NAWNOGBLLQCQHENA5CNFSM4II3CKBKYY3PNVWWK3TUL52HS4DFVREXG43VMVBW63LNMVXHJKTDN5WW2ZLOORPWSZGOD3NNNNA#issuecomment-517658292, or mute the threadhttps://github.com/notifications/unsubscribe-auth/AE2KZXUT3OS4PXDWTRVLFCLQCQHENANCNFSM4II3CKBA.

sunt05 commented 5 years ago

It turned out WRF overestimates Kdown in our London case as well; sorry I overlooked this part and thought it had been resolved.

Given the results about AOD and Kdown in Shanghai and Beijing, I suggest for our WRF-SUEWS evaluations we patch SuMin specifically for WRF-SUEWS instead of SUEWS for the following concerns:

  1. we need to apply thresholds of transmissivity at the grid/site scale;
  2. transmissivity has apparent seasonality (i.e., intra-annual variability, see below from Ao et al (2016)), suggesting specific thresholds should be applied for a short period (e.g., ~1 month);

Screenshot 2019-08-09 at 20 59 17

  1. SuMin is designed for WRF-SUEWS simulations of shorter simulation periods (e.g. ~10 days) compared with standard alone SUEWS, which is usually run for years.

As such, we can apply the thresholds in SuMin to correct the Kdown going into SUEWS kernel while keep current SUEWS as is.

However, for the long run, we may consider to introduce a chemical forcing file (e.g., hourly AOD measurements) that allows dynamic correction/adjustment of Kdown in a more responsive way.

Please let me know your thoughts.

If the SuMin approach sounds sensible, I'll start the implementation.

suegrimmond commented 5 years ago

I think we should first compare the transmissivity for the model and observations for the study periods. I.e. use the data we already have

sunt05 commented 5 years ago

results of transmissivity for London during four simulation periods

month transmissivity diff=sim-obs
Jan jan-tau jan-dif
Apr apr-tau apr-dif
Jul jul-tau jul-dif
Oct oct-tau oct-dif
suegrimmond commented 5 years ago

If we use a Kdn threshold - what is needed to remove the switch in the morning hours? (ie. from negative to positive difference)

Please could you now plot the differences normalised by time after SR (i.e. SR to SS=1) for the four periods. Please plot the four different periods (with the shading ) with one normalised time access

sunt05 commented 5 years ago

Hi @suegrimmond, I'm exploring different options for Kdown threshold but the results are not ideal. I guess your second part is to generalise diff_tau as a function of t_SR, which I'm working on now.

suegrimmond commented 5 years ago

It may be could be Zenith angle

Best wishes Sue

Prof Sue Grimmond Dept. of Meteorology, University of Reading, Reading, RG6 6BB T: 44 118 378 6248 – messages get emailed to me O:Met Building (#58 on map) rm:1U14 E: c.s.grimmond@reading.ac.ukmailto:c.s.grimmond@reading.ac.uk W: http://micromet.reading.ac.uk/

From: Ting Sun notifications@github.com Sent: 13 August 2019 15:51 To: Urban-Meteorology-Reading/WRF-SUEWS WRF-SUEWS@noreply.github.com Cc: Sue Grimmond c.s.grimmond@reading.ac.uk; Mention mention@noreply.github.com Subject: Re: [Urban-Meteorology-Reading/WRF-SUEWS] Urban haze: changing SUEWS to capture urban haze (#54)

Hi @suegrimmondhttps://github.com/suegrimmond, I'm exploring different options for Kdown threshold but the results are not ideal. I guess your second part is to generalise diff_tau as a function of t_SR, which I'm working on now.

— You are receiving this because you were mentioned. Reply to this email directly, view it on GitHubhttps://github.com/Urban-Meteorology-Reading/WRF-SUEWS/issues/54?email_source=notifications&email_token=AE2KZXR3DHYYUERCZMH2W5LQELC5LA5CNFSM4II3CKBKYY3PNVWWK3TUL52HS4DFVREXG43VMVBW63LNMVXHJKTDN5WW2ZLOORPWSZGOD4F5CZY#issuecomment-520868199, or mute the threadhttps://github.com/notifications/unsubscribe-auth/AE2KZXQ46S7TCNLW4YVQNRDQELC5LANCNFSM4II3CKBA.

sunt05 commented 5 years ago

That's a good point! re: zenith angle. My quick tentative conclusion for this is that a single cutoff might not be necessary to correct the modelled Kdown. I'll keep on posting new results that are sensible for further discussions.

suegrimmond commented 5 years ago

Maybe we can use a zenith angle as the starting time of the correction.

sunt05 commented 5 years ago

The overestimated Kdown by WRF seems to be a well-known issue that has been discussed in many papers, some of which are the following: https://dx.doi.org/10.1002/2016jd025527 https://journals.ametsoc.org/doi/10.1175/MWR-D-15-0262.1

sunt05 commented 5 years ago

London results

Note:

  1. t^~ is the time since sunrise normalised by length of day.
  2. the error bars show the confidence level (95%) of mean bootstrapped values.

Transmissivity

Untitled-tau

dif = sim - obs

Untitled

sunt05 commented 5 years ago

Swindon results

Note:

  1. t^~ is the time since sunrise normalised by length of day.
  2. the error bars show the confidence level (95%) of mean bootstrapped values.

Transmissivity

download

dif = sim - obs

download

suegrimmond commented 5 years ago

Useful to see that this is consistent with London

Best wishes Sue

Prof Sue Grimmond Meteorology, University of Reading


From: Ting Sun notifications@github.com Sent: Monday, August 19, 2019 4:35:18 AM To: Urban-Meteorology-Reading/WRF-SUEWS WRF-SUEWS@noreply.github.com Cc: Sue Grimmond c.s.grimmond@reading.ac.uk; Mention mention@noreply.github.com Subject: Re: [Urban-Meteorology-Reading/WRF-SUEWS] Urban haze: changing SUEWS to capture urban haze (#54)

Swindon results

Note:

  1. t^~ is the time since sunrise normalised by length of day.
  2. the error bars show the confidence level (95%) of mean bootstrapped values.

Transmissivity

[download]https://user-images.githubusercontent.com/1802656/63230075-e851b180-c1ff-11e9-9a9b-df515bb8c380.png

dif = sim - obs

[download]https://user-images.githubusercontent.com/1802656/63230083-f43d7380-c1ff-11e9-9e91-1c23a9ea8304.png

— You are receiving this because you were mentioned. Reply to this email directly, view it on GitHubhttps://github.com/Urban-Meteorology-Reading/WRF-SUEWS/issues/54?email_source=notifications&email_token=AE2KZXUT7YZ555HP3VKIBS3QFGXANA5CNFSM4II3CKBKYY3PNVWWK3TUL52HS4DFVREXG43VMVBW63LNMVXHJKTDN5WW2ZLOORPWSZGOD4RHUVY#issuecomment-522353239, or mute the threadhttps://github.com/notifications/unsubscribe-auth/AE2KZXRK6QXTRK4HLTUDGWLQFGXANANCNFSM4II3CKBA.

hamidrezaomidvar commented 5 years ago

I have started implementing this. However, I think we need a more sophisticated method to modify kdown rather than a threshold (for example after midday) because it leads a jump in kdown when we kick in the modification.

hamidrezaomidvar commented 5 years ago

I see this in a simple implementation I did

hamidrezaomidvar commented 5 years ago

I applied the transmissivity difference over the entire day to make sure it has a smooth profile. I did a run over April period. The run is not fully finished yet (11 to 22 is finished), but here are the results. The first two plots are original and modified Kdown.

radiation-diurnal-London

hamidrezaomidvar commented 5 years ago

The modification is done in the interface of WRF-SUEWS, so the original SUEWS code is not changed!