Search

Improving Climate Models and Projections Using Observations Workshop 2023 What: More than 70 experts from research and operational centers came together in a hybrid meeting held at the Massachusetts Institute [...]

The response to Antarctic sea ice loss within a coupled modeling framework is examined in comparison to the response to Arctic sea ice loss and within the context of general greenhouse warming. Sea ice loss responses are found to be linear (particularly in response to Antarctic or global sea ice loss) with respect to the degree of imposed perturbation and additive when perturbations are applied in hemispheres separately and concurrently. Globally, and in the tropical Pacific in particular, Antarctic sea ice loss plays a relatively larger role than Arctic sea ice loss in both the atmosphere and the ocean, within the parameters of our experiments. The pattern of response to Antarctic sea ice loss is also found to more closely resemble that of greenhouse warming, again particularly in the tropics. An extension to multiparameter pattern scaling is developed to include a scaling factor for Antarctic change in addition to those for tropical warming and Arctic sea ice loss. The decomposition is applied to the modeled response to Antarctic sea ice loss to break it down into component partial responses that scale with Antarctic, tropical, and Arctic changes. This reveals the aspects of the response that are directly related to Antarctic change, such as an equatorward intensification of tropical precipitation in the Northern Hemisphere, and those that are modified via the induced changes in the tropics and Arctic, such as Northern Hemisphere temperature change. With this, we hope to gain a deeper understanding of the role of each of these changes for the development of physical mechanisms of the response. [ABSTRACT FROM AUTHOR]

The effect of anthropogenic climate change on extreme near-surface wind speeds is uncertain. Observed trends are weak and difficult to disentangle from internal variability, and model projections disagree on the sign and magnitude of trends. Standard coarse-resolution climate models represent the fine structures of relevant physical phenomena such as extratropical cyclones (ETCs), upper-level jet streaks, surface energy fluxes, and land surface variability less skillfully than their high-resolution counterparts. Here, we use simulations with the NCAR Community Earth System Model with both uniform (110 km) resolution and the variable-resolution configuration (VR-CESM-SONT, from 110 to 7 km) to study the effect of refined spatial resolution on projections of extreme strong and weak wind speeds in the Great Lakes region under end-of-century RCP8.5 forcing. The variable-resolution configuration projects strengthening of strong-wind events in the refined region with the opposite occurring in the uniform-resolution simulation. The two configurations provide consistent changes to synoptic-scale circulations associated with high-wind events. However, only the variable-resolution configuration projects weaker static stability, enhanced turbulent vertical mixing, and consequentially enhanced surface wind speeds because boundary layer dynamics are better captured in the refined region. Both models project increased frequency of extreme weak winds, though only VR-CESM-SONT resolves the cold-season inversions and summertime high temperatures associated with stagnant wind events. The identifiable mechanism of the changes to strong winds in VR-CESM-SONT provides confidence in its projections and demonstrates the value of enhanced spatial resolution for the study of extreme winds under climate change. Significance Statement: In this study, we compare climate change projections of high and low extreme wind speeds in the Great Lakes region between a standard coarse-resolution simulation and a high-resolution simulation performed using the same climate model. The fine-resolution simulation projects strengthening high wind speeds, opposite to the coarse-resolution simulation. Both project increasing frequency of extreme weak winds, but the human-health-related impacts of stagnant winds are only captured at fine resolution. The changes in the coarse-resolution simulation are explained by changes to large-scale circulation, while the fine-resolution changes are linked to local processes the coarse model does not resolve. This helps explain the diverging projections of strong winds and gives greater credibility to the fine-resolution simulation. [ABSTRACT FROM AUTHOR]

Current global reanalyses show marked discrepancies in snow mass and snow cover extent for the Northern Hemisphere. Here, benchmark snow datasets are produced by driving a simple offline snow model, the Brown Temperature Index Model (B-TIM), with temperature and precipitation from each of three reanalyses. B-TIM offline snow performs comparably to or better than online (coupled land-atmosphere) reanalysis snow when evaluated against in situ snow measurements. Sources of discrepancy in snow climatologies, which are difficult to isolate when comparing online reanalysis snow products amongst themselves, are partially elucidated by separately bias-adjusting temperature and precipitation in B-TIM. Interannual variability in snow mass and snow spatial patterns is far more self-consistent amongst offline B-TIM snow products than amongst online reanalysis snow products, and specific artifacts related to temporal inhomogeneity in snow data assimilation are revealed in the analysis. B-TIM, released here as an open-source, self-contained Python package, provides a simple benchmarking tool for future updates to more sophisticated online and offline snow datasets. [ABSTRACT FROM AUTHOR]

The European Space Agency's Aeolus mission, launched in August 2018, provides the first global horizontal line‐of‐sight wind profile measurements. Previous studies have shown that Aeolus winds in global forecast systems improve the overall forecast skill, especially in the upper tropospheric tropics and in other data‐sparse regions. In this study, we use a series of observing system experiments with the latest version of the reprocessed Aeolus wind product (2B11) to better characterize the locations and drivers of improved skill from Aeolus with Environment and Climate Change Canada's Global Deterministic Prediction System. Observing system experiments that test the impact of Aeolus winds and the impact of all operational wind observations are carried out, covering the period summer 2019 and winter 2019–2020. Assimilation of operational winds improves the tropospheric wind forecast over the Tropics by a reduction of 8% in the forecast error, and adding the Aeolus winds to the assimilations results in an extra 0.7–0.9%. Aeolus wind impacts are improvements are 0.7–0.9% for the Arctic, and 0.4–0.6% over the Northern and Southern Hemisphere extratropics. The scale dependence of these impacts is investigated using spatial spectra (spherical harmonic decomposition). The improvement is quantified using the difference of the 250 hPa kinetic energy forecast error spectra between experiments. The operational winds largely improve the forecast of planetary scale to intermediate scale for spherical wave numbers between 1 and 20 in the short‐range forecasts. The operational wind impact decreases as the forecast lead time increases. On the other hand, the impact of Aeolus is mostly seen in the intermediate to large scale range with a peak around spherical wave number 8. The Aeolus‐related improvement around this wave number increases with forecast lead time. This analysis suggests that Aeolus winds provide estimates of the wind state that are valuable and complementary to that provided from current operational winds. [ABSTRACT FROM AUTHOR]

Extreme near-surface wind speeds in cities can have major societal impacts but are not well represented in climate models. Despite this, large-scale dynamics in the free troposphere, which models resolve better, could provide reliable constraints on local extreme winds. This study identifies synoptic circulations associated with midlatitude extreme wind events and assesses how resolution affects their representation in analysis products and a climate model framework. Composites of reanalysis (ERA5) sea level pressure and upper-tropospheric winds during observed extreme wind events reveal distinct circulation structures for each quadrant of the surface-wind rose. Enhanced resolution of the analysis product (ERA5 versus the higher-resolution ECMWF Operational Analysis) reduced wind speed biases but has little impact on capturing occurrences of wind extremes seen in station observations. Composite circulations for surface wind extremes in a climate model (CESM) skillfully reproduce circulations found in reanalysis. Regional refinement of CESM over a region centered on southern Ontario, Canada, using variable resolution (VR-CESM) improves representation of surface ageostrophic circulations and the strength of vertical coupling between upper-level and near-surface winds. We thus can distinguish situations for which regional refinement (dynamical downscaling) is necessary for realistic representation of the large-scale atmospheric circulations associated with extreme winds, from situations where the coarse resolution of standard GCMs is sufficient. [ABSTRACT FROM AUTHOR]

Evaporation adds moisture to the atmosphere, while condensation removes it. Condensation also adds thermal energy to the atmosphere, which must be removed from the atmosphere by radiative cooling. As a result of these two processes, there is a net flow of energy driven by surface evaporation adding energy and radiative cooling removing energy from the atmosphere. Here, we calculate the implied heat transport of this process to find the atmospheric heat transport in balance with the surface evaporation. In modern-day Earth-like climates, evaporation varies strongly between the equator and the poles, while the net radiative cooling in the atmosphere is nearly meridionally uniform, and as a consequence, the heat transport governed by evaporation is similar to the total poleward heat transport of the atmosphere. This analysis is free from cancellations between moist and dry static energy transports, which greatly simplifies the interpretation of atmospheric heat transport and its relationship to the diabatic heating and cooling that governs the atmospheric heat transport. We further demonstrate, using a hierarchy of models, that much of the response of atmospheric heat transport to perturbations, including increasing CO2 concentrations, can be understood from the distribution of evaporation changes. These findings suggest that meridional gradients in surface evaporation govern atmospheric heat transport and its changes. [ABSTRACT FROM AUTHOR]

Snow on Arctic sea ice has many, contrasting effects on ice thickness and extent. Furthermore, estimates of snow depth on Arctic sea ice are a key input for ice thickness estimates from satellite altimeters such as ICESat‐2. Models such as the NASA Eulerian Snow on Sea Ice Model (NESOSIM) have been recently utilized by the sea ice community to provide time‐varying basin‐wide estimates of snow depth and density on Arctic sea ice. NESOSIM is a two‐snow‐layer model with simple representations of snow accumulation, wind packing, loss due to blowing snow, and redistribution due to sea ice motion. Two free parameters in NESOSIM, which dictate the bulk effect of wind packing (densification) and blowing snow processes, lack direct observational constraints. We present an indirect calibration of these parameters using a Markov Chain Monte Carlo (MCMC) approach. NESOSIM output is calibrated to observations of snow depth from Operation IceBridge and CRREL‐Dartmouth buoys, and density from historical drifting stations. OIB measurements alone are found to more strictly constrain the blowing snow parameter, and including additional observations yields more physically reasonable density estimates. The MCMC‐calibrated model output is further used to estimate sea ice thickness and uncertainty from model parameter uncertainty using ICESat‐2 freeboard measurements. Despite visible differences in density, the change in ice thickness is minimal. We also find that the model is relatively insensitive to parameter variations, and hence, the snow model uncertainty contribution to ice thickness is small compared to the systematic uncertainty from snow in the current ICESat‐2 thickness product. Key Points: Free parameters in a snow‐on‐sea‐ice model are automatically calibrated to observations using a Markov Chain Monte Carlo (MCMC) methodThe MCMC‐optimized model configuration yields small snow depth changes, but larger snow density changes compared to the prior configurationThe impact of snow model parameter calibration on sea ice thickness and uncertainty is minimal [ABSTRACT FROM AUTHOR]

In August 2018, the European Space Agency (ESA) launched the Aeolus satellite, whose Atmospheric LAser Doppler INstrument (ALADIN) is the first space-borne Doppler wind lidar to regularly measure vertical profiles of horizontal line-of-sight (HLOS) winds with global sampling. This mission is intended to assess improvement to numerical weather prediction provided by wind observations in regions poorly constrained by atmospheric mass, such as the tropics, but also, potentially, in polar regions such as the Arctic where direct wind observations are especially sparse. There remain gaps in the evaluation of the Aeolus products over the Arctic region, which is the focus of this contribution. Here, an assessment of the Aeolus Level-2B (L2B) wind product is carried out, progressing from specific locations in the Canadian North to the pan-Arctic. In particular, Aeolus data are compared to a limited sample of coincident ground-based Ka-band radar measurements at Iqaluit, Nunavut, to a larger set of coincident radiosonde measurements over the Canadian North, to Environment and Climate Change Canada (ECCC)'s short-range forecast, and to the reanalysis product, ERA5, from the European Centre for Medium-Range Weather Forecasts (ECMWF). Periods covered include the early phase of the first laser flight model (flight model A – FM-A; September to October 2018), the early phase of the second laser flight model (flight model B – FM-B; August to September 2019), and the middle phase of FM-B (December 2019 to January 2020). The adjusted r -squared between Aeolus and other local datasets is around 0.9 except for lower values for the comparison to the Ka-band radar, reflecting limited sampling opportunities with the radar data. This consistency is degraded by about 10 % for the Rayleigh winds in the summer due to solar background noise and other possible errors. Over the pan-Arctic, consistency, with correlation greater than 0.8, is found in the Mie channel from the planetary boundary layer to the lower stratosphere (near surface to 16 km a.g.l.) and in the Rayleigh channel from the troposphere to the stratosphere (2 to 25 km a.g.l.). In all three periods, Aeolus standard deviations are found to be 5 % to 40 % greater than those from ECCC-B and ERA5. We found that the L2B estimated error product for Aeolus is coherent with the differences between Aeolus and the other datasets and can be used as a guide for expected consistency. Our work shows that the high quality of the Aeolus dataset that has been demonstrated globally applies to the sparsely sampled Arctic region. It also demonstrates the lack of available independent wind measurements in the Canadian North, lending urgency to the need to augment the observing capacity in this region to ensure suitable calibration and validation of future space-borne Doppler wind lidar (DWL) missions. [ABSTRACT FROM AUTHOR]

To assess the effect of ocean‐atmosphere coupling in the climate response to forced sea ice loss, the Polar Amplification Model Intercomparison Project protocol includes centennial coupled atmosphere‐ocean general circulation model simulations with imposed sea ice loss. The protocol, which specifies sea ice concentration and thickness distribution targets, does not prescribe a method for achieving them. Although different methods for imposing sea ice loss (or growth) in models have been documented, testing of the method‐dependence of the resulting climate responses has been limited. Achieving the targeted sea ice state has proven to be challenging using the so‐called ghost‐flux nudging method, which induces ice melt from below, as this method does not constrain the partitioning between thickness and concentration. We propose, describe and test a simple method that combines the advantages of direct sea ice nudging and ghost‐flux nudging. The hybrid nudging method better captures the partitioning between thickness and concentration while conserving total water content. We document how this novel sea ice constraining method reaches specific targets, enhances surface turbulent heat flux responses to sea ice loss, and induces tropospheric warming for both polar regions. Plain Language Summary: The Arctic is warming faster than the global average due to several processes that, once combined, lead to so‐called Arctic Amplification. Part of this anomalous polar warming comes from an intense reduction in ice cover allowing heat into the ocean, warming the Arctic ocean near the surface, hence melting more ice. A joint effort by several climate modeling groups called the Polar Amplification Model Intercomparison Project (PAMIP) aims at better understanding Arctic Amplification through a coordinated set of climate simulations. Among this ensemble of simulations is a set of centennial simulations performed with fully coupled state‐of‐the‐art climate models. In these experiments, Arctic (and Antarctic) sea ice are forced to reach specific states in order to better isolate Arctic Amplification and sea ice loss from the rest of anthropogenic global warming. In this paper, we propose a simple technique to nudge sea ice models to specific states such as prescribed by PAMIP. This new method combines advantages from existing techniques to improve the control over the extent and the thickness of the ice. We document how our novel method produces atmospheric warming and circulation responses that are more consistent with previously published results than the existing techniques. Key Points: Already established sea ice nudging method leads to issues in the partitioning between sea ice concentration (SIC) and sea ice thickness (SIT)New hybrid nudging method significantly improves capture of both SIC and SIT targets for Polar Amplification Model Intercomparison Project, in particular for Arctic sea iceHybrid scheme increases Arctic turbulent heat flux and warming [ABSTRACT FROM AUTHOR]

Stationary wave nonlinearity describes the self-interaction of stationary waves and is important in maintaining the zonally asymmetric atmospheric general circulation. However, the dynamics of stationary wave nonlinearity, which is often calculated explicitly in stationary wave models, is not well understood. Stationary wave nonlinearity is examined here in the simplified setting of the response to localized topographic forcing in quasigeostrophic barotropic dynamics in the presence and absence of transient eddies. It is shown that stationary wave nonlinearity accounts for most of the difference between the linear and full nonlinear response, particularly if the adjustment of the zonal-mean flow to the stationary waves is taken into account. The separate impact of transient eddy forcing is also quantified. Wave activity analysis shows that stationary wave nonlinearity in this setting is associated with Rossby wave critical layer reflection. A nonlinear stationary wave model, similar to those used in baroclinic stationary wave model studies, is also tested and is shown to capture the basic features of the full nonlinear stationary wave solution. DOI: 10.1175/2010JAS3332.1

Sahel summertime precipitation declined from the 1950s to 1970s and recovered from the 1970s to 2000s. Anthropogenic aerosol contributions to this evolution are typically attributed to interhemispheric gradient changes of Atlantic Ocean sea surface temperature (SST). However recent work by Hirasawa et al. indicates a more complex picture, with the response being a combination of "fast" direct atmospheric (DA) processes and "slow" ocean-mediated (OM) processes. Here, we extend this understanding using the Community Atmosphere Model 5 to determine the role of regional ocean-basin perturbations and regional aerosol emission changes in the overall aerosol-driven OM and DA responses, respectively. From the 1950s to 1970s, there was an OM Sahel wetting response due to Pacific Ocean cooling that was offset by drying due to Atlantic cooling. By contrast, from the 1970s to 2000s, Atlantic trends reversed and amplified the Pacific cooling-induced wetting. This wetting was partially offset by drying driven by Indian Ocean cooling. Thus, the OM Sahel precipitation response to aerosol crucially depends on the balance of responses to Atlantic, Pacific, and Indian Ocean SST anomalies. From the 1950s to 1970s, there is DA Sahel drying that was principally due to North American aerosol emissions, with negligible effect from European emissions. DA drying from the 1970s to 2000s was mainly due to African aerosol emissions. Thus, the shifting roles of regional OM and DA effects reveal a complex interplay of direct driving and remote teleconnections in determining the time evolution of Sahel precipitation due to aerosol forcing in the late twentieth century. Significance Statement: Studies of global climate models consistently indicate that anthropogenic aerosol emissions were a significant contributor to a severe drought that occurred in the Sahel region of Africa in the late twentieth century. The drying influence of aerosol forcing is the combined result of rapid atmospheric responses directly due to the forcing and slower responses due to forced ocean temperature changes. Using a set of simulations targeted at determining the influences from different ocean basins and different emission regions for two periods in the late twentieth century, we find there is a surprising range of mechanisms through which aerosol emissions affect the Sahel. This results in a complex interplay of at times competing and at times complementary regional influences. [ABSTRACT FROM AUTHOR]

Copyright of Atmosphere -- Ocean (Taylor & Francis Ltd) is the property of Taylor & Francis Ltd and its content may not be copied or emailed to multiple sites or posted to a listserv without the copyright holder's express written permission. However, users may print, download, or email articles for individual use. This abstract may be abridged. No warranty is given about the accuracy of the copy. Users should refer to the original published version of the material for the full abstract. (Copyright applies to all Abstracts.)

A simple atmospheric general circulation model (GCM) is used to investigate the transient response of the stratosphere-troposphere system to externally imposed pulses of lower-tropospheric planetary wave activity. The atmospheric GCM is a dry, hydrostatic, global primitive-equations model, whose circulation includes an active polar vortex and a tropospheric jet maintained by baroclinic eddies. Planetary wave activity pulses are generated by a perturbation of the solid lower boundary that grow and decay over a period of 10 days. The planetary wave pulses propagate upward and break in the stratosphere. Subsequently, a zonal-mean circulation anomaly propagates downward, often into the troposphere, at lags of 30-100 days. The evolution of the response is found to be dependent on the state of the stratosphere-troposphere system at the time the pulse is generated. In particular, on the basis of a large ensemble of these simulations, it is found that the length of time the signal takes to propagate downward from the stratosphere is controlled by initial anomalies in the zonal-mean circulation and in the zonal-mean wave drag. Criteria based on these anomaly patterns can be used, therefore, to predict the long-term surface response of the stratosphere troposphere system to a planetary wave pulse up to 90 days after the pulse is generated. In an independent test, it is verified that the initial states that most strongly satisfy these criteria respond in the expected way to the lower-tropospheric wave activity pulse.

An exceptionally strong stratospheric sudden warming (SSW) that spontaneously occurs in a very simple stratosphere--troposphere AGCM is discussed. The model is a dry, hydrostatic, primitive equation model without planetary stationary waves. Transient baroclinic wave-wave interaction in the troposphere thus provides the only source of upward-propagating wave activity into the stratosphere. The model's SSW is grossly similar to the Southern Hemisphere major SSW of 2002: it occurs after weaker warmings 'precondition' the polar vortex for breaking, it involves a split of the polar vortex, and it has a downward-propagating signature. These similarities suggest that the Southern Hemisphere SSW of 2002 might itself have been caused by transient baroclinic wave-wave interaction. The simple model used for this study also provides some insight into how often such extreme events might occur. The frequency distribution of SSWs in the model has exponential, as opposed to Gaussian, tails. This suggests that very large amplitude SSWs, though rare, might occur with higher frequency than might be naively expected.

The influence of zonally asymmetric boundary conditions on the leading modes of variability in a suite of atmospheric general circulation models is investigated. The set of experiments consists of nine model configurations, with varying degrees and types of zonal asymmetry in their boundary conditions. The structure of the leading EOF varies with the zonal asymmetry of the base state for each model configuration. In particular, a close relationship is found between the structure of the EOF and the model storm tracks. An approximately linear relationship is found to hold between the magnitude of the zonal asymmetry of the leading EOF and of the storm tracks in the models. It is shown that this linear relationship extends to the observations. One-point correlation maps centered on the regions where the EOFs reach their maximum amplitude show similar structures for all configurations. These structures consist of a north-south dipole, resembling the observed structure of the North Atlantic Oscillation (NAO). They are significantly more zonally localized than the leading EOF, but do resemble one-point correlation maps and sector EOFs calculated for a simulation with zonally symmetric boundary conditions. Thus, the leading EOF for each simulation appears to represent the longitudinal distribution of zonally localized NAO-like patterns. This longitudinal distribution appears to be tied to the distribution of high-frequency eddies, as represented by the storm tracks. A detailed momentum budget for each case confirms that high-frequency eddies play a crucial role in producing the NAO-like patterns. Other dynamical processes also play an important role, but vary with the details of the simulation.

A simple dynamical model is presented for the basic spatial and temporal structure of the large-scale modes of intraseasonal variability and associated variations in the zonal index. Such variability in the extratropical atmosphere is known to be represented by fairly well-defined patterns, and among the most prominent are the North Atlantic Oscillation (NAO) and a more zonally symmetric pattern known as an annular mode, which is most pronounced in the Southern Hemisphere. These patterns may be produced by the momentum fluxes associated with large-scale midlatitude stirring, such as that provided by baroclinic eddies. It is shown how such stirring, as represented by a simple stochastic forcing in a barotropic model, leads to a variability in the zonal flow via a variability in the eddy momentum flux convergence and to patterns similar to those observed. Typically, the leading modes of variability may be characterized as a mixture of 'wobbles' in the zonal jet position and 'pulses' in the zonal jet strength. If the stochastic forcing is statistically zonally uniform, then the resulting patterns of variability as represented by empirical orthogonal functions are almost zonally uniform and the pressure pattern is dipolar in the meridional direction, resembling an annular mode. If the forcing is enhanced in a zonally localized region, thus mimicking the effects of a storm track over the ocean, then the resulting variability pattern is zonally localized, resembling the North Atlantic Oscillation. This suggests that the North Atlantic Oscillation and annular modes are produced by the same mechanism and are manifestations of the same phenomenon. The time scale of variability of the patterns is longer than the decorrelation time scale of the stochastic forcing, because of the temporal integration of the forcing by the equations of motion limited by the effects of nonlinear dynamics and friction. For reasonable parameters these produce a decorrelation time of the order of 5-10 days. The model also produces some long-term (100 days or longer) variability, without imposing such variability via the external parameters except insofar as it is contained in the nearly white stochastic forcing.

Analyzing a multimodel ensemble of coupled climate model simulations forced with Arctic sea ice loss using a two-parameter pattern-scaling technique to remove the cross-coupling between low- and high-latitude responses, the sensitivity to high-latitude sea ice loss is isolated and contrasted to the sensitivity to low-latitude warming. Despite some differences in experimental design, the Northern Hemisphere near-surface atmospheric sensitivity to sea ice loss is found to be robust across models in the cold season; however, a larger intermodel spread is found at the surface in boreal summer, and in the free tropospheric circulation. In contrast, the sensitivity to low-latitude warming is most robust in the free troposphere and in the warm season, with more intermodel spread in the surface ocean and surface heat flux over the Northern Hemisphere. The robust signals associated with sea ice loss include upward turbulent and longwave heat fluxes where sea ice is lost, warming and freshening of the Arctic Ocean, warming of the eastern North Pacific Ocean relative to the western North Pacific with upward turbulent heat fluxes in the Kuroshio Extension, and salinification of the shallow shelf seas of the Arctic Ocean alongside freshening in the subpolar North Atlantic Ocean. In contrast, the robust signals associated with low-latitude warming include intensified ocean warming and upward latent heat fluxes near the western boundary currents, freshening of the Pacific Ocean, salinification of the North Atlantic, and downward sensible and longwave fluxes over the ocean. [ABSTRACT FROM AUTHOR]