Updates to radiation, cloud, and land processes for improved surface temperature and radiative flux biases

[rhaesung]

Changes from @RuiyuSun and @HelinWei-NOAA.

This PR addresses several key issues identified in the HRv4 and UFS/Noah-MP models, focusing on radiative flux biases, soil thermal conductivity, and evapotranspiration improvements. The goal of these updates is to enhance model performance and better align simulations with observed data.

1. Cloud/Radiation adjustments:

• To reduce the high bias in downward shortwave (DSW) radiative flux over land, the cloud number concentration over land has been increased from 100/cm³ to 150/cm³.
• The convective cloud ice radius in radiation has been reduced from 50 µm to 25 µm.
• The 2020 summer test results (HRv4Ntc100) show a reduction in the high DSW bias over land and a corresponding reduction in the low bias of upward shortwave (SW) radiative flux at the top of the atmosphere (TOA).

2. Soil thermal conductivity enhancement:

• A persistent nighttime and early morning warm surface temperature bias remains in both the GFS/Noah and UFS/Noah-MP models. This is partly due to overlooking the influence of soil organic matter (SOM) on soil thermal conductivity.
• By incorporating a lookup table developed from the PNNL soil carbon dataset and VIIRS surface type data, the effects of SOM on soil thermal properties have been accounted for.
• Tests show that including SOM results in better simulation of ground heat flux, reducing heat storage in the soil during the day, and mitigating the nighttime warm bias, particularly over the Great Plains.

3. Evapotranspiration adjustment:

• The aquifer hydraulic conductivity has been updated to use the algorithm average rather than the harmonic average. This increases the amount of water uptake by the soil layer, improving evapotranspiration modeling.

4. Bug Fixes:

• Snow aging at night: Fixed an issue to allow snow aging to occur at night, ensuring more realistic snow dynamics.
• Unit correction for potential evaporation: Corrected the unit for potential evaporation from W/m² to mm/s to ensure proper unit consistency and accuracy in the model. See this change in Updates to radiation, cloud, and land processes for improved surface temperature and radiative flux biases NOAA-EMC/ufsatm#903.
• Address Bug fix in physics/module_sf_noahmp_glacier.F90: only compute glacier albedo in snowalb_bats_glacier if cosz > 0 NCAR/ccpp-physics#1091

[RuiyuSun]

@rhaesung please use radiation_clouds.f at /lfs/h2/emc/physics/noscrub/ruiyu.sun/git/global_workflow.git.hrv4/sorc/ufs_model.fd/FV3/ccpp/physics/physics/Radiation to replace the one you are using in this PR. This file uses 25 um as the convective cloud ice radius over land in radiation instead of 50 µm.

[HelinWei-NOAA]

@rhaesung please use radiation_clouds.f at /lfs/h2/emc/physics/noscrub/ruiyu.sun/git/global_workflow.git.hrv4/sorc/ufs_model.fd/FV3/ccpp/physics/physics/Radiation to replace the one you are using in this PR. This file uses 25 um as the convective cloud ice radius over land in radiation instead of 50 µm.

@RuiyuSun A little confused here? The original value is 50 and you proposed to change it to 25, right?

[RuiyuSun]

@rhaesung please use radiation_clouds.f at /lfs/h2/emc/physics/noscrub/ruiyu.sun/git/global_workflow.git.hrv4/sorc/ufs_model.fd/FV3/ccpp/physics/physics/Radiation to replace the one you are using in this PR. This file uses 25 um as the convective cloud ice radius over land in radiation instead of 50 µm.

@RuiyuSun A little confused here? The original value is 50 and you proposed to change it to 25, right?

@HelinWei-NOAA Right. But it is only changed over land.

[rhaesung]

@RuiyuSun radiation_clouds.f has been updated.

[climbfuji]

Reviewed the module_mp_thompson and radiation_clouds parts.

[grantfirl]

@rhaesung This looks good from a CCPP point of view. There will need to be a similar change for the NoahMP component model code in order to stay in sync.

[rhaesung]

@rhaesung This looks good from a CCPP point of view. There will need to be a similar change for the NoahMP component model code in order to stay in sync.

@barlage @uturuncoglu @HelinWei-NOAA FYI

[barlage]

@HelinWei-NOAA @rhaesung @climbfuji I will add a review here shortly. In addition to the albedo changes in the non-glacier land model, similar changes should also be made to address NCAR#1091 with this PR.

[rhaesung]

@barlage The files module_sf_noahmplsm.F90 and module_sf_noahmp_glacier.F90 have been updated according to your comments. FYI, @HelinWei-NOAA and @cenlinhe.

[rhaesung]

@barlage Thanks! The code has been updated with your changes.

[jkbk2004]

Tests are done ok at ufs-community/ufs-weather-model#2548. @grantfirl @rhaesung @dustinswales Can you merge this pr?