Distinct Impacts of Land Use and Land Management on Summer Temperatures

Land use has been recognized as an important anthropogenic forcing of climate change in recent studies. However, climatic effects of land management practices have been little discussed and compared to land-use impacts. As land-atmosphere interactions via surface fluxes are particularly strong during the warm season, we investigate the impacts of historical land use and present irrigation practices on summer temperatures in the Northern Hemisphere using the most recent version of Community Earth System Model. Our results suggest that historical land use leads to an overall cooling in the middle latitudes and a warming in the tropics, and the sign and magnitude of the changes in temperature depend on the type of land cover change. On the other hand, summer irrigation leads to a significant cooling over many irrigated areas due to enhanced evapotranspiration, and the local cooling is comparable to and even stronger than the land-use effects. Both land use and irrigation can also significantly influence the intensity and frequency of hot extremes. Land use shows stronger impacts during the night through ground heat flux feedback, while irrigation shows stronger impacts during the day through latent heat flux feedback. Our comparison demonstrates the importance of irrigation in local and regional climate, highlighting the necessity of considering such land management practices in future assessments of regional climate change and climate mitigation.


INTRODUCTION
It is widely recognized that land use/land cover change (hereafter referred to simply as land use) affects the overlying atmosphere through land-atmosphere interactions, and thus modifies the local and broader-scale climate. The importance of land use in the climate system promoted the Coupled Model Intercomparison Project Phase 5 (CMIP5) to include land use forcing in its climate projections (Hurtt et al., 2011). Based on the CMIP5 simulations, many studies have been carried out to investigate the impacts of land use on climate (Brovkin et al., 2013;Kumar et al., 2013;Di Vittorio et al., 2014) and extreme events (Lejeune et al., 2017. It is found that land use can have impacts comparable to increased greenhouse gasses or sea surface temperature variations for many climate variables de Noblet-Ducoudré et al., 2012), and land-use forcing can be as important as other anthropogenic forcings in explaining the historical changes in temperature extremes over the regions with extensive land use .
Besides land use, which usually refers to conversions from one land cover to another due to human activities, land management practices and their climatic impacts have drawn increasing attention recently (Luyssaert et al., 2014;Mahmood et al., 2014). As an important agricultural practice to maintain adequate soil moisture for stable crop production, irrigation not only modifies the surface water budget, but also affects the energy balance of the land surface, thereby altering the climate (Kueppers et al., 2007). For instance, irrigation can change surface partitioning of available energy between sensible and latent heat fluxes by allowing more evapotranspiration through increased soil moisture and larger vegetation coverage, potentially lowering near-surface temperatures over the irrigated areas.
Recent modeling studies have demonstrated the significant impacts of irrigation on near-surface climate and atmospheric circulation (e.g., Boucher et al., 2004;Kueppers et al., 2007;Lobell et al., 2009;Sacks et al., 2009;Lo and Famiglietti, 2013;Wei et al., 2013;Lu and Kueppers, 2015;de Vrese et al., 2016;Huang and Ullrich, 2016;Krakauer et al., 2016), especially the strong cooling effects on daytime temperatures or hot extremes. It is proposed that irrigation can be an effective way to mitigate the regional warming of hot extremes Thiery et al., 2017). Studies based on satellite observations also indicate the substantial local cooling of irrigation (Ambika and Mishra, 2019;Chen and Dirmeyer, 2019b), and imply that the cooling effects may be underestimated by global climate models due to their coarse spatial resolutions (Sorooshian et al., 2011;Chen and Dirmeyer, 2019b).
Despite its fundamental importance in altering the regional climate, the full scope of land management has rarely been considered in recent climate assessment projects. For instance, CMIP5 did not include irrigation as a historical forcing (Kumar et al., 2013), and many of the participating Earth system models did not implement irrigation schemes in their land surface models (Singh et al., 2018). The project Land-Use and Climate, Identification of Robust Impacts (LUCID), for example, is mainly focused on the biogeophysical impacts of historical land use . There are a few studies demonstrating the siginificance of irrigation compared to other land surface changes. For instance, the effects of irrigation on water vapor flows are equally as important as deforestation (Gordon et al., 2005); the cooling of irrigation is most pronounced compared to increased leaf area index and reduced tillage (Lobell et al., 2006). However, few studies have investigated the impacts of irrigation and its relative importance compared to historical land use in modifying local and regional temperatures.
With the goal of comparing these two types of land management, this study uses a state-of-the-art Earth system model, Community Earth System Model version 2 (CESM2), to investigate the separate impacts of land use and irrigation on summer temperatures and extreme heat events. We focus on boreal summer temperature because of the strong landatmosphere interactions and high irrigation demands during the warm season, when these two land surface forcings should have the most profound climatic impacts and the greatest social implications.

Land-Use Experiments in CESM2
The land-use experiments are conducted with CESM2, which is a coupled Earth System model composed of separate climate system components for atmosphere, ocean, land, sea ice and land ice. The major focus of this study is on the interactions between the land surface and atmosphere, so only the Community Atmosphere Model (CAM6) and Community Land Model (CLM5) components are used in our simulations (no interactive ocean or sea ice), with the component set F2000climo at a spatial resolution of 0.9 • × 1.25 • . The F2000climo component set allows climatology simulations with cyclic circa-year-2000 forcing, in which there are prescribed sea surface temperature (SST) and sea ice cover with a fixed CO 2 concentration of 367.0 ppm. Monthly mean climatology of SST and sea ice cover are derived from a merged product based on the monthly mean Hadley Center Sea Ice and SST dataset, version 1 (HadISST1), and version 2 of the NOAA weekly Optimum Interpolation SST (OISST2) analysis during the period 1995-2005 (Hurrell et al., 2008). Compared with the previous version of CLM, improvements are made in CLM5 to better represent soil and plant hydrology, snow density, carbon and nitrogen cycling and coupling, the representation of crops and human land management (Lawrence et al., 2019).
In CLM5, land surface heterogeneity is represented as a nested subgrid hierarchy, in which grid cells are composed of multiple land units (such as the vegetated unit and crop unit). The vegetated unit can be composed of different plant functional types (PFTs), and the crop unit can be composed of different crop functional types (CFTs). Land use is represented as the changes in percent PFTs and CFTs within a grid cell. In this study, two land-use experiments (pre-industrial and present) are carried out with prescribed pre-industrial and present land cover conditions ( Table 1). The difference between present and pre-industrial can be considered as the impacts of historical land use. A two-tailed Student's t-test is conducted to assess the significance of the landuse induced difference. Details of the pre-industrial land cover can be found in Hurtt et al. (2006) and Lawrence et al. (2012), while the present land cover is derived from MODIS satellite data as described in Lawrence and Chase (2007).

Irrigation in CLM5
In the CLM5, there are a total of 31 managed crop types . Each crop type has rainfed and irrigated CFTs, which are placed on separate soil columns. The proportion of irrigated CFTs is based on a dataset of areas equipped for irrigation (Portmann et al., 2010), which is shown in Figure 1D. When irrigation is enabled, the model checks if the crops need irrigation in the first time step after 6 AM local time every day. If the crop leaf area is greater than zero and the available soil water is below a specified threshold, irrigation will be triggered. Irrigated water is removed from river water storage and applied directly to the ground surface. Details of the irrigation scheme in CLM5 can be found in Lawrence et al. (2018). To separate the impacts of land use and land management, an irrigation experiment, irrig, is carried out with the present land cover conditions ( Table 1).
Irrigation is not automatically activated in the pre-industrial and present experiments, so the difference between irrig and present can be considered as the effects of irrigation. A two-tailed Student's t-test is used to assess the significance of the irrigation induced difference.
Because the prognostic biogeochemical model is activated to estimate the vegetation phenology, land initial conditions for each experiment are obtained from a separate 700−year offline spin−up run, then each experiment is conducted for 60 years with output saved at three-hourly intervals.

Indices of Hot Extremes
Previous studies suggest that land use and irrigation have significant impacts on hot extremes Pitman et al., 2012;Christidis et al., 2013;Hirsch et al., 2017;Thiery et al., 2017). Therefore, we assess the changes in four warm extreme temperature indices, listed and defined in Table 2, as recommended by the CCl/CLIVAR/JCOMM Expert Team on Climate Change Detection and Indices (ETCCDI; Alexander et al., 2006). These indices are based on daily maximum and minimum temperature and are developed to assess the intensity and frequency of extreme temperature events.
Additionally, to understand the mechanism of temperature responses to changes in land use and land management, other land surface variables (such as surface fluxes and radiation) are also included in our analysis. With the focus on daily minimum and maximum temperature, we identify the timing of daily maximum (minimum) land surface temperatures (T s ) based on three-hourly model output, and extract the corresponding surface fluxes and radiation when T s reaches the daily maximum (minimum).

Historical Land Use and Irrigation Activities
Figures 1A-C shows the historical land cover change from pre-industrial to present day over the Northern Hemisphere. Deforestation has mainly occurred in the Midwest of the US, eastern Europe, India, eastern China, and in tropical areas such as Central America, Africa, and southeastern Asia. Most of the deforestation is driven by agricultural demand through cropland expansion in the middle latitudes, which has also led to grassland loss in such regions as the Great Plains, Eastern Europe,  (Karl et al., 1999;Peterson et al., 2001).
and India. Meanwhile, the conversion of cropland back into forest or grassland occurs in some regions of the northeastern US and western Europe. Grassland degradation (conversion of grassland into bare ground) is found in some regions of the Tibetan Plateau, which might be associated with over-stocking (Li et al., 2013). Areas of the irrigated cropland are shown in Figure 1D. The irrigation systems are mainly distributed in northern India, eastern China, the Middle East, Southern Europe, the Central Valley in California, the North American Great Plains, and the Mississippi Embayment, according to Portmann et al. (2010). Based on the irrigation experiment, the total amount of irrigated water can be estimated in CLM5 (Figure 2A), which shows good agreement with observationally based estimates , indicating that CESM2 is an appropriate tool to investigate the impact of irrigation (Huang and Ullrich, 2016;Thiery et al., 2017). Seven subregions with high irrigation rates are selected for regional analysis. We also identified the month with the highest irrigation rate ( Figure 2B). Summertime peak irrigation occurs mostly in the middle latitudes, such as North America, the Mediterranean, and western/eastern Asia. In the tropics and regions that are influenced by monsoons (such as the East Asian Monsoon and North American Monsoon), the highest irrigation rate is mainly found in spring, which is the dry season prior to the onset monsoon rains. However, there is still a large amount of irrigated water applied during the summer in those regions, according to the irrigation algorithm in CLM5 (Supplementary Figure S1), exerting impacts on the surface climate then as well.

Changes in Temperature
Due to possible different responses of daytime temperature to deforestation at the land surface (T s ) and 2-m air (T 2m ) reported in our previous work (Chen and Dirmeyer, 2019a), we present the results of T s and T 2m separately (Figure 3). Historical deforestation leads to a significant daytime warming at the land surface (T s ) in the tropics, India, and the central US ( Figure 3A). However, the change in daily maximum T 2m largely depends on the type of land use. Cooling effects on T 2m are found in parts of the central US, Central America and Columbia, where the major land use is the conversion of forest into cropland and grassland. For regions like India and Southeast Asia where forest and grassland are converted into cropland, there is increased daily maximum T 2m . In eastern Europe, the combined effects of deforestation and grassland conversion do  not result in significant changes to either T s or T 2m . In eastern China, although there is no significant change in daily maximum T s , land use can significantly decrease daily maximum T 2m . Conversion of grassland into cropland over the Great Plains tends to substantially decrease both daily maximum T s and T 2m . Additionally, grassland degradation in the Tibetan Plateau leads to significant warming.
During the night, there are consistent changes in T s and T 2m (Figures 3E,G), which show a similar pattern to the changes in daily maximum T s (Figure 3A). Significant cooling is found in the Great Plains and eastern Europe, while warming land surface and 2-m air are found in a majority of the deforested areas.
Irrigation activities can lead to a significant cooling at both the land surface and the 2-m air throughout the day, with greater changes in daytime temperatures. The cooling effects are more local, but comparable to (and even stronger than) the impacts of the historical land use. Table 3 shows the regional changes in daily maximum and monthly mean T 2m and T s due to historical FIGURE 3 | Changes in summer daily maximum (A-D) and minimum (E-H) land surface temperature (T s in K) and 2-m air temperature (T 2m in K) due to historical land use (left) and irrigation (right). Stippling indicates significance at the 95% confidence level.
TABLE 3 | Regional average temperature changes due to land use and irrigation in the irrigated grid cells (with total summertime irrigation greater than 30 mm) in seven sub-regions. Bold numbers indicate over half of the grid cells have statistically significant temperature changes at the 95% confidence level.  Figures 1C,D), significant change is only manifested in the irrigation impacts on daily maximum T s . There are also areas of warming away from the irrigated regions, which result from shifts in the general circulation driven by irrigation. Broader atmospheric effects are discussed further in section "Changes above the surface."

Changes in Surface Fluxes
To understand the mechanism of the temperature responses, we investigate the changes in the diurnal cycle of surface fluxes (Figure 4 for daytime and Figure 5 for nighttime). During the daytime, there is increased latent heat flux (LE) and decreased sensible heat flux (H) in the regions with the conversion of grassland to cropland (Figure 4A), such as the northern Great Plains, Eastern Europe, and northeastern China, where there is decreased T s (Figure 3A). Increased LE is also found in the afforested regions, such as the northeastern US and western Europe, where there is a slight cooling at the land surface. In most areas of the tropics, deforestation leads to significantly decreased LE, suggesting the less evapotranspiration and consequently more H (Figure 4C) over the open land than the forest, resulting in a warming effect at the land surface ( Figure 3A). However, in the United States Midwest, there is an evident decrease in H mainly due to decreased surface roughness after deforestation, which can explain the warming land surface but cooling of the air (Figures 3A,C). Compared to LE and H, ground heat flux (G) responds to land use with relatively small magnitudes during the daytime. Therefore, its change is not a major factor influencing the daytime temperature change. Generally, deforestation leads to more heat flux stored in the ground, while afforestation or grassland conversion results in a reduction of heat entering the ground because of more surface available energy being used as latent heat flux. During nighttime, there is an upward ground heat flux. In other words, the heat stored in the soil layers during the daytime can feed back to the land surface and enter the lower atmosphere. Therefore, changes in G show opposite signs during the daytime and nighttime (Figures 4E, 5E), and the nighttime T s and T 2m (Figures 3E,G) are mainly determined by the changes in G.
Compared with land use, irrigation leads to consistently increased latent heat flux and decreased sensible heat flux throughout the day in the irrigated areas, especially in the Great Plains, California, and northwestern India (Figures 5B,D, 6B,D). Reduced T s and T 2m are mainly associated with the evaporative cooling of irrigation.

Changes Above the Surface
Land use and land management do not only modify the local land surface features; they are also able to influence the atmosphere and large-scale circulation. Therefore, the temperature changes (especially for T 2m ) can also be associated with the atmospheric feedback to land use or land management. Figure 6 shows the changes in summer total cloud cover and precipitation due to land use and irrigation. Land use leads to a significant increase in cloud cover and precipitation over the northern plains of North America and parts of eastern Europe, which corresponds to the cooling effects in those regions. A significant decrease in cloud cover and precipitation is found in northern India and tropical Africa, which corresponds to the land use-induced warming effects in those regions. Irrigation generally increases cloud cover in many regions, possibly due to enhanced local convective processes, but the significant increase in precipitation is only found in northwestern India and limited regions of the Great   Plains in the United States, which exhibit significant cooling effects of irrigation.
Although the changes in large-scale circulation and remote effects of land use/land management are not the scope of this study, the evident atmospheric feedback certainly demonstrates the impacts of land surface change on temperature through land-atmosphere interactions. For instance, the conversion of grassland into cropland in the northern plains of North America enhances evapotranspiration at the land surface, which exerts cooling effects based on the surface energy budget. On the other hand, higher evapotranspiration may lead to more cloud cover, which in turn decreases surface and near-surface temperature, and also potentially enhances precipitation, which further increases soil moisture and surface evapotranspiration.

Changes in Hot Extremes
We also examine the changes in hot extremes due to land use and land management (Figure 7). Generally, deforestation slightly reduces the intensity and frequency of hot events in middle latitudes during daytime, but leads to evidently more intense and frequent hot events during nighttime. The conversion of grassland into cropland, especially in the northern plains, shows significantly reduced hot extremes. On the contrary, the impacts of irrigation on hot extremes are more evident during daytime. Significantly reduced hot extremes are found in the irrigated regions (except eastern China), where the irrigation-induced cooling can be stronger than the land-use impacts. It should be noted that maximum temperature in India occurs in late spring before the monsoon when there is a high irrigation demand. Therefore, strong cooling effects are found in this region, but are not manifested in the summertime temperature responses discussed in previous sections.

DISCUSSION AND CONCLUSION
This study presents a comparison between land use and irrigation with respect to their impacts on summer temperatures. Tropical deforestation shows an evident warming effect, especially at the land surface and during the nighttime. However, landuse impacts in the middle latitudes are dominated by the cooling effect of grassland conversion, while the extensive deforestation in the United States Midwest does not show strong impacts on local or regional temperature. A few caveats should be noted regarding the deforestation-temperature relationship.
First, uncertainties still remain in current land surface models in representing the energy partitioning between latent and sensible heat flux Cai et al., 2019). Unlike the observed decrease in ET after deforestation (Meier et al., 2018), the models suggest a slight increase in latent heat flux (Figure 4A), consequently exerting a cooling effect. Second, the model suggests a warmer land surface but cooler air mainly due to the decreased sensible heat flux after deforestation. The conversion of forest lands to agriculture reduces surface roughness and increases aerodynamic resistance, so the heat transfer becomes less efficient over the open land (Chen and Dirmeyer, 2019a). Consequently, less efficient turbulent heat exchange with the atmosphere leads to land surface warming ( Figure 3A) but cooling of the lower-atmosphere ( Figure 3C). The dominance of surface roughness (or aerodynamic resistance) in deforestationinduced biophysical effects has been documented in previous studies (Liao et al., 2018;Winckler et al., 2019a). The different responses in T s and T 2m to deforestation are also found in another climate model (Winckler et al., 2019a). However, further assessments of land-use impacts on temperature are necessarily taken using multi-model experiments, such as The Land Use Model Intercomparison Project [LUMIP, (Lawrence et al., 2016)]. Furthermore, observation-based studies suggest nighttime cooling (Schultz et al., 2017;Liao et al., 2018), which can be attributed to the decreased turbulence after deforestation and less heat aloft being brought to the surface (Lee et al., 2011). However, our results suggest nighttime warming of deforestation, which is mainly driven by the heat storage in the ground. Similar nighttime warming is also shown in other climate models (Vanden Broucke et al., 2015;Winckler et al., 2019b), suggesting ground heat flux and the nocturnal boundary layer need further attention in future model evaluation and development.
Meanwhile, a strong cooling effect is found in the irrigated areas, especially during the daytime and at the land surface. For instance, irrigation in California reduces daytime T s by over 2 K. Such a cooling effect of irrigation is even greater than the land-use impacts in many regions ( Table 2). Even though the irrigation-induced cooling is also confirmed through groundbased and satellite observations (Bonfils and Lobell, 2007;Lobell et al., 2008;Chen and Dirmeyer, 2019b), caution should be taken to interpret the simulated irrigation effects. The local cooling of irrigation can be even stronger if the model is run at a higher spatial resolution (Kueppers et al., 2007;Lobell et al., 2009;Sacks et al., 2009;Lo and Famiglietti, 2013;Lu and Kueppers, 2015;Huang and Ullrich, 2016;Chen and Dirmeyer, 2019b). For instance, there is up to 7 K of cooling due to irrigation over California based on a regional climate model (e.g., Kueppers et al., 2007), which highlights the impacts of irrigation at local or regional scales.
It also should be noted that uncertainties exist in current irrigation simulations. Observation-based studies suggest a nighttime warming in irrigated areas (Kanamaru and Kanamitsu, 2008;Chen and Jeong, 2018) because irrigation increases soil heat capacity (Kalnay and Cai, 2003) and thermal conductivity (Kanamaru and Kanamitsu, 2008), which allows more heat to be stored during daytime and released from the ground to the air during nighttime. Meanwhile, elevated moisture of the air would enhance downward longwave radiation, which would have a warming effect especially during night (Souri et al., 2020). In our results, although increased ground heat flux (implying more heat entered into the soil) is evident in the irrigated areas during daytime (Figure 4F), decreased ground heat flux (implying more heat released from the ground) is not found during nighttime (Figure 5F), and irrigation leads to a slight cooling in daily minimum temperature (Figures 3F,H). The issue with the simulated response of minimum temperature to irrigation has also been documented in other model-based studies Huang and Ullrich, 2016). As ground heat flux is calculated from the residual of the surface energy balance in order to precisely conserve energy , the model bias may be associated with the uncertainties in surface flux parameterization in land surface models. The uncertainties in the representation of irrigation characteristics are also documented in other climate models, such as the Max Planck Institute for Meteorology's Earth System model (MPI-ESM), and regional climate models. The irrigation-induced temperature change can be highly related to the land-atmosphere coupling scheme used in the model, irrigation effectiveness, timing of irrigation (de Vrese and Hagemann, 2017). Kueppers et al. (2008) also found that the simulated effects of irrigation vary among different regional climate models due to different model physics and irrigation parameterizations. These factors should to considered in future assessments of irrigation impacts, and multiple-model approaches are necessary to further evaluate the uncertainties in the representation of irrigation in different Earth system models.
Additionally, there are possibly remote impacts of land use or irrigation on climate beyond the regions where the land use and land management occur. For instance, there are evident changes in temperature (Figure 3), evapotranspiration (Figure 4), cloud cover (Figure 6) in the northeast of Caspian Sea, although irrigation is only applied in the south of this region. Such a nonlocal effect of irrigation is also shown in previous studies using different climate models (e.g., Puma and Cook, 2010). Although many studies have discussed the remote effects of land use and land management (e.g., Wei et al., 2013;Badger and Dirmeyer, 2016;de Vrese et al., 2016;Swann et al., 2018), the detected remote effects can be influenced by many factors, such as the length of simulations (Swann et al., 2018), representation of ocean circulation (Badger and Dirmeyer, 2016;Krakauer et al., 2016), and model physics (de Vrese et al., 2016), which need further investigation in future studies.
Our results mainly present the climatic impacts of land use and land management during summer because of the strong landatmosphere interactions and intensive agricultural practices in the warm season. We acknowledge land cover change also has significant impacts on temperature during winter, especially for boreal deforestation, which exerts a cooling effect through snowalbedo effects (Davin and de Noblet-Ducoudré, 2010); irrigation also occurs in winter in some regions and has the potential to alter temperature then (Chen and Jeong, 2018). Supplementary Figure S2 shows a significant cooling in the deforested areas in the United States Midwest and Northern Plains and a significant warming effect of tropical deforestation during the boreal winter. However, no evident change in surface temperature is found in the irrigation regions in this study.
As the goal of this paper is to emphasize the relative importance of land management in regional climate  and the implications in hot extreme adaptation under the background of global climate change , our results suggest that the local cooling effects of irrigation can be as strong as the land-use impacts during the boreal summer. As the cropland area is projected to continue growing in all the shared socioeconomic pathways (SSPs) in CMIP6 and about a half of the SSP scenarios suggest an increase in irrigated cropland area in the future (Lawrence et al., 2016), the importance of land management in local and regional climate will be amplified by potentially intensified irrigation activities under a warming climate. Additionally, considering the enhanced ET due to irrigation (Figure 4B), which can potentially raise wetbulb temperatures (Kang and Eltahir, 2018), humidity should also be considered in future assessments of irrigation impacts on heat waves and their social impacts.
In summary, this study highlights the importance of irrigation in local and regional climate based on the experiments using the most recent version of CESM. It should be noted that other land management practices, such no-till farming (Davin et al., 2014) or cover crops (Lombardozzi et al., 2018), may also have certain climatic consequences from the biogeophysical perspective through radiative processes. Also, there are many approaches to irrigation (e.g., surface irrigation, sprinklers including center-pivot irrigation, drip irrigation) drawing from both surface water and groundwater. Current models do not represent this variety, yet each has different potential effects on soil moisture and surface fluxes including evaporative cooling. Therefore, toward a comprehensive understanding of the climatic impacts of land management, different agricultural practices should be considered and compared when assessing regionalscale climate adaptation and mitigation (Lobell et al., 2006;Seneviratne et al., 2018).

DATA AVAILABILITY STATEMENT
The datasets generated for this study are available on request to the corresponding author.

AUTHOR CONTRIBUTIONS
LC and PD designed the study. LC performed the model simulations and the data analysis. Both authors contributed to the writing of the manuscript.

FUNDING
This study was supported by the National Science Foundation (AGS-1419445). LC was also supported by the Illinois State Water Survey, Prairie Research Institute at the University of Illinois in Urbana-Champaign.

ACKNOWLEDGMENTS
All views and opinions expressed do not necessarily reflect those of these institutions. We would like to acknowledge high-performance computing support from Cheyenne (doi: 10. 5065/D6RX99HX) provided by NCAR's Computational and Information Systems Laboratory, sponsored by the National Science Foundation. The Community Earth System Model is freely available at http://www.cesm.ucar.edu/models/cesm2/. We also are grateful to the reviewers whose insightful comments helped improve our manuscript.