This dataset contains monthly average output files from the iCAM6 simulations used in the manuscript "Enhancing understanding of the hydrological cycle via pairing of process-oriented and isotope ratio tracers," in review at the Journal of Advances in Modeling Earth Systems. A file corresponding to each of the tagged and isotopic variables used in this manuscript is included. Files are at 0.9° latitude x 1.25° longitude, and are in NetCDF format. Data from two simulations are included: 1) a simulation where the atmospheric model was "nudged" to ERA5 wind and surface pressure fields, by adding an additional tendency (see section 3.1 of associated manuscript), and 2) a simulation where the atmospheric state was allowed to freely evolve, using only boundary conditions imposed at the surface and top of atmosphere. Specific information about each of the variables provided is located in the "usage notes" section below. Associated article abstract: The hydrologic cycle couples the Earth's energy and carbon budgets through evaporation, moisture transport, and precipitation. Despite a wealth of observations and models, fundamental limitations remain in our capacity to deduce even the most basic properties of the hydrological cycle, including the spatial pattern of the residence time (RT) of water in the atmosphere and the mean distance traveled from evaporation sources to precipitation sinks. Meanwhile, geochemical tracers such as stable water isotope ratios provide a tool to probe hydrological processes, yet their interpretation remains equivocal despite several decades of use. As a result, there is a need for new mechanistic tools that link variations in water isotope ratios to underlying hydrological processes. Here we present a new suite of “process-oriented tags,” which we use to explicitly trace hydrological processes within the isotopically enabled Community Atmosphere Model, version 6 (iCAM6). Using these tags, we test the hypotheses that precipitation isotope ratios respond to parcel rainout, variations in atmospheric RT, and preserve information regarding meteorological conditions during evaporation. We present results for a historical simulation from 1980 to 2004, forced with winds from the ERA5 reanalysis. We find strong evidence that precipitation isotope ratios record information about atmospheric rainout and meteorological conditions during evaporation, but little evidence that precipitation isotope ratios vary with water vapor RT. These new tracer methods will enable more robust linkages between observations of isotope ratios in the modern hydrologic cycle or proxies of past terrestrial environments and the environmental processes underlying these observations. Details about the simulation setup can be found in section 3 of the associated open-source manuscript, "Enhancing understanding of the hydrological cycle via pairing of process‐oriented and isotope ratio tracers." In brief, we conducted two simulations of the atmosphere from 1980-2004 using the isotope-enabled version of the Community Atmosphere Model 6 (iCAM6) at 0.9x1.25° horizontal resolution, and with 30 vertical hybrid layers spanning from the surface to ~3 hPa. In the first simulation, wind and surface pressure fields were "nudged" toward the ERA5 reanalysis dataset by adding a nudging tendency, preventing the model from diverging from observed/reanalysis wind fields. In the second simulation, no additional nudging tendency was included, and the model was allowed to evolve 'freely' with only boundary conditions provided at the top (e.g., incoming solar radiation) and bottom (e.g., observed sea surface temperatures) of the model. In addition to the isotopic variables, our simulation included a suite of 'process-oriented tracers,' which we describe in section 2 of the manuscript. These variables are meant to track a property of water associated with evaporation, condensation, or atmospheric transport. Metadata are provided about each of the files below; moreover, since the attached files are NetCDF data - this information is also provided with the data files. NetCDF metadata can be accessed using standard tools (e.g., ncdump). Each file has 4 variables: the tagged quantity, and the associated coordinate variables (time, latitude, longitude). The latter three are identical across all files, only the tagged quantity changes. Twelve files are provided for the nudged simulation, and an additional three are provided for the free simulations: Nudged simulation files iCAM6_nudged_1980-2004_mon_RHevap: Mass-weighted mean evaporation source property: RH (%) with respect to surface temperature. iCAM6_nudged_1980-2004_mon_Tevap: Mass-weighted mean evaporation source property: surface temperature in Kelvin iCAM6_nudged_1980-2004_mon_Tcond: Mass-weighted mean condensation property: temperature (K) iCAM6_nudged_1980-2004_mon_columnQ: Total (vertically integrated) precipitable water (kg/m2). Not a tagged quantity, but necessary to calculate depletion times in section 4.3 (e.g., Fig. 11 and 12). iCAM6_nudged_1980-2004_mon_d18O: Precipitation d18O (‰ VSMOW) iCAM6_nudged_1980-2004_mon_d18Oevap_0: Mass-weighted mean evaporation source property - d18O of the evaporative flux (e.g., the 'initial' isotope ratio prior to condensation), (‰ VSMOW) iCAM6_nudged_1980-2004_mon_dxs: Precipitation deuterium excess (‰ VSMOW) - note that precipitation d2H can be calculated from this file and the precipitation d18O as d2H = d-excess - 8*d18O. iCAM6_nudged_1980-2004_mon_dexevap_0: Mass-weighted mean evaporation source property - deuterium excess of the evaporative flux iCAM6_nudged_1980-2004_mon_lnf: Integrated property - ln(f) calculated from the constant-fractionation d18O tracer (see section 3.2). iCAM6_nudged_1980-2004_mon_precip: Total precipitation rate in m/s. Note there is an error in the metadata in this file - it is total precipitation, not just convective precipitation. iCAM6_nudged_1980-2004_mon_residencetime: Mean atmospheric water residence time (in days). iCAM6_nudged_1980-2004_mon_transportdistance: Mean atmospheric water transport distance (in km). Free simulation files iCAM6_free_1980-2004_mon_d18O: Precipitation d18O (‰ VSMOW) iCAM6_free_1980-2004_mon_dxs: Precipitation deuterium excess (‰ VSMOW) - note that precipitation d2H can be calculated from this file and the precipitation d18O as d2H = d-excess - 8*d18O. iCAM6_free_1980-2004_mon_precip: Total precipitation rate in m/s. Note there is an error in the metadata in this file - it is total precipitation, not just convective precipitation.
more »
« less
Conservation of Dry Air, Water, and Energy in CAM and Its Potential Impact on Tropical Rainfall
Abstract For the Community Atmosphere Model version 6 (CAM6), an adjustment is needed to conserve dry air mass. This adjustment exposes an inconsistency in how CAM6’s energy budget incorporates water—in CAM6 water in the vapor phase has energy, but condensed phases of water do not. When water vapor condenses, only its latent energy is retained in the model, while its remaining internal, potential, and kinetic energy are lost. A global fixer is used in the default CAM6 model to maintain global energy conservation, but locally the energy tendency associated with water changing phase violates the divergence theorem. This error in energy tendency is intrinsically tied to the water vapor tendency, and reaches its highest values in regions of heavy rainfall, where the error can be as high as 40 W m −2 annually averaged. Several possible changes are outlined within this manuscript that would allow CAM6 to satisfy the divergence theorem locally. These fall into one of two categories: 1) modifying the surface flux to balance the local atmospheric energy tendency and 2) modifying the local atmospheric tendency to balance the surface plus top-of-atmosphere energy fluxes. To gauge which aspects of the simulated climate are most sensitive to this error, the simplest possible change—where condensed water still does not carry energy and a local energy fixer is used in place of the global one—is implemented within CAM6. Comparing this experiment with the default configuration of CAM6 reveals precipitation, particularly its variability, to be highly sensitive to the energy budget formulation. Significance Statement This study examines and explains spurious regional sources and sinks of energy in a widely used climate model. These energy errors result from not tracking energy associated with water after it transitions from the vapor phase to either liquid or ice. Instead, the model used a global fixer to offset the energy tendency related to the energy sources and sinks associated with condensed water species. We replace this global fixer with a local one to examine the model sensitivity to the regional energy error and find a large sensitivity in the simulated hydrologic cycle. This work suggests that the underlying thermodynamic assumptions in the model should be revisited to build confidence in the model-simulated regional-scale water and energy cycles.
more »
« less
- Award ID(s):
- 1912134
- NSF-PAR ID:
- 10333885
- Date Published:
- Journal Name:
- Journal of Climate
- Volume:
- 35
- Issue:
- 9
- ISSN:
- 0894-8755
- Page Range / eLocation ID:
- 2895 to 2917
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
more » « less
Abstract A comparative analysis between observational data from McMurdo Station, Antarctica and the Community Atmosphere Model version 6 (CAM6) simulation is performed focusing on cloud characteristics and their thermodynamic conditions. Ka‐band Zenith Radar (KAZR) and High Spectral Resolution Lidar (HSRL) retrievals are used as the basis of cloud fraction and cloud phase identifications. Radiosondes released at 12‐h increments provide atmospheric profiles for evaluating the simulated thermodynamic conditions. Our findings show that the CAM6 simulation consistently overestimates (underestimates) cloud fraction above (below) 3 km in four seasons of a year. Normalized by total in‐cloud samples, ice and mixed phase occurrence frequencies are underestimated and liquid phase frequency is overestimated by the model at cloud fractions above 0.6, while at cloud fractions below 0.6 ice phase frequency is overestimated and liquid‐containing phase frequency is underestimated by the model. The cloud fraction biases are closely associated with concurrent biases in relative humidity (RH), that is, high (low) RH biases above (below) 2 km. Frequencies of correctly simulating ice and liquid‐containing phase increase when the absolute biases of RH decrease. Cloud fraction biases also show a positive correlation with RH biases. Water vapor mixing ratio biases are the primary contributor to RH biases, and hence, likely a key factor controlling the cloud biases. This diagnosis of the evident shortfalls of representations of cloud characteristics in CAM6 simulation at McMurdo Station brings new insight in improving the governing model physics therein.
-
Abstract The Indian Ocean has an intriguing intertropical convergence zone (ITCZ) south of the equator year-round, which remains largely unexplored. Here we investigate this Indian Ocean ITCZ and the mechanisms for its origin. With a weak semiannual cycle, this ITCZ peaks in January–February with the strongest rainfall and southernmost location and a northeast–southwest orientation from the Maritime Continent to Madagascar, reaches a minimum around May with a zonal orientation, grows until its secondary maximum around September with a northwest–southeast orientation, weakens slightly until December, and then regains its mature phase in January. During austral summer, the Indian Ocean ITCZ exists over maximum surface moist static energy (MSE), consistent with convective quasi-equilibrium theory. This relationship breaks up during boreal summer when the surface MSE maximizes in the northern monsoon region. The position and orientation of the Indian Ocean ITCZ can be simulated well in both a linear dynamical model and the state-of-the-art Community Atmosphere Model version 6 (CAM6) when driven by observed sea surface temperature (SST). To quantify the contributions of the planetary boundary layer (PBL) and free-atmosphere processes to this ITCZ, we homogenize the free-atmosphere diabatic heating over the Indian Ocean in CAM6. In response, the ITCZ weakens significantly, owing to a weakened circulation and deep convection. Therefore, in CAM6, the SST drives the Indian Ocean ITCZ directly through PBL processes and indirectly via free-atmosphere diabatic heating. Their contributions are comparable during most seasons, except during the austral summer when the free-atmosphere diabatic heating dominates the mature-phase ITCZ. Significance Statement The intertropical convergence zone (ITCZ) is the globe-encircling band where trade winds converge and strong rainfall occurs in the tropics. Its rains provide life-supporting water to billions of people. Its associated latent heating invigorates the tropical atmospheric circulation and influences climate and weather across the planet. The ITCZ is located north of the equator in most tropical oceans, except in the Indian Ocean where it sits south of the equator year-around. In contrast to the well-known northern ITCZs, the origin of the southern ITCZ in the Indian Ocean remains unknown. This work provides the first explanation for how ocean surface temperature works together with processes in the lower and upper atmosphere to shape the unique ITCZ in the Indian Ocean.more » « less
-
more » « less
Abstract Three climate models are evaluated using in situ airborne observations from the Southern Ocean Clouds, Radiation, Aerosol Transport Experimental Study (SOCRATES) campaign. The evaluation targets cloud phases, microphysical properties, thermodynamic conditions, and aerosol indirect effects from −40°C to 0°C. Compared with 580‐s averaged observations (i.e., 100 km horizontal scale), the Community Atmosphere Model version 6 (CAM6) shows the most similar result for cloud phase frequency distribution and allows more liquid‐containing clouds below −10°C compared with its predecessor—CAM5. The Energy Exascale Earth System Model (E3SM) underestimates (overestimates) ice phase frequencies below (above) −20°C. CAM6 and E3SM show liquid and ice water contents (i.e., LWC and IWC) similar to observations from −25°C to 0°C, but higher LWC and lower IWC than observations at lower temperatures. Simulated in‐cloud RH shows higher minimum values than observations, possibly restricting ice growth during sedimentation. As number concentrations of aerosols larger than 500 nm (Na500) increase, observations show increases of LWC, IWC, liquid, and ice number concentrations (Nliq, Nice). Number concentrations of aerosols larger than 100 nm (Na100) only show positive correlations with LWC and Nliq. From −20°C to 0°C, higher aerosol number concentrations are correlated with lower glaciation ratio and higher cloud fraction. From −40°C to −20°C, large aerosols show positive correlations with glaciation ratio. CAM6 shows small increases of LWC and Nliqwith Na500and Na100. E3SM shows small increases of Nicewith Na500. Overall, CAM6 and E3SM underestimate aerosol indirect effects on ice crystals and supercooled liquid droplets over the Southern Ocean.
-
Aerosols affect cirrus formation and evolution, yet quantification of these effects remain difficult based on in-situ observations due to the complexity of nucleation mechanisms and large variabilities in ice microphysical properties. This work employed a method to distinguish five evolution phases of cirrus clouds based on in-situ aircraft-based observations from seven U.S. National Science Foundation (NSF) and five NASA flight campaigns. Both homogeneous and heterogeneous nucleation were captured in the 1-Hz aircraft observations, inferred from the distributions of relative humidity in the nucleation phase. Using linear regressions to quantify the correlations between cirrus microphysical properties and aerosol number concentrations, we found that ice water content (IWC) and ice crystal number concentration (Ni) show strong positive correlations with larger aerosols (> 500 nm) in the nucleation phase, indicating strong contributions of heterogeneous nucleation when ice crystals first start to nucleate. For the later growth phase, IWC and Ni show similar positive correlations with larger and smaller (i.e., > 100 nm) aerosols, possibly due to fewer remaining ice nucleating particles in the later growth phase that allows more homogeneous nucleation to occur. Both 200-m and 100-km observations were compared with the nudged simulations from the National Center for Atmospheric Research (NCAR) Community Atmosphere Model version 6 (CAM6). Simulated aerosol indirect effects are weaker than the observations for both larger and smaller aerosols. Observations show stronger aerosol indirect effects (i.e., positive correlations between IWC, Ni and Na) in the Southern Hemisphere (SH) compared with the Northern Hemisphere (NH), while the simulations show negative correlations in the SH. The simulations underestimate IWC by a factor of 3 – 30 in the early/later growth phase, indicating that the low bias of simulated IWC was due to insufficient ice particle growth. Such hypothesis is consistent with the model biases of lower frequencies of ice supersaturation and lower vertical velocity standard deviation in the early/later growth phases. Overall, these findings show that aircraft observations can capture the competitions between heterogeneous and homogeneous nucleation, and their contributions vary as cirrus clouds evolve. Future model development is also recommended to evaluate and improve the representation of water vapor and vertical velocity on the sub-grid scale to resolve the insufficient ice particle growth.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1175/JCLI-D-21-0512.1
