Urban-Meteorology-Reading / WRF-SUEWS

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

Paper in SUEWS parameters for vegetated surfaces #66

Closed hamidrezaomidvar closed 4 years ago

hamidrezaomidvar commented 4 years ago

Here the rough plan for starting the paper to investigate the SUEWS parameters for vegetated surfaces in order to reference it in the WRF-SUEWS paper:

Important parameters to focus on:

At the end, we like to purpose a table containing the suitable parameters for each land cover for users to use it when running SUEWS/WRF-SUEWS.

hamidrezaomidvar commented 4 years ago

Albedo for DecTr:

I tested SUEWS model in calculating albedo of DecTr land cover for 3 sites each for 4 years. Below you can see the results of these sites. During the winter, because of the snow, there are some outliers as you can see for two north sites and for the souther site (MMS), the outliers are less during the winter.

I also plotted the temperature and BaseT and BaseTe as well as calculated LAI for all the years. I think the plots give a good idea about how SUEWS is doing for albedo of DecTr and what do we need to change for better albedos I will work on EveTr and grass.

US-UMB ![US-UMB-albedo](https://user-images.githubusercontent.com/44125994/65048728-7bdfe480-d95c-11e9-91ec-9b5d29f5de3a.png)
US-MMS ![US-MMS-albedo](https://user-images.githubusercontent.com/44125994/65048757-87331000-d95c-11e9-8c17-6bf916ce071b.png)
US-Oho ![US-Oho-albedo](https://user-images.githubusercontent.com/44125994/65048774-8f8b4b00-d95c-11e9-99f6-e6c2c1cea581.png)
hamidrezaomidvar commented 4 years ago

Here are similar results for EveTr land cover. You can see a different behaviour in Albedo, and lot's of snow effects for CA-Qcu during the winter (because of the location of the site)

CA-Obs ![CA-Obs-albedo](https://user-images.githubusercontent.com/44125994/65054973-5657d880-d966-11e9-8927-c854be8eaaa0.png)
US-Blk ![US-Blk-albedo](https://user-images.githubusercontent.com/44125994/65055003-62dc3100-d966-11e9-9fb6-a6d57c61beb1.png)
CA-Qcu ![CA-Qcu-albedo](https://user-images.githubusercontent.com/44125994/65055019-6a033f00-d966-11e9-8644-15ef093b0d38.png)
hamidrezaomidvar commented 4 years ago

And here is the Grass land cover for two sites:

US-Var ![US-Var-albedo](https://user-images.githubusercontent.com/44125994/65059829-9622be00-d96e-11e9-9ad0-e40483d93a3f.png)
US-AR1 ![US-AR1-albedo](https://user-images.githubusercontent.com/44125994/65059838-99b64500-d96e-11e9-898c-5eefe303dd8f.png)
suegrimmond commented 4 years ago

These are all very good – it would be good to know how the LAI looks relative to MODIS for these sites (others vegetation types as well.) Have you talked with Mathew – to make certain you both look at the same first sites? So we can determine the next steps

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: Hamidreza Omidvar notifications@github.com Sent: 17 September 2019 17:15 To: Urban-Meteorology-Reading/WRF-SUEWS WRF-SUEWS@noreply.github.com Cc: Subscribed subscribed@noreply.github.com Subject: Re: [Urban-Meteorology-Reading/WRF-SUEWS] Paper in SUEWS parameters for vegetated surfaces (#66)

And here is the Grass land cover for two sites: US-Var

[US-Var-albedo]https://user-images.githubusercontent.com/44125994/65059829-9622be00-d96e-11e9-9ad0-e40483d93a3f.png US-AR1

[US-AR1-albedo]https://user-images.githubusercontent.com/44125994/65059838-99b64500-d96e-11e9-898c-5eefe303dd8f.png

— 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/66?email_source=notifications&email_token=AE2KZXWZOYJR2JLJFMNMUE3QKD7AVA5CNFSM4IXMHC7KYY3PNVWWK3TUL52HS4DFVREXG43VMVBW63LNMVXHJKTDN5WW2ZLOORPWSZGOD65CFBY#issuecomment-532292231, or mute the threadhttps://github.com/notifications/unsubscribe-auth/AE2KZXSP77ZP3VOE7A427VDQKD7AVANCNFSM4IXMHC7A.

suegrimmond commented 4 years ago

It would be good to mark on the GDD and SDD points as well as the BaseTs so we can see what are likely to be more appropriate values (as I think I can see how they will probably need to be changed (But should not start that until the MODIS data are being looked at) ). We can skype if you want.

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: Hamidreza Omidvar notifications@github.com Sent: 17 September 2019 17:15 To: Urban-Meteorology-Reading/WRF-SUEWS WRF-SUEWS@noreply.github.com Cc: Subscribed subscribed@noreply.github.com Subject: Re: [Urban-Meteorology-Reading/WRF-SUEWS] Paper in SUEWS parameters for vegetated surfaces (#66)

And here is the Grass land cover for two sites: US-Var

[US-Var-albedo]https://user-images.githubusercontent.com/44125994/65059829-9622be00-d96e-11e9-9ad0-e40483d93a3f.png US-AR1

[US-AR1-albedo]https://user-images.githubusercontent.com/44125994/65059838-99b64500-d96e-11e9-898c-5eefe303dd8f.png

— 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/66?email_source=notifications&email_token=AE2KZXWZOYJR2JLJFMNMUE3QKD7AVA5CNFSM4IXMHC7KYY3PNVWWK3TUL52HS4DFVREXG43VMVBW63LNMVXHJKTDN5WW2ZLOORPWSZGOD65CFBY#issuecomment-532292231, or mute the threadhttps://github.com/notifications/unsubscribe-auth/AE2KZXSP77ZP3VOE7A427VDQKD7AVANCNFSM4IXMHC7A.

hamidrezaomidvar commented 4 years ago

Ok, I will talk to Mathew to work on these sites.

hamidrezaomidvar commented 4 years ago

The following parameters control the LAI in SUEWS:

From these parameters, I found that BaseTe (because of Climate) and LaiMax (For same vegetation might be still different across sites) are site dependent (for same type of vegetation). So to calculate these, first I calculated them based on MMS-2010. Then for each other sites, I only calculated BaseTe and LaiMax and tested them for 3 other years of the same site. So basically, BaseT,GDDFUll,SDDFUll, and LAIMin are the same in all sites and tests (based on MMS-2010), but BaseTe and LaiMax varies for different sites (calculated based on 1 year of each site). Here are the results of tests for DecTr sites:

Parameters:

US-MMS ![US-MMS-LAI](https://user-images.githubusercontent.com/44125994/65268426-09633600-db0f-11e9-87b9-da2ce65361b0.png)
US-UMB ![US-UMB-LAI](https://user-images.githubusercontent.com/44125994/65268457-17b15200-db0f-11e9-96a8-cbc549aaa790.png)
US-Oho ![US-Oho-LAI](https://user-images.githubusercontent.com/44125994/65268464-1f70f680-db0f-11e9-871f-7926f4d6e57a.png)
sunt05 commented 4 years ago

Good job! @hamidrezaomidvar

hamidrezaomidvar commented 4 years ago

Here are the results for EveTr. It seems that GDDFull and SDDFull of previous sites are still working for these sites. So to be consistent, I did not change them. For one of the sites (CA-Qcu) where it is very cold during the winter, it seems it needed a lower BaseT otherwise the vegetation would not grow much during the growing season. Parameters:

US-Blk ![US-Blk-LAI](https://user-images.githubusercontent.com/44125994/65314308-ff344c80-db8d-11e9-9c6d-6f5582b8cf54.png)
CA-Obs ![CA-Obs-LAI](https://user-images.githubusercontent.com/44125994/65314316-02c7d380-db8e-11e9-9d22-1b5363966d9a.png)
CA-Qcu ![CA-Qcu-LAI](https://user-images.githubusercontent.com/44125994/65314326-078c8780-db8e-11e9-9267-b46758e9b872.png)
hamidrezaomidvar commented 4 years ago

for the Grass landcover, it seems a more variability in LAI even for one specific site which I think it might be related to mowing etc. For example, for US-AR1, I tuned on year 2010:

Here ![US-AR1-2010-LAI](https://user-images.githubusercontent.com/44125994/65316707-5b00d480-db92-11e9-9659-c378e529644e.png)

and here are the results for other two years which Laimax is different:

Here ![US-AR1-LAI](https://user-images.githubusercontent.com/44125994/65316735-681dc380-db92-11e9-86ba-15d1914efc44.png)

or for example for US-Var (California) , it seems the behaviour of LAI is completely different and the pick is not in the growing season?!

Here ![US-Var-2004-LAI](https://user-images.githubusercontent.com/44125994/65316920-ba5ee480-db92-11e9-8162-6c1753fa1731.png) ![US-Var-LAI](https://user-images.githubusercontent.com/44125994/65317084-03af3400-db93-11e9-82ee-4ed959f079bc.png)
hamidrezaomidvar commented 4 years ago

for US-Var, increasing the BaseTe as it is warmer, and decreasing GDDFull as it growing starts sooner would help to have a better representation of LAI:

Here ![US-Var-LAI](https://user-images.githubusercontent.com/44125994/65317787-5f2df180-db94-11e9-8423-16d1df6829b9.png)
hamidrezaomidvar commented 4 years ago

After getting proper parameters for LAI for different surfaces, I used them to tune (max and min of albedo for different sites) and test albedo for different vegetated surfaces. Here are the results

DecTr

US-MMS ![US-MMS-albedo](https://user-images.githubusercontent.com/44125994/65528594-95e56e00-deec-11e9-8d5f-f73c53882bb9.png)
US-UMB ![US-UMB-albedo](https://user-images.githubusercontent.com/44125994/65528612-9f6ed600-deec-11e9-9f99-1b016300d266.png)
US-Oho ![US-Oho-albedo](https://user-images.githubusercontent.com/44125994/65528645-abf32e80-deec-11e9-8522-e07b36294e21.png)

EveTr

CA-Obd ![CA-Obs-albedo](https://user-images.githubusercontent.com/44125994/65528691-c2998580-deec-11e9-86fc-263b499af9db.png)
US-Blk ![US-Blk-albedo](https://user-images.githubusercontent.com/44125994/65528733-d513bf00-deec-11e9-8247-d1c8d667f26c.png)
CA-Qcu ![CA-Qcu-albedo](https://user-images.githubusercontent.com/44125994/65528752-df35bd80-deec-11e9-810b-91ad57469027.png)

Grass

US-AR1 ![US-AR1-albedo](https://user-images.githubusercontent.com/44125994/65528795-f2e12400-deec-11e9-8170-574892e90603.png)
US-KUT ![US-KUT-albedo](https://user-images.githubusercontent.com/44125994/65528810-fb395f00-deec-11e9-88ac-9d65d7d9b389.png)
hamidrezaomidvar commented 4 years ago

Finally, I used the parameters from LAI and Albedo of previous sections to fit g1-g6 for different sites and tested them. Here are the results of the tests:

DecTr g1=1,g2=266.83,g3=0.135,g4=0.429,g5=35.2,g6=0.020

Test:

US-MMS ![diurnal_test-US-MMS](https://user-images.githubusercontent.com/44125994/65529219-a5b18200-deed-11e9-83f6-1011e929f5b3.png) ![compare_test-US-MMS](https://user-images.githubusercontent.com/44125994/65529229-a9450900-deed-11e9-8ef4-ffe62c3f0978.png)
US-UMB ![compare_test-US-UMB](https://user-images.githubusercontent.com/44125994/65529247-b3ff9e00-deed-11e9-875e-2acc24ea2980.png) ![diurnal_test-US-UMB](https://user-images.githubusercontent.com/44125994/65529248-b3ff9e00-deed-11e9-9d03-72104b25e7c2.png)
US-Oho ![compare_test-US-Oho](https://user-images.githubusercontent.com/44125994/65529267-bc57d900-deed-11e9-9248-6d0de26b6092.png) ![diurnal_test-US-Oho](https://user-images.githubusercontent.com/44125994/65529269-bc57d900-deed-11e9-9453-0690393172ea.png)

EveTr g1=1,g2=266.95,g3=0.360,g4=0.757,g5=33.81,g6=0.024

Test:

CA-Obs ![compare_test-CA-Obs](https://user-images.githubusercontent.com/44125994/65529404-ef01d180-deed-11e9-818d-202e923c19f6.png) ![diurnal_test-CA-Obs](https://user-images.githubusercontent.com/44125994/65529406-ef9a6800-deed-11e9-9aa5-b8be50365239.png)
US-Blk ![compare_test-US-Blk](https://user-images.githubusercontent.com/44125994/65529427-fc1ec080-deed-11e9-82b4-7a5541ac2c68.png) ![diurnal_test-US-Blk](https://user-images.githubusercontent.com/44125994/65529428-fc1ec080-deed-11e9-8ee2-835833ac2c98.png)
CA-Qcu ![compare_test-CA-Qcu](https://user-images.githubusercontent.com/44125994/65529453-0c36a000-deee-11e9-8a2a-0131690bf9ae.png) ![diurnal_test-CA-Qcu](https://user-images.githubusercontent.com/44125994/65529454-0ccf3680-deee-11e9-8cd8-045f6424e2f4.png)

Grass g1=1,g2=266.92,g3=0.350,g4=0.763,g5=33.9,g6=0.021

Test:

US-KUT ![compare_test-US-KUT-2007](https://user-images.githubusercontent.com/44125994/65529554-3daf6b80-deee-11e9-9255-e6715e6a5d4a.png) ![diurnal_test-US-KUT-2007](https://user-images.githubusercontent.com/44125994/65529555-3daf6b80-deee-11e9-9185-fa3bbf1a0b90.png)
US-AR1 ![compare_test-US-AR1](https://user-images.githubusercontent.com/44125994/65529575-486a0080-deee-11e9-97be-d2d72f693c70.png) ![diurnal_test-US-AR1](https://user-images.githubusercontent.com/44125994/65529576-486a0080-deee-11e9-90ae-5bc9964827ad.png)
hamidrezaomidvar commented 4 years ago

In summary, we need the following processes for each vegetated surface in order to get suitable SUEWS parameters:

SUEWS_Parameters