IMILAST: A community effort to intercompare extratropical cyclone detection and tracking algorithms: assessing method-related uncertainties

The variability of results from different automated methods of detection and tracking of extratropical cyclones is assessed in order to identify uncertainties related to the choice of method. Fifteen international teams applied their own algorithms to the same dataset - the period 1989-2009 of interim European Centre for Medium-Range Weather Forecasts (ECMWF) Re-Analysis (ERAInterim) data. This experiment is part of the community project Intercomparison of Mid Latitude Storm Diagnostics (IMILAST; see www.proclim.ch/imilast/index.html). The spread of results for cyclone frequency, intensity, life cycle, and track location is presented to illustrate the impact of using different methods. Globally, methods agree well for geographical distribution in large oceanic regions, interannual variability of cyclone numbers, geographical patterns of strong trends, and distribution shape for many life cycle characteristics. In contrast, the largest disparities exist for the total numbers of cyclones, the detection of weak cyclones, and distribution in some densely populated regions. Consistency between methods is better for strong cyclones than for shallow ones. Two case studies of relatively large, intense cyclones reveal that the identification of the most intense part of the life cycle of these events is robust between methods, but considerable differences exist during the development and the dissolution phases.

Extratropical fronts in the lower troposphere-global perspectives obtained from two automated methods

In this study two commonly used automated methods to detect atmospheric fronts in the lower troposphere are compared in various synoptic situations. The first method is a thermal approach, relying on the gradient of equivalent potential temperature (TH), while the second method is based on temporal changes in the 10 m wind (WND). For a comprehensive objective comparison of the outputs of these methods of frontal identification, both schemes are firstly applied to an idealised strong baroclinic wave simulation in the absence of topography. Then, two case-studies (one in the Northern Hemisphere (NH) and one in the Southern Hemisphere (SH)) were conducted to contrast fronts detected by the methods. Finally, we obtain global winter and summer frontal occurrence climatologies (derived from ERA-Interim for 1979–2012) and compare the structure of these. TH is able to identify cold and warm fronts in strong baroclinic cases that are in good agreement with manual analyses. WND is particularly suited for the detection of strongly elongated, meridionally oriented moving fronts, but has very limited ability to identify zonally oriented warm fronts. We note that the areas of the main TH frontal activity are shifted equatorwards compared to the WND patterns and are located upstream of regions of main WND front activity. The number of WND fronts in the NH shows more interseasonal variations than TH fronts, decreasing by more than 50% from winter to summer. In the SH there is a weaker seasonal variation of the number of observed WND fronts, however TH front activity reduces from summer (DJF) to winter (JJA). The main motivation is to give an overview of the performance of these methods, such that researchers can choose the appropriate one for their particular interest.

The role of upper-level dynamics and surface processes for the Pakistan flood of July 2010

In July and August 2010 floods of unprecedented impact afflicted Pakistan. The floods resulted from a series of intense multi-day precipitation events in July and early August. At the same time a series of blocking anticyclones dominated the upper-level flow over western Russia and breaking waves i.e. equatorward extrusions of stratospheric high potential vorticity (PV) air formed along the downstream flank of the blocks. Previous studies suggested that these extratropical upper-level breaking waves were crucial for instigating the precipitation events in Pakistan. Here a detailed analysis is provided of the extratropical forcing of the precipitation. Piecewise PV inversion is used to quantify the extratropical upper-level forcing associated with the wave breaking and trajectories are calculated to study the pathways and source regions of the moisture that precipitated over Pakistan. Limited-area model simulations are carried out to complement the Lagrangian analysis. The precipitation events over Pakistan resulted from a combination of favourable boundary conditions with strong extratropical and monsoonal forcing factors. Above-normal sea-surface temperatures in the Indian Ocean led to an elevated lower-tropospheric moisture content. Surface monsoonal depressions ensured the transport of moist air from the ocean towards northeastern Pakistan. Along this pathway the air parcel humidity increased substantially (60–90% of precipitated moisture) via evapotranspiration from the land surface. Extratropical breaking waves influenced the surface wind field substantially by enhancing the wind component directed towards the mountains which reinforced the precipitation.

Land use change and nitrogen feedbacks constrain the trajectory of the land carbon sink

Our understanding of Earth's carbon climate system depends critically upon interactions between rising atmospheric CO2, changing land use, and nitrogen limitation on vegetation growth. Using a global land model, we show how these factors interact locally to generate the global land carbon sink over the past 200 years. Nitrogen constraints were alleviated by N2 fixation in the tropics and by atmospheric nitrogen deposition in extratropical regions. Nonlinear interactions between land use change and land carbon and nitrogen cycling originated from three major mechanisms: (i) a sink foregone that would have occurred without land use conversion; (ii) an accelerated response of secondary vegetation to CO2 and nitrogen, and (iii) a compounded clearance loss from deforestation. Over time, these nonlinear effects have become increasingly important and reduce the present-day net carbon sink by ~40% or 0.4 PgC yr−1.

Warm Conveyor Belts in the ERA-Interim Dataset (1979–2010): Part I: Climatology and Potential Vorticity Evolution

A global climatology of warm conveyor belts (WCBs) is presented for the years 1979–2010, based on trajectories calculated with Interim ECMWF Re-Analysis (ERA-Interim) data. WCB trajectories are identified as strongly ascending air parcels (600 hPa in 2 days) near extratropical cyclones. Corroborating earlier studies, WCBs are more frequent during winter than summer and they ascend preferentially in the western ocean basins between 25° and 50° latitude. Before ascending, WCB trajectories typically approach from the subtropics in summer and from more midlatitude regions in winter. Considering humidity, cloud water, and potential temperature along WCBs confirms that they experience strong condensation and integrated latent heating during the ascent (typically >20 K). Liquid and ice water contents along WCBs peak at about 700 and 550 hPa, respectively. The mean potential vorticity (PV) evolution shows typical tropospheric values near 900 hPa, followed by an increase to almost 1 potential vorticity unit (PVU) at 700 hPa, and a decrease to less than 0.5 PVU at 300 hPa. These low PV values in the upper troposphere constitute significant negative anomalies with amplitudes of 1–3 PVU, which can strongly influence the downstream flow. Considering the low-level diabatic PV production, (i) WCBs starting at low latitudes (<40°) are unlikely to attain high PV (due to weak planetary vorticity) although they exhibit the strongest latent heating, and (ii) for those ascending at higher latitudes, a strong vertical heating gradient and high absolute vorticity are both important. This study therefore provides climatological insight into the cloud diabatic formation of significant positive and negative PV anomalies in the extratropical lower and upper troposphere, respectively.

A Lagrangian Analysis of the Northern Hemisphere Subtropical Jet

This study presents a 5-yr climatology of 7-day back trajectories started from the Northern Hemisphere subtropical jet. These trajectories provide insight into the seasonally and regionally varying angular momentum and potential vorticity characteristics of the air parcels that end up in the subtropical jet. The trajectories reveal preferred pathways of the air parcels that reach the subtropical jet from the tropics and the extratropics and allow estimation of the tropical and extratropical forcing of the subtropical jet. The back trajectories were calculated 7 days back in time and started every 6 h from December 2005 to November 2010 using the Interim European Centre for Medium-Range Weather Forecasts (ECMWF) Re-Analysis (ERA-Interim) dataset as a basis. The trajectories were started from the 345-K isentrope in areas where the wind speed exceeded a seasonally varying threshold and where the wind shear was confined to upper levels. During winter, the South American continent, the Indian Ocean, and the Maritime Continent are preferred areas of ascent into the upper troposphere. From these areas, air parcels follow an anticyclonic pathway into the subtropical jet. During summer, the majority of air parcels ascend over the Himalayas and Southeast Asia. Angular momentum is overall well conserved for trajectories that reach the subtropical jet from the deep tropics. In winter and spring, the hemispheric-mean angular momentum loss amounts to approximately 6%; in summer, it amounts to approximately 18%; and in fall, it amounts to approximately 13%. This seasonal variability is confirmed using an independent potential vorticity–based method to estimate tropical and extratropical forcing of the subtropical jet.

On the Co-Occurrence of Warm Conveyor Belt Outflows and PV Streamers*

The co-occurrence of warm conveyor belts (WCBs), strongly ascending moist airstreams in extratropical cyclones, and stratospheric potential vorticity (PV) streamers, indicators for breaking Rossby waves on the tropopause, is investigated for a 21-yr period in the Northern Hemisphere using Interim European Centre for Medium-Range Weather Forecasts (ECMWF) Re-Analysis (ERA-Interim) data. WCB outflows and PV streamers are respectively identified as two- and three-dimensional objects and tracked during their life cycle. PV streamers are more frequent than WCB outflows and nearly 15% of all PV streamers co-occur with WCBs during their life cycle, whereas about 60% of all WCB outflows co-occur with PV streamers. Co-occurrences are most frequent over the North Atlantic and North Pacific in spring and winter. WCB outflows are often located upstream of the PV streamers and form earlier, indicating the importance of diabatic processes for downstream Rossby wave breaking. Less frequently, PV streamers occur first, leading to the formation of new WCBs.

Tracking eolian dust with helium and thorium: Impacts of grain size and provenance

Reconstructions of the deposition rate of windblown mineral dust in ocean sediments offer an important means of tracking past climate changes and of assessing the radiative and biogeochemical impacts of dust in past climates. Dust flux estimates in ocean sediments have commonly been based on the operationally defined lithogenic fraction of sediment samples. More recently, dust fluxes have been estimated from measurements of helium and thorium, as rare isotopes of these elements (He-3 and Th-230) allow estimates of sediment flux, and the dominant isotopes (He-4 and Th-232) are uniquely associated with the lithogenic fraction of marine sediments. In order to improve the fidelity of dust flux reconstructions based on He and Th, we present a survey of He and Th concentrations in sediments from dust source areas in East Asia, Australia and South America. Our data show systematic relationships between He and Th concentrations and grain size, with He concentrations decreasing and Th concentrations increasing with decreasing grain size. We find consistent He and Th concentrations in the fine fraction (<5 μm) of samples from East Asia, Australia and Central South America (Puna-Central West Argentina), with Th concentrations averaging 14 μg/g and He concentrations averaging 2 μcc STP/g. We recommend use of these values for estimating dust fluxes in sediments where dust is dominantly fine-grained, and suggest that previous studies may have systematically overestimated Th-based dust fluxes by 30%. Source areas in Patagonia appear to have lower He and Th contents than other regions, as fine fraction concentrations average 0.8 μcc STP/g and 9 μg/g for 4He and 232Th, respectively. The impact of grain size on lithogenic He and Th concentrations should be taken into account in sediments proximal to dust sources where dust grain size may vary considerably. Our data also have important implications for the hosts of He in long-traveled dust and for the 3He/4He ratio used for terrigenous He in studies of extraterrestrial He in sediments and ice. We also investigate the use of He/Th ratios as a provenance tracer. Our results suggest differences in fine fraction He/Th ratios between East Asia, Australia, central South America and Patagonia, with ratios showing a positive relationship with the geological age of source rocks. He/Th ratios may thus provide useful provenance information, for example allowing separation of Patagonian sources from Puna-Central West Argentina or Australian dust sources. He/Th ratios in open-ocean marine sediments are similar to ratios in the fine fraction of upwind dust source areas. He/Th ratios in mid-latitude South Atlantic sediments suggest that dust in this region primarily derives from the Puna-Central West Argentina region (23–32°S) rather than Patagonia (>38°S). In the equatorial Pacific, He/Th ratios are much lower than in extratropical Pacific sediments or potential source areas measured as a part of this study (East Asia, South America, Australia) for reasons that are at present unclear, complicating their use as provenance tracers in this region.

Characterising the three-dimensional ozone distribution of a tidally locked Earth-like planet

We simulate the 3D ozone distribution of a tidally locked Earth-like exoplanet using the high-resolution, 3D chemistry climate model CESM1(WACCM) and study how the ozone layer of a tidally locked Earth (TLE) (ΩTLE = 1/365 days) differs from that of our present-day Earth (PDE) (ΩPDE = 1/1 day). The middle atmosphere reaches a steady state a symptotically within the first 80 days of the simulation. An upwelling, centred on the subsolar point, is present on the day side while a downwelling, centred on the antisolar point, is present on the night side. In the mesosphere, we find similar global ozone distributions for the TLE and the PDE, with decreased ozone on the day side and enhanced ozone on the night side. In the lower mesosphere, a jet stream transitions into a large-scale vortex around a low-pressure system, located at low latitudes of the TLE night side. In the middle stratosphere, the concentration of odd oxygen is approximately equal to that of the ozone [(Ox) ≈ (O3)]. At these altitudes, the lifetime of odd oxygen is ~16 h and the transport processes significantly contribute to the global distribution of stratospheric ozone. Compared to the PDE, where the strong Coriolis force acts as a mixing barrier between low and high latitudes, the transport processes of the TLE are governed by jet streams variable in the zonal and meridional directions. In the middle stratosphere of the TLE, we find high ozone values on the day side, due to the increased production of atomic oxygen on the day side, where it immediately recombines with molecular oxygen to form ozone. In contrast, the ozone is depleted on the night side, due to changes in the solar radiation distribution and the presence of a downwelling. As a result of the reduced Coriolis force, the tropical and extratropical air masses are well mixed and the global temperature distribution of the TLE stratosphere has smaller horizontal gradients than the PDE. Compared to the PDE, the total ozone column global mean is reduced by ~19.3 %. The day side and the night side total ozone column means are reduced by 23.21 and 15.52 %, respectively. Finally, we present the total ozone column (TOC) maps as viewed by a remote observer for four phases of the TLE during its revolution around the star. The mean TOC values of the four phases of the TLE vary by up to 23 %.