Analysis Scheme in the Ensemble Kalman Filter

This paper discusses an important issue related to the implementation and interpretation of the analysis scheme in the ensemble Kalman filter . I t i s shown that the obser vations must be treated as random variables at the analysis steps. That is, one should add random perturbations with the correct statistics to the obser vations and generate an ensemble of obser vations that then is used in updating the ensemble of model states. T raditionally , this has not been done in previous applications of the ensemble Kalman filter and, as will be shown, this has resulted in an updated ensemble with a variance that is too low . This simple modification of the analysis scheme results in a completely consistent approach if the covariance of the ensemble of model states is interpreted as the prediction error covariance, and there are no further requirements on the ensemble Kalman filter method, except for the use of an ensemble of sufficient size. Thus, there is a unique correspondence between the error statistics from the ensemble Kalman filter and the standard Kalman filter approach

Application of data assimilation to ocean and climate prediction

Ocean prediction systems are now able to analyse and predict temperature, salinity and velocity structures within the ocean by assimilating measurements of the ocean’s temperature and salinity into physically based ocean models. Data assimilation combines current estimates of state variables, such as temperature and salinity, from a computational model with measurements of the ocean and atmosphere in order to improve forecasts and reduce uncertainty in the forecast accuracy. Data assimilation generally works well with ocean models away from the equator but has been found to induce vigorous and unrealistic overturning circulations near the equator. A pressure correction method was developed at the University of Reading and the Met Office to control these circulations using ideas from control theory and an understanding of equatorial dynamics. The method has been used for the last 10 years in seasonal forecasting and ocean prediction systems at the Met Office and European Center for Medium-range Weather Forecasting (ECMWF). It has been an important element in recent re-analyses of the ocean heat uptake that mitigates climate change.

Lightning as a space-weather hazard: UK thunderstorm activity modulated by the passage of the heliospheric current sheet

Lightning flash rates, RL, are modulated by corotating interaction regions (CIRs) and the polarity of the heliospheric magnetic field (HMF) in near-Earth space. As the HMF polarity reverses at the heliospheric current sheet (HCS), typically within a CIR, these phenomena are likely related. In this study, RL is found to be significantly enhanced at the HCS and at 27 days prior/after. The strength of the enhancement depends on the polarity of the HMF reversal at the HCS. Near-Earth solar and galactic energetic particle fluxes are also ordered by HMF polarity, though the variations qualitatively differ from RL, with the main increase occurring prior to the HCS crossing. Thus, the CIR effect on lightning is either the result of compression/amplification of the HMF (and its subsequent interaction with the terrestrial system) or that energetic particle preconditioning of the Earth system prior to the HMF polarity change is central to solar wind lightning coupling mechanism.

The impact of a thermodynamic sea-ice module in the COSMO numerical weather prediction model on simulations for the Laptev Sea, Siberian Arctic

Previous versions of the Consortium for Small-scale Modelling (COSMO) numerical weather prediction model have used a constant sea-ice surface temperature, but observations show a high degree of variability on sub-daily timescales. To account for this, we have implemented a thermodynamic sea-ice module in COSMO and performed simulations at a resolution of 15 km and 5 km for the Laptev Sea area in April 2008. Temporal and spatial variability of surface and 2-m air temperature are verified by four automatic weather stations deployed along the edge of the western New Siberian polynya during the Transdrift XIII-2 expedition and by surface temperature charts derived from Moderate Resolution Imaging Spectroradiometer (MODIS) satellite data. A remarkable agreement between the new model results and these observations demonstrates that the implemented sea-ice module can be applied for short-range simulations. Prescribing the polynya areas daily, our COSMO simulations provide a high-resolution and high-quality atmospheric data set for the Laptev Sea for the period 14-30 April 2008. Based on this data set, we derive a mean total sea-ice production rate of 0.53 km3/day for all Laptev Sea polynyas under the assumption that the polynyas are ice-free and a rate of 0.30 km3/day if a 10-cm-thin ice layer is assumed. Our results indicate that ice production in Laptev Sea polynyas has been overestimated in previous studies.

Coordinated weather balloon solar radiation measurements during a solar eclipse

Solar eclipses provide a rapidly changing solar radiation environment. These changes can be studied using simple photodiode sensors, if the radiation reaching the sensors is unaffected by cloud.Transporting the sensors aloft using standard meteorological instrument packages modified to carry extra sensors, provides one promising but hitherto unexploited possibility for making solar eclipse radiation measurements. For the 20th March 2015 solar eclipse, a coordinated campaign of balloon-carried solar radiation measurements was undertaken from Reading (51.44N, 0.94W), Lerwick (60.15N, 1.13W) and Reykjavik (64.13N, 21.90W), straddling the path of the eclipse.The balloons reached sufficient altitude at the eclipse time for eclipse-induced variations in solar radiation and solar limb darkening to be measured above cloud. Because the sensor platforms were free to swing, techniques have been evaluated to correct the measurements for their changing orientation. In the swing-averaged technique, the mean value across a set of swings was used to approximate the radiation falling on a horizontal surface; in the swing-maximum technique, the direct beam was estimated by assuming the sensing surface becomes normal to the solar beam direction at a maximum swing. Both approaches, essentially independent,give values that agree with theoretical expectations for the eclipse-induced radiation changes.

El Niño 2015/2016 impact analysis, monthly outlook November 2015

During the summer and autumn 2015, El Niño conditions in the east and central Pacific have strengthened, disrupting weather patterns throughout the tropics and into the mid-latitudes. For example, rainfall during this summer’s Indian monsoon was approximately 15% below normal. The continued strong El Niño conditions have the potential to trigger damaging impacts (e.g., droughts, famines, floods), particularly in less-developed tropical countries, which would require a swift and effective humanitarian response to mitigate damage to life and property (e.g., health, migration, infrastructure). This analysis uses key climatic variables (temperature, soil moisture and precipitation) as measures to monitor the ongoing risk of these potentially damaging impacts. The previous 2015-2016 El Niño Impact Analysis was based on observations over the past 35 years and produced Impact Tables showing the likelihood and severity of the impacts on temperature and rainfall by season. The current report is an extension of this work providing information from seasonal forecast models to give a more detailed monthly outlook of the potential near-term impacts of the current El Niño conditions by region. This information has been added to the Impact Tables in the form of a monthly outlook column. This monthly outlook is an indication of the average likely conditions for that month and region and is not a definite prediction of weather impacts.

El Niño 2015/2016 impact analysis, monthly outlook December 2015

During the summer and autumn 2015, El Niño conditions in the east and central Pacific have strengthened, disrupting weather patterns throughout the tropics and into the mid-latitudes. For example, rainfall during this summer’s Indian monsoon was approximately 15% below normal. The continued strong El Niño conditions have the potential to trigger damaging impacts (e.g. droughts, famines, floods), particularly in less-developed tropical countries, which would require a swift and effective humanitarian response to mitigate damage to life and property (e.g. health, migration, infrastructure). This analysis uses key climatic variables (temperature, soil moisture and precipitation) as measures to monitor the ongoing risk of these potentially damaging impacts. The previous 2015-2016 El Niño Impact Analysis was based on observations over the past 35 years and produced Impact Tables showing the likelihood and severity of the impacts on temperature and rainfall by season. The current report is an extension of this work providing information from observations and seasonal forecast models to give a more detailed outlook of the potential near-term impacts of the current El Niño conditions by region. This information has been added to the Impact Tables in the form of an ‘Observations and Outlook’ row. This consists of observational information for the past seasons of JJA 2015 and SON 2015, a detailed monthly outlook from 5 modeling centres for Dec 2015 and then longer-term seasonal forecast information from 2 modeling centres for the future seasons of JF 2016 and MAM 2016. The seasonal outlook information is an indication of the average likely conditions for that coming month (or season) and region and is not a definite prediction of weather impacts.

Exploiting existing ground-based remote sensing networks to improve high-resolution weather forecasts

A new generation of high-resolution (1 km) forecast models promises to revolutionize the prediction of hazardous weather such as windstorms, flash floods, and poor air quality. To realize this promise, a dense observing network, focusing on the lower few kilometers of the atmosphere, is required to verify these new forecast models with the ultimate goal of assimilating the data. At present there are insufficient systematic observations of the vertical profiles of water vapor, temperature, wind, and aerosols; a major constraint is the absence of funding to install new networks. A recent research program financed by the European Union, tasked with addressing this lack of observations, demonstrated that the assimilation of observations from an existing wind profiler network reduces forecast errors, provided that the individual instruments are strategically located and properly maintained. Additionally, it identified three further existing European networks of instruments that are currently underexploited, but with minimal expense they could deliver quality-controlled data to national weather services in near–real time, so the data could be assimilated into forecast models. Specifically, 1) several hundred automatic lidars and ceilometers can provide backscatter profiles associated with aerosol and cloud properties and structures with 30-m vertical resolution every minute; 2) more than 20 Doppler lidars, a fairly new technology, can measure vertical and horizontal winds in the lower atmosphere with a vertical resolution of 30 m every 5 min; and 3) about 30 microwave profilers can estimate profiles of temperature and humidity in the lower few kilometers every 10 min. Examples of potential benefits from these instruments are presented.

Can site and landscape-scale environmental attributes buffer bird populations against weather events?

Projected impacts of climate change on the populations and distributions of species pose a challenge for conservationists. In response, a number of adaptation strategies to enable species to persist in a changing climate have been proposed. Management to maximise the quality of habitat at existing sites may reduce the magnitude or frequency of climate-driven population declines. In addition large-scale management of landscapes could potentially improve the resilience of populations by facilitating inter-population movements. A reduction in the obstacles to species’ range expansion, may also allow species to track changing conditions better through shifts to new locations, either regionally or locally. However, despite a strong theoretical base, there is limited empirical evidence to support these management interventions. This makes it difficult for conservationists to decide on the most appropriate strategy for different circumstances. Here extensive data from long-term monitoring of woodland birds at individual sites are used to examine the two-way interactions between habitat and both weather and population count in the previous year. This tests the extent to which site-scale and landscape-scale habitat attributes may buffer populations against variation in winter weather (a key driver of woodland bird population size) and facilitate subsequent population growth. Our results provide some support for the prediction that landscape-scale attributes (patch isolation and area of woodland habitat) may influence the ability of some woodland bird species to withstand weather-mediated population declines. These effects were most apparent among generalist woodland species. There was also evidence that several, primarily specialist, woodland species are more likely to increase following population decline where there is more woodland at both site and landscape scales. These results provide empirical support for the concept that landscape-scale conservation efforts may make the populations of some woodland bird species more resilient to climate change. However in isolation, management is unlikely to provide a universal benefit to all species.

Understanding climate change projections for precipitation over western Europe with a weather typing approach

Precipitation over western Europe (WE) is projected to increase (decrease) roughly northward (equatorward) of 50°N during the 21st century. These changes are generally attributed to alterations in the regional large-scale circulation, e.g., jet stream, cyclone activity, and blocking frequencies. A novel weather typing within the sector (30°W–10°E, 25–70°N) is used for a more comprehensive dynamical interpretation of precipitation changes. A k-means clustering on daily mean sea level pressure was undertaken for ERA-Interim reanalysis (1979–2014). Eight weather types are identified: S1, S2, S3 (summertime types), W1, W2, W3 (wintertime types), B1, and B2 (blocking-like types). Their distinctive dynamical characteristics allow identifying the main large-scale precipitation-driving mechanisms. Simulations with 22 Coupled Model Intercomparison Project 5 models for recent climate conditions show biases in reproducing the observed seasonality of weather types. In particular, an overestimation of weather type frequencies associated with zonal airflow is identified. Considering projections following the (Representative Concentration Pathways) RCP8.5 scenario over 2071–2100, the frequencies of the three driest types (S1, B2, and W3) are projected to increase (mainly S1, +4%) in detriment of the rainiest types, particularly W1 (−3%). These changes explain most of the precipitation projections over WE. However, a weather type-independent background signal is identified (increase/decrease in precipitation over northern/southern WE), suggesting modifications in precipitation-generating processes and/or model inability to accurately simulate these processes. Despite these caveats in the precipitation scenarios for WE, which must be duly taken into account, our approach permits a better understanding of the projected trends for precipitation over WE.

El Niño 2015/2016 impact analysis, monthly outlook January 2016

During the summer and autumn of 2015, El Niño conditions in the east and central Pacific have strengthened, disrupting weather patterns throughout the tropics and into the mid-latitudes. For example, rainfall during this summer’s Indian monsoon was approximately 15% below normal. The continued strong El Niño conditions have the potential to trigger damaging impacts (e.g., droughts, famines, floods), particularly in less-developed tropical countries, which would require a swift and effective humanitarian response to mitigate damage to life and property (e.g., health, migration, infrastructure). This analysis uses key climatic variables (temperature, soil moisture and precipitation) as measures to monitor the ongoing risk of these potentially damaging impacts. The previous 2015-2016 El Niño Impact Analysis was based on observations over the past 35 years and produced Impact Tables showing the likelihood and severity of the impacts on temperature and rainfall by season. The current report is an extension of this work providing information from observations and seasonal forecast models to give a more detailed outlook of the potential near-term impacts of the current El Niño conditions by region. This information has been added to the Impact Tables in the form of an ‘Observations and Outlook’ row. This consists of observational information for the past seasons of JJA 2015, SON 2015 and Dec 2015, a detailed monthly outlook from 4 modeling centres for Jan 2016 and then longer-term seasonal forecast information from 2 modeling centres for the future seasons of Feb 2016, MAM 2016 and Jun 2016. The seasonal outlook information is an indication of the average likely conditions for that coming month (or season) and region and is not a definite prediction of weather impacts.

El Niño 2015/2016 impact analysis, monthly outlook February 2016

During the summer and autumn of 2015, El Niño conditions in the east and central Pacific strengthened, disrupting weather patterns throughout the tropics and into the mid-latitudes. For example, rainfall during the summer’s Indian monsoon was approximately 15% below normal. The continued strong El Niño conditions have the potential to trigger damaging impacts (e.g., droughts, famines, floods), particularly in less-developed tropical countries, which would require a swift and effective humanitarian response to mitigate damage to life and property (e.g., health, migration, infrastructure). This analysis uses key climatic variables (temperature, soil moisture and precipitation) as measures to monitor the ongoing risk of these potentially damaging impacts. The previous 2015-2016 El Niño Impact Analysis was based on observations over the past 35 years and produced Impact Tables showing the likelihood and severity of the impacts on temperature and rainfall by season. The current report is an extension of this work, providing information from observations and seasonal forecast models to give a more detailed outlook of the potential near-term impacts of the current El Niño conditions by region. This information has been added to the Impact Tables in the form of an ‘Observations and Outlook’ row. This consists of observational information for the past seasons of JJA 2015, SON 2015 and DJ 2015/2016, a detailed monthly outlook from 5 modeling centres for Feb 2016 and then longer-term seasonal forecast information from 2 modeling centres for the future seasons of MAM 2016 and JJ 2016. The seasonal outlook information is an indication of the average likely conditions for that coming month (or season) and region and is not a definite prediction of weather impacts. This report has been produced by University of Reading for Evidence on Demand with the assistance of the UK Department for International Development (DFID) contracted through the Climate, Environment, Infrastructure and Livelihoods Professional Evidence and Applied Knowledge Services (CEIL PEAKS) programme, jointly managed by DAI (which incorporates HTSPE Limited) and IMC Worldwide Limited.

Understanding advances in the simulation of the Madden-Julian Oscillation in a numerical weather prediction model.

The Madden-Julian Oscillation (MJO) is the dominant mode of intraseasonal variability in the Trop- ics. It can be characterised as a planetary-scale coupling between the atmospheric circulation and organised deep convection that propagates east through the equatorial Indo-Pacific region. The MJO interacts with weather and climate systems on a near-global scale and is a crucial source of predictability for weather forecasts on medium to seasonal timescales. Despite its global signifi- cance, accurately representing the MJO in numerical weather prediction (NWP) and climate models remains a challenge. This thesis focuses on the representation of the MJO in the Integrated Forecasting System (IFS) at the European Centre for Medium-Range Weather Forecasting (ECMWF), a state-of-the-art NWP model. Recent modifications to the model physics in Cycle 32r3 (Cy32r3) of the IFS led to ad- vances in the simulation of the MJO; for the first time the observed amplitude of the MJO was maintained throughout the integration period. A set of hindcast experiments, which differ only in their formulation of convection, have been performed between May 2008 and April 2009 to asses the sensitivity of MJO simulation in the IFS to the Cy32r3 convective parameterization. Unique to this thesis is the attribution of the advances in MJO simulation in Cy32r3 to the mod- ified convective parameterization, specifically, the relative-humidity-dependent formulation for or- ganised deep entrainment. Increasing the sensitivity of the deep convection scheme to environmen- tal moisture is shown to modify the relationship between precipitation and moisture in the model. Through dry-air entrainment, convective plumes ascending in low-humidity environments terminate lower in the atmosphere. As a result, there is an increase in the occurrence of cumulus congestus, which acts to moisten the mid-troposphere. Due to the modified precipitation-moisture relationship more moisture is able to build up which effectively preconditions the tropical atmosphere for the transition to deep convection. Results from this thesis suggest that a tropospheric moisture control on convection is key to simulating the interaction between the physics and large-scale circulation associated with the MJO.

Exploiting satellite-based rainfall for weather index insurance: the challenges of spatial and temporal aggregation

Lack of access to insurance exacerbates the impact of climate variability on smallholder famers in Africa. Unlike traditional insurance, which compensates proven agricultural losses, weather index insurance (WII) pays out in the event that a weather index is breached. In principle, WII could be provided to farmers throughout Africa. There are two data-related hurdles to this. First, most farmers do not live close enough to a rain gauge with sufficiently long record of observations. Second, mismatches between weather indices and yield may expose farmers to uncompensated losses, and insurers to unfair payouts – a phenomenon known as basis risk. In essence, basis risk results from complexities in the progression from meteorological drought (rainfall deficit) to agricultural drought (low soil moisture). In this study, we use a land-surface model to describe the transition from meteorological to agricultural drought. We demonstrate that spatial and temporal aggregation of rainfall results in a clearer link with soil moisture, and hence a reduction in basis risk. We then use an advanced statistical method to show how optimal aggregation of satellite-based rainfall estimates can reduce basis risk, enabling remotely sensed data to be utilized robustly for WII.