Segregation of migration by feeding ground origin in North Atlantic humpback whales (Megaptera novaeangliae)

Results from a large-scale, capture–recapture study of humpback whales Megaptera novaeangliae in the North Atlantic show that migration timing is influenced by feeding ground origin. No significant differences were observed in the number of individuals from any feeding area that were re-sighted in the common breeding area in the West Indies. However, there was a relationship between the proportion (logit transformed) of West Indies sightings and longitude (r2 = 0.97, F1,3 = 98.27, P = 0.0022) suggesting that individuals feeding farther to the east are less likely to winter in the West Indies. A relationship was also detected between sighting date in the West Indies and feeding area. Mean sighting dates in the West Indies for individuals identified in the Gulf of Maine and eastern Canada were significantly earlier than those for animals identified in Greenland, Iceland and Norway (9.97 days, t179 = 3.53, P = 0.00054). There was also evidence for sexual segregation in migration; males were seen earlier on the breeding ground than were females (6.63 days, t105 = 1.98, P = 0.050). This pattern was consistently observed for animals from all feeding areas; a combined model showed a significant effect for both sex (F1 = 5.942, P = 0.017) and feeding area (F3 =4.756, P=0.0038). The temporal difference in occupancy of the West Indies between individuals from different feeding areas, coupled with sexual differences in migratory patterns, presents the possibility that there are reduced mating opportunities between individuals from different high latitude areas.

Notes on Microtine Rodents from the Brooks Range, Arctic Alaska

In connection with parasitological studies carried on during 1949 for the U.S. Public Health Service, the writer collected a considerable number of mammals from the northern edge of the "Endicott" section of the Brooks Range, in Arctic Alaska. The mammalian fauna of this region is poorly known, since apparently no previous collecting has been done here. About 200 microtine rodents were collected, mostly near Tolugak Lake (latitude 68° 24' N, longitude 151° 26' W), near the head of the Anaktuvuk River Valley. A few specimens were also taken at Umiat, on the Colville River, about 80 miles north of Tolugak Lake (latitude 69° 23' N, longitude 152° 10' W). Five species are represented: Clethrionomys rutilis dawsoni (Merriam), Microtus oeconomus macfarlani Merriam, Microtus miurus paneaki, n. subsp., Lemmus trimucronatus alascensis Merriam and Dicrostonyx groenlandicus rubricatus (Richardson). More complete details concerning their ecology and reproduction will be presented in a later paper, at which time the other mammals obtained will also be considered. The specimens have been deposited in the U.S. National Museum.

Validation of the MAPGEO2010 and comparison with recent geopotential models

This paper presents the offorts to calculate the geoid model for Brazil. It is limited by 6 degrees N and 35 degrees S in latitude and 30 degrees W and 75 degrees W in longitude. The terrestrial gravity data for the continent have been updated by means of the most recent surveys in Brazil and in the neighbour countries. The complete Bouguer and Helmert gravity anomalies have been derived through the Canadian package SHGEO. The short wavelength component was estimated via FFT. The geopotential model EGM2008 was used as a reference field restricted to degree and order 150. The model was validated over 844 GPS observations on Bench Marks of the spirit leveling network. The height anomalies plus a topography dependent correction term derived from EGM2008 (degree 2190 and order 2159), GO_CONS_GCF_2_DIR_R2 (degree and order 240), GOCO02S (degree and order 250), EIGEN 51C (degree and order 359) and EIGEN 6C (degree and order 1420), geoidal height derived from MAPGEO2004 (old official geoid model in Brazil) have also been compared to the GPS points on Bench Marks.

Genotype by environment interaction for birth and weaning weights of composite beef cattle in different regions of Brazil

The objectives of the present study were to characterize and define homogenous production environments of composite beef cattle in Brazil in terms of climatic and geographic variables using multivariate exploratory techniques and to use them to assess the presence of G x E for birth weight (BW) and weaning weight (WW). Data from animals born between 1995 and 2008 on 36 farms located in 27 municipalities of the Brazilian states were used. Fifteen years of climate observations (mean minimum and maximum annual temperature and mean annual rainfall) and geographic (latitude, longitude and altitude) data were obtained for each municipality where the farms were located for characterization of the production environments. Hierarchical and nonhierarchical cluster analysis was used to group farms located in regions with similar environmental variables into clusters. Six clusters of farms were formed. The effect of sire-cluster interaction was tested by single-trait analysis using deviance information criterion (DIC). Genetic parameters were estimated by multi-trait analysis considering the same trait to be different in each cluster. According to the values of DIC, the inclusion of sire-cluster effect did not improve the fit of the genetic evaluation model for BW and WW. Estimates of genetic correlations among clusters ranged from -0.02 to 0.92. The low genetic correlation among the most studied regions permits us to suggest that a separate genetic evaluation for some regions should be undertaken. (C) 2012 Elsevier B.V. All rights reserved.

Estabilidade temporal e variabilidade espacial da distribuição da armazenagem de água no solo numa sucessão feijão/aveia-preta

A armazenagem de água no solo é muito variável no tempo e no espaço devido à influência de vários fatores ambientais e topográficos. Objetivou-se neste estudo: quantificar a armazenagem de água no solo; avaliar sua estabilidade temporal e sua variabilidade espacial em um local próximo e outro distante do sistema radicular numa sucessão feijão/aveia-preta; e constatar sua variabilidade espacial em função do relevo. Em área experimental de 1.500 m², situada em Piracicaba-SP (latitude de 22° 42' 30" S, longitude de 47° 38' 00" W e 546 m de altitude), estabeleceram-se 60 pontos de amostragem, distanciados entre si de 5 m, numa grade de 10 por 6 pontos (50 x 30 m). Os valores da armazenagem de água no solo apresentaram comportamento-padrão para o solo estudado, com valores maiores em profundidade do que na camada mais superficial. Houve maior estabilidade temporal da armazenagem de água no solo durante o período chuvoso, porém na fase de secagem do solo a estabilidade temporal também foi constatada, mas com valores de coeficiente de correlação mais elevados na camada de 0,0-0,80 m, mostrando que essa estabilidade foi claramente devida à posição topográfica desses pontos, os quais estão localizados na porção mais baixa do relevo. Assim, o ponto 52 foi escolhido como representativo da média na fase de recarga de água no solo em ambas as camadas estudadas, e o ponto 46, na camada superior durante a fase de secagem do solo. A dependência espacial da variabilidade da armazenagem de água no solo foi detectada em ambas as camadas de solo, porém com maior alcance na camada de 0,0-0,40 m, quando houve ocorrência de precipitação. Quando as chuvas cessaram, o alcance foi maior na camada de 0,0-0,80 m. Durante o período chuvoso, o padrão de variabilidade espacial foi muito semelhante em ambas as camadas de solo, com armazenagens maiores na faixa de maior inclinação e depressões do terreno.

Tidal perturbations of the rotation of the moon: the creep tide approach

In this communication we report results from the application to the study of the rotation of the Moon of the creeping tide theory just proposed (Ferraz-Mello, Cel. Mech. Dyn. Astron., submitted. ArXiv astro-ph 1204.3957). The choice of the Moon for the first application of this new theory is motivated by the fact that the Moon is one of the best observed celestial bodies and the comparison of the theoretical predictions of the theory with observations i may validate the theory or point out the need of further improvements. Particularly, the tidal perturbations of the rotation of the Moon - the physical libration of the Moon - have been detected in the Lunar Laser Ranging measurements (Williams et al. JGR 106, 27933, 2001). The major difficulty in this application comes from the fact that tidal torques in a planet-satellite system are very sensitive to the distance between the two-bodies, which is strongly affected by Solar perturbations. In the case of the Moon, the main solar perturbations - the Evection and the Variation - are more important than most of the Keplerian oscillations, being smaller only than the first Keplerian harmonic (equation of the centre). Besides, two of the three components of the Moon's libration in longitude whose tidal contributions were determined by LLR are related to these perturbations. The results may allow us to determine the main parameter of a possible Moon's creeping tide. The preliminary results point to a relaxation factor (gamma) 2 to 4 times smaller than the one predicted from the often cited values of thr Moon's quality factor Q (between 30 and 40), and points to larger Q values.

Estimating the monthly pCO2 distribution in the North Atlantic using a self-organizing neural network

[EN] Here we present monthly, basin-wide maps of the partial pressure of carbon dioxide (pCO2) for the North Atlantic on a latitude by longitude grid for years 2004 through 2006 inclusive. The maps have been computed using a neural network technique which reconstructs the non-linear relationships between three biogeochemical parameters and marine pCO2. A self organizing map (SOM) neural network has been trained using 389 000 triplets of the SeaWiFSMODIS chlorophyll-a concentration, the NCEP/NCAR reanalysis sea surface temperature, and the FOAM mixed layer depth. The trained SOM was labelled with 137 000 underway pCO2 measurements collected in situ during 2004, 2005 and 2006 in the North Atlantic, spanning the range of 208 to 437atm. The root mean square error (RMSE) of the neural network fit to the data is 11.6?atm, which equals to just above 3 per cent of an average pCO2 value in the in situ dataset. The seasonal pCO2 cycle as well as estimates of the interannual variability in the major biogeochemical provinces are presented and discussed. High resolution combined with basin-wide coverage makes the maps a useful tool for several applications such as the monitoring of basin-wide air-sea CO2 fluxes or improvement of seasonal and interannual marine CO2 cycles in future model predictions. The method itself is a valuable alternative to traditional statistical modelling techniques used in geosciences.

Automatische objektive Identifikation und Berechnung der Eigenschaften von Rossbywellenzügen

In dieser Arbeit wird eine Methode entwickelt, die Rossbywellenzüge automatisch identifiziert und und deren Eigenschaften quantifiziert. Mit dieser Methode wird der Wellenzug als eine Einheit in Raum und Zeit interpretiert. Einheit im Raum heißt, dass nicht die einzelnen Tröge und Rücken eines Wellenzugs betrachtet werden, sondern deren Einhüllende. Einheit in der Zeit bedeutet, dass der Wellenzug nicht nur zu einem Zeitpunkt betrachtet wird, sondern über seine gesamte Lebensdauer hinweg. Um den Wellenzug als räumliche und zeitliche Einheit zu erhalten, werden die Einhüllenden der Wellenzüge in Längengrad-Zeit Diagrammen, sogenannten Hovmöllerdiagrammen, betrachtet. Dort werden zusammenhängende Regionen als Objekte, die jeweils einen Wellenzug repräsentieren, identifiziert. Deren Eigenschaften werden dann automatisch berechnet. Diese Eigenschaften können nun direkt dem zugrunde liegenden Rossbywellenzug zugeordnet werden.rnDie neue Methode wird in zwei verschiedenen Szenarien angewendet: erstens zur Beurteilung der Vorhersagequalität eines einzelnen Rossbywellenzugs und zweitens für die klimatologische Betrachtung von Rossbywellenzügen im ERA-40 Reanalysedatensatz. Sie wurde weiterhin mit bisher verwendeten Methoden zur Identifikation und Quantifizierung von Rossbywellenzügen verglichen.rnDie Untersuchung der Vorhersagequalität ergab, dass in dem betrachteten Fall die Übereinstimmung der Vorhersage mit der Analyse des Wellenzugs gering war, sofern das Modell initialisiert wurde, bevor der Rossbywellenzug eingesetzt hatte. Im Gegensatz dazu nahm die Vorhersagequalität deutlich zu, wenn der Wellenzug bereits in den Vorhersagedaten enthalten war. Dies deutet darauf hin, dass es in dem vorliegenden Fall problematisch ist, mit dem Modell den Auslösemechanismus korrekt voherzusagen. Für die weitere Untersuchung der Vorhersagequalität wurde eine spezielle Art der Darstellung der Daten verwendet, mit deren Hilfe deutlich wurde, dass das verwendete Modell in der Lage ist, diesen Wellenzug ungefähr sechs Tage im Voraus vorherzusagen. Diese Zeitspanne ist deutlich kürzer als die Lebensdauer des Wellenzugs, die etwa 10 Tage beträgt.rnIm Rahmen der klimatologischen Studie ergab sich eine positive Korrelation zwischen der Lebensdauer eines Rossbywellenzugs und des Bereichs den dieser Wellenzug während seiner gesamten Existenz in zonaler Richtung überstreicht. Für Wellenzüge mit einer kurzen Lebensdauer ergab sich eine ebenfalls positive Korrelation zwischen der mittleren Amplitude und der Dauer des Wellenzugs. Für eine längere Lebensdauer geht diese Korrelation aber in eine Sättigung über und die mittlere Amplitude steigt nicht mehr weiter an. Als eine mögliche Erklärung für dieses Verhalten wird angeführt, dass eine gewisse Stärke der Amplitude benötigt wird um stromabwärtige Entwicklung zu erhalten aber zu große Amplituden im Allgemeinen zum Brechen der Welle führen. Das Brechen leitet den letzten Abschnitt im Lebenszyklus eines Rossbywellenzuges ein, welcher im Anschluss meist zerfällt. Ein weiteres Ergebnis der klimatologischen Untersuchung ist das Auffinden bevorzugter Regionen der Entstehung und des Abklingens von Rossbywellenzügen. Diese Regionen unterscheiden sich erheblich für Rossbywellenzüge unterschiedlicher minimaler Lebensdauer. Langlebige Rossbywellenzüge entstehen demnach hauptsächlich über Ostasien und dem Westpazifik und vergehen dann über Europa.rnSchließlich wurde die entwickelte Methode in einen systematischen Vergleich anderer Methoden zur Identifikation und Quantifizierung von Rossbywellenzügen eingereiht. Die betrachteten Methoden beinhalten verschiedene Trog-und-Rücken Hovmöllerdiagramme des Meridionalwindes, Methoden die Rossbywellenzüge als eine Einheit identifizieren und Methoden die den Beitrag verschiedener physikalischer Aspekte zu der Entwicklung von Rossbywellenenzügen quantifizieren. Der Vergleich macht deutlich, dass jede Methode ihre individuellen Stärken und Schwächen hat. Dies bedeutet insbesondere, dass die Eignung der Methode von dem Stadium des Lebenszyklus, in dem sich der Rossbywellenzug befindet und dem Fokus, den man bei der Betrachtung hat, abhängt. Ideal ist eine Kombination mehrerer Methoden, da dies ein vollständigeres Bild eines Rossbywellenzuges ergibt als einzelne Methoden es zu liefern vermögen. Obwohl alle Methoden für die Anwendungen, für die sie jeweils konzipiert wurden, geeignet sind, ergeben sich bei der Diagnose der Rossbywellenzüge beträchtliche Unterschiede. Letztendlich stellt sich heraus, dass sogar die Defintion eines Rossbywellenzugs bis zu einem gewissen Grad von der zu seiner Identifizierung verwendeten Methode abhängt.

Polychaete fauna of the Northwest Portuguese Coastal Shelf: ecology, diversity and distribution

Polychaetes are one of the larger groups of macroinvertebrates with more than 9000 species recognised, distributed worldwide. Thanks to the broad ecological adaptability and high abundaces, this taxon plays a leading role and is considered an important component of all benthic assemblages. Our knowledge about the West Iberian Coast polychaete fauna are scarce, and the only studies are recent. In this sense, the aim of this work was to investigate the composition and the spatial distribution of the polychaete fauna along the NW Portuguese Coastal Shelf, focusing on their relationship to environmental factors (depth, grain size, longitude and latitude) and to add new data to the existing biological dataset. A total of 39 sites were analysed, collected in an area of about 5665 km², between 20 and 150 m depth, distributed in a way to cover the overall grain size gradient. A total of 9352 specimens belonging to 41 families were found, and the analysis based on the abundance of polychaete species revealed five affinity groups: (a) nearshore medium sand characterised by Pisione parapari and Hesionura elongata; (b) very coarse sand that showed the highest abundance of Syllidae and was characterised by Protodorvillea kefersteini and Syllis garciai; (c) fine sand dominated by Spiophanes bombyx and Glycera tridactyla; (d) very fine sand with Nepthys assimilis and Amage sp. and the highest abundance of Paraonidae; (d) mud characterised by Labioleanira yhleni and Ampharete finmarchica. The combination of the environmental variables and the biological data, done with BIOENV routine, demonstrated that depth, grain size and fine contents were the best related with the biological data (rho=0.598). In general, the results agree with the composition and the spatial distribution of the polychaete fauna in other parts of the world; further polychaete assemblages related to mud sediments were firstly recorded in the Northwestern Portuguese Coastal Shelf.

Using Micro-Gravity Techniques to Map Alluvium Thickness and Pleistocene Location of the West Branch of the Susquehanna River Near Muncy, Pennsylvania

Laurentide glaciation during the early Pleistocene (~970 ka) dammed the southeast-flowing West Branch of the Susquehanna River (WBSR), scouring bedrock and creating 100-km-long glacial Lake Lesley near the Great Bend at Muncy, Pennsylvania (Ramage et al., 1998). Local drill logs and well data indicate that subsequent paleo-outwash floods and modern fluvial processes have deposited as much as 30 meters of alluvium in this area, but little is known about the valley fill architecture and the bedrock-alluvium interface. By gaining a greater understanding of the bedrock-alluvium interface the project will not only supplement existing depth to bedrock information, but also provide information pertinent to the evolution of the Muncy Valley landscape. This project determined if variations in the thickness of the valley fill were detectable using micro-gravity techniques to map the bedrock-alluvium interface. The gravity method was deemed appropriate due to scale of the study area (~30 km2), ease of operation by a single person, and the available geophysical equipment. A LaCoste and Romberg Gravitron unit was used to collect gravitational field readings at 49 locations over 5 transects across the Muncy Creek and Susquehanna River valleys (approximately 30 km2), with at least two gravity base stations per transect. Precise latitude, longitude and ground surface elevation at each location were measured using an OPUS corrected Trimble RTK-GPS unit. Base stations were chosen based on ease of access due to the necessity of repeat measurements. Gravity measurement locations were selected and marked to provide easy access and repeat measurements. The gravimeter was returned to a base station within every two hours and a looping procedure was used to determine drift and maximize confidence in the gravity measurements. A two-minute calibration reading at each station was used to minimize any tares in the data. The Gravitron digitally recorded finite impulse response filtered gravity measurements every 20 seconds at each station. A measurement period of 15 minutes was used for each base station occupation and a minimum of 5 minutes at all other locations. Longer or multiple measurements were utilized at some sites if drift or other externalities (i.e. train or truck traffic) were effecting readings. Average, median, standard deviation and 95% confidence interval were calculated for each station. Tidal, drift, latitude, free-air, Bouguer and terrain corrections were then applied. The results show that the gravitational field decreases as alluvium thickness increases across the axes of the Susquehanna River and Muncy Creek valleys. However, the location of the gravity low does not correspond with the present-day location of the West Branch of the Susquehanna River (WBSR), suggesting that the WBSR may have been constrained along Bald Eagle Mountain by a glacial lobe originating from the Muncy Creek Valley to the northeast. Using a 3-D inversion model, the topography of the bedrock-alluvium interface was determined over the extent of the study area using a density contrast of -0.8 g/cm3. Our results are consistent with the bedrock geometry of the area, and provide a low-cost, non-invasive and efficient method for exploring the subsurface and for supplementing existing well data.

Social Stratification of the Alans of North Caucasus (IV-IX Cents. A.D.)

Dmitry Korobov. Social Stratification of the Alans of the North Caucasus (4th-9th Centuries A.D.) There were two stages to this project: studies of the differences between local groups of catacomb cemeteries in the North Caucasus, and an analysis of social ranks among the Alans on the basis of the large cemeteries. The first stage involved a cluster analysis of 100 cemeteries on the basis of the variables of longitude and latitude, which made it possible to identify 11 local groups of cemeteries. This was followed by a study of the difference between burial rites charcterising flat graves and barrow burials and between catacombs of different periods within these two main groups. Korobov then looked at the difference between the catacombs of various local groups. The second step was an investigation of the social structure of the Alanic tribes. Here Korobov first selected gender-specific grave goods of adult burials on the basis of large cemeteries, and compared single male and female burials on the basis of the associated grave goods. He then used a cluster analysis to draw up a subdivision of male and female burials according to their grave goods, and interpreted the results obtained with the help of ethnographic and written sources.

SOLAR RADIATION PRESSURE AND LOCAL INTERSTELLAR MEDIUM FLOW PARAMETERS FROM INTERSTELLAR BOUNDARY EXPLORER LOW ENERGY HYDROGEN MEASUREMENTS

Neutral hydrogen atoms that travel into the heliosphere from the local interstellar medium (LISM) experience strong effects due to charge exchange and radiation pressure from resonant absorption and re-emission of Lyα. The radiation pressure roughly compensates for the solar gravity. As a result, interstellar hydrogen atoms move along trajectories that are quite different than those of heavier interstellar species such as helium and oxygen, which experience relatively weak radiation pressure. Charge exchange leads to the loss of primary neutrals from the LISM and the addition of new secondary neutrals from the heliosheath. IBEX observations show clear effects of radiation pressure in a large longitudinal shift in the peak of interstellar hydrogen compared with that of interstellar helium. Here, we compare results from the Lee et al. interstellar neutral model with IBEX-Lo hydrogen observations to describe the distribution of hydrogen near 1 AU and provide new estimates of the solar radiation pressure. We find over the period analyzed from 2009 to 2011 that radiation pressure divided by the gravitational force (μ) has increased slightly from μ = 0.94 ± 0.04 in 2009 to μ = 1.01 ± 0.05 in 2011. We have also derived the speed, temperature, source longitude, and latitude of the neutral H atoms and find that these parameters are roughly consistent with those of interstellar He, particularly when considering the filtration effects that act on H in the outer heliosheath. Thus, our analysis shows that over the period from 2009 to 2011, we observe signatures of neutral H consistent with the primary distribution of atoms from the LISM and a radiation pressure that increases in the early rise of solar activity.

Volcanic forcing for climate modeling: a new microphysics-based data set covering years 1600–present

As the understanding and representation of the impacts of volcanic eruptions on climate have improved in the last decades, uncertainties in the stratospheric aerosol forcing from large eruptions are now linked not only to visible optical depth estimates on a global scale but also to details on the size, latitude and altitude distributions of the stratospheric aerosols. Based on our understanding of these uncertainties, we propose a new model-based approach to generating a volcanic forcing for general circulation model (GCM) and chemistry–climate model (CCM) simulations. This new volcanic forcing, covering the 1600–present period, uses an aerosol microphysical model to provide a realistic, physically consistent treatment of the stratospheric sulfate aerosols. Twenty-six eruptions were modeled individually using the latest available ice cores aerosol mass estimates and historical data on the latitude and date of eruptions. The evolution of aerosol spatial and size distribution after the sulfur dioxide discharge are hence characterized for each volcanic eruption. Large variations are seen in hemispheric partitioning and size distributions in relation to location/date of eruptions and injected SO2 masses. Results for recent eruptions show reasonable agreement with observations. By providing these new estimates of spatial distributions of shortwave and long-wave radiative perturbations, this volcanic forcing may help to better constrain the climate model responses to volcanic eruptions in the 1600–present period. The final data set consists of 3-D values (with constant longitude) of spectrally resolved extinction coefficients, single scattering albedos and asymmetry factors calculated for different wavelength bands upon request. Surface area densities for heterogeneous chemistry are also provided.

Estimating Spatial Distribution of American Lobster Homarus Americanus Using Habitat Variables

The spatial distribution of the American lobster Homarus americanus is influenced by many factors, which are often difficult to quantify. We implemented a modeling approach for quantifying season-, size-, and sex-specific lobster spatial distribution in the Gulf of Maine with respect to environmental and spatial variables including bottom temperature, bottom salinity, latitude, longitude, depth, distance offshore, and 2 substratum features. Lobster distribution was strongly associated with temperature and depth, and differed seasonally by sex. In offshore waters in the fall, females were dominant at higher latitudes and males at lower latitudes. This segregation was not apparent in the spring although females were still dominant at higher latitudes in offshore waters. Juveniles and adults were also distributed differently; juveniles were more abundant at the lower latitudes in inshore waters, while adults were more widespread along the entire coast. These patterns are consistent with the ecology of the American lobster. This study provides a tool to evaluate changes in lobster spatial distribution with respect to changes in key habitat and other environmental variables, and consequently could be of value for the management of the American lobster.

Migration and population expansion of Abies, Fagus, Picea, and Quercus since 15000 years in and across the Alps, based on pollen-percentage threshold values

Aims: The aim of this study is to explore the migration (colonization of new areas) and subsequent population expansion (within an area) since 15 ka cal BP of Abies, Fagus, Picea, and Quercus into and through the Alps solely on the basis of high-quality pollen data. Methods: Chronologies of 101 pollen sequences are improved or created. Data from the area delimited by 45.5–48.1°N and 6–14°E are summarized in three ways: (1) in a selection of pollen-percentage threshold maps (thresholds 0.5%, 1%, 2%, 4%, 8%, 16%, and 32% of land pollen); (2) in graphic summaries of 250-year time slices and geographic segments (lengthwise and transverse in relation to the main axis of the Alps) as pollen-percentage curves, pollen-percentage difference curves, and pollen-percentage threshold ages cal BP graphed against both the length and the transverse Alpine axes; and (3) in tables showing statistical relationships of either pollen-percentage threshold ages cal BP or pollen expansion durations (=time lapse between different pollen-percentage threshold ages cal BP) with latitude, longitude, and elevation; to establish these relationships we used both simple linear regression and multiple linear regression after stepwise-forward selection. Results: The statistical results indicate that (a) the use of pollen-percentage thresholds between 0.5% and 8% yield mostly similar directions of tree migration, so the method is fairly robust, (b) Abies migrated northward, Fagus southward, Picea westward, and Quercus northward; more detail does not emerge due to an extreme scarcity of high-quality data especially along the southern foothills of the Alps and in the eastern Alps. This scarcity allows the reconstruction of one immigration route only of Abies into the southern Alps. The speed of population expansion (following arrival) of Abies increased and of Picea decreased during the Holocene, of Fagus it decreased especially during the later Holocene, and of Quercus it increased especially at the start of the Holocene.