Search

The supersaturation equation for a vertically moving adiabatic cloud parcel is analyzed. The effects of turbulent updrafts are incorporated in the shape of a stochastic Lagrangian model, with spatial and time correlations expressed in terms of turbulent kinetic energy. Using the Fokker–Planck equation, the steady-state probability distributions of supersaturation are analytically computed for a number of approximations involving the time-scale separation between updraft fluctuations and phase relaxation, and droplet or ice particle size fluctuations. While the analytical results are presented in general for single-phase clouds, the calculated distributions are used to compute mixed-phase cloud properties—mixed fraction and mean liquid water content in an initially icy cloud—and are argued to be useful for generalizing and constructing new parameterization schemes. Significance Statement: Supersaturation is the fuel for the development of clouds in the atmosphere. In this paper, our goal is to better understand the supersaturation budget of clouds embedded in a turbulent environment by analyzing the basic equations of cloud microphysics. It is found that the turbulent characteristics of an air parcel substantially affect the cloud's supersaturation budget and hence its life cycle. This is also shown in the context of mixed-phase clouds where, depending on the turbulent regime, different liquid-to-ice ratios are found. Consequently, the theoretical approach of this paper is crucial to develop tools to parameterize small-scale atmospheric features, like clouds, into global circulation models to improve climate projections for the future. [ABSTRACT FROM AUTHOR]

The potential role of the New England seamount chain (NESC) on the Gulf Stream pathway and variability has been long recognized, and the series of numerical experiments presented in this paper further emphasize the importance of properly resolving the NESC when modeling the Gulf Stream. The NESC has a strong impact on the Gulf Stream pathway and variability, as demonstrated by comparison experiments with and without the NESC. With the NESC removed from the model bathymetry, the Gulf Stream remains a stable coherent jet much farther east than in the experiment with the NESC. The NESC is the leading factor destabilizing the Gulf Stream and, when it is not properly resolved by the model's grid, its impact on the Gulf Stream's pathway and variability is surprisingly large. A high-resolution bathymetry, which better resolves the New England seamounts (i.e., narrower and rising higher in the water column), leads to a tighter Gulf Stream mean path that better agrees with the observed path and a sea surface height variability distribution that is in excellent agreement with the observations. [ABSTRACT FROM AUTHOR]

The Indo-Pacific Ocean appears exponentially stratified between 1- and 3-km depth with a decay scale on the order of 1 km. In his celebrated paper "Abyssal recipes," W. Munk proposed a theoretical explanation of these observations by suggesting a pointwise buoyancy balance between the upwelling of cold water and the downward diffusion of heat. Assuming a constant upwelling velocity w and turbulent diffusivity κ, the model yields an exponential stratification whose decay scale is consistent with observations if κ ∼ 10−4 m2 s−1. Over time, much effort has been made to reconcile Munk's ideas with evidence of vertical variability in κ, but comparably little emphasis has been placed on the even stronger evidence that w decays toward the surface. In particular, the basin-averaged w nearly vanishes at 1-km depth in the Indo-Pacific. In light of this evidence, we consider a variable-coefficient, basin-averaged analog of Munk's budget, which we verify against a hierarchy of numerical models ranging from an idealized basin-and-channel configuration to a coarse global ocean simulation. Study of the budget reveals that the decay of basin-averaged w requires a concurrent decay in basin-averaged κ to produce an exponential-like stratification. As such, the frequently cited value of 10−4 m2 s−1 is representative only of the bottom of the middepths, whereas κ must be much smaller above. The decay of mixing in the vertical is as important to the stratification as its magnitude. Significance Statement: Using a combination of theory and numerical simulations, it is argued that the observed magnitude and shape of the global ocean stratification and overturning circulation appear to demand that turbulent mixing increases quasi-exponentially toward the ocean bottom. Climate models must therefore prescribe such a vertical profile of turbulent mixing in order to properly represent the heat and carbon uptake accomplished by the global overturning circulation on centennial and longer time scales. [ABSTRACT FROM AUTHOR]

This paper describes the downscaling of an ensemble of 12 general circulation models (GCMs) using the Weather Research and Forecasting (WRF) Model at 12-km grid spacing over the period 1970–2099, examining the mesoscale impacts of global warming as well as the uncertainties in its mesoscale expression. The RCP8.5 emissions scenario was used to drive both global and regional climate models. The regional climate modeling system reduced bias and improved realism for a historical period, in contrast to substantial errors for the GCM simulations driven by lack of resolution. The regional climate ensemble indicated several mesoscale responses to global warming that were not apparent in the global model simulations, such as enhanced continental interior warming during both winter and summer as well as increasing winter precipitation trends over the windward slopes of regional terrain, with declining trends to the lee of major barriers. During summer there is general drying, except to the east of the Cascades. The 1 April snowpack declines are large over the lower-to-middle slopes of regional terrain, with small snowpack increases over the lower elevations of the interior. Snow-albedo feedbacks are very different between GCM and RCM projections, with the GCMs producing large, unphysical areas of snowpack loss and enhanced warming. Daily average winds change little under global warming, but maximum easterly winds decline modestly, driven by a preferential sea level pressure decline over the continental interior. Although temperatures warm continuously over the domain after approximately 2010, with slight acceleration over time, occurrences of temperature extremes increase rapidly during the second half of the twenty-first century. Significance Statement: This paper provides a unique high-resolution view of projected climate change over the Pacific Northwest and does so using an ensemble of regional climate models, affording a look at the uncertainties in local impacts of global warming. The paper examines regional meteorological processes influenced by global warming and provides guidance for adaptation and preparation. [ABSTRACT FROM AUTHOR]

Large-amplitude [±100 Sv (1 Sv ≡ 106 m3 s−1)], high-frequency oscillations in the Pacific Ocean's meridional overturning circulation within 10° of the equator have been found in integrations of the NEMO ocean general circulation model. Part I of this paper showed that these oscillations are dominated by two bands of frequencies with periods close to 4 and 10 days and that they are driven by the winds within about 10° of the equator. This part shows that the oscillations can be well simulated by small-amplitude, wind-driven motions on a horizontally uniform, stably stratified state of rest. Its main novelty is that, by focusing on the zonally integrated linearized equations, it presents solutions for the motions in a basin with sloping side boundaries. The solutions are found using vertical normal modes and equatorial meridional modes representing Yanai and inertia–gravity waves. Simulations of 16-day-long segments of the time series for the Pacific of each of the first three meridional and vertical modes (nine modes in all) capture between 85% and 95% of the variance of matching time series segments diagnosed from the NEMO integrations. The best agreement is obtained by driving the solutions with the full wind forcing and the full pressure forces on the bathymetry. Similar results are obtained for the corresponding modes in the Atlantic and Indian Oceans. Slower variations in the same meridional and vertical modes of the MOC are also shown to be well simulated by a quasi-stationary solution driven by zonal wind and pressure forces. [ABSTRACT FROM AUTHOR]

Modern general circulation models continue to require parameterizations of subgrid transport due to planetary boundary layer (PBL) turbulence and convection. Some schemes that unify these processes rely on assumed joint probability distributions of vertical velocity and moist conserved thermodynamic variables to predict the subgrid-scale contribution to the mean state of the atmosphere. The multivariate double-Gaussian mixture has been proposed as an appropriate model for PBL turbulence and shallow convection, but it is unable to reproduce important features of shallow cumulus convection. In this study, a novel unified PBL turbulence–convection–cloud macrophysics scheme is presented based on the eddy-diffusivity/mass-flux framework. The new scheme augments the double-Gaussian representation of subgrid variability with multiple stochastic mass-flux plumes at minimal added computational cost. Improved results for steady-state maritime and transient continental shallow convection from a single-column model implementation of the new scheme are shown with respect to reference large-eddy simulations. Improvements are seen in the cloud layer due to mass-flux plumes occupying the extreme moist, low liquid-water potential temperature tail of the joint temperature–moisture distribution. Significance Statement: Computer models of the atmosphere used to predict future climate are unable to directly represent air motion at small spatial scales because it would take too long to run the model over the entire planet. Instead, models typically use coarse model grid spacing and a simplified statistical representation of the physical processes that cause small-scale motions. This paper improves a particular simplified representation by adding a mechanism to represent statistically rare events of strong small-scale air motion that coherently transport air from near the surface to higher in the atmosphere. This increased transport also improves the representation of clouds, a particularly difficult phenomenon to simulate in models. [ABSTRACT FROM AUTHOR]

Concentrated poleward flows near the eastern boundaries between 2- and 4-km depth have been observed repeatedly, particularly in the Southern Hemisphere. These deep eastern boundary currents (DEBCs) play an important role in setting the large-scale tracer distribution and have nonnegligible contribution to global transports of mass, heat, and tracers, but their dynamics are not well understood. In this paper, we first demonstrate the significant role of DEBCs in the southeastern Atlantic, Indian, and Pacific Oceans, using the Southern Ocean State Estimate (SOSE) data assimilating product, and using high-resolution regional general circulation model configurations. The vorticity balances of these DEBCs reveal that, over most of the width of such currents, they are in an interior-like vorticity budget, with the meridional advection of planetary vorticity by and vortex stretching fwz being the largest two terms, and with contributions of nonlinearity and friction that are of smaller spatial scale. The stretching is shown, using a temperature budget, to be largely forced by resolved or parameterized eddy temperature transport. Strongly decaying signals from the eastern boundary in friction and stretching form the dominant balance in a sublayer close to the eastern boundary. The temporal variability of DEBCs is then examined, to help to interpret observations that tend to be sporadic in both time and space. The probability distribution functions of northward velocity inDEBCregions are broad, implying that flow reversals are common. Although the regions of the simulated DEBCs are generally local minima of eddy kinetic energy, they are still constantly releasing westward-propagating Rossby waves. [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]

Several statistical downscaling models have been developed in the past couple of decades to assess the hydrologic impacts of climate change by projecting the station-scale hydrological variables from large-scale atmospheric variables simulated by general circulation models (GCMs). This paper presents and compares different statistical downscaling models that use multiple linear regression (MLR), positive coefficient regression (PCR), stepwise regression (SR), and support vector machine (SVM) techniques for estimating monthly rainfall amounts in the state of Florida. Mean sea level pressure, air temperature, geopotential height, specific humidity, U wind, and V wind are used as the explanatory variables/predictors in the downscaling models. Data for these variables are obtained from the National Centers for Environmental Prediction-National Center for Atmospheric Research (NCEP-NCAR) reanalysis dataset and the Canadian Centre for Climate Modelling and Analysis (CCCma) Coupled Global Climate Model, version 3 (CGCM3) GCM simulations. The principal component analysis (PCA) and fuzzy c-means clustering method (FCM) are used as part of downscaling model to reduce the dimensionality of the dataset and identify the clusters in the data, respectively. Evaluation of the performances of the models using different error and statistical measures indicates that the SVM-based model performed better than all the other models in reproducing most monthly rainfall statistics at 18 sites. Output from the third-generation CGCM3 GCM for the A1B scenario was used for future projections. For the projection period 2001-10, MLR was used to relate variables at the GCM and NCEP grid scales. Use of MLR in linking the predictor variables at the GCM and NCEP grid scales yielded better reproduction of monthly rainfall statistics at most of the stations (12 out of 18) compared to those by spatial interpolation technique used in earlier studies. [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]

Previous studies have shown that Kelvin–Rossby instability is a viable mechanism for producing equatorial superrotation in small and/or slowly rotating planets. It is shown in this paper that this mechanism can also produce superrotation with terrestrial parameters when the baroclinic forcing moves to low latitudes, explaining previous results by Williams. The transition between superrotating and subrotating flow occurs abruptly as the baroclinic forcing moves poleward. Although Kelvin–Rossby instability weakens when the baroclinic zone moves away from the equator, the key factor explaining the abrupt transition is the change in the baroclinic eddies. When differential heating is contained within the tropics, baroclinic eddies do not decelerate the subtropical jet and the upper-tropospheric flow approximately conserves angular momentum, providing conditions favorable for Kelvin–Rossby instability. In contrast, when baroclinic eddies are generated in the extratropics, they decelerate the subtropical jet and prevent the Kelvin–Rossby coupling. Due to this sensitivity to baroclinic eddies, the system exhibits hysteresis: near the transition parameter, extratropical eddies can prevent superrotation when they are initially present. [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]

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]

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]

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]

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]

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]

The westerly jet (WJ) is an important component of atmospheric circulation, which is characterized by prominent seasonal variations in intensity and position. However, the response of the WJ over Asia during the Last Glacial Maximum (LGM) is still not clear. Using general circulation model experiments, the seasonal behaviors of the WJ over central Asia and Japan are analyzed in this paper. The results show that, compared to the present day (PD), the WJ presents a complicated response during the LGM, both in intensity and position. Over central Asia, it becomes weaker in both summer and winter. But over Japan, it is enhanced in summer but becomes diminished in winter. In terms of position, the WJ over central Asia shifts southward in both summer and winter, whereas the WJ over Japan moves southward in summer but does not change obviously relative to PD in winter. Such WJ changes are well explained by meridional temperature gradients in high troposphere, which is closely linked to seasonal thermal anomalies over the Tibetan Plateau (TP). Despite cooler LGM conditions, the anomalous warming center over the TP becomes stronger in summer. Derived from the heat budget equation, the stronger heating center is mainly caused by the weaker adiabatic cooling generated from ascending motion over the area south of the TP. In winter, the cooling over the TP is also strengthened, mostly owing to the subsidence-induced weaker adiabatic heating. Due to the importance of the WJ, the potential role of TP thermal effects should be a focus when explaining the East Asian climate change during the LGM. [ABSTRACT FROM AUTHOR]

This paper concerns the simulation of the water table elevation in shallow unconfined aquifers where infiltration is assumed as the main mechanism of recharge. The main aim is to provide a reliable tool for groundwater management that satisfies water supply managers. Such a tool is a candidate as a physically based alternative to the use of empirical methods or general circulation models. It is based on the use of two widely available sets of data: the water table elevation measurements and soil moisture time series. In fact, the former are usually provided by government agencies on public websites whereas the latter are included in the atmospheric global datasets (reanalysis). It is notable that data from reanalysis are accessible to any citizen and organization around the world on an open-access basis (e.g., Copernicus). In the proposed method, the measured water table elevations are correlated quantitatively with the water fluxes toward the aquifer evaluated using the soil moisture data from ERA5 reanalysis (provided by ECMWF) within a Richards equation-based approach. The analysis is executed using data from the Umbria region (Italy) on both a daily and monthly scale. In fact, these are the time intervals of interest for a proper management of groundwater resources. The proposed relationships include both a logarithmic and linear term and point out the possible different regimes of the shallow aquifers with regard to the recharge due to infiltration. These different mechanisms reflect in the different role played by the water fluxes toward the aquifer in terms of water table elevation changes according to the considered time scale. [ABSTRACT FROM AUTHOR]

When record-breaking climate and weather extremes occur, decision-makers and planners want to know whether they are random natural events with historical levels of reoccurrence or are reflective of an altered frequency or intensity as a result of climate change. This paper describes a method to attribute extreme weather and climate events to observed increases in atmospheric CO2 using an initialized subseasonal to seasonal coupled global climate prediction system. Application of this method provides quantitative estimates of the contribution arising from increases in the level of atmospheric CO2 to individual weather and climate extreme events. Using a coupled subseasonal to seasonal forecast system differs from other methods because it has the merit of being initialized with the observed conditions and subsequently reproducing the observed events and their mechanisms. This can aid understanding when the reforecasts with and without enhanced CO2 are compared and communicated to a general audience. Atmosphere–ocean interactions are accounted for. To illustrate the method, we attribute the record Australian heat event of October 2015. We find that about half of the October 2015 Australia-wide temperature anomaly is due to the increase in atmospheric CO2 since 1960. This method has the potential to provide attribution statements for forecast events within an outlook period (i.e., before they occur). This will allow for informed messaging to be available as required when an extreme event occurs, which is of particular use to weather and climate services. [ABSTRACT FROM AUTHOR]

Ocean lee waves occur on length scales that are smaller than the grid scale of global circulation models (GCMs). Therefore, such models must parameterize the drag associated with launching lee waves. This paper compares the lee wave drag predicted by existing parameterizations with the drag measured in highresolution nonhydrostatic numerical simulations of a lee wave over periodic sinusoidal bathymetry. The simulations afford a time-varying glimpse at the nonlinear and nonhydrostatic oceanic lee wave spinup process and identify a characteristic time scale to reach steady state. The maximum instantaneous lee wave drag observed during the spinup period is found to be well predicted by linear lee wave theory for all hill heights. In steady state, the simulations demonstrate the applicability of parameterizing the drag based on applying linear theory to the lowest overtopping streamline of the flow (LOTS), as is currently employed in GCMs. However, because existing parameterizations are based only on the height of the LOTS, they implicitly assume hydrostatic flow. For hills tall enough to trap water in their valleys, the simulations identify a set of nonhydrostatic processes that can result in a reduction of the lee wave drag from that given by hydrostatic parameterizations. The simulations suggest implementing a time-dependent nonhydrostatic version of the LOTS-based parameterization of lee wave drag and demonstrate the remarkable applicability of linear lee wave theory to oceanic lee waves. [ABSTRACT FROM AUTHOR]

Atmospheric stochastic forcing associated with the North Atlantic Oscillation (NAO) and intrinsic ocean modes associated with the large-scale baroclinic instability of the North Atlantic Current (NAC) are recognized as two strong paradigms for the existence of the Atlantic multidecadal oscillation (AMO). The degree to which each of these factors contribute to the low-frequency variability of the North Atlantic is the central question in this paper. This issue is addressed here using an ocean general circulation model run under a wide range of background conditions extending from a supercritical regime where the oceanic variability spontaneously develops in the absence of any atmospheric noise forcing to a damped regime where the variability requires some noise to appear. The answer to the question is captured by a single dimensionless number Γ measuring the ratio between the oceanic and atmospheric contributions, as inferred from the buoyancy variance budget of the western subpolar region. Using this diagnostic, about two-thirds of the sea surface temperature (SST) variance in the damped regime is shown to originate from atmospheric stochastic forcing whereas heat content is dominated by internal ocean dynamics. Stochastic wind stress forcing is shown to substantially increase the role played by damped ocean modes in the variability. The thermal structure of the variability is shown to differ fundamentally between the supercritical and damped regimes, with abrupt modifications around the transition between the two regimes. Ocean circulation changes are further shown to be unimportant for setting the pattern of SST variability in the damped regime but are fundamental for a preferred time scale to emerge. [ABSTRACT FROM AUTHOR]

In this paper, it is shown that westward-propagating monsoon low pressure system–like disturbances in the South Asian monsoon region can be simulated in an idealized moist general circulation model through the addition of a simplified parameterization of land. Land is parameterized as having one-tenth the heat capacity of the surrounding slab ocean, with evaporation limited by a bucket hydrology model. In this model, the prominent topography of the Tibetan Plateau does not appear to be necessary for these storm systems to form or propagate; therefore, focus is placed on the simulation with land but no topography. The properties of the simulated storms are elucidated using regression analysis and compared to results from composites of storms from comprehensive GCMs in prior literature and reanalysis. The storms share a similar vertical profile in anomalous Ertel potential vorticity to those in reanalysis. Propagation, however, does not seem to be strongly dictated by beta drift. Rather, it seems to be more closely consistent with linear moisture vortex instability theory, with the exception of the importance of the vertical advection term in the Ertel potential vorticity budget toward the growth and maintenance of disturbances. The results presented here suggest that a simplified GCM configuration might be able to be used to gain a clearer understanding of the sensitivity of monsoon low pressure systems to changes in the mean state climate. [ABSTRACT FROM AUTHOR]

Coupled data assimilation (DA) provides a consistent framework for assimilating satellite observations that are sensitive to several components of the Earth system. In this paper, we focus on low-peaking infrared satellite channels that are sensitive to the lower atmosphere and Earth surface temperature (EST) over both ocean and land. Our atmospheric hybrid-4DVAR system [the Navy Global Environmental Model (NAVGEM)] is extended to include the following: 1) variability in the sea surface temperature (both diurnal variability and climatological perturbations to the ensemble members), 2) the coupled Jacobians of the radiative transfer model for the infrared sensors, and 3) the coupled covariances between the EST and the atmosphere. Our coupling approach is found to improve forecast accuracy and to provide corrections to the EST that are in balance with the atmospheric analysis. The largest impact of the coupling is found on near-surface atmospheric temperature and humidity in the tropics, but the impact extends all the way to the stratosphere. The role of each coupling element on the performance of the global atmospheric circulation model is investigated. Inclusion of variability in the sea surface temperature has the strongest positive impact on the forecast quality. Additional inclusion of the coupled Jacobian and ensemble-based coupled covariances led to further improvements in scores and to modification of the corrections to the ocean boundary layer. Coupled DA had significant impact on latent and sensible heat fluxes over land, locations of western boundary currents, and along the ice edge. [ABSTRACT FROM AUTHOR]

In aiming for better access to climate change information and for providing climate service, it is important to obtain reliable high-resolution temperature simulations. Systematic comparisons are still deficient between statistical and dynamic downscaling techniques because of their inherent unavoidable uncertainties. In this paper, 20 global climate models (GCMs) and one regional climate model [Providing Regional Climates to Impact Studies (PRECIS)] are employed to evaluate their capabilities in reproducing average trends of mean temperature (Tm), maximum temperature (Tmax), minimum temperature (Tmin), diurnal temperature range (DTR), and extreme events represented by frost days (FD) and heat-wave days (HD) across China. It is shown generally that bias of temperatures from GCMs relative to observations is over ±1°C across more than one-half of mainland China. PRECIS demonstrates better representation of temperatures (except for HD) relative to GCMs. There is relatively better performance in Huanghuai, Jianghuai, Jianghan, south Yangzi River, and South China, whereas estimation is not as good in Xinjiang, the eastern part of northwest China, and the Tibetan Plateau. Bias-correction spatial disaggregation is used to downscale GCMs outputs, and bias correction is applied for PRECIS outputs, which demonstrate better improvement to a bias within ±0.2°C for Tm, Tmax, Tmin, and DTR and ±2 days for FD and HD. Furthermore, such improvement is also verified by the evidence of increased spatial correlation coefficient and symmetrical uncertainty, decreased root-mean-square error, and lower standard deviation for reproductions. It is seen from comprehensive ranking metrics that different downscaled models show the most improvement across different climatic regions, implying that optional ensembles of models should be adopted to provide sufficient high-quality climate information. [ABSTRACT FROM AUTHOR]

Africa is one of the three key regions of deep convection in the global tropics. There is a wealth of information on the intensity, variability, and change of convection and associated rainfall in regions across the continent but almost all of this literature is regionally focused and confined to specific seasons. This fragmented approach precludes a continent-wide view of deep convection leaving the following key issues unanswered: When is deep convection the most widespread across Africa? Where on the continent is deep convection most active? Where does widespread convection have the most interannual variability? This paper confronts these questions using a satellite-derived integral of deep convection. At the continental scale, March exhibits the most extensive deep convection whereas the West African monsoon during June-July exhibits the least. El Niño generally suppresses pan-African convective activity while La Niña enhances this activity. These pan-African signals are largely determined by regional hotspots: the eastern Congo hosts the most persistent widespread deep convection, southeastern southern Africa displays the highest interannual variability, and regional highlands maintain local convective activity hotspots. Furthermore, pan-African annual mean convective activity has increased ;10% between 1983 and 2015 with increases of .20% recorded in local hotspots. Results in this study provide a climatological baseline for both observational and model-based studies of African climates and offer insights into when African convection has the greatest potential impact on the general circulation. [ABSTRACT FROM AUTHOR]

This study evaluates the feasibility of applying the Estimating the Circulation and Climate of the Ocean (ECCO) data assimilation (DA) framework to a submesoscale-resolving model [Massachusetts Institute of Technology General Circulation Model (MITgcm), at 1-km grid resolution]. This is in preparation for future studies to understand and assimilate the novel swath measurements of sea surface height from the Surface Water and Ocean Topography (SWOT) satellite mission. The model domain is centered at the SWOT calibration/validation (CalVal) site, located 300 km offshore of Monterey Bay, California. We assimilate vertical profiles of temperature and salinity from a linear array of three SWOT prelaunch CalVal moorings, between September and December 2019. Two model solutions are analyzed: 1) a nonassimilating forward simulation termed "first-guess" and 2) an optimized solution that assimilates the in situ observations. Both runs are nested within the global 1/12° Hybrid Coordinate Ocean Model that uses Navy Coupled Ocean Data Assimilation (HYCOM/NCODA) analysis. We evaluate the performance by comparing the two model solutions against assimilated and withheld in situ observations. We show that by assimilating hydrographic data, the model performance over the first-guess solution is improved with an error-to-observation reduction of 30%–45% for temperature, 14%–41% for salinity, and 44%–56% for steric height. Over the study period, the average steric height error at the three assimilated moorings was 1.27 cm for the optimized solution versus 2.6 cm for the first-guess solution. A comparison to withheld glider observations shows that the optimized solution outperforms the first-guess solution with an error-to-observation reduction of approximately 38% in steric height (1.5 versus 2.4 cm). This study indicates, for the first time, that the MITgcm–ECCO framework can be successfully applied to the reconstruction of submesoscale ocean variability, via the nesting of a high-resolution regional domain into a global outer domain. [ABSTRACT FROM AUTHOR]

Using a range of Detection and Attribution Model Intercomparison Project (DAMIP) simulations from phase 6 of Coupled Model Intercomparison Project (CMIP6), we study the response of dynamic sea level (DSL) to external anthropogenic climate forcing [greenhouse gases (GHGs), aerosols, and stratospheric ozone] with a focus on the differences over the twentieth and twenty-first century. In the second half of the twentieth century, the DSL nonuniformity in the Northern Hemisphere (NH) was relatively small due to a cancellation between the effects of increasing GHGs and aerosols. In contrast, the DSL signal in the Southern Hemisphere (SH) over this period was large because stratospheric ozone depletion reinforced the effects of increasing GHGs. In the twenty-first century, the DSL response has been intensified in the NH because the warming effects of diminishing aerosols have acted to reinforce the effects of increasing GHGs. Meanwhile, the distribution of SH DSL has also become uneven although stratospheric ozone recovery has partially offset the effects of rising GHGs. Using a global ocean circulation model, we decompose the changes in the twenty-first century DSL into distinct responses to surface forcings including sea surface temperature, salinity, and wind stress. Our results show that the dipole-like pattern of DSL in the North Pacific can be attributed largely to sea surface warming, while the dipole-like pattern in the North Atlantic is attributed to subpolar surface salinity freshening. The belted pattern of DSL changes in the Southern Ocean is induced by both surface warming and intensifying/poleward-shifting westerly winds. [ABSTRACT FROM AUTHOR]

The time at which climate change signal can be clearly distinguished from noise is known as time of emergence (ToE) and is typically detected by a general circulation model (GCM) signal-to-noise ratio exceeding a certain threshold. ToE is commonly estimated at large scales from GCMs, although management decisions and adaptation strategies are implemented locally. This paper proposes a methodology to estimate ToE for both precipitation and temperature at local scales (i.e., river basin). The methodology considers local climatic conditions and unbiased GCM projections to estimate ToE by using the statistical power to find when the climate significantly differs from the historical one. The method suggests that ToE for temperature already occurred in three Chilean basins (Limarí, Maipo, and Maule). However, in terms of precipitation, an earlier ToE is clearly identified for the Maule basin, indicating that risk assessment and adaptation measures should be implemented first in this basin. [ABSTRACT FROM AUTHOR]

Low-frequency intraseasonal oscillations in the tropical atmosphere in general circulation models (GCMs) were studied extensively in many previous studies. However, the simulation of the quasi-biweekly oscillation (QBWO), which is an important component of the intraseasonal oscillations, in GCMs has not received much attention. This paper evaluates the QBWO features over the South China Sea in early [May–June (MJ)] and late [August–September (AS)] summer in the National Center for Atmospheric Research (NCAR) Community Atmosphere Model, version 5.3 (CAM5), using observations and reanalysis data. Results show that the major features of the spatial distribution of the QBWO in both MJ and AS are simulated reasonably well by the model, although the amplitude of the variation is overestimated. CAM5 captures the local oscillation in MJ and the westward propagation in AS of the QBWO. Although there are important biases in geographical location and intensity in MJ, the model represents the QBWO horizontal and vertical structure qualitatively well in AS. The diagnosis of the eddy vorticity budget is conducted to better understand the QBWO activities in the model. Both horizontal advection of relative vorticity and that of planetary vorticity (Coriolis parameter) are important for the local evolution of the QBWO in MJ in observations as well as model simulation, whereas advection of planetary vorticity contributes to the westward propagation of QBWO vorticity anomalies in AS. Since the Coriolis parameter f only changes with latitude, this suggests that the correct simulation of anomalous meridional wind is a key factor in the realistic simulation of the QBWO in the model. [ABSTRACT FROM AUTHOR]

The evidence for, baroclinic energetics of, and geographic distribution of parametric subharmonic instability (PSI) arising from both diurnal and semidiurnal tides in a global ocean general circulation model is investigated using 1/12.5° and 1/25° simulations that are forced by both atmospheric analysis fields and the astronomical tidal potential. The paper examines whether PSI occurs in the model, and whether it accounts for a significant fraction of the tidal baroclinic energy loss. Using energy transfer calculations and bispectral analyses, evidence is found for PSI around the critical latitudes of the tides. The intensity of both diurnal and semidiurnal PSI in the simulations is greatest in the upper ocean, consistent with previous results from idealized simulations, and quickly drops off about 5° from the critical latitudes. The sign of energy transfer depends on location; the transfer is positive (from the tides to subharmonic waves) in some locations and negative in others. The net globally integrated energy transfer is positive in all simulations and is 0.5%–10% of the amount of energy required to close the baroclinic energy budget in the model. The net amount of energy transfer is about an order of magnitude larger in the 1/25° semidiurnal simulation than the 1/12.5° one, implying the dependence of the rate of energy transfer on model resolution. [ABSTRACT FROM AUTHOR]

Monthly Weather Review, one of the oldest continuously published meteorological journals in the world, publishes its 150th volume this year. In January 1873, the U.S. War Department's Army Signal Service began producing this monthly report summarizing weather across the United States. Its first issue consisted of a one-page narrative of weather conditions and storms and one chart depicting tracks of low pressure centers during that month. In 1891, Monthly Weather Review continued as a government publication with the transfer of the nation's weather service from the military to the newly established U.S. Weather Bureau. Over time and sometimes erratically, it grew into a scientific journal. In 1974, Monthly Weather Review was transferred to the American Meteorological Society, who continues to publish it to this day (although a 2003 proposal might have ended it). This Historical Review discusses some of the journal's history and impact, as well as its legacy. This review also compiles for the first time a complete list of Monthly Weather Review editors. The research published within Monthly Weather Review has included highly cited, ground-breaking articles on weather and climate phenomena (e.g., extratropical and tropical cyclones, El Niño–Southern Oscillation, Madden–Julian oscillation), general circulation modeling, and numerical weather prediction. The data published in the early issues have been used—and continue to be used to this day—for a variety of applied research and historical analysis purposes. The composition and content of Monthly Weather Review have changed over the past century and a half, continuing to evolve with the modern publishing landscape, with color figures at no additional cost, open-access articles, open data, and, in the near future, embedded figure animations. [ABSTRACT FROM AUTHOR]

We assess the representation of mesoscale stirring in a suite of models against an estimate derived from microstructure data collected during the North Atlantic Tracer Release Experiment (NATRE). We draw heavily from the approximate temperature variance budget framework of Ferrari and Polzin. This framework assumes two sources of temperature variance away from boundaries: first, the vertical stirring of the large-scale mean vertical gradient by small-scale turbulence; and second, the lateral stirring of large-scale mean along-isopycnal gradients by mesoscale eddies. Temperature variance so produced is transformed and on average transferred down scales for ultimate dissipation at the microscale at a rate χ estimated using microstructure observations. Ocean models represent these pathways by a vertical mixing parameterization, and an along-isopycnal lateral mixing parameterization (if needed). We assess the rate of variance production by the latter as a residual from the NATRE dataset and compare against the parameterized representations in a suite of model simulations. We find that variance production due to lateral stirring in a Parallel Ocean Program version 2 (POP2) 1/10° simulation agrees well, to within the estimated error bars, with that inferred from the NATRE estimate. A POP2 1° simulation and the Estimating the Circulation and Climate of the Ocean Version 4 release 4 (ECCOV4r4) simulation appear to dissipate an order of magnitude too much variance by applying a lateral diffusivity, when compared to the NATRE estimate, particularly below 1250 m. The ECCOV4r4-adjusted lateral diffusivities are elevated where the microstructure suggests elevated χ sourced from mesoscale stirring. Such elevated values are absent in other diffusivity estimates suggesting the possibility of compensating errors and caution in interpreting ECCOV4r4's adjusted lateral diffusivities. Significance Statement: We look at whether microstructure turbulence observations can provide a useful metric for judging the fidelity of representation of mesoscale stirring in a suite of models. We focus on the region of the North Atlantic Tracer Release Experiment (NATRE), the site of a major ocean turbulence observation campaign, and use an approximate variance budget framework for the region with observational estimates from Ferrari and Polzin (2005). The approach provides a novel framework to evaluate the approximate representation of mesoscale stirring in a variety of models. [ABSTRACT FROM AUTHOR]

Robust and accurate Gulf Stream separation remains an unsolved problem in general circulation modeling whose resolution will positively impact the ocean and climate modeling communities. Oceanographic literature does not face a shortage of plausible hypotheses that attempt to explain the dynamics of the Gulf Stream separation, yet a single theory that the community agrees on is missing. In this paper, the authors investigate the impact of the deep western boundary current (DWBC), coastline curvature, and continental shelf steepening on the Gulf Stream separation within regional configurations of theMassachusetts Institute of Technology General CirculationModel. Artificial modifications to the regional bathymetry are introduced to investigate the sensitivity of the separation to each of these factors. Metrics for subsurface separation detection confirm the direct link between flow separation and the surface expression of theGulf Stream in the Mid-Atlantic Bight. It is shown that the Gulf Stream separation and mean surface position are most sensitive to the continental slope steepening, consistent with a theory proposed by Melvin Stern in 1998. In contrast, the Gulf Stream separation exhibits minimal sensitivity to the presence of the DWBC and coastline curvature. The implications of these results to the development of a "separation recipe" for oceanmodeling are discussed. This study concludes adequate topographic resolution is a necessary, but not sufficient, condition for proper Gulf Stream separation. [ABSTRACT FROM AUTHOR]

The paper demonstrates that sudden stratospheric warmings (SSWs) can be simulated in an ensemble of dry dynamical cores that miss the typical SSW forcing mechanisms like moist processes, land-sea contrasts, or topography. These idealized general circulation model (GCM) simulations are driven by a simple Held-Suarez-Williamson (HSW) temperature relaxation and low-level Rayleigh friction. In particular, the four dynamical cores of NCAR's Community Atmosphere Model, version 5 (CAM5), are used, which are the semi-Lagrangian (SLD) and Eulerian (EUL) spectral-transform models and the finite-volume (FV) and the spectral element (SE) models. Three research themes are discussed. First, it is shown that SSW events in such idealized simulations have very realistic flow characteristics that are analyzed via the SLD model. A single vortex-split event is highlighted that is driven by wavenumber-1 and -2 wave-mean flow interactions. Second, the SLD simulations are compared to the EUL, FV, and SE dynamical cores, which sheds light on the impact of the numerical schemes on the circulation. Only SLD produces major SSWs, while others only exhibit minor stratospheric warmings. These differences are caused by SLD's more vigorous wave-mean flow interactions in addition to a warm pole bias, which leads to relatively weak polar jets in SLD. Third, it is shown that tropical quasi-biennial oscillation (QBO)-like oscillations and SSWs can coexist in such idealized HSW simulations. They are present in the SLD dynamical core that is used to analyze the QBO-SSW interactions via a transformed Eulerian-mean (TEM) analysis. The TEM results provide support for the Holton-Tan effect. [ABSTRACT FROM AUTHOR]

Satellite observation of sea surface height (SSH) may soon have sufficient accuracy and resolution to map geostrophic currents in Lake Superior. A dynamic atmosphere correction will be needed to remove SSH variance due to basinwide seiching. Here, the dynamics of rotating barotropic gravity modes are examined using numerical models and lake-level gauges. Gravity modes explain 94% of SSH variance in a general circulation model and evolve as forced, damped oscillators. These modes have significant SSH, but negligible kinetic energy (2 J m−2) and dissipation rates (0.01 mW m−2) relative to other motions in Lake Superior. Removing gravity modes from instantaneous SSH allows geostrophic currents to be accurately computed. Complex empirical orthogonal functions (CEOFs) from 50 years of data at eight lake-level gauges show patterns consistent with the first two gravity modes. The frequency spectra of these CEOFs are consistent with forced, damped oscillators with natural frequencies of 3.05 and 4.91 cycles per day (cpd) and decay time scales of 4.5 and 1.0 days. Modal amplitudes from the general circulation model and lake-level gauges are 80% coherent at 1 cpd, but only 50% coherent at 3 cpd, indicating that the atmospheric reanalysis used to force the general circulation model is not accurate at the high natural frequencies of the gravity modes. The results indicate that a dynamic atmosphere correction should combine modeled gravity modes below 1 cpd and observed mode-1 and mode-2 amplitudes (from lake-level gauges) at higher frequencies. An inverted barometer correction is also recommended to account for low-frequency atmospheric pressure gradients that do not project onto gravity modes. [ABSTRACT FROM AUTHOR]

One year of Spread spectrum Interferometric Multistatic meteor radar Observing Network (SIMONe) measurements are analyzed and compared for the first time between two low-latitude locations in Peru: Jicamarca (12°S, 77°W) and Piura (5°S, 80°W). Investigation of the mean horizontal winds and tides reveals that mesosphere and lower thermosphere (MLT) planetary-scale dynamics are similar between these two locations, although differences can be seen in some tidal components, e.g., the diurnal tide. On the other hand, 28-day median values of the momentum fluxes obtained with 4-h, 4-km time–altitude bins indicate that the mesoscale dynamics differ significantly between Jicamarca and Piura, places separated by approximately 850 km. From the middle of July until October 2021, a strong acceleration of the background zonal wind by westward-propagating gravity waves (GWs) is observed above ∼90 km at both locations, although with larger amplitudes over Jicamarca. From the middle of January until April 2022, a second strong acceleration of the background zonal wind, again by westward-propagating GWs, is observed, but this time with larger amplitudes over Piura. The latter is further supported by the dominance of negative vertical gradients of the zonal momentum flux above 89 km of altitude. Thus, these results observationally confirm the previous studies based on general circulation model simulations indicating that the directions of the zonal GW drag and the zonal background wind coincide in the low-latitude MLT. The weak correlations between the horizontal wind gradients over Jicamarca and Piura reinforce the fact that the mesoscale dynamics are different at these two locations. [ABSTRACT FROM AUTHOR]

The impact of melt ponds on sea ice albedo has been observed and documented. In general circulation models, ponds are now accounted for through indirect diagnostic treatments ("implicit" schemes) or prognostic melt-pond parameterizations ("explicit" schemes). However, there has been a lack of studies showing the impacts of these schemes on simulated Arctic climate. We focus here on rectifying this using the general circulation model HadGEM3, one of the few models with a detailed explicit pond scheme. We identify the impact of melt ponds on the sea ice and climate, and associated ice–ocean–atmosphere interactions. We run a set of constant forcing simulations for three different periods and show, for the first time, that using mechanistically different pond schemes can lead to very significantly different sea ice and climate states. Under near-future conditions, an implicit scheme never yields an ice-free summer Arctic, while an explicit scheme yields an ice-free Arctic in 35% of years and raises autumn Arctic air temperatures by 5° to 8°C. We find that impacts on climate and sea ice depend on the ice state: under near-future and last-interglacial conditions, the thin sea ice is very sensitive to pond formation and parameterization, whereas during the preindustrial period the thicker sea ice is less sensitive to the pond scheme choice. Both of these two commonly used parameterizations of sea ice albedo yield similar results under preindustrial conditions but in warmer climates lead to very different Arctic sea ice and ocean and atmospheric temperatures. Thus, changes to physical parameterizations in the sea ice model can have large impacts on simulated sea ice, ocean, and atmosphere. Significance Statement: This study investigates the impacts of melt ponds on Arctic sea ice under different climate conditions, using the HadGEM3-GC3.1-LL general circulation model (GCM). Additionally, we study the impact of changing the type of pond scheme used. We find that changing the pond scheme causes large differences to how a GCM simulates Arctic sea ice, the ocean, and the atmosphere, for both near-future and warmer paleoclimate conditions. These large differences have not been found previously, because this is one of the first GCM studies of this type. Our results demonstrate the importance of melt ponds, and their wider impacts on ocean and atmosphere. Furthermore, they suggest that better evaluation of the representation of sea ice processes is vital for the robust projection of future climate change. [ABSTRACT FROM AUTHOR]

The simple scaling relation for internal-tide generation proposed by Jayne and St. Laurent is widely used for parameterizing turbulent mixing induced by breaking of internal tides. Based on the internal-tide generation derived from a 0.1° ocean general circulation model, we show that depending on which stratification is used, this relation produces different vertical distributions of internal-tide generation. When using the buoyancy frequency at the seafloor, which is a common practice, the scaling relation produces, relative to the model, too-strong internal-tide generation in the upper 2000 m and too-weak internal-tide generation in the lower 2000 m. Moreover, the different vertical distributions in the different ocean basins, characterized by a generally decreasing internal tide generation with increasing depth in the Indo-Pacific but not-decreasing or even increasing internal tide generation with increasing depth in the upper 3000 m of the Atlantic, cannot be captured when using bottom stratification. These unsatisfactory features can be easily removed by replacing the buoyancy frequency at the seafloor by a buoyancy frequency averaged over a large part of the water column. To our knowledge, this sensitivity to stratification has not been explicitly quantified for the global ocean. Because of this sensitivity, the scaling relation of Jayne and St. Laurent should be used with an averaged stratification to ensure a more adequate representation of turbulent diffusivity due to tidal mixing and water mass transformation in the deep oceans. [ABSTRACT FROM AUTHOR]

To understand the intensification process of tropical cyclones (TCs), we analyzed the relationship between the TC intensification rate (I) and environmental variables along TC tracks during the time from TC genesis (tG) to maximum TC strength (tX), hereinafter τGX ≡ tX − tG, using a state-of-the-art general circulation model (GCM), observed TC tracks, and ERA5 data. During τGX, strong TCs with high I (sTCs) consume more convective available potential energy (CAPE) than weak TCs with low I (wTCs) and bring more CAPE from the equator to sustain sTCs. Compared to wTCs, sTCs prefer an unstable atmosphere with higher sea surface temperature (SST), stronger grid-mean upward flow at 500 hPa (ω500), more moisture convergence (MC), and weaker wind shear (Vs). Our GCM simulation shows that MC and CAPE have a single regression slope with I applicable both within and across climate regimes. Using machine learning, we found that the best combination of environmental variables (V6) for predicting I consists of ω500, MC, SST, midtropospheric stability (MTS), Vs, and latitude (|f|). Machine learning with V6 reproduces well the spatial distribution and interclimate changes of I: TCs are intensified in regions of stronger upward ω500, more MC, warmer SST, weaker MTS, smaller Vs, and larger |f|; TCs in a warmer climate have higher I than TCs in a colder climate due to more MC, warmer SST, but stronger MTS. These results are consistent with the conceptual understanding that TCs are intensified by the release of latent heat. [ABSTRACT FROM AUTHOR]

Radiative feedbacks over interannual time scales can be potentially useful for global warming estimation. However, the diversity of the lead–lag relationships in global mean surface temperature (GMST) and net radiation flux at the top of the atmosphere (GMTOA) create uncertainty during the estimation of radiative feedbacks. In this study, key physical processes controlling lead–lag relationships were elucidated by categorizing preindustrial control simulations of CMIP6 into three groups based on cross correlation values of GMTOA against GMST at lag 0 and lag +1 year. The diversity in the lead–lag relationships was primarily caused by the climatological state difference of the atmosphere over the equatorial Pacific, which modulated the strength of convective activity and sensitivity of low-level clouds. Diminished atmospheric stability caused enhanced convective activity, more efficient energy release, and smaller lags. In addition, enhanced stability in the lower atmosphere rendered the low-level clouds more sensitive to sea surface temperature changes and considerably delayed the radiative response. The climatological state difference of the atmosphere resulted from model-inherent atmospheric conditions. These findings suggest that the diversity of lead–lag relationships of GMST and GMTOA over interannual time scales could represent the characteristics of general atmospheric circulation models and possible solutions of the actual atmosphere, which could also affect long-term feedback features. [ABSTRACT FROM AUTHOR]

An idealized case of gradual oceanic transition from shallow to deep convection is simulated at three different horizontal resolutions: one that resolves most of the turbulent eddies, one typical of cloud-resolving models, and one typical of general circulation models. The former serves as a reference and allows the identification of clouds as individual objects, distinguishing shallow cumulus, congestus, and cumulonimbus. At coarser resolutions, parameterizations of convection are included and assessed, with a particular focus on congestus clouds and precipitation associated with shallow convective clouds. Congestus clouds are found to contribute the most to turbulent transport during the transition, while occupying a volume comparable to shallow cumulus and cumulonimbus. Kilometer-scale horizontal resolutions prove to be insufficient to resolve congestus, and parameterization schemes of shallow and deep convection are not necessarily appropriate to represent those intermediate clouds. The representation of rainfall in the shallow convection scheme plays a key role in the transition. Sensitivity experiments show that enhanced rainfall inhibits convection in single-column simulations, while it favors resolved convection and spatial heterogeneities in three-dimensional simulations with kilometer-scale resolution. Results highlight the need for an appropriate parameterization of congestus in both kilometer-scale and large-scale models. The case study and the methods presented here are proposed as a useful framework to evaluate models and their parameterizations in a shallow-to-deep convection transition context. [ABSTRACT FROM AUTHOR]

Throughout 2020, a persistent marine heatwave (MHW) event occurred in the northeast (NE) Pacific, which exhibited a record-breaking intensity and long duration. Three peaks of sea surface temperature (SST) anomaly are identified in April, July, and November 2020, respectively, all of which are caused by the surface latent heat flux (LH) term. A Taylor expansion analysis for the LH term reveals that positive LH anomaly (LHA) played the decisive role in increasing the latent heating effect, whereas the mixed layer depth anomaly was not important in the 2020 MHW. The positive LHA was primarily caused by the southerly wind anomalies, which advected more humid air from lower latitudes to the NE Pacific and thus reduced the difference between the saturated specific humidity at sea surface and actual surface specific humidity. This reduced sea–air humidity difference was conducive to less evaporation, leading to positive LHA. Using an atmospheric general circulation model forced by observed SSTs, we find that although the atmospheric circulation anomalies associated with the 2020 MHW can be partly constrained by both tropical Pacific SST and local SST over the North Pacific, the role of internal atmospheric variability in the 2020 MHW's formation cannot be overlooked. The prediction skill for the 2020 MHW in the North American Multimodel Ensemble models is limited to one month, indicating the challenge of accurately predicting MHW in the NE Pacific. The finding about the contribution of sea–air humidity difference to the LHA provides a new insight into the formation mechanism of MHW. Significance Statement: An exceptionally strong and persistent marine heatwave occurred in the northeast Pacific throughout 2020, with three peaks in April, July, and November, causing considerable damage to local fisheries and ecosystems. This study found that the evolution of this marine heatwave was due to the less-than-normal release of latent heat at the ocean surface. This was mainly caused by the reduced sea–air humidity difference, which was attributed to the increased near-surface moisture transported by southerly wind anomalies. Our model experiments demonstrated that the associated anomaly fields in the atmosphere partly stemmed from stochastic processes. The state-of-the-art dynamic models can only predict this marine heatwave one month in advance, indicating the necessity of improving prediction skill for marine heatwave's evolution in current models. [ABSTRACT FROM AUTHOR]

This paper discusses whether and how a baroclinic atmosphere can equilibrate with very small bottom friction in a dry primitive equation general circulation model. The model is forced by a Newtonian relaxation of temperature to a prescribed temperature profile, and it is damped by a linear friction near the lower boundary. When friction is decreased by four orders of magnitude, kinetic energy dissipation by friction gradually becomes negligible, while 'energy recycling' becomes dominant. In this limit kinetic energy is converted back into potential energy at the largest scales, thus closing the energy cycle without significant frictional dissipation. The momentum fluxes are of opposite sign in the upper and lower atmosphere: in the upper atmosphere, eddies converge momentum into the westerly jets; however, in the lower atmosphere, the eddies diverge momentum out of the westerly jets. The secondary circulation driven by the meridional eddy momentum fluxes thus acts to increase the baroclinicity of the westerly jet. This regime may be relevant for the Jovian atmosphere, where the frictional time scale may be much larger than the radiative damping time scale. [ABSTRACT FROM AUTHOR]

Cumulus parameterizations in general circulation models (GCMs) frequently apply mass-flux schemes in their description of tropical convection. Mass flux constitutes the product of the fractional area covered by cumulus clouds in a model grid box and the vertical velocity within the cumulus clouds. The cumulus area fraction profiles can be derived from precipitating radar reflectivity volumes. However, the vertical velocities are difficult to observe, making the evaluation of mass-flux schemes difficult. In this paper, the authors develop and evaluate a parameterization of vertical velocity in convective (cumulus) clouds using only radar reflectivities collected by a C-band polarimetric research radar (CPOL), operating at Darwin, Australia. The parameterization is trained using vertical velocity retrievals from a dual-frequency wind profiler pair located within the field of view of CPOL. The parametric model uses two inputs derived from CPOL reflectivities: the 0-dB Z echo-top height (0-dB Z ETH) and a height-weighted column reflectivity index ( ZHWT). The 0-dB Z ETH determines the shape of the vertical velocity profile, while ZHWT determines its strength. The evaluation of these parameterized vertical velocities using (i) the training dataset, (ii) an independent wind-profiler-based dataset, and (iii) 1 month of dual-Doppler vertical velocity retrievals indicates that the statistical representation of vertical velocity is reasonably accurate up to the 75th percentile. However, the parametric model underestimates the extreme velocities. The method allows for the derivation of cumulus mass flux and its variability on current GCM scales based only on reflectivities from precipitating radar, which could be valuable to modelers. [ABSTRACT FROM AUTHOR]

Conventional altimetry measures a one-dimensional profile of sea surface height (SSH) along the satellite track. Two-dimensional SSH can be reconstructed using mapping techniques; however, the spatial resolution is quite coarse even when data from several altimeters are analyzed. A new satellite mission based on radar interferometry is scheduled to be launched in 2020. This mission, called Surface Water and Ocean Topography (SWOT), will measure SSH at high resolution along a wide swath, thus providing two-dimensional images of the ocean surface topography. This new capability will provide a large amount of data even though they are contaminated with instrument noise and geophysical errors. This paper presents a tool that simulates synthetic observations of SSH from the future SWOT mission using SSH from any ocean general circulation model (OGCM). SWOT-like data have been generated from a high-resolution model and analyzed to investigate the sampling and accuracy characteristics of the future SWOT data. This tool will help explore new ideas and methods for optimizing the retrieval of information from future SWOT missions. [ABSTRACT FROM AUTHOR]

This study has developed a rigorous and efficient maximum likelihood method for estimating the parameters in stochastic energy balance models (with any k > 0 number of boxes) given time series of surface temperature and top-of-the-atmosphere net downward radiative flux. The method works by finding a state-space representation of the linear dynamic system and evaluating the likelihood recursively via the Kalman filter. Confidence intervals for estimated parameters are straightforward to construct in the maximum likelihood framework, and information criteria may be used to choose an optimal number of boxes for parsimonious k-box emulation of atmosphere–ocean general circulation models (AOGCMs). In addition to estimating model parameters the method enables hidden state estimation for the unobservable boxes corresponding to the deep ocean, and also enables noise filtering for observations of surface temperature. The feasibility, reliability, and performance of the proposed method are demonstrated in a simulation study. To obtain a set of optimal k-box emulators, models are fitted to the 4 × CO2 step responses of 16 AOGCMs in CMIP5. It is found that for all 16 AOGCMs three boxes are required for optimal k-box emulation. The number of boxes k is found to influence, sometimes strongly, the impulse responses of the fitted models. [ABSTRACT FROM AUTHOR]