城市冠层模型UCM与城市土地利用的修改

Kecheng

大家好，我在使用wrf-ucm的过程中遇到了这样的问题：

我使用的是modis21类土地资料，运行geogrid完成后，我使用我写的python代码根据人口密度将geo_em里LU_INDEX中的城市类别13替换成了31低密度居民区、32高密度居民区和33商业区，并且运行metgrid之后生成的met_em文件里LU_INDEX也是31、32、33。

可是我发现在real完成之后，发现wrfinput里的LU_INDEX又变成了13，其中我在namelist里的设置num_land_cat =21，并且运行wrf的过程中我发现有这样一句话：
USING DEFAULT URBAN MORPHOLOGY
我不确定它是否使用了我设置的URBAPARM.TBL。

如果我将wrfinput里的LU_INDEX手动把13改成31、32、33，程序也会报错无法运行。

因此我有几个问题想要寻求帮助：

第一，是否我这样修改土地利用的做法是错误的？应该怎么把城市类别替换成31、32、33呢？是否需要修改LANDUSE.TBL等文件？
第二，URBAPARM.TBL的开头是这样：
# The parameters in this table may vary greatly from city to city.
# The default values are probably not appropriate for any given city.
# Users should adapt these values based on the city they are working
# with.

# Urban Parameters depending on Urban type
# USGS

Number of urban categories: 3

#
#  Where there are multiple columns of values, the values refer, in
#  order, to: 1) Low density residential, 2) High density residential,
#  and 3) Commercial:  I.e.:
#
#  Index:     1           2              3
#  Type:  Low-dens Res, Hi-dens Res, Commercial
#
我该怎么让31、32、33的地区采用对应的URBAPARM.TBL的设置呢？还是它会默认识别31、32和33？

谢谢大家的帮助！

tbag

是这样，问下为什么你要改到31, 32, 33?

如果我没记错的话，默认MODIS或者USGS的urban land use是只有一种: high-density residential
来源在此: https://forum.mmm.ucar.edu/threa ... -urban-class.20451/

31-33的landuse 是来自于NLCD数据，这个NLCD数据只包含了北美
你可以看看这个链接: https://ral.ucar.edu/sites/defau ... l/WRF-LSM-Urban.pdf

tbag

能用LCZ就用LCZ吧，LCZ数据使用起来也简单。

Kecheng

tbag 发表于 2025-7-26 10:54
是这样，问下为什么你要改到31, 32, 33?

如果我没记错的话，默认MODIS或者USGS的urban land use是只有一 ...

我想着改到31、32、33，再使用对应的URBAPARM.TBL里的数值，模拟得应该会更准确一些？

Kecheng

tbag 发表于 2025-7-26 10:57
能用LCZ就用LCZ吧，LCZ数据使用起来也简单。

好的，谢谢您的建议，我去了解一下LCZ的用法。方便的话能请问一下您有什么推荐的资料吗？

tbag

Kecheng 发表于 2025-7-26 18:10
我想着改到31、32、33，再使用对应的URBAPARM.TBL里的数值，模拟得应该会更准确一些？

按道理来说没错
主要是看你模拟的范围，如果你模拟范围在北美，那用NLCD没啥问题。

LCZ可以研究下这里: https://ral.ucar.edu/sites/defau ... enlinhe20230113.pdf

tbag

数据下载在这里: https://www2.mmm.ucar.edu/wrf/us ... urces_wps_geog.html

Kecheng

tbag 发表于 2025-7-26 22:35
按道理来说没错
主要是看你模拟的范围，如果你模拟范围在北美，那用NLCD没啥问题。

我模拟的区域主要是香港，因为想和之前没开城市冠层的结果进行比对，所以想继续用modis数据，然后把13换成31、32、33再运行，但是就遇到了上面的问题就很奇怪。
LCZ我昨天试了一下确实会方便很多，不过感觉默认的URBPARM_LCZ.TBL还需要修改，想请问一下您知道目前有已经做好的（类似我另外一个帖子里的那些URBAPAR.TBL）的URBPARM_LCZ.TBL吗？

pigzero527

LANDUSE.TBL里要进行相应的更改，比如找到modis的类型，然后给对应的31、32、33分配上low,middle,high
31,    10.,   .10,   .97,   80.,    3.,  1.67, 18.9e5, 'Low Intensity Residential'
32,    10.,   .10,   .97,   80.,    3.,  1.67, 18.9e5, 'High Intensity Residential'
33,    10.,   .10,   .97,   80.,    3.,  1.67, 18.9e5, 'Industrial or Commercial'

另外VEGPARM.TBL里也要加上31-33的类型
LOW_DENSITY_RESIDENTIAL
31
HIGH_DENSITY_RESIDENTIAL
32
HIGH_INTENSITY_INDUSTRIAL
33

Kecheng

pigzero527 发表于 2025-7-27 22:34
LANDUSE.TBL里要进行相应的更改，比如找到modis的类型，然后给对应的31、32、33分配上low,middle,high
31, ...

非常感谢！我明天去尝试一下！