From mid-May to late June, the eastward propagation of some mid-tropospheric cyclones (MTCs) originating in the Bay of Bengal (BOB)-Burma (BUM) region across northern Indochina and southwest China to East Asia are observed. From 1979 to 2016, 137 MTCs were identified for the BOB-BUM region. Eighty-eight MTCs were trapped over northern Indochina, but 49 were propagated by the southwesterly flow of the BOB mid-tropospheric trough to reach East Asia. Thirty-one moved across Taiwan, but 18 did not. Eventually, all these MTCs moved to Okinawa. The former group of rainstorms produced rainfall >= 100 mm.day(-1) over Taiwan, while 44 rainstorms generated rainfall >= 50 mm.day(- 1) over Okinawa. The 49 eastward propagating through-Indochina MTCs exhibited the following salient features: (a) initial formation mostly occurred in the mid-troposphere ahead/within the BOB trough in the morning aver the BOB, but in the afternoon-evening over BUM. (b) The eastward propagating BOB-BUM MTCs underwent two descending processes: descending one mandatory level 6 hr after formation, and descending again from 700 to 925 hPa over southeastern China to interact with a low-level cyclonic shear flow. (c) Intensified (weakened) by dipole anomalous circulation cells, the confluent monsoon southwesterly flow west (east) of Taiwan facilitated the merger of the MTC with a local front. (d) Rainfall produced by rainstorms over Taiwan/ Okinawa was maintained by the intensified convergence of water vapour flux toward Taiwan/Okinawa. This hydrological process was stronger toward Taiwan than Okinawa and generated more rainfall in the former region. The impact of the BOB-BUM MTC on East Asian weather is confirmed by this study.

The ability to find and recognize patterns in high-dimensional geophysical data is fundamental to climate science and critical for meaningful interpretation of weather and climate processes. Archetypal analysis (AA) is one technique that has recently gained traction in the geophysical science community for its ability to find patterns based on extreme conditions. While traditional empirical orthogonal function (EOF) analysis can reveal patterns based on data covariance, AA seeks patterns from the points located at the edges of the data distribution. The utility of any objective pattern method depends on the properties of the data to which it is applied and the choices made in implementing the method. Given the relative novelty of the application of AA in geophysics it is important to develop experience in applying the method. We provide an assessment of the method, implementation, sensitivity, and interpretation of AA with respect to geophysical data. As an example for demonstration, we apply AA to a 39-yr sea surface temperature (SST) reanalysis dataset. We show that the decisions made to implement AA can significantly affect the interpretation of results, but also, in the case of SST, that the analysis is exceptionally robust under both spatial and temporal coarse graining. Significance Statement Archetypal analysis (AA), when applied to geophysical fields, is a technique designed to find typical configurations or modes in underlying data. This technique is relatively new to the geophysical science community and has been shown to be beneficial to the interpretation of extreme modes of the climate system. The identification of extreme modes of variability and their expression in day-to-day weather or state of the climate at longer time scales may help in elucidating the interplay between major teleconnection drivers and their evolution in a changing climate. The purpose of this work is to bring together a comprehensive report of the AA methodology using an SST reanalysis for demonstration. It is shown that the AA results are significantly affected by each implementation decision, but also can be resilient to spatiotemporal averaging. Any application of AA should provide a clear documentation of the choices made in applying the method.

Abstract Wildfire can cause significant adverse impacts to society and the environment. Weather and climate play an important role in modulating wildfire activity. We explore the joint occurrence of global fire weather and meteorological drought using a compound events framework. We show that, for much of the globe, burned area increases when periods of heightened fire weather compound with dry antecedent conditions. Regions associated with wildfire disasters, such as southern Australia and the western USA, are prone to experiencing years of compound drought and fire weather. Such compound events have increased in frequency for much of the globe, driven primarily by increases in fire weather rather than changes in precipitation. El Ni $$\tilde{{{{\rm{n}}}}}$$ n ̃ o Southern Oscillation is associated with widespread, spatially compounding drought and fire weather. In the Northern Hemisphere, a La Ni $$\tilde{{{{\rm{n}}}}}$$ n ̃ a signature is evident, whereas El Ni $$\tilde{{{{\rm{n}}}}}$$ n ̃ o is associated with such events in the tropics and, to a lesser degree, the Southern Hemisphere. Other climate modes and regional patterns of atmospheric circulation are also important, depending on the region. We show that the lengths of the fire weather seasons in eastern Australia and western North America have increased substantially since 2000, raising the likelihood of overlapping fire weather events in these regions. These cross-hemispheric events may be linked to the occurrence of El Ni $$\tilde{{{{\rm{n}}}}}$$ n ̃ o, although the sea-surface temperature magnitudes are small. Instead, it is likely that anthropogenic climate change is the primary driver of these changes.

With increased interest in climate forecasts and projections, it is important to understand more about their sources and levels of skill. A starting point here is to describe the nature of the skill associated with forecasts and projections. Climate forecasts and projections typically both include time varying forcing of the climate, but only forecasts have initial conditions set close to the observed climate state. Climate forecasts therefore derive skill from both initial conditions and from forcing. The character of the initial condition skill and forcing skill is different. Skill from initial conditions results in a narrowing of expectations relative to a climatological distribution and points toward a more favoured part of the distribution. Forcing skill could result from a shift in the preferred parts of the climatological distribution in response to forcing, or it could result from a shift in the entire distribution, or both. Assessments of forcing skill require time averages of the target variable that are long enough so that the contributions from internal variations are small compared to the forced response. The assessment of skill of climate forecasts and projections is inherently partial because of the small number of repeated trials possible on typical climate time scales but is nonetheless the only direct measure of their performance.

Large-scale cloud features referred to as cloudbands are known to be related to widespread and heavy rain via the transport of tropical heat and moisture to higher latitudes. The Australian northwest cloudband is such a feature that has been identified in simple searches of satellite imagery but with limited investigation of its atmospheric dynamical support. An accurate, long-term climatology of northwest cloudbands is key to robustly assessing these events. A dynamically based search algorithm has been developed that is guided by the presence and orientation of the subtropical jet stream. This jet stream is the large-scale atmospheric feature that determines the development and alignment of a cloudband. Using a new 40-yr dataset of cloudband events compiled by this search algorithm, composite atmospheric and ocean surface conditions over the period 1979–2018 have been assessed. Composite cloudband upper-level flow revealed a tilted low pressure trough embedded in a Rossby wave train. Composites of vertically integrated water vapor transport centered around the jet maximum during northwest cloudband events reveal a distinct atmospheric river supplying tropical moisture for cloudband rainfall. Parcel backtracking indicated multiple regions of moisture support for cloudbands. A thermal wind anomaly orientated with respect to an enhanced sea surface temperature gradient over the Indian Ocean was also a key composite cloudband feature. A total of 300 years of a freely coupled control simulation of the ACCESS-D system was assessed for its ability to simulate northwest cloudbands. Composite analysis of model cloudbands compared reasonably well to reanalysis despite some differences in seasonality and frequency of occurrence.

The CSIRO Climate retrospective Analysis and Forecast Ensemble system: version 1 (CAFE60v1) provides a large (96 member) ensemble retrospective analysis of the global climate system from 1960 to present with sufficiently many realizations and at spatio-temporal resolutions suitable to enable probabilistic climate studies. Using a variant of the ensemble Kalman filter, 96 climate state estimates are generated over the most recent six decades. These state estimates are constrained by monthly mean ocean, atmosphere and sea ice observations such that their trajectories track the observed state while enabling estimation of the uncertainties in the approximations to the retrospective mean climate over recent decades. For the atmosphere, we evaluate CAFE60v1 in comparison to empirical indices of the major climate teleconnections and blocking with various reanalysis products. Estimates of the large scale ocean structure, transports and biogeochemistry are compared to those derived from gridded observational products and climate model projections (CMIP). Sea ice (extent, concentration and variability) and land surface (precipitation and surface air temperatures) are also compared to a variety of model and observational products. Our results show that CAFE60v1 is a useful, comprehensive and unique data resource for studying internal climate variability and predictability, including the recent climate response to anthropogenic forcing on multi-year to decadal time scales.

From time to time atmospheric flows become organized and form coherent long-lived structures. Such structures could be propagating, quasi-stationary, or recur in place. We investigate the ability of Principal Components Analysis (PCA) and Archetypal Analysis (AA) to identify long-lived events, excluding propagating forms. Our analysis is carried out on the Southern Hemisphere mid-tropospheric flow represented by geopotential height at 500hPa ( Z 500 ). The leading basis patterns of Z 500 for PCA and AA are similar and describe structures representing (or similar to) the Southern Annular Mode (SAM) and Pacific South American (PSA) pattern. Long-lived events are identified here from sequences of 8 days or longer where the same basis pattern dominates for PCA or AA. AA identifies more long-lived events than PCA using this approach. The most commonly occurring long-lived event for both AA and PCA is the annular SAM-like pattern. The second most commonly occurring event is the PSA-like Pacific wavetrain for both AA and PCA. For AA the flow at any given time is approximated as weighted contributions from each basis pattern, which lends itself to metrics for discriminating among basis patterns. These show that the longest long-lived events are in general better expressed than shorter events. Case studies of long-lived events featuring a blocking structure and an annular structure show that both PCA and AA can identify and discriminate the dominant basis pattern that most closely resembles the flow event.

Subseasonal forecast skill is not homogeneous in time, and prior assessment of the likely forecast skill would be valuable for end-users. We propose a method for identifying periods of high forecast confidence using atmospheric circulation patterns, with an application to southern Australia precipitation. In particular, we use archetypal analysis to derive six patterns, called archetypes, of daily 500-hPa geopotential height (Z(500)) fields over Australia. We assign Z(500) reanalysis fields to the closest-matching archetype and subsequently link the archetypes to precipitation for three key regions in the Australian agriculture and energy sectors: the Murray Basin, southwest Western Australia, and western Tasmania. Using a 20-yr hindcast dataset from the European Centre for Medium-Range Weather Forecasts subseasonal-to-seasonal prediction system, we identify periods of high confidence as when hindcast Z(500) fields closely match an archetype according to a distance criterion. We compare the precipitation hindcast accuracy during these confident periods compared to normal. Considering all archetypes, we show that there is greater skill during confident periods for lead times of less than 10 days in the Murray Basin and western Tasmania, and for greater than 6 days in southwest Western Australia, although these conclusions are subject to substantial uncertainty. By breaking down the skill results for each archetype individually, we highlight how skill tends to be greater than normal for those archetypes associated with drier-than-average conditions.

Between June 2019 and March 2020, thousands of wildfires spread devastation across Australia at the tragic cost of many lives, vast areas of burnt forest, and estimated economic losses upward of AU$100 billion. Exceptionally hot and dry weather conditions, and preceding years of severe drought across Australia, contributed to the severity of the wildfires. Here we present analysis of a very large ensemble of initialized climate simulations to assess the likelihood of the concurrent drought and fire-weather conditions experienced at that time. We focus on a large region in southeast Australia where these fires were most widespread and define two indices to quantify the susceptibility to fire from drought and fire weather. Both indices were unprecedented in the observed record in 2019. We find that the likelihood of experiencing such extreme susceptibility to fire in the current climate was 0.5%, equivalent to a 200 year return period. The conditional probability is many times higher than this when we account for the states of key climate modes that impact Australian weather and climate. Drought and fire-weather conditions more extreme than those experienced in 2019 are also possible in the current climate.

Assessments of climate forecast skill depend on choices made by the assessor. In this perspective, we use forecasts of the El Niño-Southern-Oscillation to outline the impact of bias-correction on skill. Many assessments of skill from hindcasts (past forecasts) are probably overestimates of attainable forecast skill because the hindcasts are informed by observations over the period assessed that would not be available to real forecasts. Differences between hindcast and forecast skill result from changes in model biases from the period used to form forecast anomalies to the period over which the forecast is made. The relative skill rankings of models can change between hindcast and forecast systems because different models have different changes in bias across periods.

There is evidence suggesting that fire activity in southeast Australia has increased in the past two decades. There is also anecdotal evidence that pyrocumulonimbus events have become more common, although the observed record is short. We explore the extent to which possible changes in fire and pyrocumulonimbus occurrence in southeast Australia can be explained by changes in the number of extreme (above the 95th percentile) days per year of two fire-weather indices, the Forest Fire Danger Index (FFDI) and the continuous Haines index (C-Haines). For the period 1958 through 2020, we show that there is a dependence between the number of extreme FFDI days per year and burned area for forested parts of southeast Australia. To a lesser extent, there is a relationship between the FFDI and pyrocumulonimbus occurrences. We find an increase in the occurrence of extreme FFDI days from the late 1990s, with up to 24 more extreme days per year. In the southwest of the state of Victoria, there are concomitant increases in temperature and wind speed, and decreases in relative humidity. We show that the structure of C-Haines has not changed, but extreme days more frequently coincide with extreme FFDI days. This suggests that possible observed increases in pyrocumulonimbus occurrences may partly be explained as a consequence of the increases in co-occurring extreme FFDI and C-Haines days, rather than due to changes in the atmospheric stability. We also find that changes in the atmospheric circulation are consistent with those seen in the fire weather indices. Increased surface and mid-tropospheric pressure since 1998 are consistent with a poleward shift of the storm track. In particular, an increase in anticyclonic anomalies in the Tasman Sea is indicative of enhanced blocking activity. Although blocking frequency in these longitudes has decreased overall, the proportion of extreme FFDI days that are considered blocked has increased.

Attribution of an extreme magnitude 1-day rainfall event in Hobart is inhibited by small sample size. For moderate magnitude Hobart daily rainfall extremes, models suggest that the associated extratropical lows will deliver more rainfall with weaker pressure anomalies in a warmer world.