The ability of CMIP5 models to simulate North Atlantic extratropical cyclones

The ability of the climate models participating in phase 5 of the Coupled Model Intercomparison Project (CMIP5) to simulate North Atlantic extratropical cyclones in winter [December–February (DJF)] and summer [June–August (JJA)] is investigated in detail. Cyclones are identified as maxima in T42 vorticity at 850 hPa and their propagation is tracked using an objective feature-tracking algorithm. By comparing the historical CMIP5 simulations (1976–2005) and the ECMWF Interim Re-Analysis (ERA-Interim; 1979–2008), the authors find that systematic biases affect the number and intensity of North Atlantic cyclones in CMIP5 models. In DJF, the North Atlantic storm track tends to be either too zonal or displaced southward, thus leading to too few and weak cyclones over the Norwegian Sea and too many cyclones in central Europe. In JJA, the position of the North Atlantic storm track is generally well captured but some CMIP5 models underestimate the total number of cyclones. The dynamical intensity of cyclones, as measured by either T42 vorticity at 850 hPa or mean sea level pressure, is too weak in both DJF and JJA. The intensity bias has a hemispheric character, and it cannot be simply attributed to the representation of the North Atlantic large- scale atmospheric state. Despite these biases, the representation of Northern Hemisphere (NH) storm tracks has improved since CMIP3 and some CMIP5 models are able of representing well both the number and the intensity of North Atlantic cyclones. In particular, some of the higher-atmospheric-resolution models tend to have a better representation of the tilt of the North Atlantic storm track and of the intensity of cyclones in DJF.

A multimodel assessment of future projections of North Atlantic and European extratropical cyclones in the CMIP5 climate models

The response of North Atlantic and European extratropical cyclones to climate change is investigated in the climate models participating in phase 5 of the Coupled Model Intercomparison Project (CMIP5). In contrast to previous multimodel studies, a feature-tracking algorithm is here applied to separately quantify the re- sponses in the number, the wind intensity, and the precipitation intensity of extratropical cyclones. Moreover, a statistical framework is employed to formally assess the uncertainties in the multimodel projections. Under the midrange representative concentration pathway (RCP4.5) emission scenario, the December–February (DJF) response is characterized by a tripolar pattern over Europe, with an increase in the number of cyclones in central Europe and a decreased number in the Norwegian and Mediterranean Seas. The June–August (JJA) response is characterized by a reduction in the number of North Atlantic cyclones along the southern flank of the storm track. The total number of cyclones decreases in both DJF (24%) and JJA (22%). Classifying cyclones according to their intensity indicates a slight basinwide reduction in the number of cy- clones associated with strong winds, but an increase in those associated with strong precipitation. However, in DJF, a slight increase in the number and intensity of cyclones associated with strong wind speeds is found over the United Kingdom and central Europe. The results are confirmed under the high-emission RCP8.5 scenario, where the signals tend to be larger. The sources of uncertainty in these projections are discussed.

Impact of atmospheric forcing on Antarctic continental shelf water masses

The Antarctic continental shelf seas feature a bimodal distribution of water mass temperature, with the Amundsen and Bellingshausen Seas flooded by Circumpolar Deep Water that is several degrees Celsius warmer than the cold shelf waters prevalent in the Weddell and Ross Seas. This bimodal distribution could be caused by differences in atmospheric forcing, ocean dynamics, ocean and ice feedbacks, or some combination of these factors. In this study, a highly simplified coupled sea ice–mixed layer model is developed to investigate the physical processes controlling this situation. Under regional atmospheric forcings and parameter choices the 10-yr simulations demonstrate a complete destratification of the Weddell Sea water column in winter, forming cold, relatively saline shelf waters, while the Amundsen Sea winter mixed layer remains shallower, allowing a layer of deep warm water to persist. Applying the Weddell atmospheric forcing to the Amundsen Sea model destratifies the water column after two years, and applying the Amundsen forcing to the Weddell Sea model results in a shallower steady-state winter mixed layer that no longer destratifies the water column. This suggests that the regional difference in atmospheric forcings alone is sufficient to account for the bimodal distribution in Antarctic shelf-sea temperatures. The model prediction of mixed layer depth is most sensitive to the air temperature forcing, but a switch in all forcings is required to prevent destratification of the Weddell Sea water column.

Rheology of discrete failure regimes of anisotropic sea ice

A rheological model of sea ice is presented that incorporates the orientational distribution of ice thickness in leads embedded in isotropic floe ice. Sea ice internal stress is determined by coulombic, ridging and tensile failure at orientations where corresponding failure criteria are satisfied at minimum stresses. Because sea ice traction increases in thinner leads and cohesion is finite, such failure line angles are determined by the orientational distribution of sea ice thickness relative to the imposed stresses. In contrast to the isotropic case, sea ice thickness anisotropy results in these failure lines becoming dependent on the stress magnitude. Although generally a given failure criteria type can be satisfied at many directions, only two at most are considered. The strain rate is determined by shearing along slip lines accompanied by dilatancy and closing or opening across orientations affected by ridging or tensile failure. The rheology is illustrated by a yield curve determined by combining coulombic and ridging failure for the case of two pairs of isotropically formed leads of different thicknesses rotated with regard to each other, which models two events of coulombic failure followed by dilatancy and refreezing. The yield curve consists of linear segments describing coulombic and ridging yield as failure switches from one lead to another as the stress grows. Because sliding along slip lines is accompanied by dilatancy, at typical Arctic sea ice deformation rates a one-day-long deformation event produces enough open water that these freshly formed slip lines are preferential places of ridging failure.

[Global climate] global aerosols [in state of the climate in 2008]

The global mean temperature in 2008 was slightly cooler than that in 2007; however, it still ranks within the 10 warmest years on record. Annual mean temperatures were generally well above average in South America, northern and southern Africa, Iceland, Europe, Russia, South Asia, and Australia. In contrast, an exceptional cold outbreak occurred during January across Eurasia and over southern European Russia and southern western Siberia. There has been a general increase in land-surface temperatures and in permafrost temperatures during the last several decades throughout the Arctic region, including increases of 1° to 2°C in the last 30 to 35 years in Russia. Record setting warm summer (JJA) air temperatures were observed throughout Greenland. The year 2008 was also characterized by heavy precipitation in a number of regions of northern South America, Africa, and South Asia. In contrast, a prolonged and intense drought occurred during most of 2008 in northern Argentina, Paraguay, Uruguay, and southern Brazil, causing severe impacts to agriculture and affecting many communities. The year began with a strong La Niña episode that ended in June. Eastward surface current anomalies in the tropical Pacific Ocean in early 2008 played a major role in adjusting the basin from strong La Niña conditions to ENSO-neutral conditions by July–August, followed by a return to La Niña conditions late in December. The La Niña conditions resulted in far-reaching anomalies such as a cooling in the central tropical Pacific, Arctic Ocean, and the regions extending from the Gulf of Alaska to the west coast of North America; changes in the sea surface salinity and heat content anomalies in the tropics; and total column water vapor, cloud cover, tropospheric temperature, and precipitation patterns typical of a La Niña. Anomalously salty ocean surface salinity values in climatologically drier locations and anomalously fresh values in rainier locations observed in recent years generally persisted in 2008, suggesting an increase in the hydrological cycle. The 2008 Atlantic hurricane season was the 14th busiest on record and the only season ever recorded with major hurricanes each month from July through November. Conversely, activity in the northwest Pacific was considerably below normal during 2008. While activity in the north Indian Ocean was only slightly above average, the season was punctuated by Cyclone Nargis, which killed over 145,000 people; in addition, it was the seventh-strongest cyclone ever in the basin and the most devastating to hit Asia since 1991. Greenhouse gas concentrations continued to rise, increasing by more than expected based on with CO2 the 1979 to 2007 trend. In the oceans, the global mean uptake for 2007 is estimated to be 1.67 Pg-C, about CO2 0.07 Pg-C lower than the long-term average, making it the third-largest anomaly determined with this method since 1983, with the largest uptake of carbon over the past decade coming from the eastern Indian Ocean. Global phytoplankton chlorophyll concentrations were slightly elevated in 2008 relative to 2007, but regional changes were substantial (ranging to about 50%) and followed long-term patterns of net decreases in chlorophyll with increasing sea surface temperature. Ozone-depleting gas concentrations continued to fall globally to about 4% below the peak levels of the 2000–02 period. Total column ozone concentrations remain well below pre-1980, levels and the 2008 ozone hole was unusually large (sixth worst on record) and persistent, with low ozone values extending into the late December period. In fact the polar vortex in 2008 persisted longer than for any previous year since 1979. Northern Hemisphere snow cover extent for the year was well below average due in large part to the record-low ice extent in March and despite the record-maximum coverage in January and the shortest snow cover duration on record (which started in 1966) in the North American Arctic. Limited preliminary data imply that in 2008 glaciers continued to lose mass, and full data for 2007 show it was the 17th consecutive year of loss. The northern region of Greenland and adjacent areas of Arctic Canada experienced a particularly intense melt season, even though there was an abnormally cold winter across Greenland's southern half. One of the most dramatic signals of the general warming trend was the continued significant reduction in the extent of the summer sea-ice cover and, importantly, the decrease in the amount of relatively older, thicker ice. The extent of the 2008 summer sea-ice cover was the second-lowest value of the satellite record (which started in 1979) and 36% below the 1979–2000 average. Significant losses in the mass of ice sheets and the area of ice shelves continued, with several fjords on the northern coast of Ellesmere Island being ice free for the first time in 3,000–5,500 years. In Antarctica, the positive phase of the SAM led to record-high total sea ice extent for much of early 2008 through enhanced equatorward Ekman transport. With colder continental temperatures at this time, the 2007–08 austral summer snowmelt season was dramatically weakened, making it the second shortest melt season since 1978 (when the record began). There was strong warming and increased precipitation along the Antarctic Peninsula and west Antarctica in 2008, and also pockets of warming along coastal east Antarctica, in concert with continued declines in sea-ice concentration in the Amundsen/Bellingshausen Seas. One significant event indicative of this warming was the disintegration and retreat of the Wilkins Ice Shelf in the southwest peninsula area of Antarctica.

Creation of the WATCH Forcing Data and Its Use to Assess Global and Regional Reference Crop Evaporation over Land during the Twentieth Century

The Water and Global Change (WATCH) project evaluation of the terrestrial water cycle involves using land surface models and general hydrological models to assess hydrologically important variables including evaporation, soil moisture, and runoff. Such models require meteorological forcing data, and this paper describes the creation of the WATCH Forcing Data for 1958–2001 based on the 40-yr ECMWF Re-Analysis (ERA-40) and for 1901–57 based on reordered reanalysis data. It also discusses and analyses modelindependent estimates of reference crop evaporation. Global average annual cumulative reference crop evaporation was selected as a widely adopted measure of potential evapotranspiration. It exhibits no significant trend from 1979 to 2001 although there are significant long-term increases in global average vapor pressure deficit and concurrent significant decreases in global average net radiation and wind speed. The near-constant global average of annual reference crop evaporation in the late twentieth century masks significant decreases in some regions (e.g., the Murray–Darling basin) with significant increases in others.

Generation of a buoyancy-driven coastal current by an Antarctic Polynya

Descent and spreading of high salinity water generated by salt rejection during sea ice formation in an Antarctic coastal polynya is studied using a hydrostatic, primitive equation three-dimensional ocean model called the Proudman Oceanographic Laboratory Coastal Ocean Modeling System (POLCOMS). The shape of the polynya is assumed to be a rectangle 100 km long and 30 km wide, and the salinity flux into the polynya at its surface is constant. The model has been run at high horizontal spatial resolution (500 m), and numerical simulations reveal a buoyancy-driven coastal current. The coastal current is a robust feature and appears in a range of simulations designed to investigate the influence of a sloping bottom, variable bottom drag, variable vertical turbulent diffusivities, higher salinity flux, and an offshore position of the polynya. It is shown that bottom drag is the main factor determining the current width. This coastal current has not been produced with other numerical models of polynyas, which may be because these models were run at coarser resolutions. The coastal current becomes unstable upstream of its front when the polynya is adjacent to the coast. When the polynya is situated offshore, an unstable current is produced from its outset owing to the capture of cyclonic eddies. The effect of a coastal protrusion and a canyon on the current motion is investigated. In particular, due to the convex shape of the coastal protrusion, the current sheds a dipolar eddy.

Contribution of land use changes to near-surface air temperatures during recent summer extreme heat events in the Phoenix Metropolitan Area

The impact of 1973–2005 land use–land cover (LULC) changes on near-surface air temperatures during four recent summer extreme heat events (EHEs) are investigated for the arid Phoenix, Arizona, metropolitan area using the Weather Research and Forecasting Model (WRF) in conjunction with the Noah Urban Canopy Model. WRF simulations were carried out for each EHE using LULC for the years 1973, 1985, 1998, and 2005. Comparison of measured near-surface air temperatures and wind speeds for 18 surface stations in the region show a good agreement between observed and simulated data for all simulation periods. The results indicate consistent significant contributions of urban development and accompanying LULC changes to extreme temperatures for the four EHEs. Simulations suggest new urban developments caused an intensification and expansion of the area experiencing extreme temperatures but mainly influenced nighttime temperatures with an increase of up to 10 K. Nighttime temperatures in the existing urban core showed changes of up to 2 K with the ongoing LULC changes. Daytime temperatures were not significantly affected where urban development replaced desert land (increase by 1 K); however, maximum temperatures increased by 2–4 K when irrigated agricultural land was converted to suburban development. According to the model simulations, urban landscaping irrigation contributed to cooling by 0.5–1 K in maximum daytime as well as minimum nighttime 2-m air temperatures in most parts of the urban region. Furthermore, urban development led to a reduction of the already relatively weak nighttime winds and therefore a reduction in advection of cooler air into the city.

The effects of rotation and ice shelf topography on frazil-laden ice shelf water plumes

A model of the dynamics and thermodynamics of a plume of meltwater at the base of an ice shelf is presented. Such ice shelf water plumes may become supercooled and deposit marine ice if they rise (because of the pressure decrease in the in situ freezing temperature), so the model incorporates both melting and freezing at the ice shelf base and a multiple-size-class model of frazil ice dynamics and deposition. The plume is considered in two horizontal dimensions, so the influence of Coriolis forces is incorporated for the first time. It is found that rotation is extremely influential, with simulated plumes flowing in near-geostrophy because of the low friction at a smooth ice shelf base. As a result, an ice shelf water plume will only rise and become supercooled (and thus deposit marine ice) if it is constrained to flow upslope by topography. This result agrees with the observed distribution of marine ice under Filchner–Ronne Ice Shelf, Antarctica. In addition, it is found that the model only produces reasonable marine ice formation rates when an accurate ice shelf draft is used, implying that the characteristics of real ice shelf water plumes can only be captured using models with both rotation and a realistic topography.

Optimization of a sea ice model using basinwide observations of Arctic sea ice thickness, extent, and velocity

A stand-alone sea ice model is tuned and validated using satellite-derived, basinwide observations of sea ice thickness, extent, and velocity from the years 1993 to 2001. This is the first time that basin-scale measurements of sea ice thickness have been used for this purpose. The model is based on the CICE sea ice model code developed at the Los Alamos National Laboratory, with some minor modifications, and forcing consists of 40-yr ECMWF Re-Analysis (ERA-40) and Polar Exchange at the Sea Surface (POLES) data. Three parameters are varied in the tuning process: Ca, the air–ice drag coefficient; P*, the ice strength parameter; and α, the broadband albedo of cold bare ice, with the aim being to determine the subset of this three-dimensional parameter space that gives the best simultaneous agreement with observations with this forcing set. It is found that observations of sea ice extent and velocity alone are not sufficient to unambiguously tune the model, and that sea ice thickness measurements are necessary to locate a unique subset of parameter space in which simultaneous agreement is achieved with all three observational datasets.

Interpretation of ageostrophic winds and implications for jet stream maintenance

The use of ageostrophic flow to infer the presence of vertical circulations in the entrances and exits of the climatological jet streams is questioned. Problems of interpretation arise because of the use of different definitions of geostrophy in theoretical studies and in analyses of atmospheric data. The nature and role of the ageostrophic flow based on constant and variable Coriolis parameter definitions of geostrophy vary. In the latter the geostrophic divergence cannot be neglected, so the vertical motion is not associated solely with the ageostrophic flow. Evidence is presented suggesting that ageostrophic flow in the climatological jet streams is primarily determined by the kinematic requirements of wave retrogression rather than by a forcing process. These requirements are largely met by the rotational flow, with the divergent circulations present being geostrophically forced, and so playing a secondary, restoring role.

Toward seamless prediction: calibration of climate change projections using seasonal forecasts

The CWRF is developed as a climate extension of the Weather Research and Forecasting model (WRF) by incorporating numerous improvements in the representation of physical processes and integration of external (top, surface, lateral) forcings that are crucial to climate scales, including interactions between land, atmosphere, and ocean; convection and microphysics; and cloud, aerosol, and radiation; and system consistency throughout all process modules. This extension inherits all WRF functionalities for numerical weather prediction while enhancing the capability for climate modeling. As such, CWRF can be applied seamlessly to weather forecast and climate prediction. The CWRF is built with a comprehensive ensemble of alternative parameterization schemes for each of the key physical processes, including surface (land, ocean), planetary boundary layer, cumulus (deep, shallow), microphysics, cloud, aerosol, and radiation, and their interactions. This facilitates the use of an optimized physics ensemble approach to improve weather or climate prediction along with a reliable uncertainty estimate. The CWRF also emphasizes the societal service capability to provide impactrelevant information by coupling with detailed models of terrestrial hydrology, coastal ocean, crop growth, air quality, and a recently expanded interactive water quality and ecosystem model. This study provides a general CWRF description and basic skill evaluation based on a continuous integration for the period 1979– 2009 as compared with that of WRF, using a 30-km grid spacing over a domain that includes the contiguous United States plus southern Canada and northern Mexico. In addition to advantages of greater application capability, CWRF improves performance in radiation and terrestrial hydrology over WRF and other regional models. Precipitation simulation, however, remains a challenge for all of the tested models.

Dependence of sea ice yield-curve shape on ice thickness

In this note, the authors discuss the contribution that frictional sliding of ice floes (or floe aggregates) past each other and pressure ridging make to the plastic yield curve of sea ice. Using results from a previous study that explicitly modeled the amount of sliding and ridging that occurs for a given global strain rate, it is noted that the relative contribution of sliding and ridging to ice stress depends upon ice thickness. The implication is that the shape and size of the plastic yield curve is dependent upon ice thickness. The yield-curve shape dependence is in addition to plastic hardening/weakening that relates the size of the yield curve to ice thickness. In most sea ice dynamics models the yield-curve shape is taken to be independent of ice thickness. The authors show that the change of the yield curve due to a change in the ice thickness can be taken into account by a weighted sum of two thickness-independent rheologies describing ridging and sliding effects separately. It would be straightforward to implement the thickness-dependent yield-curve shape described here into sea ice models used for global or regional ice prediction.

Influence of mass diffusion in sea ice dynamical models

The mixing of floes of different thickness caused by repeated deformation of the ice cover is modeled as diffusion, and the mass balance equation for sea ice accounting for mass diffusion is developed. The effect of deformational diffusion on the ice thickness balance is shown to reach 1% of the divergence effect, which describes ridging and lead formation. This means that with the same accuracy the mass balance equation can be written in terms of mean velocity rather than mean mass-weighted velocity, which one should correctly use for a multicomponent fluid such as sea ice with components identified by floe thickness. Mixing (diffusion) of sea ice also occurs because of turbulent variations in wind and ocean drags that are unresolved in models. Estimates of the importance of turbulent mass diffusion on the dynamic redistribution of ice thickness are determined using empirical data for the turbulent diffusivity. For long-time-scale prediction (≫5 days), where unresolved atmospheric motion may have a length scale on the order of the Arctic basin and the time scale is larger than the synoptic time scale of atmospheric events, turbulent mass diffusion can exceed 10% of the divergence effect. However, for short-time-scale prediction, for example, 5 days, the unresolved scales are on the order of 100 km, and turbulent diffusion is about 0.1% of the divergence effect. Because inertial effects are small in the dynamics of the sea ice pack, diffusive momentum transfer can be disregarded.

The effect of a new drag-law parameterization on ice shelf water plume dynamics

A drag law accounting for Ekman rotation adjacent to a flat, horizontal bou ndary is proposed for use in a plume model that is written in terms of the depth-mean velocity. The drag l aw contains a variable turning angle between the mean velocity and the drag imposed by the turbulent bound ary layer. The effect of the variable turning angle in the drag law is studied for a plume of ice shelf wat er (ISW) ascending and turning beneath an Antarctic ice shelf with draft decreasing away from the groundi ng line. As the ISW plume ascends the sloping ice shelf–ocean boundary, it can melt the ice shelf, wh ich alters the buoyancy forcing driving the plume motion. Under these conditions, the typical turning ang le is of order 10° over most of the plume area for a range of drag coefficients (the minus sign arises for th e Southern Hemisphere). The rotation of the drag with respect to the mean velocity is found to be signifi cant if the drag coefficient exceeds 0.003; in this case the plume body propagates farther along and across the b ase of the ice shelf than a plume with the standard quadratic drag law with no turning angle.