Latitudinal gradients in butterfly population variability are influenced by landscape heterogeneity

The variability of populations over time is positively associated with their risk of local extinction. Previous work has shown that populations at the high-latitude boundary of species’ ranges show higher inter-annual variability, consistent with increased sensitivity and exposure to adverse climatic conditions. However, patterns of population variability at both high- and low-latitude species range boundaries have not yet been concurrently examined. Here, we assess the inter-annual population variability of 28 butterfly species between 1994 and 2009 at 351 and 18 sites in the United Kingdom and Catalonia, Spain, respectively. Local population variability is examined with respect to the position of the species’ bioclimatic envelopes (i.e. whether the population falls within areas of the ‘core’ climatic suitability or is a climatically ‘marginal’ population), and in relation to local landscape heterogeneity, which may influence these range location – population dynamic relationships. We found that butterfly species consistently show latitudinal gradients in population variability, with increased variability in the more northerly UK. This pattern is even more marked for southerly distributed species with ‘marginal’ climatic suitability in the UK but ‘core’ climatic suitability in Catalonia. In addition, local landscape heterogeneity did influence these range location – population dynamic relationships. Habitat heterogeneity was associated with dampened population dynamics, especially for populations in the UK. Our results suggest that promoting habitat heterogeneity may promote the persistence of populations at high-latitude range boundaries, which may potentially aid northwards expansion under climate warming. We did not find evidence that population variability increases towards southern range boundaries. Sample sizes for this region were low, but there was tentative evidence, in line with previous ecological theory, that local landscape heterogeneity may promote persistence in these retracting low-latitude range boundary populations.

Spatial and temporal variability of the concentration field from localized releases in a regular building array

Spatial and temporal fluctuations in the concentration field from an ensemble of continuous point-source releases in a regular building array are analyzed from data generated by direct numerical simulations. The release is of a passive scalar under conditions of neutral stability. Results are related to the underlying flow structure by contrasting data for an imposed wind direction of 0 deg and 45 deg relative to the buildings. Furthermore, the effects of distance from the source and vicinity to the plume centreline on the spatial and temporal variability are documented. The general picture that emerges is that this particular geometry splits the flow domain into segments (e.g. “streets” and “intersections”) in each of which the air is, to a first approximation, well mixed. Notable exceptions to this general rule include regions close to the source, near the plume edge, and in unobstructed channels when the flow is aligned. In the oblique (45 deg) case the strongly three-dimensional nature of the flow enhances mixing of a scalar within the canopy leading to reduced temporal and spatial concentration fluctuations within the plume core. These fluctuations are in general larger for the parallel flow (0 deg) case, especially so in the long unobstructed channels. Due to the more complex flow structure in the canyon-type streets behind buildings, fluctuations are lower than in the open channels, though still substantially larger than for oblique flow. These results are relevant to the formulation of simple models for dispersion in urban areas and to the quantification of the uncertainties in their predictions.

Eddies and variability in the Mozambique Channel

Between 1995 and 2000, on average 4 eddies per year are observed from satellite altimetry to propagate southward through the Mozambique Channel, into the upstream Agulhas region. Further south, these eddies have been found to control the timing and frequenc yof Agulhas ring shedding. Within the Mozambique Channel, anomalous SSH amplitudes rise to 30 cm ; in agreement with in situ measured velocities. Comparison of an observed velocit ysection with GCM model results shows that the Mozambique Channel eddies in these models are too surface intensified. Also, the number of eddies formed in the models is in disagreement with our observational analysis. Moored current meter measurements observing the passage of three eddies in 2000 are extended to a 5-year time series b yreferencing the anomalous surface currents estimated from altimeter data to a s ynoptic LADCP velocit y measurement. The results show intermittent edd ypassage at the mooring location. A statistical analysis of SSH observations in different parts of the Mozambique Channel shows a southward decrease of the dominant frequency of the variability, going from 7 per year in the extension of the South Equatorial Current north of Madagascar to 4 per year south of Madagascar. The observations suggest that frequency reduction is related to the Rossb ywaves coming in from the east

On the ability of statistical wind-wave models to capture the variability and long-term trends of the North Atlantic winter wave climate

A dynamical wind-wave climate simulation covering the North Atlantic Ocean and spanning the whole 21st century under the A1B scenario has been compared with a set of statistical projections using atmospheric variables or large scale climate indices as predictors. As a first step, the performance of all statistical models has been evaluated for the present-day climate; namely they have been compared with a dynamical wind-wave hindcast in terms of winter Significant Wave Height (SWH) trends and variance as well as with altimetry data. For the projections, it has been found that statistical models that use wind speed as independent variable predictor are able to capture a larger fraction of the winter SWH inter-annual variability (68% on average) and of the long term changes projected by the dynamical simulation. Conversely, regression models using climate indices, sea level pressure and/or pressure gradient as predictors, account for a smaller SWH variance (from 2.8% to 33%) and do not reproduce the dynamically projected long term trends over the North Atlantic. Investigating the wind-sea and swell components separately, we have found that the combination of two regression models, one for wind-sea waves and another one for the swell component, can improve significantly the wave field projections obtained from single regression models over the North Atlantic.

Simulating decadal variability in the North Atlantic Ocean

Observations and climate models suggest significant decadal variability within the North Atlantic subpolar gyre (NA SPG), though observations are sparse and models disagree on the details of this variability. Therefore, it is important to understand 1) the mechanisms of simulated decadal variability, 2) which parts of simulated variability are more faithful representations of reality, and 3) the implications for climate predictions. Here, we investigate the decadal variability in the NA SPG in the state-of-the-art, high resolution (0.25◦ ocean resolution), climate model ‘HadGEM3’. We find a decadal mode with a period of 17 years that explains 30% of the annual variance in related indices. The mode arises due to the advection of heat content anomalies, and shows asymmetries in the timescale of phase reversal between positive and negative phases. A negative feedback from temperature-driven density anomalies in the Labrador Sea (LS) allows for the phase reversal. The North Atlantic Oscillation (NAO), which exhibits the same periodicity, amplifies the mode. The atmosphere-ocean coupling is stronger during positive rather than negative NAO states, explaining the asymmetry. Within the NA SPG, there is potential predictability arising partly from this mode for up to 5 years. There are important similarities between observed and simulated variability, such as the apparent role for the propagation of heat content anomalies. However, observations suggest interannual LS density anomalies are salinity-driven. Salinity control of density would change the temperature feedback to the south, possibly limiting real-world predictive skill in the southern NA SPG with this model. Finally, to understand the diversity of behaviours, we analyse 42 present-generation climate models. Temperature and salinity biases are found to systematically influence the driver of density variability in the LS. Resolution is a good predictor of the biases. The dependence of variability on the background state has important implications for decadal predictions.

Annular modes and apparent eddy feedbacks in the Southern Hemisphere

Lagged correlation analysis is often used to infer intraseasonal dynamical effects but is known to be affected by non-stationarity. We highlight a pronounced quasi-two-year peak in the anomalous zonal wind and eddy momentum flux convergence power spectra in the Southern Hemisphere, which is prima facie evidence for non-stationarity. We then investigate the consequences of this non-stationarity for the Southern Annular Mode and for eddy momentum flux convergence. We argue that positive lagged correlations previously attributed to the existence of an eddy feedback are more plausibly attributed to non-stationary interannual variability external to any potential feedback process in the mid-latitude troposphere. The findings have implications for the diagnosis of feedbacks in both models and re-analysis data as well as for understanding the mechanisms underlying variations in the zonal wind.

Climate variability and socio-environmental changes in the northern Aegean (NE Mediterranean) during the last 1500 years

We provide new evidence on sea surface temperature (SST) variations and paleoceanographic/paleoenvironmental changes over the past 1500 years for the north Aegean Sea (NE Mediterranean). The reconstructions are based on multiproxy analyses, obtained from the high resolution (decadal to multidecadal) marine record M2 retrieved from the Athos basin. Reconstructed SSTs show an increase from ca. 850 to 950 AD and from ca. 1100 to 1300 AD. A cooling phase of almost 1.5 �C is observed from ca. 1600 AD to 1700 AD. This seems to have been the starting point of a continuous SST warming trend until the end of the reconstructed period, interrupted by two prominent cooling events at 1832 ± 15 AD and 1995 ± 1 AD. Application of an adaptive Kernel smoothing suggests that the current warming in the reconstructed SSTs of the north Aegean might be unprecedented in the context of the past 1500 years. Internal variability in atmospheric/oceanic circulations systems as well as external forcing as solar radiation and volcanic activity could have affected temperature variations in the north Aegean Sea over the past 1500 years. The marked temperature drop of approximately ~2 �C at 1832 ± 15 yr AD could be related to the 1809 АD ‘unknown’ and the 1815 AD Tambora volcanic eruptions. Paleoenvironmental proxy-indices of the M2 record show enhanced riverine/continental inputs in the northern Aegean after ca. 1450 AD. The paleoclimatic evidence derived from the M2 record is combined with a socio-environmental study of the history of the north Aegean region. We show that the cultivation of temperature-sensitive crops, i.e. walnut, vine and olive, co-occurred with stable and warmer temperatures, while its end coincided with a significant episode of cooler temperatures. Periods of agricultural growth in Macedonia coincide with periods of warmer and more stable SSTs, but further exploration is required in order to identify the causal links behind the observed phenomena. The Black Death likely caused major changes in agricultural activity in the north Aegean region, as reflected in the pollen data from land sites of Macedonia and the M2 proxy-reconstructions. Finally, we conclude that the early modern peaks in mountain vegetation in the Rhodope and Macedonia highlands, visible also in the M2 record, were very likely climate-driven.

Heart rate variability is associated with amygdala functional connectivity with MPFC across younger and older adults

The ability to regulate emotion is crucial to promote well-being. Evidence suggests that the medial prefrontal cortex (mPFC) and adjacent anterior cingulate (ACC) modulate amygdala activity during emotion regulation. Yet less is known about whether the amygdala-mPFC circuit is linked with regulation of the autonomic nervous system and whether the relationship differs across the adult lifespan. The current study tested the hypothesis that heart rate variability (HRV) reflects the strength of mPFC-amygdala interaction across younger and older adults. We recorded participants’ heart rates at baseline and examined whether baseline HRV was associated with amygdala-mPFC functional connectivity during rest. We found that higher HRV was associated with stronger functional connectivity between the amygdala and the mPFC during rest across younger and older adults. In addition to this age-invariant pattern, there was an age-related change, such that greater HRV was linked with stronger functional connectivity between amygdala and ventrolateral PFC (vlPFC) in younger than in older adults. These results are in line with past evidence that vlPFC is involved in emotion regulation especially in younger adults. Taken together, our results support the neurovisceral integration model and suggest that higher heart rate variability is associated with neural mechanisms that support successful emotional regulation across the adult lifespan.

The variability of baroreflex sensitivity in juvenile, spontaneously hypertensive rats

In this study the baroreflex sensitivity of conscious, juvenile, spontaneously hypertensive rats (SHRs) was compared. The study population consisted of 19 eight-week-old male SHRs. The baroreflex sensitivity was quantified as the derivative of the variation in heart rate (HR) and the variation of mean arterial pressure (baroreflex sensitivity = Delta HR/Delta MAP). MAP was manipulated with sodium nitroprusside (SNP) and phenylephrine (PHE), administered via an inserted cannula in the right femoral vein. The SHRs were divided into four groups: (1) low bradycardic baroreflex (LB) where the baroreflex gain (BG) was between 0 and 1 bpm/mmHg with PHE; (2) high bradycardic baroreflex (HB), where the BG was < -1 bpm/mmHg with PHE; (3) low tachycardic baroreflex (LT) where the BC was between 0 and 3 bpm/mmHg with SNP; (4) high tachycardic baroreflex (HT) where the BG was > 3 bpm/mmHg with SNP. We noted that 36.8% of the rats presented with an increased bradycardic reflex, while 27.8% demonstrated an attenuated tachycardic reflex. No significant alterations were noted regarding the basal MAP and HR. There were significant differences in the baroreflex sensitivity between SHRs in the same laboratory. One should be careful when interpreting studies employing the SHR as a research model.

Climatic simulations of the eastern Andes low-level jet and its dependency on convective parameterizations

The South American low level jet (SALLJ) of the Eastern Andes is investigated with Regional Climate Model version 3 (RegCM3) simulations during the 2002-2003 austral summer using two convective parameterizations (Grell and Emanuel). The simulated SALLJ is compared with the special observations of SALLJEX (SALLJ Experiment). Both the Grell and Emanuel schemes adequately simulate the low level flow over South America. However, there are some intensity differences. Due to the larger (smaller) convective activity, the Emanuel (Grell) scheme simulates more intense (weaker) low level wind than analysis in the tropics and subtropics. The objectives criteria of Sugahara (SJ) and Bonner (BJ) were used for LLJ identification. When applied to the observations, both criteria suggest a larger frequency of the SALLJ in Santa Cruz, followed by Mariscal, Trinidad and Asuncin. In Mariscal and Asuncin, the diurnal cycle indicates that SJ occurs mainly at 12 UTCs (morning), while the BJ criterion presents the SALLJ as more homogenously distributed. The concentration into two of the four-times-a-day observations does not allow conclusions about the diurnal cycle in Santa Cruz and Trinidad. The simulated wind profiles result in a lower than observed frequency of SALLJ using both the SJ and BJ criteria, with fewer events obtained with the BJ. Due to the stronger simulated winds, the Emanuel scheme produces an equal or greater relative frequency of SALLJ than the Grell scheme. However, the Grell scheme using the SJ criterion simulates the SALLJ diurnal cycle closer to the observed one. Although some discrepancies between observed and simulated mean vertical profiles of the horizontal wind are noted, there is large agreement between the composites of the vertical structure of the SALLJ, especially when the SJ criterion is used with the Grell scheme. On an intraseasonal scale, a larger southward displacement of SALLJ in February and December when compared with January has been noted. The Grell and Emanuel schemes simulated this observed oscillation in the low-level flow. However, the spatial pattern and intensity of rainfall and circulation anomalies simulated by the Grell scheme are closer to the analyses than those obtained with the Emanuel scheme.

The Impacts of Inter-El Nino Variability on the Tropical Atlantic and Northeast Brazil Climate

In this study, observations and numerical simulations are used to investigate how different El Nino events affect the development of SST anomalies in the Atlantic and how this relates to the Brazilian northeast (NE) precipitation. The results show that different types of El Nino have different impacts on the SST anomalies of the equatorial and tropical South Atlantic but a similar SST response in the tropical North Atlantic. Strong and long (weak and short) El Ninos with the main heating source located in the eastern (central) Pacific generate cold (warm) anomalies in the cold tongue and Benguela upwelling regions during boreal winter and spring. When the SST anomalies in the eastern equatorial and tropical South Atlantic are cold (warm), the meridional SST gradient across the equator is positive (negative) and the ITCZ is not allowed (allowed) to move southward during the boreal spring; as a consequence, the precipitation is below (above) the average over the NE. Thus, strong and long (weak and short) El Ninos are followed by dry (wet) conditions in the NE. During strong and long El Ninos, changes in the Walker circulation over the Atlantic and in the Pacific-South Atlantic (PSA) wave train cause easterly wind anomalies in the western equatorial Atlantic, which in turn activate the Bjerknes mechanism, establishing the cold tongue in boreal spring and summer. These easterly anomalies are also responsible for the Benguela upwelling. During short and weak El Ninos, westerly wind anomalies are present in the western equatorial Atlantic accompanied by warm anomalies in the eastern equatorial and tropical South Atlantic; a positive phase of the South Atlantic dipole develops during boreal winter. The simulations highlight the importance of ocean dynamics in establishing the correct slope of the equatorial thermocline and SST anomalies, which in turn determine the correct rainfall response over the NE.

Origin of Convectively Coupled Kelvin Waves over South America

Convectively coupled Kelvin waves over the South American continent are examined through the use of temporal and spatial filtering of reanalysis, satellite, and gridded rainfall data. They are most prominent from November to April, the season analyzed herein. The following two types of events are isolated: those that result from preexisting Kelvin waves over the eastern Pacific Ocean propagating into the continent, and those that apparently originate over Amazonia, forced by disturbances propagating equatorward from central and southern South America. The events with precursors in the Pacific are mainly upper-level disturbances, with almost no signal at the surface. Those events with precursors over South America, on the other hand, originate as upper-level synoptic wave trains that pass over the continent and resemble the ""cold surges`` documented by Garreaud and Wallace. As the wave train propagates over the Andes, it induces a southerly low-level wind that advects cold air to the north. Precipitation associated with a cold front reaches the equator a few days later and subsequently propagates eastward with the characteristics of a Kelvin wave. The structures of those waves originating over the Pacific are quite similar to those originating over South America as they propagate to eastern South America and into the Atlantic. South America Kelvin waves that originate over neither the Pacific nor the midlatitudes of South America can also be identified. In a composite sense, these form over the eastern slope of the Andes Mountains, close to the equator. There are also cases of cold surges that reach the equator yet do not form Kelvin waves. The interannual variability of the Pacific-originating events is related to sea surface temperatures in the central-eastern Pacific Ocean. When equatorial oceanic conditions are warm, there tends to be an increase in the number of disturbances that reach South America from the Pacific.

The role of the South Indian and Pacific oceans in South American monsoon variability

This work has investigated the impact of three different low-frequency sea surface temperature (SST) variability modes located in the Indian and the Pacific Oceans on the interannual variability of the South American Monsoon System (SAMS) using observed and numerical data. Rotated Empirical Orthogonal Function (REOF) analysis and numerical simulations with a General Circulation Model (GCM) were used. One of the three SST variability modes is located close to southeastern Africa. According to the composites, warmer waters over this region are associated with enhanced austral summer precipitation over the sub-tropics. The GCM is able to reproduce this anomalous precipitation pattern, simulating a wave train emanating from the Indian Ocean towards South America (SA). A second SST variability mode was located in the western Pacific Ocean. REOF analysis indicates that warmer waters are associated with drought conditions over the South Atlantic Convergence Zone (SACZ) and enhanced precipitation over the sub-tropics. The GCM indicates that the warmer waters over Indonesia generate drought conditions over tropical SA through a Pacific South America-like (PSA) wave pattern emanating from the western Pacific. Finally, the third SST variability mode is located over the southwestern South Pacific. The composites indicate that warmer waters are associated with enhanced precipitation over the SACZ and drought conditions over the sub-tropics. There is a PSA-like wave train emanating from Indonesia towards SA, and another crossing the Southern Hemisphere in the extra-tropics, probably associated with transient activity. The GCM is able to reproduce the anomalous precipitation pattern, although it is weaker than observed. The PSA-like pattern is simulated, but the model fails in reproducing the extra-tropical wave activity.

Precipitation variability in Sao Paulo State, Brazil

The State of Sao Paulo is the richest in Brazil, responsible for over 30% of the Brazilian gross rate. It has a population of around 30 million and its economy is based on agriculture and industrial products. Any change in climate can have a profound influence on the socio-economics of the State. In order to determine changes in total and extreme rainfall over Sao Paulo State, climate change indices derived from daily precipitation data were calculated using specially designed software. Maps of trends for a subset of 59 rain gauge stations were analysed for the period 1950-1999 and also for a subset of this period, 1990-1999, representing more recent climate. A non-parametric Mann-Kendall test was applied to the time series. Maps of trends for six annual precipitation indices (annual total precipitation (PRCPTOT), very heavy precipitation days (R20mm), events greater than the 95th percentile (R95p), maximum five days precipitation total (RX5day), the length of the largest wet spell (CWD) and the length of the largest dry spell (CDD)) were analysed for the entire period. These exhibited statistically significant trends associated with a wetter climate. A significant increase in PRCPTOT, associated with very heavy precipitation days, were observed at more than 45% of the rain gauge stations. The Mann-Kendall test identified that the positive trend in PRCPTOT is possibly related to the increase in the R95p and R20mm indices. Therefore, the results suggest that there has been a change in precipitation intensity. In contrast, the indices for the more recent shorter time series are significantly different to the longer term indices. The results indicate that intense precipitation is becoming concentrated in a few days and spread over the period when the CDD and R20mm indices show positive trends, while negative ones are seen in the RX5day index. The trends found could be related to many anthropogenic aspects such as biomass burning aerosols and land use.