Search

In this paper, we investigate the relationship between the El Niño–Southern Oscillation (ENSO) spring persistence barrier (PB) and predictability barrier (PD) and apply it to explain the interdecadal modulation of ENSO prediction skill using the anomaly correlation coefficient (ACC). Previous studies showed that a longer persistence (i.e., autocorrelation) tends to produce a higher prediction skill. Using the recharge oscillator model of ENSO, both analytical and numerical solutions suggest that the predictability (i.e., ACC) is related to the persistence of sea surface temperature (SST) and cross correlation between SST and subsurface ocean heat content in the tropical Pacific. In particular, a larger damping rate in SST anomalies will lead to a lower persistence and ACC and a stronger PD. However, a shortened ENSO period, which controls the cross correlation, will lead to a lower persistence but a higher ACC associated with a weaker PD. Finally, we apply our solutions to observations and suggest that a higher ACC associated with a weaker PD after 1960 is caused by the shortened ENSO period. [ABSTRACT FROM AUTHOR]

This paper examines the processes that drive Arctic anomalous surface warming and sea ice loss during winter-season tropospheric energy flux events, synoptic periods of increased tropospheric energy flux convergence (Ftrop), using the NASA MERRA-2 reanalysis. During an event, a poleward anomaly in Ftrop initially increases the sensible and latent energy of the Arctic troposphere; as the warm and moist troposphere loses heat, the anomalous energy source is balanced by a flux upward across the tropopause and a downward net surface flux. A new metric for the Arctic surface heating efficiency (Etrop) is defined, which measures the fraction of the energy source that reaches the surface. Composites of high-, medium-, and low-efficiency events help identify key physical factors, including the vertical structure of Ftrop and Arctic surface preconditioning. In high-efficiency events (Etrop ≥ 0.63), a bottom-heavy poleward Ftrop occurs in the presence of an anomalously warm and unstratified Arctic—a consequence of decreased sea ice—resulting in increased vertical mixing, enhanced near-surface warming and moistening, and further sea ice loss. Smaller Etrop, and thus weaker surface impacts, are found in events with anomalously large initial sea ice extent and more vertically uniform Ftrop. These differences in Etrop are manifested primarily through turbulent heat fluxes rather than downward longwave radiation. The frequency of high-efficiency events has increased from the period 1980–99 to the period 2000–19, contributing to Arctic surface warming and sea ice decline. [ABSTRACT FROM AUTHOR]

An interesting aspect of the El Ni˜no-Southern Oscillation (ENSO) phenomenon is the asymmetry between its two phases. This paper evaluates the simulations of this property of ENSO by the CoupledModel Intercomparison Project phase 6 (CMIP6) models. Both the surface and subsurface signals of ENSO are examined for this purpose. The results show that the models still underestimate ENSO asymmetry as shown in the SST field, but do a better job in the subsurface. A much weaker negative feedback from the net surface heat flux during La Ni˜na in the models is identified as a factor causing the degradation of the ENSO asymmetry at the surface. The simulated asymmetry in the subsurface is still weaker than the observations owing to a weaker dynamic coupling between the atmosphere and ocean. Consistent with the finding of a weaker dynamic coupling strength, the precipitation response to the SST changes is also found to be weaker in the models. The results underscore that a more objective assessment of the simulation of ENSO by climate models may have to involve the examination of the subsurface signals. Future improvements in simulating ENSO will likely require a better simulation of the surface heat flux feedback from the atmosphere as well as the dynamical coupling strength between the atmosphere and ocean. [ABSTRACT FROM AUTHOR]

Understanding the contribution of ocean circulation to glacial-interglacial climate change is a major focus of paleoceanography. Specifically, many have tried to determine whether the volumes and depths of Antarctic- and North Atlantic-sourced waters in the deep ocean changed at the Last Glacial Maximum (LGM; ∼22-18 kyr BP) when atmospheric pCO2 concentrations were 100 ppm lower than the preindustrial. Measurements of sedimentary geochemical proxies are the primary way that these deep ocean structural changes have been reconstructed. However, the main proxies used to reconstruct LGM Atlantic water mass geometry provide conflicting results as to whether North Atlantic-sourced waters shoaled during the LGM. Despite this, a number of idealized modeling studies have been advanced to describe the physical processes resulting in shoaled North Atlantic waters. This paper aims to critically assess the approaches used to determine LGMAtlantic circulation geometry and lay out best practices for future work. We first compile existing proxy data and paleoclimate model output to deduce the processes responsible for setting the ocean distributions of geochemical proxies in the LGM Atlantic Ocean. We highlight how small-scale mixing processes in the ocean interior can decouple tracer distributions from the large-scale circulation, complicating the straightforward interpretation of geochemical tracers as proxies for water mass structure. Finally, we outline promising paths toward ascertaining the LGMcirculation structure more clearly and deeply. [ABSTRACT FROM AUTHOR]

This paper examines the physical controls of extratropical humidity and clouds by isolating the effects of cloud physics factors in an idealized model. The Held–Suarez dynamical core is used with the addition of passive water vapor and cloud tracers, allowing cloud processes to be explored cleanly. Separate saturation adjustment and full cloud scheme controls are used to consider the strength of advection–condensation theory. Three sets of perturbations to the cloud scheme are designed to test the model’s sensitivity to the physics of condensation, sedimentation, and precipitation formation. The condensation and sedimentation perturbations isolate two key differences between the control cases. First, the sub-grid-scale relative humidity distribution assumed for the cloud macrophysics influences the location and magnitude of the extratropical cloud maxima, which interrupt the isentropic transport of moisture to the polar troposphere. Second, within the model’s explicit treatment of cloud microphysics, re-evaporation of hydrometeors moistens and increases clouds in the lower troposphere. In contrast, microphysical processes of precipitation formation (specifically, the ratio of accretion to autoconversion) have negligible effects on humidity, cloudiness, and precipitation apart from the strength of the largescale condensation and formation cycle. In addition, counterintuitive relationships}such as cloud condensate and cloud fraction responding in opposing directions}emphasize the need for careful dissection of physical mechanisms. In keeping with advection–condensation theory, circulation sets the patterns of humidity, clouds, and precipitation to first order, with factors explored herein providing secondary controls. The results substantiate the utility of such idealized modeling and highlight key cloud processes to constrain. [ABSTRACT FROM AUTHOR]

Global climate models, like any in silico numerical experiments, are affected by different types of bias. Uncertainty quantification remains a challenge in any climate detection and attribution analysis. A fundamental methodological question is to determine which statistical summaries, while bringing relevant signals, can be robust with respect to multimodel errors. In this paper, we propose a simple statistical framework that significantly improves signal detection in climate attribution studies. We show that the complex bias correction step can be entirely bypassed for models for which bias between the simulated and unobserved counterfactual worlds is the same as between the simulated and unobserved factual worlds. To illustrate our approach, we infer emergence times in precipitation from the CMIP5 and CMIP6 archives. The detected anthropogenic signal in yearly maxima of daily precipitation clearly emerges at the beginning of the twenty-first century. In addition, no CMIP model seems to outperform the others and a weighted linear combination of all improves the estimation of emergence times. Significance Statement: We show that the bias in multimodel global climate simulations can be efficiently handled when the appropriate metric is chosen. This metric leads to an easy-to-implement statistical procedure based on a checkable assumption. This allows us to demonstrate that optimal convex combinations of CMIP outputs can improve the signal strength in finding emergence times. Our data analysis procedure is applied to yearly maximum of precipitation from CMIP5 and CMIP6 databases. The attribution of the anthropogenic forcing clearly emerges in extreme precipitation at the beginning of the twenty-first century. [ABSTRACT FROM AUTHOR]

Although tropical interbasin interactions at interannual time scales are presently receiving much attention, their controlling factors and variations on longer time scales are under debate. Tropical convection plays a crucial role in the occurrence and nonstationarity of them. In this paper, we investigate the dependence of interannual tropical Atlantic–Pacific basin interactions on convection-related features of the tropical oceans' climatology, especially the ITCZ position. We contrast a CGCM control simulation with an experiment in which tropical convection is modified by an artificial perturbation outside the tropics that reduces the incident shortwave radiation in a region of the South Atlantic. Based on previous work, this modification is expected to shift in latitude the climatological position of the simulated ITCZ. The experiment shows altered Walker circulations, stronger interannual variability over the tropical oceans, a westward extension of the Atlantic Niño pattern and of convection, and shallower thermocline in the Pacific, making the basin more sensitive to both local and remote perturbations. As a consequence, the experiment shows enhanced interannual Atlantic–Pacific basin interactions at the equator, and weaker teleconnections between the north tropical Atlantic and the equatorial Pacific. The latter seems to occur because the impact of the warm Atlantic SST anomalies is offset by the presence of warm SST anomalies in El Niño region. Despite the uncertainties raised because the simulations are relatively short, we conclude that this work presents a potential explanation for the long-term changes in the tropical basin interactions and offers a novel and useful methodology for their analysis. [ABSTRACT FROM AUTHOR]

In southern Africa, models from phase 5 of the Coupled Model Intercomparison Project (CMIP5) predict robust future drying associated with a delayed rainy-season onset in the austral spring and a range of wetting and drying patterns in the austral summer. This paper relates these rainfall changes to dynamical shifts in two classes of weather systems: the Congo Air Boundary (CAB) and tropical lows. Objective algorithms are used to track these features in CMIP5 model output. It is then established that the climatological locations and frequencies of these systems are reasonably well represented in the CMIP5 models. RCP8.5 end-of-twenty-first-century projections are compared with historical end-of-twentieth-century simulations. Future projections in tropical-low locations and frequencies diverge, but indicate an overall average decrease of 15% and in some cases a northward shift. The projected spatial change in the tropical-low frequency distribution is weakly positively correlated to the projected spatial change in the austral summer rainfall distribution. Meanwhile, future projections indicate a 13% increase in CAB frequency from October to December. This is associated with the gradual climatological CAB breakdown occurring half a month later on average in end-of-twenty-first-century RCP8.5 projections. A delay in the gradual seasonal decline of the CAB prevents rainfall to the south of the CAB’s mean position, most of which is shown to occur on CAB breakdown days, hence creating the austral spring drying signal and delayed wet-season onset. Intermodel variability in the magnitude of CAB frequency increase is able to explain intermodel variability in the projected drying. [ABSTRACT FROM AUTHOR]

Cyclones can be identified from gridded pressure data at different levels of the troposphere, with vertical structure known to influence the temporal development and impacts of midlatitude cyclones. However, studies of midlatitude cyclones typically focus on cyclones identified at a single atmospheric level. This paper examines how the frequency of vertically organized or deep cyclones varies around the world, with a focus on southeastern Australia. About 50% of global cyclones identified from mean sea level pressure show a coherent vertical structure extending to at least 500 hPa, based on ERA-Interim reanalysis data, and shallow cyclones are most common in the global midlatitudes. Using a combination of reanalysis data and satellite-based rainfall and lightning, we show that in southeast Australia deep cyclones have higher intensities, longer durations, and more severe winds and rainfall than either shallow surface cyclones or upper-level cyclones with no surface low, motivating a three-dimensional approach for future cyclone analyses. [ABSTRACT FROM AUTHOR]

Better understanding of how moisture, clouds, and precipitation covary under climate warming lacks a comprehensive observational view. This paper analyzes the tropical atmospheric water cycle's evolution with sea surface temperature (SST), using for the first time, the synergistic dataset of instantaneous observations of the relative humidity profile from the Megha-Tropiques satellite, clouds from the CALIPSO satellite, and near-surface precipitation from the CloudSat satellite, and quantifies their rates of change with SST warming. The dataset is partitioned into three vertical velocity regimes, with cloudy grid boxes categorized by phase (ice or liquid), opacity (opaque or thin), and the presence of near-surface precipitation. Opaque cloud cover is always larger in the presence of near-surface precipitation (high ice clouds especially). Low liquid water clouds in the descending regime dominate for SSTs < 299.25 K, where the free troposphere is dry (~20%), and opaque liquid water cloud cover decreases with SST warming (−8% K−1) and thin liquid water cloud cover stays constant (~20%). High ice clouds dominate the ascending regime in which, for 299.25 < SST < 301.75 K, humidity increases with SST in the lower free troposphere and peaks around 302 K. Over the warm SST range (>301.75 K), in the ascending regime, opaque high ice cloud cover decreases with SST (−13% K−1), while thin ice cloud cover increases (+6% K−1). Over the warm SST range, total cloudiness decreases with warming in all regimes. This paper characterizes fundamental relationships between aspects of the tropical atmospheric water cycle and SST. [ABSTRACT FROM AUTHOR]

Evidence from previous climate model simulations has suggested a potentially low efficacy of contrails to force global mean surface temperature changes. In this paper, a climate model with a state-of-the-art contrail cirrus representation is used for fixed sea surface temperature simulations in order to determine the effective radiative forcing (ERF) from contrail cirrus. ERF is expected to be a good metric for intercomparing the quantitative importance of different contributions to surface temperature and climate impact. Substantial upscaling of aviation density is necessary to ensure statistically significant results from our simulations. The contrail cirrus ERF is found to be less than 50% of the respective instantaneous or stratosphere adjusted radiative forcings, with a best estimate of roughly 35%. The reduction of ERF is much more substantial for contrail cirrus than it is for a CO2 increase when both stratosphere adjusted forcings are of similar magnitude. Analysis of all rapid radiative adjustments contributing to the ERF indicates that the reduction is mainly induced by a compensating effect of natural clouds that provide a negative feedback. Compared to the CO2 reference case, a less positive combined water vapor and lapse rate adjustment also contributes to a more distinct reduction of contrail cirrus ERF, but not as much as the natural cloud adjustment. Based on the experience gained in this paper, respective contrail cirrus simulations with interactive ocean will be performed as the next step toward establishing contrail cirrus efficacy. ERF results of contrail cirrus from other climate models equipped with suitable parameterizations are regarded as highly desirable. [ABSTRACT FROM AUTHOR]

A key question in climate modeling is to what extent sea surface temperature and upper-ocean heat content are driven passively by air-sea heat fluxes, as opposed to forcing by ocean dynamics. This paper investigates the question using a climate model at different resolutions, and observations, for monthly variability. At the grid scale in a high-resolution climate model with resolved mesoscale ocean eddies, ocean dynamics (i.e., ocean heat flux convergence) dominates upper 50 m heat content variability over most of the globe. For deeper depths of integration to 400 m, the heat content variability at the grid scale is almost totally controlled by ocean heat flux convergence. However, a strong dependence on spatial scale is found--for the upper 50 m of ocean, after smoothing the data to around 7°, air-sea heat fluxes, augmented by Ekman heat transports, dominate. For deeper depths of integration to 400 m, the transition scale becomes larger and is above 10° in western boundary currents. Comparison of climate model results with observations show that the small-scale influence of ocean intrinsic variability is well captured by the high-resolution model but is missing from a comparable model with parameterized ocean-eddy effects. In the deep tropics, ocean dynamics dominates in all cases and all scales. In the subtropical gyres at large scales, air-sea heat fluxes play the biggest role. In the midlatitudes, at large scales >10°, atmosphere-driven air-sea heat fluxes and Ekman heat transport variability are the dominant processes except in the western boundary currents for the 400 m heat content. [ABSTRACT FROM AUTHOR]

Precipitation sustains life and supports human activities, making its prediction one of the most societally relevant challenges in weather and climate modeling. Limitations in modeling precipitation underscore the need for diagnostics and metrics to evaluate precipitation in simulations and predictions. While routine use of basic metrics is important for documenting model skill, more sophisticated diagnostics and metrics aimed at connecting model biases to their sources and revealing precipitation characteristics relevant to how model precipitation is used are critical for improving models and their uses. This paper illustrates examples of exploratory diagnostics and metrics including 1) spatiotemporal characteristics metrics such as diurnal variability, probability of extremes, duration of dry spells, spectral characteristics, and spatiotemporal coherence of precipitation; 2) process-oriented metrics based on the rainfall–moisture coupling and temperature–water vapor environments of precipitation; and 3) phenomena-based metrics focusing on precipitation associated with weather phenomena including low pressure systems, mesoscale convective systems, frontal systems, and atmospheric rivers. Together, these diagnostics and metrics delineate the multifaceted and multiscale nature of precipitation, its relations with the environments, and its generation mechanisms. The metrics are applied to historical simulations from phases 5 and 6 of the Coupled Model Intercomparison Project. Models exhibit diverse skill as measured by the suite of metrics, with very few models consistently ranked as top or bottom performers compared to other models in multiple metrics. Analysis of model skill across metrics and models suggests possible relationships among subsets of metrics, motivating the need for more systematic analysis to understand model biases for informing model development. [ABSTRACT FROM AUTHOR]

In a recent paper, we argued that ocean dynamics increase the variability of midlatitude sea surface temperatures (SSTs) on monthly to interannual time scales, but act to damp lower-frequency SST variability over broad midlatitude regions. Here, we use two configurations of a simple stochastic climate model to provide new insights into this important aspect of climate variability. The simplest configuration includes the forcing and damping of SST variability by observed surface heat fluxes only, and the more complex configuration includes forcing and damping by ocean processes, which are estimated indirectly from monthly observations. It is found that the simple model driven only by the observed surface heat fluxes generally produces midlatitude SST power spectra that are too red compared to observations. Including ocean processes in the model reduces this discrepancy by whitening the midlatitude SST spectra. In particular, ocean processes generally increase the SST variance on <2-yr time scales and decrease it on >2-yr time scales. This happens because oceanic forcing increases the midlatitude SST variance across many time scales, but oceanic damping outweighs oceanic forcing on >2-yr time scales, particularly away from the western boundary currents. The whitening of midlatitude SST variability by ocean processes also operates in NCAR's Community Earth System Model (CESM). That is, midlatitude SST spectra are generally redder when the same atmospheric model is coupled to a slab rather than dynamically active ocean model. Overall, the results suggest that forcing and damping by ocean processes play essential roles in driving midlatitude SST variability. [ABSTRACT FROM AUTHOR]

The literature offers conflicting findings about which regions contribute most to increases in the global radiative feedback after a forcing increase. This paper explains the disagreement by discriminating between two common definitions of the local feedback, which use either local temperature or global temperature as their basis. Although the two definitions of feedback have been previously compared in aquaplanet models with slab oceans, here the definitions are compared for the first time in an atmosphere-ocean general circulation model (MPI-ESM1.2) integrated over four doublings of atmospheric CO2 concentrations. Large differences between the definitions can be seen in all feedbacks, but especially in the temperature and water vapor feedbacks. Differences of up to 10 W m-2 K-1 over the Southern Ocean can be explained by the pattern of surface warming, which weights the local feedbacks and reduces their contribution to the global mean. This finding is, however, dependent on the resolution of analysis, because the local-temperature definition is mathematically inconsistent across spatial scales. Furthermore, attempts to estimate the effect of "pattern weighting" by separating local feedbacks and warming patterns at the gridcell level fail, because the radiative change in key tropical regions is also determined by tropospheric stability via the global circulation. These findings indicate that studies of regional feedback change are more sensitive to methodological choices than previously thought, and that the tropics most likely dominate regional contributions to global radiative feedback change on decadal to centennial time scales. [ABSTRACT FROM AUTHOR]

Climate models project an intensification of the wintertime North Atlantic Ocean storm track, over its downstream region, by the end of this century. Previous studies have suggested that ocean–atmosphere coupling plays a key role in this intensification, but the precise role of the different components of the coupling has not been explored and quantified. In this paper, using a hierarchy of ocean coupling experiments, we isolate and quantify the respective roles of thermodynamic (changes in surface heat fluxes) and dynamic (changes in ocean heat flux convergence) ocean coupling in the projected intensification of North Atlantic transient eddy kinetic energy (TEKE). We show that dynamic coupling accounts for nearly all of the future TEKE strengthening as it overcomes the much smaller effect of surface heat flux changes to weaken the TEKE. We further show that by reducing the Arctic amplification in the North Atlantic, ocean heat flux convergence increases the meridional temperature gradient aloft, causing a larger eddy growth rate and resulting in the strengthening of North Atlantic TEKE. Our results stress the importance of better monitoring and investigating the changes in ocean heat transport, for improving climate change adaptation strategies. Significance Statement: By the end of this century, the North Atlantic Ocean storm track is projected to intensify on its eastward flank. Such intensification will have large societal impacts, mostly over western Europe. Thus, it is critical to better understand the mechanism underlying the intensification of the storm track. Here we investigate the role of ocean coupling in the future intensification of the North Atlantic storm track and find that ocean heat transport processes are responsible for the strengthening of the storm track. Our results suggest that better monitoring the changes in ocean heat transport will hopefully improve climate change adaption strategies. [ABSTRACT FROM AUTHOR]

Projections of future West African monsoon (WAM) precipitation change in response to increasing greenhouse gases are uncertain, and an improved understanding of the drivers of WAM precipitation change is needed to help aid model development and better inform adaptation policies in the region. This paper addresses one of those drivers, the direct radiative effect of increased CO2 (i.e., the impact of increased CO2 in the absence of SST warming and changes in plant physiology). An atmosphere-only model is used to examine both the equilibrium response and the evolution of the change over the days following the instantaneous CO2 increase. In response to the direct radiative effect, WAM precipitation increases due to a weakening of the shallow meridional circulation over North Africa, advecting less dry air into the convective column associated with the monsoon. Changes in the shallow circulation are associated with atmospheric and surface warming patterns over North Africa. A large-scale atmospheric warming pattern, whereby North Africa warms more than the monsoon region, leads to a northward shift in the Saharan heat low. In response to increased precipitation in the Sahel, local soil moisture feedbacks play a key role in determining the low-level circulation change and the location of the intertropical discontinuity. The large-scale warming patterns over North Africa result from differing levels of constraint applied by convective quasi-equilibrium. While this constraint acts strongly in the equatorial WAM region, preventing the region from warming in response to the direct radiative effect, North Africa is not strongly constrained and is therefore able to warm. [ABSTRACT FROM AUTHOR]

Realistic ocean subsurface simulations of thermal structure and variation are critically important to success in climate prediction and projection; currently, substantial systematic subsurface biases still exist in the state-of-the-art ocean and climate models. In this paper, subsurface biases in the tropical Atlantic Ocean (TA) are investigated by analyzing simulations from the Ocean Model Intercomparison Project (OMIP) and conducting ocean-only experiments that are based on the Parallel Ocean Program, version 2 (POP2). The subsurface biases are prominent in almost all OMIP simulations, characterized by two warm-bias patches off the equator. By conducting two groups of POP2-based ocean-only experiments, two potential origins of the biases are explored, including uncertainties in wind forcing and vertical mixing parameterization, respectively. It is illustrated that the warm bias near 10°N can be slightly reduced by modulating the prescribed wind field, and the warm biases over the entire basin are significantly reduced by reducing background diffusivity in the ocean interior in ways to match observations. By conducting a heat-budget analysis, it is found that the improved subsurface simulations are attributed to the enhanced cooling effect by constraining the vertical mixing diffusivity in terms of the observational estimate, implying that overestimation of vertical mixing is primarily responsible for the subsurface warm biases in the TA. Since the climate simulation is very sensitive to the vertical mixing parameterization, more accurate representations of ocean vertical mixing are clearly needed in ocean and climate models. Significance Statement: The purpose of our study is to analyze the characteristics of subsurface temperature biases in the tropical Atlantic Ocean and to investigate the causes for the biases. This is important because subsurface biases greatly reduce the reliability of models in climate prediction and projection. It is found that significant subsurface warm biases arise in 100–150 m over the entire tropical Atlantic basin and the biases are mainly attributed to overestimated ocean vertical mixing. Our work highlights that subsurface ocean simulations are highly sensitive to vertical mixing parameterization, and further research is necessary for more accurate representations of ocean vertical mixing in ocean and climate modeling. [ABSTRACT FROM AUTHOR]

Large meridional excursions of a jet stream are conducive to blocking and related midlatitude weather extremes, yet the physical mechanism of jet meandering is not well understood. This paper examines the mechanisms of jet meandering in boreal winter through the lens of a potential vorticity (PV)-like tracer advected by reanalysis winds in an advection–diffusion model. As the geometric structure of the tracer displays a compact relationship with PV in observations and permits a linear mapping from tracer to PV at each latitude, jet meandering can be understood by the geometric structure of tracer field that is only a function of prescribed advecting velocities. This one-way dependence of tracer field on advecting velocities provides a new modeling framework to quantify the effects of time mean flow versus transient eddies on the spatiotemporal variability of jet meandering. It is shown that the mapped tracer wave activity resembles the observed spatial pattern and magnitude of PV wave activity for the winter climatology, interannual variability, and blocking-like wave events. The anomalous increase in tracer wave activity for the composite over interannual variability or blocking-like wave events is attributed to weakened composite mean winds, indicating that the low-frequency winds are the leading factor for the overall distributions of wave activity. It is also found that the tracer model underestimates extreme wave activity, likely due to the lack of feedback mechanisms. The implications for the mechanisms of jet meandering in a changing climate are also discussed. [ABSTRACT FROM AUTHOR]

Recent rapid melting of the summer Greenland ice sheet (GrIS) and its impact on Earth's climate has attracted much attention. In this paper, we establish a connection between the melting of GrIS and the variability of summer sea surface temperature (SST) anomalies over the North Atlantic on interannual to interdecadal time scales through changes in subseasonal Greenland blocking (GB). It is found that the latitude and width of GB are important for the spatial patterns of the GrIS melting. The melting of the GrIS on interdecadal time scales is most prominent for the positive Atlantic multidecadal oscillation phase (AMO+) because the high-latitude GB and its large width, long lifetime, and slow decay are favored. However, the North Atlantic mid-high latitude warm–cold–warm (cold–warm–cold) tripole, referred to as the NAT+ (NAT−) pattern, on interannual time scales tends to strengthen (weaken) the role of AMO+ in the GrIS melting, especially on the northern or northeastern periphery of Greenland, by promoting (inhibiting) high-latitude GB and increasing (decreasing) its width. It is further revealed that AMO+ (NAT+) favors the persistence and width of GB mainly through producing weak summer zonal winds and a small summer meridional potential vorticity gradient (PVy) in the North Atlantic mid-high latitudes at 55°–70°N (55°–65°N) compared to the role of negative AMO (NAT−). The event frequency and zonal width of GB events and their poleward shift are favored by the combination of NAT+ with AMO+. In contrast, the combination of NAT− and AMO+ tends to suppress reduced summer zonal winds and PVy, thus inhibiting the event frequency of GB events and their poleward shift and zonal width. Significance Statement: Rapid melting of the summer Greenland ice sheet (GrIS) has been observed to frequently occur, especially after the year 2000, leading to a rise in sea level and other effects on Earth's climate. The physical cause of the rapid melting of the GrIS is an important area of research. We establish a connection between the summer melting of the GrIS and different sea surface temperature (SST) modes in the North Atlantic via changes in Greenland blocking. Although the positive Atlantic multidecadal oscillation (AMO+) phase favors the overall melting of GrIS, the phase of the North Atlantic tripole (NAT) SST pattern modulates the strength and location of the GrIS melting. The positive NAT phase (NAT+) with a warm–cold–warm tripole structure in the North Atlantic mid-high latitudes and AMO+ combine to result in a strong warm SST anomaly in the high latitudes of the North Atlantic north of 60°N, which promotes the melting of GrIS on the western, northern, and northeastern peripheries of Greenland via high-latitude Greenland blocking with an increased zonal width. The combination of the negative NAT phase (NAT−) with a cold–warm–cold tripole structure and AMO+ tends to suppress this effect. Thus, our results provide a new understanding of why the melting of GrIS shows a strong variability in strength and region. [ABSTRACT FROM AUTHOR]

Interdecadal variations of the relationship between El Niño–Southern Oscillation (ENSO) and the Indo-China Peninsula (ICP) surface air temperature (SAT) in winter are investigated in the study. Generally, there exists a positive correlation between them during 1958–2015 because the ENSO-induced anomalous western North Pacific anticyclone (WNPAC) is conducive to pronounced temperature advection anomalies over the ICP. However, such correlation is unstable in time, having experienced a high-to-low transition around the mid-1970s and a recovery since the early 1990s. This oscillating relationship is owing to the anomalous WNPAC intensity in different decades. During the epoch of high correlation, the anomalous WNPAC and associated southwesterly winds over the ICP are stronger, which brings amounts of warm temperature advection and markedly heats the ICP. In contrast, a weaker WNPAC anomaly and insignificant ICP SAT anomalies are the circumstances for the epoch of low correlation. It is also found that substantial southwesterly wind anomalies over the ICP related to the anomalous WNPAC occur only when large sea surface temperature (SST) anomalies over the northwest Indian Ocean (NWIO) coincide with ENSO (viz., when the ENSO–NWIO SST connection is strong). The NWIO SST anomalies are capable of driving favorable atmospheric circulation that effectively alters ICP SAT and efficiently modulates the ENSO–ICP SAT correlation, which is further supported by numerical simulations utilizing the Community Atmospheric Model, version 4 (CAM4). This paper emphasizes the non-stationarity of the ENSO–ICP SAT relationship and also uncovers the underlying modulation factors, which has important implications for the seasonal prediction of the ICP temperature. [ABSTRACT FROM AUTHOR]

In this paper, the different effects of the eastern equatorial Pacific (EP) and central equatorial Pacific (CP) Ocean El Niño–Southern Oscillation (ENSO) events on interannual variation in the diurnal sea surface temperature (SST) are explored in both the Niño-3 and Niño-4 regions. In the Niño-3 region, the diurnal SST anomaly (DSSTA) is negative during both EP and CP El Niño events and becomes positive during both EP and CP La Niña events. However, the DSSTA in the Niño-4 region is positive in El Niño years and negative in La Niña years, which is opposite to that in the Niño-3 region. Further analysis indicates that the incident shortwave radiation (SWR), wind stress (WS), and upward latent heat flux (LHF) are the main factors causing the different interannual variations in the DSST. In the Niño-3 region, decreased SWR and increased LHF lead to a negative DSSTA in EP El Niño years, and enhanced WS and increased LHF cause a negative DSSTA in CP El Niño years. Conversely, in that same region, increased SWR and decreased LHF lead to a positive DSSTA in EP La Niña years, and reduced WS and decreased LHF cause a positive DSSTA in CP La Niña years. In the Niño-4 region, the reduced trade wind plays a key role in producing the positive DSSTA, whereas the decreased SWR has an opposite effect that reduces the range of the DSSTA during both EP and CP El Niño events, and conversely the enhanced trade wind plays a key role in producing the negative DSSTA, whereas the increased SWR has an opposite effect that increases the range of the DSSTA during both EP and CP La Niña events. [ABSTRACT FROM AUTHOR]

A traditional view is that the ocean outside of the tropics responds passively to atmosphere forcing, which implies that air–sea heat fluxes are mainly driven by atmosphere variability. This paper tests this viewpoint using state-of-the-art air–sea turbulent heat flux observational analyses and a climate model run at different resolutions. It is found that in midlatitude ocean frontal zones the variability of air–sea heat fluxes is not predominantly driven by the atmosphere variations but instead is forced by sea surface temperature (SST) variations arising from intrinsic oceanic variability. Meanwhile in most of the tropics and subtropics wind is the dominant driver of heat flux variability, and atmosphere humidity is mainly important in higher latitudes. The predominance of ocean forcing of heat fluxes found in frontal regions occurs on scales of around 700 km or less. Spatially smoothing the data to larger scales results in the traditional atmosphere-driving case, while filtering to retain only small scales of 5° or less leads to ocean forcing of heat fluxes over most of the globe. All observational analyses examined (1° OAFlux; 0.25° J-OFURO3; 0.25° SeaFlux) show this general behavior. A standard resolution (1°) climate model fails to reproduce the midlatitude, small-scale ocean forcing of heat flux: refining the ocean grid to resolve eddies (0.1°) gives a more realistic representation of ocean forcing but the variability of both SST and of heat flux is too high compared to observational analyses. [ABSTRACT FROM AUTHOR]

In this paper and Part II a comprehensive picture of the annual cycle of the Northern Hemisphere storm tracks is presented and discussed for the first time. It is based on both feature tracking and Eulerian-based diagnostics, applied to vorticity and meridional wind in the upper and lower troposphere. Here, the storm tracks, as diagnosed using both variables and both diagnostic techniques, are presented for the four seasons for each of the two levels. The oceanic storm tracks retain much of their winter mean intensity in spring with only a small change in their latitude. In the summer they are much weaker, particularly in the Pacific and are generally farther poleward. In autumn the intensities are larger again, comparable with those in spring, but the latitude is still nearer to that of summer. However, in the lower troposphere in the eastern ocean basins the tracking metrics show northern and southern tracks that change little with latitude through the year. The Pacific midwinter minimum is seen in upper-troposphere standard deviation diagnostics, but a richer picture is obtained using tracking. In winter there are high intensities over a wide range of latitudes in the central and eastern Pacific, and the western Pacific has high track density but weak intensity. In the lower troposphere all the diagnostics show that the strength of the Pacific and Atlantic storm tracks are generally quite uniform over the autumn–winter–spring period. There is a close relationship between the upper-tropospheric storm track, particularly that based on vorticity, and tropopause-level winds and temperature gradients. In the lower troposphere, in winter the oceanic storm tracks are in the region of the strong meridional SST gradients, but in summer they are located in regions of small or even reversed SST gradients. However, over North America the lower-tropospheric baroclinicity and the upstream portion of the Atlantic storm track stay together throughout the year. [ABSTRACT FROM AUTHOR]

A frequency–wavenumber power () spectrum was constructed using satellite-derived outgoing longwave radiation (OLR) and brightness temperature for the tropical latitudes. Since the two datasets overlap for over 34 years with nonintersecting sources in error and compare relatively well with each other, it is possible to diagnose trends in the tropical wave activity from the two datasets with confidence. The results suggest a weakening trend in characterized by high interannual variability for wave activity occurring in the low-frequency part of the spectrum and a steady increase in with relatively low interannual variability for wave activity occurring in the high-frequency part of the spectrum. The results show the parts of the spectrum representing the Madden–Julian oscillation and equatorial Rossby wave losing and other parts of the spectrum representing Kelvin waves, mixed Rossby–gravity waves, and tropical disturbance–like wave activity gaining . Similar results were obtained when trends in variance corresponding to the first principal component were produced using spectrally filtered OLR data representative of atmospheric equatorial waves. Spatial trends in the active phase of wave events and the mean duration of events are also shown for the different wave types. Linear trends in for the entire spectrum and regional means in the spectrum corresponding to the abovementioned five wave types with confidence intervals are also presented in the paper. Finally, we demonstrate that El Niño–Southern Oscillation variability does not appear to control the overall spatial patterns and trends observed in the spectrum. [ABSTRACT FROM AUTHOR]

The path of the Gulf Stream as it leaves the continental shelf near Cape Hatteras is marked by a sharp gradient in ocean temperature known as the North Wall. Previous work in the literature has considered processes related to the North Atlantic Oscillation (NAO) in triggering latitudinal displacements of the North Wall position. This paper presents evidence that the Atlantic meridional mode (AMM) also impacts interannual variations of the North Wall position. The AMM signal from the tropics propagates to the Gulf Stream near the 200-m depth, and there are two time scales for this interaction. Anomalous Ekman suction induced by AMM cools the tropical Atlantic. The cold water in the Caribbean Sea is entrained into the currents feeding the Gulf Stream, and this cooling signal reaches the North Wall within a year. A second mechanism involves cold anomalies in the western tropical Atlantic, which initially propagate westward as baroclinic planetary waves, reaching the Gulf Stream and resulting in a southward shift in the North Wall position after a delay of about one year. In an analysis for the period 1961–2015, AMM's signal dominates North Wall fluctuations in the upper 300 m, while NAO is the major influence below ~500 m; the influence of both the teleconnections is seen between 300 and 500 m. The relationship between the Atlantic meridional overturning circulation (AMOC) and the North Wall is investigated for the 2005–15 period and found to be statistically significant only at the sea surface in one of the three North Wall indices used. [ABSTRACT FROM AUTHOR]

Global warming projections of dynamics are less robust than projections of thermodynamics. However, robust aspects of the thermodynamics can be used to constrain some dynamical aspects. This paper argues that tropospheric expansion under global warming (a thermodynamical process) explains changes in the amplitude of the quasi-biennial oscillation (QBO) in the lower and middle stratosphere (a dynamical process). A theoretical scaling for tropospheric expansion of approximately 6 hPa K−1 is derived, which agrees well with global climate model (GCM) experiments. Using this theoretical scaling, the response of QBO amplitude to global warming is predicted by shifting the climatological QBO amplitude profile upward by 6 hPa per kelvin of global warming. In global warming simulations, QBO amplitude in the lower to middle stratosphere shifts upward as predicted by tropospheric expansion. Applied to observations, the tropospheric expansion framework suggests a historical weakening of QBO amplitude at 70 hPa of 3% decade−1 from 1953 to 2020. This expected weakening trend is half of the 6% decade−1 from 1953 to 2012 detected and attributed to global warming in a recent study. The previously reported trend was reinforced by record low QBO amplitudes during the mid-2000s, from which the QBO has since recovered. Given the modest weakening expected on physical grounds, past decadal modulations of QBO amplitude are reinterpreted as a hitherto unrecognized source of internal variability. This large internal variability dominates over the global warming signal such that, despite 65 years of observations, there is not yet a statistically significant weakening trend. [ABSTRACT FROM AUTHOR]

This paper studies the influences of the Madden–Julian oscillation (MJO) on East Asian (EA) winter rainfall using the singular value decomposition (SVD) approach. This method uses two-dimensional instead of latitudinally averaged variables in the commonly used real-time multivariate MJO (RMM) index. A comparison of the two approaches is made using the same OLR and zonal wind data over 37 boreal winter seasons of December–March. The SVD composite reveals a more conspicuous and coherent variation throughout the MJO cycle, while the RMM composite is more ambiguous. In particular, the SVD analysis identifies the convection anomalies over the Maritime Continent and the subtropical western Pacific (MCWP) as a major cause of enhanced rainfall in EA at RMM phases 8 and 1. This is at least one-eighth of a cycle earlier than the phases of convection development over the Indian Ocean (IO) that were emphasized by previous studies. A linearized global baroclinic model is used to demonstrate the mechanism of MJO forcing on EA rainfall during various phases, with a focus on the MCWP cooling. The result shows that the anomalous MCWP cooling and the resultant low-level anticyclonic flow interact with the East Asian jet, leading to an overall weakened EA winter monsoon circulation. The associated anomalous overturning circulation, with ascending motion and low-level horizontal moisture convergence in EA, contributes to the enhanced rainfall. This model result supports the interpretation of the SVD analysis, in that the MCWP cooling induced anomalous meridional circulation is a more direct cause of enhanced EA rainfall than the IO heating (or the IO–MCWP heating dipole) induced Rossby wave teleconnection. [ABSTRACT FROM AUTHOR]

In this paper, we investigate the relationship between upper ocean heat content (OHC) and El Niño–Southern Oscillation (ENSO) sea surface temperature (SST) anomalies mainly using the neutral recharge oscillator (NRO) model both analytically and numerically. Previous studies showed that spring OHC, which leads SST by 6–12 months, represents a major source of predictability for ENSO. It is suggested that this seasonality is caused by the seasonally varying growth rate in SST anomalies. Moreover, a shortened ENSO period will lead to a reduced SST predictability from OHC, with the most significant decrease occurring in the latter half of the calendar year. The cross-correlation relationship between OHC and ENSO SST anomalies is further identified in the damped and self-excited version of the recharge oscillator model. Finally, we suggest that the seasonal growth rate of ENSO anomalies is the cause of the seasonality in the effectiveness of OHC as a predictor in ENSO forecasting, particularly as it relates to the boreal spring persistence barrier and associated spring predictability barrier. We also explain the shorter lead time between spring OHC and ENSO SST anomalies after the turn of the twenty-first century in terms of the apparent higher frequency of the ENSO period. [ABSTRACT FROM AUTHOR]

This paper introduces a set of descriptors applied to weather regimes that allow for a detailed monitoring of the location and intensity of their atmospheric centers of action (e.g., troughs and ridges) and the gradients between them, when applicable. Descriptors are designed to document the effect of climate variability and change in modulating the character of daily weather regimes, rather than merely their occurrence statistics. As a case study, the methodology is applied to Aotearoa New Zealand (ANZ), using ERA5 ensemble reanalysis data for the period 1979-2019. Here, we analyze teleconnections between the regimes and their descriptors and large-scale climate variability. Results show a significant modulation of centers of action by the phase of the southern annular mode, with a strong relationship identified with the latitude of atmospheric ridges. Significant associations with El Niño-Southern Oscillation are also identified. Modes of large-scale variability have a stronger influence on the regimes' intrinsic features than their occurrence. This demonstrates the usefulness of such descriptors, which help explain the relationship betweenmidlatitude transient perturbations and large-scale modes of climate variability. In future research, this methodological framework will be applied to analyze (i) low-frequency changes in weather regimes under climate change, in line with the southward shift of storm tracks, and (ii) regional-scale effects on the climate of ANZ, resulting from interaction with its topography. [ABSTRACT FROM AUTHOR]

Previous studies have recognized the societal relevance of climatic extremes on the seasonal time scale and examined physical processes leading to individual high-impact extreme seasons (e.g., extremely wet orwarm seasons). However, these findings have not yet been generalized beyond case studies since at any specific location only very few seasonal events of such rarity occurred in the observational record. In this concept paper, a pragmatic approach to pool seasonal extremes across space is developed and applied to investigate hot summers and cold winters in ERA-Interim and the Community Earth System Model Large Ensemble (CESM-LENS). We identify spatial extreme season objects as contiguous regions of extreme seasonal mean temperatures based on statistical modeling. Regional pooling of extreme season objects in CESM-LENS then yields considerable samples of analogs to even themost extreme ERA-Interimevents. This approach offers numerous opportunities for systematically analyzing large samples of extreme seasons, and several such analyses are illustrated. We reveal a striking co-occurrence of ElNiño to La Niña transitions and the largest ERA-Interim midlatitude extreme summer events. Moreover, we perform a climate model evaluation with regard to extreme season size and intensitymeasures and estimate how often an extreme winter like the coldNorth American 2013/14 winter is expected anywhere inmidlatitude regions. Furthermore, we present a large set of simulated analogs to this event, whichmakes it possible to studycommonalities and differences of their underlying physical processes. Finally, substantial but spatially varying climatological differences in the size of extreme summer and extreme winter objects are identified. [ABSTRACT FROM AUTHOR]

A number of studies have suggested that heat waves will increase in frequency, intensity, and duration under anthropogenic global warming. However, changes to heat extremes and temperature distributions in general are less understood in regions of complex terrain and substantial land-water contrasts. Such surface variations are important in the northwest United States, where synoptic- and mesoscale circulations interact with local topography and land-water interfaces to produce heat waves that have substantial impacts on fire weather, air quality, wind energy, and the population in general. Thus it is crucial to identify how the synoptic- and mesoscale circulations that drive such local extremes will change in a warming world. This paper analyzes changes to the conditions associated with heat extremes over the northwest United States by utilizing global and regional climate models. A maximum in 700-hPa warming is projected over the northwest United States and southwest Canada, with low-level warming attenuated near the coast. Soil moisture declines are projected over the region, which further enhances future extremes. It is found that low-level zonal wind distributions over the northwest United States become narrower, leading to fewer days with strong offshore flow. This is important since extreme warming events are associated with the strongest offshore/downslope flow. Historical and future regional temperature distributions are described and it is shown that although CMIP5 models predict increases in heat extremes for west Oregon and Washington, these increases are less than for inland areas. [ABSTRACT FROM AUTHOR]

The impact of two different coupled cirrus microphysics-radiation parameterizations on the zonally averaged temperature and humidity biases in the tropical tropopause layer (TTL) of a Met Office climate model configuration is assessed. One parameterization is based on a linear coupling between a model prognostic variable, the ice mass mixing ratio q i, and the integral optical properties. The second is based on the integral optical properties being parameterized as functions of q i and temperature, T c, where the mass coefficients (i.e., scattering and extinction) are parameterized as nonlinear functions of the ratio between q i and T c. The cirrus microphysics parameterization is based on a moment estimation parameterization of the particle size distribution (PSD), which relates the mass moment (i.e., second moment if mass is proportional to size raised to the power of 2) of the PSD to all other PSD moments through the magnitude of the second moment and T c. This same microphysics PSD parameterization is applied to calculate the integral optical properties used in both radiation parameterizations and, thus, ensures PSD and mass consistency between the cirrus microphysics and radiation schemes. In this paper, the temperature-non-dependent and temperature-dependent parameterizations are shown to increase and decrease the zonally averaged temperature biases in the TTL by about 1 K, respectively. The temperature-dependent radiation parameterization is further demonstrated to have a positive impact on the specific humidity biases in the TTL, as well as decreasing the shortwave and longwave biases in the cloudy radiative effect. The temperature-dependent radiation parameterization is shown to be more consistent with TTL and global radiation observations. [ABSTRACT FROM AUTHOR]

This paper shows that joint temperature-precipitation information over a global domain provides a more accurate estimate of aerosol forced responses in climate models than does any other combination of temperature, precipitation, or sea level pressure. This fact is demonstrated using a new quantity called potential detectability, which measures the extent to which a forced response can be detected in a model. In particular, this measure can be evaluated independently of observations and therefore permits efficient exploration of a large number of variable combinations before performing optimal fingerprinting on observations. This paper also shows that the response to anthropogenic aerosol forcing can be separated from that of other forcings using only spatial structure alone, leaving the time variation of the response to be inferred from data, thereby demonstrating that temporal information is not necessary for detection. The spatial structure of the forced response is derived by maximizing the signal-to-noise ratio. For single variables, the north-south hemispheric gradient and equator-to-pole latitudinal gradient are important spatial structures for detecting anthropogenic aerosols in some models but not all. Sea level pressure is not an independent detection variable because it is derived partly from surface temperature. In no case does sea level pressure significantly enhance potential detectability beyond that already possible using surface temperature. Including seasonal or land-sea contrast information does not significantly enhance detectability of anthropogenic aerosol responses relative to annual means over global domains. [ABSTRACT FROM AUTHOR]

A correction to the article "Toward consistent diagnostics of the coupled atmosphere and ocean energy budgets," by Michael Mayer and colleagues, that was published in a 2017 issue is presented.

The author talks about the discontinuation of Expedited Contributions (ECs) article type by the American Meteorological Society (AMS). Topics discussed include the aim of the ECs to provide short and concise contributions to the journals of AMS, the use of digital object identifiers (DOIs) in the online articles, and the reduced time of submission of research papers for publication.

In this paper we develop a method to quantify the accuracy of different pattern extraction techniques for the additive space–time modes often assumed to be present in climate data. It has previously been shown that the standard technique of principal component analysis (PCA; also known as empirical orthogonal functions) may extract patterns that are not physically meaningful. Here we analyze two modern pattern extraction methods, namely dynamical mode decomposition (DMD) and slow feature analysis (SFA), in comparison with PCA. We develop a Monte Carlo method to generate synthetic additive modes that mimic the properties of climate modes described in the literature. The datasets composed of these generated modes do not satisfy the assumptions of any pattern extraction method presented. We find that both alternative methods significantly outperform PCA in extracting local and global modes in the synthetic data. These techniques had a higher mean accuracy across modes in 60 out of 60 mixed synthetic climates, with SFA slightly outperforming DMD. We show that in the majority of simple cases PCA extracts modes that are not significantly better than a random guess. Finally, when applied to real climate data these alternative techniques extract a more coherent and less noisy global warming signal, as well as an El Niño signal with a clearer spectral peak in the time series, and more a physically plausible spatial pattern. [ABSTRACT FROM AUTHOR]

Observations of warm past climates and projections of future climate change show that the Arctic warms more than the global mean, particularly during winter months. Previous work has attributed this reduced Arctic land seasonality to the effects of sea ice or clouds. In this paper, we show that the reduced Arctic land seasonality is a robust consequence of the relatively small surface heat capacity of land and the nonlinearity of the temperature dependence of surface longwave emission, without recourse to other processes or feedbacks. We use a general circulation model (GCM) with no clouds or sea ice and a simple representation of land. In the annual mean, the equator-to-pole surface temperature gradient falls with increasing CO2, but this is only a near-surface phenomenon and is not caused by the change in total meridional heat transport, which is virtually unaltered. The high-latitude land has about twice as much warming in winter than in summer, whereas high-latitude ocean has very little seasonality in warming. A surface energy balance model shows how the combination of the smaller surface heat capacity of land and the nonlinearity of the temperature dependence of surface longwave emission gives rise to the reduced seasonality of the land surface. The increase in evaporation over land also leads to winter amplification of warming over land, although amplification still occurs without it. While changes in clouds, sea ice, and ocean heat transport undoubtedly play a role in high-latitude warming, these results show that enhanced land surface temperature warming in winter can happen in their absence for robust reasons. [ABSTRACT FROM AUTHOR]

In this paper, we investigate the potential factors that control the relationship between the El Niño–Southern Oscillation (ENSO) persistence barriers (PBs) in sea surface temperature (SST) and ocean heat content (OHC) and apply them to explain observational ENSO PBs. With the addition of seasonal growth rate in SST in the neutral recharge oscillator (NRO) model, approximate analytical solutions of autocorrelation functions for SST and OHC suggest strictly that the timing of PB for OHC leads that of SST by half a year and the strengths of the two PBs are the same. The numerical solutions of the NRO model also show a similar relationship. The role of ENSO growth rate with regard to PBs in SST and OHC is then identified in the damped and unstable ENSO regime. Therefore, it is suggested that for the observational ENSO, the seasonally varying ENSO growth rate in SST controls PBs in SST and OHC simultaneously. [ABSTRACT FROM AUTHOR]

The subpolar North Atlantic (SPNA) experienced extreme cold during 2015, an event often called the "cold blob." The evolution of this event in the Community Earth System Model version 1 Decadal Prediction Large Ensemble (CESM1-DPLE) hindcast initialized in November 2014 is compared to observations. This CESM1-DPLE hindcast failed to predict cold conditions during 2015 despite already cold SPNA initial conditions and despite having high sea surface temperature skill in the SPNA in all other years. The goal of this paper is to understand what led to this prediction failure in order to provide insight for future decadal prediction efforts. Our analysis shows that strongly positive North Atlantic Oscillation (NAO) conditions during winter and spring 2015 likely sustained the cold blob but were not simulated in any CESM1-DPLE members. We examine the rarity of the 2015 event using the CESM1-DPLE's uninitialized counterpart, the CESM1 Large Ensemble (CESM1-LE). Results from the CESM1-LE indicate that the exceptional state of the observed NAO in the winter of 2015 is at least part of the explanation for why this event was not encompassed in the CESM1-DPLE spread. To test another possibility—namely, that deficiencies in the initial conditions degraded the prediction—we performed additional hindcasts using the CESM1-DPLE protocol but different initial conditions. Altering the initial conditions did not improve the simulation of the 2015 cold blob, and in some cases, degraded it. Given the difficulty of predicting this event, this case could be a useful test bed for future prediction system development. [ABSTRACT FROM AUTHOR]

The reliability of general circulation models (GCMs) is commonly associated with their ability to reproduce relevant aspects of observed climate, and thus the evaluation of GCM performance has become a standard practice for climate change studies. As such, there is an ever-growing literature that focuses on developing and evaluating metrics to assess GCM performance. In this paper it is shown that some commonly applied metrics provide little information for discriminating GCMs based on their performance, once uncertainty is included. A new methodology is proposed that differs from common approaches in that it focuses on evaluating GCMs' abilities to reproduce the observed response of surface temperature to changes in external radiative forcing (RF), while controlling for observed and simulated variability. It uses formal statistical tests to evaluate two aspects of the warming trend that are central for climate change studies: 1) if the response to RF produced by a particular GCM is compatible with observations and 2) if the magnitudes of the observed and simulated rates of warming are statistically similar. We illustrate the proposed methodology by evaluating the ability of 21 GCMs to reproduce the observed warming trend at the global scale and for eight subcontinental land domains. Results show that most of the GCMs provide an adequate representation of the observed warming trend for the global scale and for domains located in the Southern Hemisphere. However, GCMs tend to overestimate the warming rate for domains in the Northern Hemisphere, particularly since the mid-1990s. [ABSTRACT FROM AUTHOR]

Studies detecting trends in climate elements typically concentrate on their local significance, ignoring the question of whether the significant local trends may or may not have occurred as a result of chance. This paper fills this gap by examining several approaches to detecting statistical significance of trends defined on a grid (i.e., on a regional scale). To this end, we introduce a novel simple procedure of significance testing that is based on counting signs of local trends (sign test), and we compare it with five other approaches to testing collective significance of trends: counting, extended Mann–Kendall, Walker, false detection rate (FDR), and regression tests. Synthetic data are used to construct null distributions of trend statistics, to determine critical values of the tests, and to assess the performance of tests in terms of type-II error. For lower values of spatial and temporal autocorrelations, the sign test and extended Mann–Kendall test perform slightly better than the counting test; these three tests outperform the Walker, FDR, and regression tests by a wide margin. For high autocorrelations, which is a more realistic case, all tests become similar in their performance, with the exception of the regression test, which performs somewhat worse. Some tests cannot be used under specific conditions because of their construction: the Walker and FDR tests for high temporal autocorrelations, and the sign test under high spatial autocorrelations. [ABSTRACT FROM AUTHOR]

The present paper is the second part of a two-part study on empirical modeling and prediction of climate variability. This paper deals with spatially distributed data, as opposed to the univariate data of Part I. The choice of a basis for effective data compression becomes of the essence. In many applications, it is the set of spatial empirical orthogonal functions that provides the uncorrelated time series of principal components (PCs) used in the learning set. In this paper, the basis of the learning set is obtained instead by applying multichannel singular-spectrum analysis to climatic time series and using the leading spatiotemporal PCs to construct a reduced stochastic model. The effectiveness of this approach is illustrated by predicting the behavior of the Jin-Neelin-Ghil (JNG) hybrid seasonally forced coupled ocean-atmosphere model of El Niño-Southern Oscillation. The JNG model produces spatially distributed and weakly nonstationary time series to which the model reduction and prediction methodology is applied. Critical transitions in the hybrid periodically forced coupled model are successfully predicted on time scales that are substantially longer than the duration of the learning sample. [ABSTRACT FROM AUTHOR]

Cape Farewell, Greenland's southernmost point, is characterized by a number of low-level jets that are the result of topographic flow distortion associated with passing extratropical cyclones. The heavy seas associated with these wind events are a hazard to maritime traffic in the region. In addition, the air-sea heat flux associated with these weather systems plays an important role in the climate system by contributing to the forcing of the lower limb of the Atlantic meridional overturning circulation. In this paper, the North American Regional Reanalysis will be used to generate a higher-resolution climatology of these mesoscale jets as compared to previous studies. Through the use of a diagnostic that partitions the occurrence frequency of high-speed wind events by wind direction, the author shows that there are four different types of Cape Farewell tips jets that are characterized as having either northwesterly, southwesterly, northeasterly, or southeasterly wind direction. All four types have distinct regions in the vicinity of Cape Farewell where their respective occurrence frequencies and air-sea heat fluxes are at a maximum. The southwesterly and northeasterly jets closely resemble the wind systems previously identified as being westerly and easterly tip jets. There are also instances where one type evolves into another and so it is possible to view westerly tip jets as a continuum with the northwesterly and southwesterly events identified in this paper representing end members with a similar picture for easterly tip jets. The position of a particular event along these continua will determine its impact on local weather and the coupled climate system. [ABSTRACT FROM AUTHOR]

Future projections of precipitation at regional scales are vital to inform climate change adaptation activities. Therefore, is it important to quantify projected changes and associated uncertainty, and understand model processes responsible. This paper addresses these challenges for southern Africa and the adjacent Indian Ocean focusing on the local wet season. Precipitation projections for the end of the twenty-first century indicate a pronounced dipole pattern in the CMIP5 multimodel mean. The dipole indicates future wetting (drying) to the north (south) of the climatological axis of maximum rainfall, implying a northward shift of the ITCZ and south Indian Ocean convergence zone that is not consistent with a simple "wet get wetter" pattern. This pattern is most pronounced in early austral summer, suggesting a later and shorter wet season over much of southern Africa. Using a decomposition method we determine physical mechanisms underlying this dipole pattern of projected change, and the associated intermodel uncertainty. The projected dipole pattern is largely associated with the dynamical component of change indicative of shifts in the location of convection. Over the Indian Ocean, this apparent northward shift in the ITCZ may reflect the response to changes in the north-south SST gradient over the Indian Ocean, consistent with a "warmest get wetter" mechanism. Over land subtropical drying is relatively robust, particularly in the early wet season. This has contributions from dynamical shifts in the location of convection, which may be related to regional SST structures in the southern Indian Ocean, and the thermodynamic decline in relative humidity. Implications for understanding and potentially constraining uncertainty in projections are discussed. [ABSTRACT FROM AUTHOR]

Although El Niño-Southern Oscillation (ENSO) influences the Indian summer monsoon, its impact on moisture transport toward the southern Tibetan Plateau (TP) remains poorly understood. Precipitation stable isotopes are useful indices for climate change in the TP. Classical interpretations of variations of precipitation stable isotopes focus on the local surface air temperature or precipitation amount. However, several of the latest studies suggested they may correlate with large-scale modes of variability, such as ENSO. This paper presents a detailed study of ENSO's effect on annual variations of the oxygen stable isotopic composition of precipitation (δ18Op) at Lhasa in the southern TP for up to 10 years. The stable isotopic composition of water vapor from satellite data [Tropospheric Emission Spectrometer (TES)] and simulations from an isotopically enabled atmospheric general circulation model (zoomed LMDZiso) are used to explore the mechanism that leads to variations of δ18Op at Lhasa. Statistically significant correlations between δ18Op and ENSO indices [Southern Oscillation index (SOI) and Niño-3.4 sea surface temperature index (Niño-3.4)] are observed. This paper shows that ENSO's effects on the location and intensity of convection over the Arabian Sea, the Bay of Bengal, and the tropical Indian Ocean, along moisture transport paths toward Lhasa, further impact convection from the northern Tibetan Plateau. The changing of this convection results in lower δ18Op at Lhasa in 2007, a La Niña year, and higher δ18Op in 2006, an El Niño year. The study presented here confirms that the regional upstream convection related to ENSO teleconnections plays an important role in variations of δ18Op at the interannual scale and that the more depleted oxygen stable isotopic composition of vapor (δ18Oυ) from the northwestern region of India during a La Niña year intensifies the lower δ18Op at Lhasa in a La Niña year. The study's results have implications for the interpretation of past variations of archives with precipitation stable isotopes, such as ice cores, tree rings, lake sediments, and speleothems, in this region. [ABSTRACT FROM AUTHOR]

A high-resolution regional reanalysis of the Indian Monsoon Data Assimilation and Analysis (IMDAA) project is made available to researchers for deeper understanding of the Indian monsoon and its variability. This 12-km resolution reanalysis covering the satellite era from 1979 to 2018 using a 4D-Var data assimilation method and the U.K. Met Office Unified Model is presently the highest resolution atmospheric reanalysis carried out for the Indian monsoon region. Conventional and satellite observations from different sources are used, including Indian surface and upper air observations, of which some had not been used in any previous reanalyses. Various aspects of this reanalysis, including quality control and bias correction of observations, data assimilation system, land surface analysis, and verification of reanalysis products, are presented in this paper. Representation of important weather phenomena of each season over India in the IMDAA reanalysis verifies reasonably well against India Meteorological Department (IMD) observations and compares closely with ERA5. Salient features of the Indian summer monsoon are found to be well represented in the IMDAA reanalysis. Characteristics of major semipermanent summer monsoon features (e.g., low-level jet and tropical easterly jet) in IMDAA reanalysis are consistent with ERA5. The IMDAA reanalysis has captured the mean, interannual, and intraseasonal variability of summer monsoon rainfall fairly well. IMDAA produces a slightly cooler winter and a hotter summer than the observations; the reverse is true for ERA5. IMDAA captured the fine-scale features associated with a notable heavy rainfall episode over complex terrain. In this study, the fine grid spacing nature of IMDAA is compromised due to the lack of comparable resolution observations for verification. [ABSTRACT FROM AUTHOR]

This paper attempts to enhance our understanding of the causes of Atlantic multidecadal variability (AMV). Following the literature, we define the AMV as the SST averaged over the North Atlantic basin, linearly detrended and low-pass filtered. There is an ongoing debate about the drivers of the AMV, which include internal variability generated from the ocean or atmosphere (or both) and external radiative forcing. We test the role of these factors in explaining the time history, variance, and spatial pattern of the AMV using a 41-member ensemble from a fully coupled version of CESM and a 10-member ensemble of the CESM atmosphere coupled to a slab ocean. The large ensemble allows us to isolate the role of external forcing versus internal variability, and the model differences allow us to isolate the role of coupled ocean circulation. Both with and without coupled ocean circulation, external forcing explains more than half of the variance of the observed AMV time series, indicating its important role in simulating the twentieth-century AMV phases. In this model the net effect of ocean processes is to reduce the variance of the AMV. Dynamical ocean coupling also reduces the ability of the model to simulate the characteristic spatial pattern of the AMV, but forcing has little impact on the pattern. Historical forcing improves the time history and variance of the AMV simulation, while the more realistic ocean representation reduces the variance below that observed and lowers the correlation with observations. [ABSTRACT FROM AUTHOR]

A recent paper by Kato and Rose reports a negative correlation between the annual mean entropy production rate of the climate and the absorption of solar radiation in the CERES SYN1deg dataset, using the simplifying assumption that the system is steady in time. It is shown here, however, that when the nonsteady interannual storage of entropy is accounted for, the dataset instead implies a positive correlation; that is, global entropy production rates increase with solar absorption. Furthermore, this increase is consistent with the response demonstrated by an energy balance model and a radiative–convective model. To motivate this updated analysis, a detailed discussion of the conceptual relationship between entropy production, entropy storage, and entropy flows is provided. The storage-corrected estimate for the mean global rate of entropy production in the CERES dataset from all irreversible transfer processes is 81.9 mW m−2 K−1 and from only nonradiative processes is 55.2 mW m−2 K−1 (observations from March 2000 to February 2018). [ABSTRACT FROM AUTHOR]

El Niño–Southern Oscillation (ENSO) diversity has a significant impact on global climate and seasonal prediction. However, it is still a challenging problem for present-day global climate models to simulate different types of ENSO events with realistic features simultaneously. In this paper, a tropical stochastic skeleton model for the interactions among wind bursts and the Madden–Julian oscillation (MJO), the El Niño, and the Walker circulation is developed to reproduce both dynamical and statistical features of the ENSO diversity. In this model, the intraseasonal component with state-dependent noise captures general features of wind bursts and the MJO, both of which play important roles in triggering the El Niño. The thermocline feedback is the dominant mechanism for generating the eastern Pacific (EP) El Niño, while a nonlinear zonal advection is incorporated into the model that contributes to the central Pacific (CP) El Niño. Besides, a simple but effective stochastic process describing the multidecadal variation of the background Walker circulation modulates the spatial patterns and occurrence frequency of the EP and CP El Niño. This model succeeds in simulating the quasi-regular moderate EP El Niño, the super El Niño, and the CP El Niño as well as the La Niña simultaneously. It also captures the observed non-Gaussian characteristics of sea surface temperature anomalies in different Niño regions. Individual case studies highlight the outstanding skill of the model in reproducing the observed El Niño episodes and their underlying mechanisms. [ABSTRACT FROM AUTHOR]