Abstract. The Energy Exascale Earth System Model (E3SM) Project is an ongoing, state-of-the-science Earth system modeling, simulation, and prediction project developed by the U.S. Department of Energy (DOE). With an emphasis on supporting DOE's energy mission, understanding and quantifying how well the model simulates water cycle processes is of particular importance. Here, we evaluate E3SM version v1.0 for its ability to represent atmospheric rivers (ARs), which play significant roles in water vapor transport and precipitation. The characteristics and precipitation associated with global ARs in E3SM at standard resolution (1° × 1°) are compared to the Modern-Era Retrospective analysis for Research and Applications, Version 2 (MERRA2). Global pattern of AR frequencies in E3SM show high degrees of correlation (>= 0.97) with MERRA2 and low mean absolute errors (< 1 %) annually, seasonally, and across different ensemble members. However, some large-scale condition biases exist leading to AR biases - most significant of which are: the double-ITCZ, a stronger and/or equatorward shifted subtropical jet during boreal and austral winter, and enhanced northern hemisphere westerlies during summer. By comparing atmosphere-only and fully coupled simulations, we attribute the sources of the biases to the atmospheric component or to a coupling response. Using relationships revealed in Dong et al. (2021), we provide evidence showing the stronger north Pacific jet in winter and enhanced northern hemisphere westerlies during summer associated with E3SM's double-ITCZ and related weaker AMOC, respectively, are the sources of much of the AR biases found in the coupled simulations.