Registro de excentricidad de alta frecuencia en los sedimentos continentales de la Cuenca de Guadix (Cordillera Bética, España)

Los sedimentos continentales (Plioceno-Cuaternario) que afloran en el sector central de la Cuenca de Guadix (Sur de España) muestran una ciclicidad de 100 ka consistente en la alternancia de depósitos de abanicos aluviales y sedimentos fluvio-lacustres. Durante el Plioceno y el Pleistoceno la Cuenca de Guadix era endorreica, y se caracterizaba por la existencia de un sistema axial fluvial y una orla marginal de abanicos aluviales transversales. En la zona de estudio, estos sistemas estaban relacionados lateralmente, ocupando de forma alterna el valle axial en el sector central de la Cuenca de Guadix. La edad estimada para la alternancia, ca. 100 ka, cae en la banda de excentricidad de alta frecuencia de Milankovitch. Estas fases podrían interpretarse como el resultado de máximos de excentricidad (inviernos más largos y fríos, con mayor volumen de precipitaciones que favorecerían las progradaciones de los abanicos) o de excentricidad mínima (períodos más secos y fríos, con una cubierta vegetal más escasa en las áreas fuente y, por tanto, un mayor aporte de sedimento por precipitaciones muy concentradas en el tiempo al sistema aluvial). Se muestra cómo los datos paleomagnéticos no son lo suficientemente precisos para proporcionar una buena correlación de las fases de progradación con la curva de excentricidad de Laskar et al. (2004), por lo que se pone en duda su precisión a la hora de determinar el significado climático de la ciclicidad.

Astrophysical parameters and orbital solution of the peculiar X-ray transient IGR J00370+6122

Context. BD + 60° 73 is the optical counterpart of the X-ray source IGR J00370+6122, a probable accretion-powered X-ray pulsar. The X-ray light curve of this binary system shows clear periodicity at 15.7 d, which has been interpreted as repeated outbursts around the periastron of an eccentric orbit. Aims. We aim to characterise the binary system IGR J00370+6122 by deriving its orbital and physical parameters. Methods. We obtained high-resolution spectra of BD + 60° 73 at different epochs. We used the fastwind code to generate a stellar atmosphere model to fit the observed spectrum and obtain physical magnitudes. The synthetic spectrum was used as a template for cross-correlation with the observed spectra to measure radial velocities. The radial velocity curve provided an orbital solution for the system. We also analysed the RXTE/ASM and Swift/BAT light curves to confirm the stability of the periodicity. Results. BD + 60° 73 is a BN0.7 Ib low-luminosity supergiant located at a distance ~3.1 kpc, in the Cas OB4 association. We derive Teff = 24 000 K and log gc = 3.0, and chemical abundances consistent with a moderately high level of evolution. The spectroscopic and evolutionary masses are consistent at the 1-σ level with a mass M∗ ≈ 15 M⊙. The recurrence time of the X-ray flares is the orbital period of the system. The neutron star is in a high-eccentricity (e = 0.56 ± 0.07) orbit, and the X-ray emission is strongly peaked around orbital phase φ = 0.2, though the observations are consistent with some level of X-ray activity happening at all orbital phases. Conclusions. The X-ray behaviour of IGR J00370+6122 is reminiscent of “intermediate” supergiant X-ray transients, though its peak luminosity is rather low. The orbit is somewhat wider than those of classical persistent supergiant X-ray binaries, which when combined with the low luminosity of the mass donor, explains the low X-ray luminosity. IGR J00370+6122 will very likely evolve towards a persistent supergiant system, highlighting the evolutionary connection between different classes of wind-accreting X-ray sources.

Measuring the stellar wind parameters in IGR J17544-2619 and Vela X-1 constrains the accretion physics in supergiant fast X-ray transient and classical supergiant X-ray binaries

Context. Classical supergiant X-ray binaries (SGXBs) and supergiant fast X-ray transients (SFXTs) are two types of high-mass X-ray binaries (HMXBs) that present similar donors but, at the same time, show very different behavior in the X-rays. The reason for this dichotomy of wind-fed HMXBs is still a matter of debate. Among the several explanations that have been proposed, some of them invoke specific stellar wind properties of the donor stars. Only dedicated empiric analysis of the donors’ stellar wind can provide the required information to accomplish an adequate test of these theories. However, such analyses are scarce. Aims. To close this gap, we perform a comparative analysis of the optical companion in two important systems: IGR J17544-2619 (SFXT) and Vela X-1 (SGXB). We analyze the spectra of each star in detail and derive their stellar and wind properties. As a next step, we compare the wind parameters, giving us an excellent chance of recognizing key differences between donor winds in SFXTs and SGXBs. Methods. We use archival infrared, optical and ultraviolet observations, and analyze them with the non-local thermodynamic equilibrium (NLTE) Potsdam Wolf-Rayet model atmosphere code. We derive the physical properties of the stars and their stellar winds, accounting for the influence of X-rays on the stellar winds. Results. We find that the stellar parameters derived from the analysis generally agree well with the spectral types of the two donors: O9I (IGR J17544-2619) and B0.5Iae (Vela X-1). The distance to the sources have been revised and also agree well with the estimations already available in the literature. In IGR J17544-2619 we are able to narrow the uncertainty to d = 3.0 ± 0.2 kpc. From the stellar radius of the donor and its X-ray behavior, the eccentricity of IGR J17544-2619 is constrained to e< 0.25. The derived chemical abundances point to certain mixing during the lifetime of the donors. An important difference between the stellar winds of the two stars is their terminal velocities (ν∞ = 1500 km s-1 in IGR J17544-2619 and ν∞ = 700 km s-1 in Vela X-1), which have important consequences on the X-ray luminosity of these sources. Conclusions. The donors of IGR J17544-2619 and Vela X-1 have similar spectral types as well as similar parameters that physically characterize them and their spectra. In addition, the orbital parameters of the systems are similar too, with a nearly circular orbit and short orbital period. However, they show moderate differences in their stellar wind velocity and the spin period of their neutron star which has a strong impact on the X-ray luminosity of the sources. This specific combination of wind speed and pulsar spin favors an accretion regime with a persistently high luminosity in Vela X-1, while it favors an inhibiting accretion mechanism in IGR J17544-2619. Our study demonstrates that the relative wind velocity is critical in class determination for the HMXBs hosting a supergiant donor, given that it may shift the accretion mechanism from direct accretion to propeller regimes when combined with other parameters.

Effects of Interspecific Gene Flow on the Phenotypic Variance-Covariance Matrix in Lake Victoria Cichlids

Quantitative genetics theory predicts adaptive evolution to be constrained along evolutionary lines of least resistance. In theory, hybridization and subsequent interspecific gene flow may however rapidly change the evolutionary constraints of a population and eventually change its evolutionary potential, but empirical evidence is still scarce. Using closely related species pairs of Lake Victoria cichlids sampled from four different islands with different levels of interspecific gene flow, we tested for potential effects of introgressive hybridization on phenotypic evolution in wild populations. We found that these effects differed among our study species. Constraints measured as the eccentricity of phenotypic variance-covariance matrices declined significantly with increasing gene flow in the less abundant species for matrices that have a diverged line of least resistance. In contrast we find no such decline for the more abundant species. Overall our results suggest that hybridization can change the underlying phenotypic variance-covariance matrix, potentially increasing the adaptive potential of such populations.

Phase relationships between orbital forcing and the composition of air trapped in Antarctic ice cores

Orbital tuning is central for ice core chronologies beyond annual layer counting, available back to 60 ka (i.e. thousands of years before 1950) for Greenland ice cores. While several complementary orbital tuning tools have recently been developed using δ¹⁸Oatm, δO₂⁄N₂ and air content with different orbital targets, quantifying their uncertainties remains a challenge. Indeed, the exact processes linking variations of these parameters, measured in the air trapped in ice, to their orbital targets are not yet fully understood. Here, we provide new series of δO₂∕N₂ and δ¹⁸Oatm data encompassing Marine Isotopic Stage (MIS) 5 (between 100 and 160 ka) and the oldest part (340–800 ka) of the East Antarctic EPICA Dome C (EDC) ice core. For the first time, the measurements over MIS 5 allow an inter-comparison of δO₂∕N₂ and δ¹⁸Oatm records from three East Antarctic ice core sites (EDC, Vostok and Dome F). This comparison highlights some site-specific δO₂∕N₂ variations. Such an observation, the evidence of a 100 ka periodicity in the δO₂∕N₂ signal and the difficulty to identify extrema and mid-slopes in δO2∕N2 increase the uncertainty associated with the use of δO₂∕N₂ as an orbital tuning tool, now calculated to be 3–4 ka. When combining records of δ¹⁸Oatm and δO₂∕N₂ from Vostok and EDC, we find a loss of orbital signature for these two parameters during periods of minimum eccentricity (∼ 400 ka, ∼ 720–800 ka). Our data set reveals a time-varying offset between δO₂∕N₂ and δ¹⁸Oatm records over the last 800 ka that we interpret as variations in the lagged response of δ¹⁸Oatm to precession. The largest offsets are identified during Terminations II, MIS 8 and MIS 16, corresponding to periods of destabilization of the Northern polar ice sheets. We therefore suggest that the occurrence of Heinrich–like events influences the response of δ¹⁸Oatm to precession.

Geochemistry of Eocene/Oligocene sediments of ODP Hole 154-925A

Paleoproductivity, nutrient burial, and carbon cycling were investigated across the Eocene/Oligocene (E/O) boundary (begin to end; 36.9-32.7 Ma at ~40 kyr resolution, timescale of Shackleton et al. (1999, doi:10.1098/rsta.1999.0407) at Ocean Drilling Program Site 925 on the Ceara Rise in the western equatorial Atlantic (3040 m present water depth; 748.26-850.70 mbsf). Downcore bulk sediment records of biogenic barium, total reactive phosphorus, biogenic silica, and calcium carbonate are interpreted to represent export production, net nutrient burial, biogenic opal production, and inorganic carbon burial, respectively. The global positive excursion in d13C subsequent to the E/O boundary is recorded at Site 925. Export production appears to have been externally forced by orbital parameters at eccentricity frequencies during the study interval, based on spectral analysis of the biogenic barium and reactive phosphorus records. Biogenic silica production or preservation increased after the Eocene/Oligocene boundary to a higher baseline, although overall productivity and nutrient burial did not increase, based on barium and reactive phosphorus records. Thus, although absolute production did not increase at this site, a shift in relative abundance of siliceous versus carbonate productivity may have resulted in a change in relative organic carbon burial. This may have contributed to the positive excursion in global oceanic d13C subsequent to the Eocene/Oligocene boundary, although the silica maximum persists after the carbon isotope excursion ends.

Planktonic foraminiferal preservation in Plio/Pleistocene sediments from the western equatorial Atlantic and Caribbean

Records of mean sortable silt and planktonic foraminiferal preservation from the Ceará Rise (western equatorial Atlantic) and from the Caribbean are presented to analyze the Pliocene (3.5-2.2 Ma) to Pleistocene (1.6-0.3 Ma) evolution of near-bottom current strength and the carbonate corrosiveness of deep water. During the mid-Pleistocene climate transition (~1 Ma) a drastic decrease in glacial bottom current strength and an increase in carbonate corrosiveness is registered, demonstrating a substantial decrease in the glacial contribution of the Lower North Atlantic Deep Water (LNADW) to the Atlantic Ocean. Also, an increased sensitivity to eccentricity orbital forcing is registered after the MPT. By contrast, carbonate preservation increases considerably in the deep Caribbean in response to a strong and persistent stable contribution of Upper North Atlantic Deep Water (UNADW). We found evidence for the strongest and most stable circulation within the LNADW cell during the Northern Hemisphere cooling period between ~3.2 and 2.75 Ma. This is in agreement with the 'superconveyor model' which postulates that the highest NADW production took place prior to ~2.7 Ma. A considerable decrease in bottom current strength and planktonic foraminiferal preservation is observed synchronous with the first occurrence of large-scale continental ice sheets in the Northern Hemisphere. This documents the final termination of the 'superconveyor' at ca. 2.75 Ma. However, our data do not support a 'superconveyor' in the interval between 3.5 and 3.2 Ma when high-amplitude fluctuations in bottom current flow and preservation in planktonic foraminifera are observed. Because of the great sensitivity of NADW production to changes in surface water salinity, we assume that the high-amplitude fluctuations of LNADW circulation prior to ~3.2 Ma are linked to changes in the Atlantic salinity budget. After 2.75 Ma they are primarily controlled by ice-sheet forcing. In contrast to the stepwise deterioration of planktonic foraminiferal preservation in the western deep Atlantic, a trend toward better preservation from the Pliocene to Pleistocene is observed in the deep Caribbean. This indicates a long-term increase in the contribution of UNADW to the Atlantic Ocean.

Tropical Atlantic SST during the mid-Pleistocene transition (MPT) of ODP Hole 175-1077B

We compare a new mid-Pleistocene sea surface temperature (SST) record from the eastern tropical Atlantic to changes in continental ice volume, orbital insolation, Atlantic deepwater ventilation, and Southern Ocean front positions to resolve forcing mechanisms of tropical Atlantic SST during the mid-Pleistocene transition (MPT). At the onset of the MPT, a strong tropical cooling occurred. The change from a obliquity- to a eccentricity-dominated cyclicity in the tropical SST took place at about 650 kyr BP. In orbital cycles, tropical SST changes significantly preceded continental ice-volume changes but were in phase with movements of Southern Ocean fronts. After the onset of large-amplitude 100-kyr variations, additional late glacial warming in the eastern tropical Atlantic was caused by enhanced return flow of warm waters from the western Atlantic driven by strong trade winds. Pronounced 80-kyr variations in tropical SST occurred during the MPT, in phase with and likely directly forced by transitional continental ice-volume variations. During the MPT, a prominent anomalous long-term tropical warming occurred, likely generated by extremely northward displaced Southern Ocean fronts. While the overall pattern of global climate variability during the MPT was determined by changes in mean state and frequency of continental ice volume variations, tropical Atlantic SST variations were primarily driven by early changes in Subantarctic sea-ice extent and coupled Southern Ocean frontal positions.

Parameters of the Earths orbit for the last 5 Million years in 1 kyr resolution

New values for the astronomical parameters of the Earth's orbit and rotation (eccentricity, obliquity and precession) are proposed for paleoclimatic research related to the Late Miocene, the Pliocene and the Quaternary. They have been obtained from a numerical solution of the Lagrangian system of the planetary point masses and from an analytical solution of the Poisson equations of the Earth-Moon system. The analytical expansion developed in this paper allows the direct determination of the main frequencies with their phase and amplitude. Numerical and analytical comparisons with the former astronomical solution BER78 are performed so that the accuracy and the interval of time over which the new solution is valid can be estimated. The corresponding insolation values have also been computed and compared to the former ones. This analysis leads to the conclusion that the new values are expected to be reliable over the last 5 Ma in the time domain and at least over the last 10 Ma in the frequency domain.

Astronomical time scale for Paleocene and early Eocene sediments

The astronomical-tuned time scale is rapidly extended into the Paleogene but, due to the existence of an Eocene gap, different tuning options had to be presented for the Paleocene. These options differ both in number and tuning of ~405-kyr eccentricity related cycles and are only partially consistent with recalculated 40Ar/39Ar constraints for the Cretaceous/Paleogene (K/Pg) and Paleocene/Eocene (P/E) boundaries. In this paper, we evaluate the cyclostratigraphic interpretation of records from ODP Leg 198 and 208 sites, and the Zumaia section to solve the problem of the different tuning options. We found that the interval between the K/Pg boundary and the early Late Paleocene biotic event (ELPE) comprises 17 instead of 16 * ~405-kyr eccentricity related cycles as previously proposed, while the entire Paleocene contains 25 * ~405-kyr cycles. Starting from 40Ar/39Ar age constraints for the K/Pg boundary, a new tuning to 405-kyr eccentricity is presented for the Paleocene and earliest Eocene, which results in ages of ~66.0 and ~ 56.0 Ma for the K/Pg and P/E boundaries, respectively. This tuning introduces considerable differences in age for a number of nannofossil events at ODP Sites 1209 and 1262 in the interval between 61 and 63 Ma, but eliminates large and abrupt changes in the seafloor spreading rate. The tuning seems further consistent with recalculated 40Ar/39Ar ages for ash layer -17 of early Eocene age. However, despite this apparent consistency with existing radio-isotopic constraints, an alternative 405-kyr younger or, less likely, older tuning cannot be excluded at this stage.

Short-term spectral analysis of ODP Hole 114-704B

Short-term spectral analysis was carried out on geochemical logging data from ODP Site 704. The FFT was used to compute the amplitude spectra of short-term overlapping segments to produce depth-period-amplitude spectrograms of the logging data. The spectrograms provided a means of evaluating the significance of the observed periodic components. The periodic components that were consistently present and prominent across a given record interval were considered to be significant. Changes in the spectrogram characteristics seem to reflect changes in either lithology, sedimentation rates, or hiatuses and may therefore provide useful information to aid in stratigraphic and paleoenvironmental studies. The dominant periodicity during the late Pleistocene and Brunhes Chron (0.97 to 0.47 Ma) was determined to be > 100,000 yr whereas the upper Matuyama Chron was dominated by the 41,000-yr periodicity. These periodicities suggest that the sedimentation patterns within the upper Matuyama Chron (0.98-1.78 Ma) were influenced by the Milankovitch obliquity cycle and those within the latest Matuyama-Brunhes Chron (<0.98 Ma) by the eccentricity cycle. The Brunhes/Matuyama boundary therefore represents a major discontinuity. Periodicities observed within the lower Matuyama and the upper Gauss Chron did not correlate with any of the periodicities within the Milankovitch frequency bands.

Age model of sediment core V20-120 (Table 2)

Time series analyses of atmospheric and oceanic variables in a late Pleistocene record from the northwest Pacific show the complex relationship of the response of various segments of the climate system to changes in the earth's orbit. Most variance spectra of time series from this subarctic record contain frequency peaks with periods corresponding to at least one of the major orbital components of eccentricity, obliquity, or precession. Although the radiolarian faunal (water mass) assemblages have prominent spectral peaks with 41,000-year periods which are coherent with obliquity at this frequency, only the Transitional faunal assemblage contains variance focused at a frequency corresponding to the 100,000-year period of eccentricity. Three of these faunal time series also show variance concentrated at a frequency with a 20,000-year period. These three time series are not coherent at a 20,000-year frequency with either of the dominant spectral peaks of precession. They are coherent, however, with variations in the second harmonic of the obliquity cycle. Changes in obliquity apparently affect siliceous faunal abundances in the northwest Pacific region of this high-latitude site more than variations in eccentricity or precession. Maxima in the time series of quartz abundance occur during low values of eccentricity and high glacial ice volume. Because atmospheric winds serve as the major source of supply of quartz to the sediments at this site, these high quartz values reflect increased aridity at the source region. Except for short periods during interglacials, the characteristics of the surface waters in this region of the subarctic Pacific during much of the last 460,000 years were similar to those which exist today in the Sea of Okhotsk. The spectrum of winter sea surface temperature estimates, derived from siliceous microfaunal abundances, contains dominant frequency peaks at periods of 100,000, 41,000, and 23,000 years which are coherent with eccentricity, obliquity, and precession, respectively. Based upon the relationship of the Subarctic Front with the dominance of specific faunal asemblages, the front was positioned south of its present-day location throughout much of the late Pleistocene.

Milankovitch frequency peaks of ODP Site 113-693 (Table 3)

The first well logs collected below the Antarctic circle were obtained during Leg 113 at Site 693 on the Dronning Maud Land Margin (Antarctica) in the Weddell Sea. Gamma-ray, resistivity, and sonic logs were collected between 108.0 and 439.0 mbsf. The downhole logs show good agreement with the data collected from cores and provide a continuous measurement of the sedimentary record. These continuous log records show that the rather uniform Tertiary lithology seen in cores is characterized by high-frequency variability in the log data. Several thin hard streaks are identified, the largest of which coincides with a major Miocene hiatus. Associated with this hiatus is a change to lower illite content (and correspondingly lower gamma-ray counts) and to a significant increase in diatom content. Spectral analysis of the logs was performed on the lower Pliocene through upper Oligocene interval (108.0-343.0 mbsf). Between 108.0 and 245.0 mbsf, average sedimentation rates (50 and 26 m/m.y.) are high enough to show that variance is present in the orbital eccentricity (~95 k.y.) and obliquity (~41 k.y.) bands. Between 253.0 and 343.0 mbsf, the sedimentation rate (8 m/m.y.) is too low to resolve high frequency variations. The Milankovitch frequencies are best developed in the resistivity logs. Resistivity is responding to changes in porosity, which in these sediments is controlled by the abundance of biosiliceous sediments, particularly diatoms. The orbital forcing suggested by the Milankovitch frequencies may be influencing diatom productivity by inducing oscillations in upwelling, ice coverage, pack ice, and/or polynya. Although variations in diatom abundance were observed in the cores, they were not attributed to a Milankovitch signal, and therefore in this environment, downhole logs are an important contribution to the detection and understanding of orbitally influenced changes in sedimentation.

Massic magnetic susceptibility and CaCO3 content at the P. ohmi and T. hugii zones from the Rio Argos X.Ag-1 section, southeastern Spain

An orbital floating time scale of the Hauterivian-Barremian transition (Early Cretaceous) is proposed using high-resolution magnetic susceptibility measurements. Orbital tuning was performed on the Río Argos section (southeast Spain), the candidate for a Global boundary Stratotype Section and Point (GSSP) for the Hauterivian-Barremian transition. Spectral analyses of MS variations, coupled with the frequency ratio method, allow the recognition of precession, obliquity and eccentricity frequency bands. Orbitally-tuned magnetic susceptibility provides minimum durations for ammonite biozones. The durations of well-constrained ammonite zones are assessed at 0.78 myr for Pseudothurmannia ohmi (Late Hauterivian) and 0.57 myr for Taveraidiscus hugii (Early Barremian). These results are consistent with previous estimates from the other reference section (Angles, southeast France) and tend to show that the Río Argos section displays a complete succession for this time interval. They differ significantly from those proposed in the Geologic Time Scale 2008 and may help to improve the next compilation. The Faraoni Oceanic Anoxic Event, a key Early Cretaceous oceanographic perturbation occurring at the P. ohmi/P. catulloi subzone boundary has a duration estimated at 0.10-0.15 myr, which is similar to previous assessments.

Age models and sedimentation rates of sites in the tropical Atlantic

Late Pleistocene signals of calcium carbonate, organic carbon, and opaline silica concentration and accumulation are documented in a series of cores from a zonal/meridional/depth transect in the equatorial Atlantic Ocean to reconstruct the regional sedimentary history. Spectral analysis reveals that maxima and minima in biogenous sedimentation occur with glacial-interglacial cyclicity as a function of both (1) primary production at the sea surface modulated by orbitally forced variation in trade wind zonality and (2) destruction at the seafloor by variation in the chemical character of advected intermediate and deep water from high latitudes modulated by high-latitude ice volume. From these results a pattern emerges in which the relative proportion of signal variance from the productivity signal centered on the precessional (23 kyr) band decreases while that of the destruction signal centered on the obliquity (41 kyr) and eccentricity (100 kyr) periods increases below ~3600-m ocean depth.