An extreme case of plant-insect co-diversification: figs and fig-pollinating wasps

It is thought that speciation in phytophagous insects is often due to colonization of novel host plants, because radiations of plant and insect lineages are typically asynchronous. Recent phylogenetic comparisons have supported this model of diversification for both insect herbivores and specialized pollinators. An exceptional case where contemporaneous plant insect diversification might be expected is the obligate mutualism between fig trees (Ficus species, Moraceae) and their pollinating wasps (Agaonidae, Hymenoptera). The ubiquity and ecological significance of this mutualism in tropical and subtropical ecosystems has long intrigued biologists, but the systematic challenge posed by >750 interacting species pairs has hindered progress toward understanding its evolutionary history. In particular, taxon sampling and analytical tools have been insufficient for large-scale co-phylogenetic analyses. Here, we sampled nearly 200 interacting pairs of fig and wasp species from across the globe. Two supermatrices were assembled: on average, wasps had sequences from 77% of six genes (5.6kb), figs had sequences from 60% of five genes (5.5 kb), and overall 850 new DNA sequences were generated for this study. We also developed a new analytical tool, Jane 2, for event-based phylogenetic reconciliation analysis of very large data sets. Separate Bayesian phylogenetic analyses for figs and fig wasps under relaxed molecular clock assumptions indicate Cretaceous diversification of crown groups and contemporaneous divergence for nearly half of all fig and pollinator lineages. Event-based co-phylogenetic analyses further support the co-diversification hypothesis. Biogeographic analyses indicate that the presentday distribution of fig and pollinator lineages is consistent with an Eurasian origin and subsequent dispersal, rather than with Gondwanan vicariance. Overall, our findings indicate that the fig-pollinator mutualism represents an extreme case among plant-insect interactions of coordinated dispersal and long-term co-diversification.

Factors contributing to the development of extreme North Atlantic cyclones and their relationship with the NAO

The occurrence of extreme cyclones is analysed in terms of their relationship to the NAO phase and the dominating environmental variables controlling their intensification. These are latent energy (equivalent potential temperature 850 hPa is used as an indicator), upper-air baroclinicity, horizontal divergence and jet stream strength. Cyclones over the North Atlantic are identified and tracked using a numerical algorithm, permitting a detailed analysis of their life cycles. Extreme cyclones are selected as the 10% most severe in terms of intensity. Investigations focus on the main strengthening phase of each cyclone. The environmental factors are related to the NAO, which affects the location and orientation of the cyclone tracks, thus explaining why extreme cyclones occur more (less) frequently during strong positive (negative) NAO phases. The enhanced number of extreme cyclones in positive NAO phases can be explained by the larger area with suitable growth conditions, which is better aligned with the cyclone tracks and is associated with increased cyclone life time and intensity. Moreover, strong intensification of cyclones is frequently linked to the occurrence of extreme values of growth factors in the immediate vicinity of the cyclone centre. Similar results are found for ECHAM5/OM1 for present day conditions, demonstrating that relationships between the environment factors and cyclones are also valid in the GCM. For future climate conditions (following the SRES A1B scenario), the results are similar, but a small increase of the frequency of extreme values is detected near the cyclone cores. On the other hand, total cyclone numbers decrease by 10% over the North Atlantic. An exception is the region near the British Isles, which features increased track density and intensity of extreme cyclones irrespective of the NAO phase. These changes are associated with an intensified jet stream close to Europe. Moreover, an enhanced frequency of explosive developments over the British Isles is found, leading to more frequent windstorms affecting Europe.

Analysis of frequency and intensity of European winter storm events from a multi-model perspective, at synoptic and regional scales

This study focuses on the analysis of winter (October-November-December-January-February-March; ONDJFM) storm events and their changes due to increased anthropogenic greenhouse gas concentrations over Europe. In order to assess uncertainties that are due to model formulation, 4 regional climate models (RCMs) with 5 high resolution experiments, and 4 global general circulation models (GCMs) are considered. Firstly, cyclone systems as synoptic scale processes in winter are investigated, as they are a principal cause of the occurrence of extreme, damage-causing wind speeds. This is achieved by use of an objective cyclone identification and tracking algorithm applied to GCMs. Secondly, changes in extreme near-surface wind speeds are analysed. Based on percentile thresholds, the studied extreme wind speed indices allow a consistent analysis over Europe that takes systematic deviations of the models into account. Relative changes in both intensity and frequency of extreme winds and their related uncertainties are assessed and related to changing patterns of extreme cyclones. A common feature of all investigated GCMs is a reduced track density over central Europe under climate change conditions, if all systems are considered. If only extreme (i.e. the strongest 5%) cyclones are taken into account, an increasing cyclone activity for western parts of central Europe is apparent; however, the climate change signal reveals a reduced spatial coherency when compared to all systems, which exposes partially contrary results. With respect to extreme wind speeds, significant positive changes in intensity and frequency are obtained over at least 3 and 20% of the European domain under study (35–72°N and 15°W–43°E), respectively. Location and extension of the affected areas (up to 60 and 50% of the domain for intensity and frequency, respectively), as well as levels of changes (up to +15 and +200% for intensity and frequency, respectively) are shown to be highly dependent on the driving GCM, whereas differences between RCMs when driven by the same GCM are relatively small.

A genome-wide association meta-analysis identifies new childhood obesity loci

Multiple genetic variants have been associated with adult obesity and a few with severe obesity in childhood; however, less progress has been made in establishing genetic influences on common early-onset obesity. We performed a North American, Australian and European collaborative meta-analysis of 14 studies consisting of 5,530 cases (≥95th percentile of body mass index (BMI)) and 8,318 controls (<50th percentile of BMI) of European ancestry. Taking forward the eight newly discovered signals yielding association with P < 5 × 10(-6) in nine independent data sets (2,818 cases and 4,083 controls), we observed two loci that yielded genome-wide significant combined P values near OLFM4 at 13q14 (rs9568856; P = 1.82 × 10(-9); odds ratio (OR) = 1.22) and within HOXB5 at 17q21 (rs9299; P = 3.54 × 10(-9); OR = 1.14). Both loci continued to show association when two extreme childhood obesity cohorts were included (2,214 cases and 2,674 controls). These two loci also yielded directionally consistent associations in a previous meta-analysis of adult BMI(1).

Meta-analysis of the INSIG2 association with obesity including 74,345 individuals: does heterogeneity of estimates relate to study design?

The INSIG2 rs7566605 polymorphism was identified for obesity (BMI> or =30 kg/m(2)) in one of the first genome-wide association studies, but replications were inconsistent. We collected statistics from 34 studies (n = 74,345), including general population (GP) studies, population-based studies with subjects selected for conditions related to a better health status ('healthy population', HP), and obesity studies (OB). We tested five hypotheses to explore potential sources of heterogeneity. The meta-analysis of 27 studies on Caucasian adults (n = 66,213) combining the different study designs did not support overall association of the CC-genotype with obesity, yielding an odds ratio (OR) of 1.05 (p-value = 0.27). The I(2) measure of 41% (p-value = 0.015) indicated between-study heterogeneity. Restricting to GP studies resulted in a declined I(2) measure of 11% (p-value = 0.33) and an OR of 1.10 (p-value = 0.015). Regarding the five hypotheses, our data showed (a) some difference between GP and HP studies (p-value = 0.012) and (b) an association in extreme comparisons (BMI> or =32.5, 35.0, 37.5, 40.0 kg/m(2) versus BMI<25 kg/m(2)) yielding ORs of 1.16, 1.18, 1.22, or 1.27 (p-values 0.001 to 0.003), which was also underscored by significantly increased CC-genotype frequencies across BMI categories (10.4% to 12.5%, p-value for trend = 0.0002). We did not find evidence for differential ORs (c) among studies with higher than average obesity prevalence compared to lower, (d) among studies with BMI assessment after the year 2000 compared to those before, or (e) among studies from older populations compared to younger. Analysis of non-Caucasian adults (n = 4889) or children (n = 3243) yielded ORs of 1.01 (p-value = 0.94) or 1.15 (p-value = 0.22), respectively. There was no evidence for overall association of the rs7566605 polymorphism with obesity. Our data suggested an association with extreme degrees of obesity, and consequently heterogeneous effects from different study designs may mask an underlying association when unaccounted for. The importance of study design might be under-recognized in gene discovery and association replication so far.

Controls on inorganic nitrogen leaching from Finnish catchments assessed using a sensitivity and uncertainty analysis of the INCA-N model

The semi-distributed, dynamic INCA-N model was used to simulate the behaviour of dissolved inorganic nitrogen (DIN) in two Finnish research catchments. Parameter sensitivity and model structural uncertainty were analysed using generalized sensitivity analysis. The Mustajoki catchment is a forested upstream catchment, while the Savijoki catchment represents intensively cultivated lowlands. In general, there were more influential parameters in Savijoki than Mustajoki. Model results were sensitive to N-transformation rates, vegetation dynamics, and soil and river hydrology. Values of the sensitive parameters were based on long-term measurements covering both warm and cold years. The highest measured DIN concentrations fell between minimum and maximum values estimated during the uncertainty analysis. The lowest measured concentrations fell outside these bounds, suggesting that some retention processes may be missing from the current model structure. The lowest concentrations occurred mainly during low flow periods; so effects on total loads were small.

Extreme rainfall variability in Australia: Patterns, drivers and predictability

Leading patterns of observed monthly extreme rainfall variability in Australia are examined using an Empirical Orthogonal Teleconnection (EOT) method. Extreme rainfall variability is more closely related to mean rainfall variability during austral summer than in winter. The leading EOT patterns of extreme rainfall explain less variance in Australia-wide extreme rainfall than is the case for mean rainfall EOTs. We illustrate that, as with mean rainfall, the El Niño-Southern Oscillation (ENSO) has the strongest association with warm-season extreme rainfall variability, while in the cool-season the primary drivers are atmospheric blocking and the subtropical ridge. The Indian Ocean Dipole and Southern Annular Mode also have significant relationships with patterns of variability during austral winter and spring. Leading patterns of summer extreme rainfall variability have predictability several months ahead from Pacific sea surface temperatures (SSTs) and as much as a year in advance from Indian Ocean SSTs. Predictability from the Pacific is greater for wetter than average summer months than for months that are drier than average, whereas for the Indian Ocean the relationship has greater linearity. Several cool-season EOTs are associated with mid-latitude synoptic-scale patterns along the south and east coasts. These patterns have common atmospheric signatures denoting moist onshore flow and strong cyclonic anomalies often to the north of a blocking anti-cyclone. Tropical cyclone activity is observed to have significant relationships with some warm season EOTs. This analysis shows that extreme rainfall variability in Australia can be related to remote drivers and local synoptic-scale patterns throughout the year.

Importance of resolution and model configuration when downscaling extreme precipitation

Dynamical downscaling is frequently used to investigate the dynamical variables of extra-tropical cyclones, for example, precipitation, using very high-resolution models nested within coarser resolution models to understand the processes that lead to intense precipitation. It is also used in climate change studies, using long timeseries to investigate trends in precipitation, or to look at the small-scale dynamical processes for specific case studies. This study investigates some of the problems associated with dynamical downscaling and looks at the optimum configuration to obtain the distribution and intensity of a precipitation field to match observations. This study uses the Met Office Unified Model run in limited area mode with grid spacings of 12, 4 and 1.5 km, driven by boundary conditions provided by the ECMWF Operational Analysis to produce high-resolution simulations for the Summer of 2007 UK flooding events. The numerical weather prediction model is initiated at varying times before the peak precipitation is observed to test the importance of the initialisation and boundary conditions, and how long the simulation can be run for. The results are compared to raingauge data as verification and show that the model intensities are most similar to observations when the model is initialised 12 hours before the peak precipitation is observed. It was also shown that using non-gridded datasets makes verification more difficult, with the density of observations also affecting the intensities observed. It is concluded that the simulations are able to produce realistic precipitation intensities when driven by the coarser resolution data.

Can reanalyses represent extreme precipitation over England and Wales?

A new generation of reanalysis products is currently being produced that provides global gridded atmospheric data spanning more than a century. Such data may be useful for characterising the observed long-term variability of extreme precipitation events, particularly in regions where spatial coverage of surface observations is limited, and in the pre-satellite era. An analysis of extreme precipitation events is performed over England and Wales, investigating the ability of Twentieth Century Reanalysis and ERA-Interim to represent extreme precipitation accumulations as recorded in the England and Wales Precipitation dataset on accumulation time-scales from 1 to 7 days. Significant correlations are found between daily precipitation accumulation observations and both reanalysis products. A hit-rate analysis indicates that the reanalyses have hit rates (as defined by an event above the 98th percentile) of approximately 40–65% for extreme events in both summer (JJA) and winter (DJF). This suggests that both ERA-Interim and Twentieth Century Reanalysis are difficult to use for representing individual extreme precipitation events over England and Wales.

Atmospheric rivers do not explain UK summer extreme rainfall

Extreme rainfall events continue to be one of the largest natural hazards in the UK. In winter, heavy precipitation and floods have been linked with intense moisture transport events associated with atmospheric rivers (ARs), yet no large-scale atmospheric precursors have been linked to summer flooding in the UK. This study investigates the link between ARs and extreme rainfall from two perspectives: 1) Given an extreme rainfall event, is there an associated AR? 2) Given an AR, is there an associated extreme rainfall event? We identify extreme rainfall events using the UK Met Office daily rain-gauge dataset and link these to ARs using two different horizontal resolution atmospheric datasets (ERA-Interim and 20th Century Re-analysis). The results show that less than 35% of winter ARs and less than 15% of summer ARs are associated with an extreme rainfall event. Consistent with previous studies, at least 50% of extreme winter rainfall events are associated with an AR. However, less than 20% of the identified summer extreme rainfall events are associated with an AR. The dependence of the water vapor transport intensity threshold used to define an AR on the years included in the study, and on the length of the season, is also examined. Including a longer period (1900-2012) compared to previous studies (1979-2005) reduces the water vapor transport intensity threshold used to define an AR.

Theoretical insight into diagnosing observation error correlations using observation-minus-background and observation-minus-analysis statistics

To improve the quantity and impact of observations used in data assimilation it is necessary to take into account the full, potentially correlated, observation error statistics. A number of methods for estimating correlated observation errors exist, but a popular method is a diagnostic that makes use of statistical averages of observation-minus-background and observation-minus-analysis residuals. The accuracy of the results it yields is unknown as the diagnostic is sensitive to the difference between the exact background and exact observation error covariances and those that are chosen for use within the assimilation. It has often been stated in the literature that the results using this diagnostic are only valid when the background and observation error correlation length scales are well separated. Here we develop new theory relating to the diagnostic. For observations on a 1D periodic domain we are able to the show the effect of changes in the assumed error statistics used in the assimilation on the estimated observation error covariance matrix. We also provide bounds for the estimated observation error variance and eigenvalues of the estimated observation error correlation matrix. We demonstrate that it is still possible to obtain useful results from the diagnostic when the background and observation error length scales are similar. In general, our results suggest that when correlated observation errors are treated as uncorrelated in the assimilation, the diagnostic will underestimate the correlation length scale. We support our theoretical results with simple illustrative examples. These results have potential use for interpreting the derived covariances estimated using an operational system.

A priori error analysis of space–time Trefftz discontinuous Galerkin methods for wave problems

We present and analyse a space–time discontinuous Galerkin method for wave propagation problems. The special feature of the scheme is that it is a Trefftz method, namely that trial and test functions are solution of the partial differential equation to be discretised in each element of the (space–time) mesh. The method considered is a modification of the discontinuous Galerkin schemes of Kretzschmar et al. (2014) and of Monk & Richter (2005). For Maxwell’s equations in one space dimension, we prove stability of the method, quasi-optimality, best approximation estimates for polynomial Trefftz spaces and (fully explicit) error bounds with high order in the meshwidth and in the polynomial degree. The analysis framework also applies to scalar wave problems and Maxwell’s equations in higher space dimensions. Some numerical experiments demonstrate the theoretical results proved and the faster convergence compared to the non-Trefftz version of the scheme.

Regional projections of extreme apparent temperature days in Africa and the related potential risk to human health

Regional climate modelling was used to produce high resolution climate projections for Africa, under a “business as usual scenario”, that were translated into potential health impacts utilizing a heat index that relates apparent temperature to health impacts. The continent is projected to see increases in the number of days when health may be adversely affected by increasing maximum apparent temperatures (AT) due to climate change. Additionally, climate projections indicate that the increases in AT results in a moving of days from the less severe to the more severe Symptom Bands. The analysis of the rate of increasing temperatures assisted in identifying areas, such as the East African highlands, where health may be at increasing risk due to both large increases in the absolute number of hot days, and due to the high rate of increase. The projections described here can be used by health stakeholders in Africa to assist in the development of appropriate public health interventions to mitigate the potential health impacts from climate change.

Frequency of extreme rainfall events in Southern Brazil modulated by interannual and interdecadal variability

The frequency of extreme rainfall events in Southern Brazil is impacted by Ell Nino - Southern Oscillation (ENSO) episodes, especially in austral spring. There are two areas in which this impact is more significant: one is on the coast, where extreme events are more frequent during El Nino (EN) and the other one extends inland, where extreme events increase during EN and decrease during La Nina (LN). Atmospheric circulation patterns associated with severe rainfall in those areas are similar (opposite) to anomalous patterns characteristic of EN (LN) episodes, indicating why increase (decrease) of extreme events in EN (LN) episodes is favoured. The most recurrent precipitation patterns during extreme rainfall events in each of these areas are disclosed by Principal Component Analysis (PCA) and evidence the separation between extreme events in these areas: a severe precipitation event generally does not occur simultaneously in the coast and inland, although they may Occur inland and in the coastal region in sequence. Although EN predominantly enhances extreme rainfall, there are EN years in which fewer severe events occur than the average of neutral years, and also the enhancement of extreme rainfall is not uniform for different EN episodes, because the interdecadal non-ENSO variability also modulates significantly the frequency of extreme events in Southern Brazil. The inland region, which is more affected, shows increase (decrease) of extreme rainfall in association with the negative (positive) phase of the Atlantic Multidecadal Variability, with the negative (positive) phase of the Pacific Multidecadal Variability and with the positive (negative) phase of the Pacific Interdecadal Variability. Copyright (C) 2008 Royal Meteorological Society