Ensemble prediction for nowcasting with a convection-permitting model—I: description of the system and the impact of radar-derived surface precipitation rates

A key strategy to improve the skill of quantitative predictions of precipitation, as well as hazardous weather such as severe thunderstorms and flash floods is to exploit the use of observations of convective activity (e.g. from radar). In this paper, a convection-permitting ensemble prediction system (EPS) aimed at addressing the problems of forecasting localized weather events with relatively short predictability time scale and based on a 1.5 km grid-length version of the Met Office Unified Model is presented. Particular attention is given to the impact of using predicted observations of radar-derived precipitation intensity in the ensemble transform Kalman filter (ETKF) used within the EPS. Our initial results based on the use of a 24-member ensemble of forecasts for two summer case studies show that the convective-scale EPS produces fairly reliable forecasts of temperature, horizontal winds and relative humidity at 1 h lead time, as evident from the inspection of rank histograms. On the other hand, the rank histograms seem also to show that the EPS generates too much spread for forecasts of (i) surface pressure and (ii) surface precipitation intensity. These may indicate that for (i) the value of surface pressure observation error standard deviation used to generate surface pressure rank histograms is too large and for (ii) may be the result of non-Gaussian precipitation observation errors. However, further investigations are needed to better understand these findings. Finally, the inclusion of predicted observations of precipitation from radar in the 24-member EPS considered in this paper does not seem to improve the 1-h lead time forecast skill.

Quantitative estimates of glacial and Holocene temperature and precipitation change in lowland Amazonian Bolivia

Quantitative estimates of temperature and precipitation change during the late Pleistocene and Holocene have been difficult to obtain for much of the lowland Neotropics. Using two published lacustrine pollen records and a climate-vegetation model based on the modern abundance distributions of 154 Neotropical plant families, we demonstrate how family-level counts of fossil pollen can be used to quantitatively reconstruct tropical paleoclimate and provide needed information on historic patterns of climatic change. With this family-level analysis, we show that one area of the lowland tropics, northeastern Bolivia, experienced cooling (1–3 °C) and drying (400 mm/yr), relative to present, during the late Pleistocene (50,000–12,000 calendar years before present [cal. yr B.P.]). Immediately prior to the Last Glacial Maximum (LGM, ca. 21,000 cal. yr B.P.), we observe a distinct transition from cooler temperatures and variable precipitation to a period of warmer temperatures and relative dryness that extends to the middle Holocene (5000–3000 cal. yr B.P.). This prolonged reduction in precipitation occurs against the backdrop of increasing atmospheric CO2 concentrations, indicating that the presence of mixed savanna and dry-forest communities in northeastern Bolivia durng the LGM was not solely the result of low CO2 levels, as suggested previously, but also lower precipitation. The results of our analysis demonstrate the potential for using the distribution and abundance structure of modern Neotropical plant families to infer paleoclimate from the fossil pollen record.

The COnvective Precipitation Experiment (COPE): investigating the origins of heavy precipitation in the southwestern UK

A recent field campaign in southwest England used numerical modeling integrated with aircraft and radar observations to investigate the dynamic and microphysical interactions that can result in heavy convective precipitation. The COnvective Precipitation Experiment (COPE) was a joint UK-US field campaign held during the summer of 2013 in the southwest peninsula of England, designed to study convective clouds that produce heavy rain leading to flash floods. The clouds form along convergence lines that develop regularly due to the topography. Major flash floods have occurred in the past, most famously at Boscastle in 2004. It has been suggested that much of the rain was produced by warm rain processes, similar to some flash floods that have occurred in the US. The overarching goal of COPE is to improve quantitative convective precipitation forecasting by understanding the interactions of the cloud microphysics and dynamics and thereby to improve NWP model skill for forecasts of flash floods. Two research aircraft, the University of Wyoming King Air and the UK BAe 146, obtained detailed in situ and remote sensing measurements in, around, and below storms on several days. A new fast-scanning X-band dual-polarization Doppler radar made 360-deg volume scans over 10 elevation angles approximately every 5 minutes, and was augmented by two UK Met Office C-band radars and the Chilbolton S-band radar. Detailed aerosol measurements were made on the aircraft and on the ground. This paper: (i) provides an overview of the COPE field campaign and the resulting dataset; (ii) presents examples of heavy convective rainfall in clouds containing ice and also in relatively shallow clouds through the warm rain process alone; and (iii) explains how COPE data will be used to improve high-resolution NWP models for operational use.

A quantitative analysis of the phytoplankton of the Gulf of Panama III. General ecological conditions and the phytoplankton dynamics at 8°45'N, 79°23'W from November 1954 to May 1957

Bi-weekly phytoplankton samples were collected at 0, 10, and 20 m and enumerated by the Utermöhl sedimentation technique; 14C productivity measurements at 10 m, oblique zooplankton tows, and routine hydrographic observations were also made. Northerly winds induce upwelling during December-April, followed by a rainy season; a slight resurgence in upwelling may occur during July and/or August. Annual variations in upwelling intensity and rainfall occur. During upwelling, the upper 50 m, about 30 per cent of the total volume of the Gulf of Panama, is replaced with water 5 to 10 C colder than the more stratified, turbid and nutrient impoverished watermass present during the rainy season. The mean annual runoff accompanying an average annual precipitation of 2731 mm is estimated to equal a layer of fresh water 3.2 m thick. About 10 per cent of the phytoplankton phosphate and inorganic nitrogen requirements during the rainy season are accreted. (PDF contains 260 pages.)

Downscaling of surface moisture flux and precipitation in the Ebro Valley (Spain) using analogues and analogues followed by random forests and multiple linear regression

In this paper, reanalysis fields from the ECMWF have been statistically downscaled to predict from large-scale atmospheric fields, surface moisture flux and daily precipitation at two observatories (Zaragoza and Tortosa, Ebro Valley, Spain) during the 1961-2001 period. Three types of downscaling models have been built: (i) analogues, (ii) analogues followed by random forests and (iii) analogues followed by multiple linear regression. The inputs consist of data (predictor fields) taken from the ERA-40 reanalysis. The predicted fields are precipitation and surface moisture flux as measured at the two observatories. With the aim to reduce the dimensionality of the problem, the ERA-40 fields have been decomposed using empirical orthogonal functions. Available daily data has been divided into two parts: a training period used to find a group of about 300 analogues to build the downscaling model (1961-1996) and a test period (19972001), where models' performance has been assessed using independent data. In the case of surface moisture flux, the models based on analogues followed by random forests do not clearly outperform those built on analogues plus multiple linear regression, while simple averages calculated from the nearest analogues found in the training period, yielded only slightly worse results. In the case of precipitation, the three types of model performed equally. These results suggest that most of the models' downscaling capabilities can be attributed to the analogues-calculation stage.

Wave forecasts associated with hurricane landfalls in the vicinity of Cape Lookout, North Carolina

Hurricanes can cause extensive damage to the coastline and coastal communities due to wind-generated waves and storm surge. While extensive modeling efforts have been conducted regarding storm surge, there is far less information about the effects of waves on these communities and ecosystems as storms make landfall. This report describes a preliminary use of NCCOS’ WEMo (Wave Exposure Model; Fonseca and Malhotra 2010) to compute the wind wave exposure within an area of approximately 25 miles radius from Beaufort, North Carolina for estuarine waters encompassing Bogue Sound, Back Sound and Core Sound during three hurricane landfall scenarios. The wind wave heights and energy of a site was a computation based on wind speed, direction, fetch and local bathymetry. We used our local area (Beaufort, North Carolina) as a test bed for this product because it is frequently impacted by hurricanes and we had confidence in the bathymetry data. Our test bed conditions were based on two recent Hurricanes that strongly affected this area. First, we used hurricane Isabel which made landfall near Beaufort in September 2003. Two hurricane simulations were run first by passing hurricane Isabel along its actual path (east of Beaufort) and second by passing the same storm to the west of Beaufort to show the potential effect of the reversed wind field. We then simulated impacts by a hurricane (Ophelia) with a different landfall track, which occurred in September of 2005. The simulations produced a geographic description of wave heights revealing the changing wind and wave exposure of the region as a consequence of landfall location and storm intensity. This highly conservative simulation (water levels were that of low tide) revealed that many inhabited and developed shorelines would receive wind waves for prolonged periods of time at heights far above that found during even the top few percent of non-hurricane events. The simulations also provided a sense for how rapidly conditions could transition from moderate to highly threatening; wave heights were shown to far exceed normal conditions often long before the main body of the storm arrived and importantly, at many locations that could impede and endanger late-fleeing vessels seeking safe harbor. When joined with other factors, such as storm surge and event duration, we anticipate that the WEMo forecasting tool will have significant use by local emergency agencies and the public to anticipate the relative exposure of their property arising as a function of storm location and may also be used by resource managers to examine the effects of storms in a quantitative fashion on local living marine resources.

Quantitative Determination of Microcystins in Rat Plasma by LC-ESI Tandem MS

A rapid and sensitive method was developed and validated for the determination of MCYST (microcystin)-RR, -LR, and [Dha(7)] MCYST-LR in rat plasma by liquid chromatography-tandem mass spectrometry. The analytes were extracted from rat plasma by protein precipitation, followed by solid-phase extraction. Liquid chromatography with electrospray ionization mass spectrometry, operating in selected reaction monitoring (SRM) mode, was used to quantify MCYST-RR, -LR, and [Dha(7)] MCYST-LR in rat plasma. The recoveries for each analyte in rat plasma ranged from 70.8 to 88.7%. The calibration curve was linear within the range from 0.005 to 1.25 mu g mL(-1). The limit of detection were 1.4, 1.0, 0.6 ng mL(-1) for MCYST-RR, -LR, and [Dha(7)] MCYST-LR. The overall precision was determined on three different days. The values for within- and between-day precision in rat plasma were within 15%. This method was applied to the identification and quantification of microcystins in rat plasma with acute exposure of microcystins via intravenous injection.

Electro-assisted precipitation of electrolytes in poly(3,4-ethylenedioxythiophene) film

The electrolyte, NaBF4, can be enriched into the matrix of poly(3,4-ethylenedioxythiophene) (PEDOT) film during the p-doping potential cycling between 0.6 and -0.9 V. It has been demonstrated that this enrichment is originated from the mixed ion transfer between doping and dedoping, i.e. BF4- anion migrate into the PEDOT film during the oxidation process, the Ne cation insert into the film during the reduction process, and then, the electrolyte is accumulated into the film matrix after the multiple CV cycling. The quantitative analysis of energy-dispersive X-ray spectroscopy (EDX) confirmed the enrichment of NaBF4 in the PEDOT film.

Water vapor variability in the tropics and its links to dynamics and precipitation

The distribution and variability of water vapor and its links with radiative cooling and latent heating via precipitation are crucial to understanding feedbacks and processes operating within the climate system. Column-integrated water vapor (CWV) and additional variables from the European Centre for Medium-Range Weather Forecasts (ECMWF) 40-year reanalysis (ERA40) are utilized to quantify the spatial and temporal variability in tropical water vapor over the period 1979–2001. The moisture variability is partitioned between dynamical and thermodynamic influences and compared with variations in precipitation provided by the Climate Prediction Center Merged Analysis of Precipitation (CMAP) and the Global Precipitation Climatology Project (GPCP). The spatial distribution of CWV is strongly determined by thermodynamic constraints. Spatial variability in CWV is dominated by changes in the large-scale dynamics, in particular associated with the El Niño–Southern Oscillation (ENSO). Trends in CWV are also dominated by dynamics rather than thermodynamics over the period considered. However, increases in CWV associated with changes in temperature are significant over the equatorial east Pacific when analyzing interannual variability and over the north and northwest Pacific when analyzing trends. Significant positive trends in CWV tend to predominate over the oceans while negative trends in CWV are found over equatorial Africa and Brazil. Links between changes in CWV and vertical motion fields are identified over these regions and also the equatorial Atlantic. However, trends in precipitation are generally incoherent and show little association with the CWV trends. This may in part reflect the inadequacies of the precipitation data sets and reanalysis products when analyzing decadal variability. Though the dynamic component of CWV is a major factor in determining precipitation variability in the tropics, in some regions/seasons the thermodynamic component cancels its effect on precipitation variability.

Global oil production: forecasts and methodologies

A range of forecasts of global oil production made between 1956 and the present day are listed. For the majority of these the methodology used to generate the forecast is described. The paper distinguishes between three types of forecast: group 1-quantitative analyses which predict that global oil production will reach a resource-limited peak in the near term, and certainly before the year 2020; group 2-forecasts that use quantitative methods, but which see no production peak within the forecast's time horizon (typically 2020 or 2030); group 3-nonquantitative analyses that rule out a resource-limited oil peak within the foreseeable future. The paper analyses these forecast types and suggests that group 1 forecasts are the most realistic.

A new method for evaluating regional air quality forecasts

A Kriging interpolation method is combined with an object-based evaluation measure to assess the ability of the UK Met Office's dispersion and weather prediction models to predict the evolution of a plume of tracer as it was transported across Europe. The object-based evaluation method, SAL, considers aspects of the Structure, Amplitude and Location of the pollutant field. The SAL method is able to quantify errors in the predicted size and shape of the pollutant plume, through the structure component, the over- or under-prediction of the pollutant concentrations, through the amplitude component, and the position of the pollutant plume, through the location component. The quantitative results of the SAL evaluation are similar for both models and close to a subjective visual inspection of the predictions. A negative structure component for both models, throughout the entire 60 hour plume dispersion simulation, indicates that the modelled plumes are too small and/or too peaked compared to the observed plume at all times. The amplitude component for both models is strongly positive at the start of the simulation, indicating that surface concentrations are over-predicted by both models for the first 24 hours, but modelled concentrations are within a factor of 2 of the observations at later times. Finally, for both models, the location component is small for the first 48 hours after the start of the tracer release, indicating that the modelled plumes are situated close to the observed plume early on in the simulation, but this plume location error grows at later times. The SAL methodology has also been used to identify differences in the transport of pollution in the dispersion and weather prediction models. The convection scheme in the weather prediction model is found to transport more pollution vertically out of the boundary layer into the free troposphere than the dispersion model convection scheme resulting in lower pollutant concentrations near the surface and hence a better forecast for this case study.

3D-Var assimilation of insect-derived Doppler radar radial winds in convective cases using a high resolution model

The assimilation of Doppler radar radial winds for high resolution NWP may improve short term forecasts of convective weather. Using insects as the radar target, it is possible to provide wind observations during convective development. This study aims to explore the potential of these new observations, with three case studies. Radial winds from insects detected by 4 operational weather radars were assimilated using 3D-Var into a 1.5 km resolution version of the Met Office Unified Model, using a southern UK domain and no convective parameterization. The effect on the analysis wind was small, with changes in direction and speed up to 45° and 2 m s−1 respectively. The forecast precipitation was perturbed in space and time but not substantially modified. Radial wind observations from insects show the potential to provide small corrections to the location and timing of showers but not to completely relocate convergence lines. Overall, quantitative analysis indicated the observation impact in the three case studies was small and neutral. However, the small sample size and possible ground clutter contamination issues preclude unequivocal impact estimation. The study shows the potential positive impact of insect winds; future operational systems using dual polarization radars which are better able to discriminate between insects and clutter returns should provided a much greater impact on forecasts.