Modeling Land‐Atmosphere Coupling at Cloud‐Resolving Scale Within the Multiple Atmosphere Multiple Land (MAML) Framework in SP‐E3SM

Representing subgrid variabilities of land surface processes and their upscaled effects is crucial for global climate modeling. Here, we implement a multiple atmosphere multiple land (MAML) framework in the superparamaterized version of E3SM (SP‐E3SM) to explicitly simulate the subgrid variabilities of land states and fluxes at cloud‐resolving scale and their interactions with atmosphere. Comparing to the standard SP‐E3SM in which all the atmospheric columns of the cloud resolving model embedded within the global atmospheric model grid interact with the same land surface (i.e., multiple atmosphere single land (MASL)), the impact of MAML on the strength of land‐atmosphere coupling is limited, partly because the current implementation mainly facilitates one‐way coupling between the cloud‐resolving model and the land surface model. Despite such limitation, MAML increases the surface latent heat flux at the expense of sensible heat flux, and increases precipitation in India, Amazon, and Central Africa, reducing the model dry bias compared to the standard SP‐E3SM. By employing a normalized gross moist stability (NGMS) diagnostic framework, we find that the increase in precipitation minus evaporation (P‐E) is primarily driven by the change in large‐scale moisture convergence, particularly by the increase of water vapor in the lower atmosphere, while the local effect of total surface energy flux plays a minor role in the P‐E change. More specifically, MAML changes the surface energy partitioning (evaporative fraction), increases the atmosphere water vapor, and further increases P‐E by decreasing the NGMS. Finally, future development in the MAML framework is discussed. Plain Language Summary Land‐atmosphere interactions such as soil moisture‐precipitation feedback occur at a wide range of spatial scales. For example, spatial variability of surface fluxes such as radiation and precipitation can influence surface latent and sensible heat fluxes, which further influence turbulent mixing processes in the boundary layer and cloud and precipitation. Current global climate models (GCMs) have difficulties in representing land‐atmosphere interactions at scales smaller than the GCM grid (subgrid scales). To meet this challenge, we employ a novel framework, called multiple atmosphere multiple land (MAML), in a superparamaterized version of E3SM (SP‐E3SM) in which a cloud resolving model is embedded within each GCM grid to better resolve clouds and convection. MAML explicitly