Climate variability effects on field crops in Iberian Peninsula. Predictability

La presente Tesis constituye un avance en el conocimiento de los efectos de la variabilidad climática en los cultivos en la Península Ibérica (PI). Es bien conocido que la temperatura del océano, particularmente de la región tropical, es una de las variables más convenientes para ser utilizado como predictor climático. Los océanos son considerados como la principal fuente de almacenamiento de calor del planeta debido a la alta capacidad calorífica del agua. Cuando se libera esta energía, altera los regímenes globales de circulación atmosférica por mecanismos de teleconexión. Estos cambios en la circulación general de la atmósfera afectan a la temperatura, precipitación, humedad, viento, etc., a escala regional, los cuales afectan al crecimiento, desarrollo y rendimiento de los cultivos. Para el caso de Europa, esto implica que la variabilidad atmosférica en una región específica se asocia con la variabilidad de otras regiones adyacentes y/o remotas, como consecuencia Europa está siendo afectada por los patrones de circulaciones globales, que a su vez, se ven afectados por patrones oceánicos. El objetivo general de esta tesis es analizar la variabilidad del rendimiento de los cultivos y su relación con la variabilidad climática y teleconexiones, así como evaluar su predictibilidad. Además, esta Tesis tiene como objetivo establecer una metodología para estudiar la predictibilidad de las anomalías del rendimiento de los cultivos. El análisis se centra en trigo y maíz como referencia para otros cultivos de la PI, cultivos de invierno en secano y cultivos de verano en regadío respectivamente. Experimentos de simulación de cultivos utilizando una metodología en cadena de modelos (clima + cultivos) son diseñados para evaluar los impactos de los patrones de variabilidad climática en el rendimiento y su predictibilidad. La presente Tesis se estructura en dos partes: La primera se centra en el análisis de la variabilidad del clima y la segunda es una aplicación de predicción cuantitativa de cosechas. La primera parte está dividida en 3 capítulos y la segundo en un capitulo cubriendo los objetivos específicos del presente trabajo de investigación. Parte I. Análisis de variabilidad climática El primer capítulo muestra un análisis de la variabilidad del rendimiento potencial en una localidad como indicador bioclimático de las teleconexiones de El Niño con Europa, mostrando su importancia en la mejora de predictibilidad tanto en clima como en agricultura. Además, se presenta la metodología elegida para relacionar el rendimiento con las variables atmosféricas y oceánicas. El rendimiento de los cultivos es parcialmente determinado por la variabilidad climática atmosférica, que a su vez depende de los cambios en la temperatura de la superficie del mar (TSM). El Niño es el principal modo de variabilidad interanual de la TSM, y sus efectos se extienden en todo el mundo. Sin embargo, la predictibilidad de estos impactos es controversial, especialmente aquellos asociados con la variabilidad climática Europea, que se ha encontrado que es no estacionaria y no lineal. Este estudio mostró cómo el rendimiento potencial de los cultivos obtenidos a partir de datos de reanálisis y modelos de cultivos sirve como un índice alternativo y más eficaz de las teleconexiones de El Niño, ya que integra las no linealidades entre las variables climáticas en una única serie temporal. Las relaciones entre El Niño y las anomalías de rendimiento de los cultivos son más significativas que las contribuciones individuales de cada una de las variables atmosféricas utilizadas como entrada en el modelo de cultivo. Además, la no estacionariedad entre El Niño y la variabilidad climática europea se detectan con mayor claridad cuando se analiza la variabilidad de los rendimiento de los cultivos. La comprensión de esta relación permite una cierta predictibilidad hasta un año antes de la cosecha del cultivo. Esta predictibilidad no es constante, sino que depende tanto la modulación de la alta y baja frecuencia. En el segundo capítulo se identifica los patrones oceánicos y atmosféricos de variabilidad climática que afectan a los cultivos de verano en la PI. Además, se presentan hipótesis acerca del mecanismo eco-fisiológico a través del cual el cultivo responde. Este estudio se centra en el análisis de la variabilidad del rendimiento de maíz en la PI para todo el siglo veinte, usando un modelo de cultivo calibrado en 5 localidades españolas y datos climáticos de reanálisis para obtener series temporales largas de rendimiento potencial. Este estudio evalúa el uso de datos de reanálisis para obtener series de rendimiento de cultivos que dependen solo del clima, y utilizar estos rendimientos para analizar la influencia de los patrones oceánicos y atmosféricos. Los resultados muestran una gran fiabilidad de los datos de reanálisis. La distribución espacial asociada a la primera componente principal de la variabilidad del rendimiento muestra un comportamiento similar en todos los lugares estudiados de la PI. Se observa una alta correlación lineal entre el índice de El Niño y el rendimiento, pero no es estacionaria en el tiempo. Sin embargo, la relación entre la temperatura del aire y el rendimiento se mantiene constante a lo largo del tiempo, siendo los meses de mayor influencia durante el período de llenado del grano. En cuanto a los patrones atmosféricos, el patrón Escandinavia presentó una influencia significativa en el rendimiento en PI. En el tercer capítulo se identifica los patrones oceánicos y atmosféricos de variabilidad climática que afectan a los cultivos de invierno en la PI. Además, se presentan hipótesis acerca del mecanismo eco-fisiológico a través del cual el cultivo responde. Este estudio se centra en el análisis de la variabilidad del rendimiento de trigo en secano del Noreste (NE) de la PI. La variabilidad climática es el principal motor de los cambios en el crecimiento, desarrollo y rendimiento de los cultivos, especialmente en los sistemas de producción en secano. En la PI, los rendimientos de trigo son fuertemente dependientes de la cantidad de precipitación estacional y la distribución temporal de las mismas durante el periodo de crecimiento del cultivo. La principal fuente de variabilidad interanual de la precipitación en la PI es la Oscilación del Atlántico Norte (NAO), que se ha relacionado, en parte, con los cambios en la temperatura de la superficie del mar en el Pacífico Tropical (El Niño) y el Atlántico Tropical (TNA). La existencia de cierta predictibilidad nos ha animado a analizar la posible predicción de los rendimientos de trigo en la PI utilizando anomalías de TSM como predictor. Para ello, se ha utilizado un modelo de cultivo (calibrado en dos localidades del NE de la PI) y datos climáticos de reanálisis para obtener series temporales largas de rendimiento de trigo alcanzable y relacionar su variabilidad con anomalías de la TSM. Los resultados muestran que El Niño y la TNA influyen en el desarrollo y rendimiento del trigo en el NE de la PI, y estos impactos depende del estado concurrente de la NAO. Aunque la relación cultivo-TSM no es igual durante todo el periodo analizado, se puede explicar por un mecanismo eco-fisiológico estacionario. Durante la segunda mitad del siglo veinte, el calentamiento (enfriamiento) en la superficie del Atlántico tropical se asocia a una fase negativa (positiva) de la NAO, que ejerce una influencia positiva (negativa) en la temperatura mínima y precipitación durante el invierno y, por lo tanto, aumenta (disminuye) el rendimiento de trigo en la PI. En relación con El Niño, la correlación más alta se observó en el período 1981 -2001. En estas décadas, los altos (bajos) rendimientos se asocian con una transición El Niño - La Niña (La Niña - El Niño) o con eventos de El Niño (La Niña) que están finalizando. Para estos eventos, el patrón atmosférica asociada se asemeja a la NAO, que también influye directamente en la temperatura máxima y precipitación experimentadas por el cultivo durante la floración y llenado de grano. Los co- efectos de los dos patrones de teleconexión oceánicos ayudan a aumentar (disminuir) la precipitación y a disminuir (aumentar) la temperatura máxima en PI, por lo tanto el rendimiento de trigo aumenta (disminuye). Parte II. Predicción de cultivos. En el último capítulo se analiza los beneficios potenciales del uso de predicciones climáticas estacionales (por ejemplo de precipitación) en las predicciones de rendimientos de trigo y maíz, y explora métodos para aplicar dichos pronósticos climáticos en modelos de cultivo. Las predicciones climáticas estacionales tienen un gran potencial en las predicciones de cultivos, contribuyendo de esta manera a una mayor eficiencia de la gestión agrícola, seguridad alimentaria y de subsistencia. Los pronósticos climáticos se expresan en diferentes formas, sin embargo todos ellos son probabilísticos. Para ello, se evalúan y aplican dos métodos para desagregar las predicciones climáticas estacionales en datos diarios: 1) un generador climático estocástico condicionado (predictWTD) y 2) un simple re-muestreador basado en las probabilidades del pronóstico (FResampler1). Los dos métodos se evaluaron en un caso de estudio en el que se analizaron los impactos de tres escenarios de predicciones de precipitación estacional (predicción seco, medio y lluvioso) en el rendimiento de trigo en secano, sobre las necesidades de riego y rendimiento de maíz en la PI. Además, se estimó el margen bruto y los riesgos de la producción asociada con las predicciones de precipitación estacional extremas (seca y lluviosa). Los métodos predWTD y FResampler1 usados para desagregar los pronósticos de precipitación estacional en datos diarios, que serán usados como inputs en los modelos de cultivos, proporcionan una predicción comparable. Por lo tanto, ambos métodos parecen opciones factibles/viables para la vinculación de los pronósticos estacionales con modelos de simulación de cultivos para establecer predicciones de rendimiento o las necesidades de riego en el caso de maíz. El análisis del impacto en el margen bruto de los precios del grano de los dos cultivos (trigo y maíz) y el coste de riego (maíz) sugieren que la combinación de los precios de mercado previstos y la predicción climática estacional pueden ser una buena herramienta en la toma de decisiones de los agricultores, especialmente en predicciones secas y/o localidades con baja precipitación anual. Estos métodos permiten cuantificar los beneficios y riesgos de los agricultores ante una predicción climática estacional en la PI. Por lo tanto, seríamos capaces de establecer sistemas de alerta temprana y diseñar estrategias de adaptación del manejo del cultivo para aprovechar las condiciones favorables o reducir los efectos de condiciones adversas. La utilidad potencial de esta Tesis es la aplicación de las relaciones encontradas para predicción de cosechas de la próxima campaña agrícola. Una correcta predicción de los rendimientos podría ayudar a los agricultores a planear con antelación sus prácticas agronómicas y todos los demás aspectos relacionados con el manejo de los cultivos. Esta metodología se puede utilizar también para la predicción de las tendencias futuras de la variabilidad del rendimiento en la PI. Tanto los sectores públicos (mejora de la planificación agrícola) como privados (agricultores, compañías de seguros agrarios) pueden beneficiarse de esta mejora en la predicción de cosechas. ABSTRACT The present thesis constitutes a step forward in advancing of knowledge of the effects of climate variability on crops in the Iberian Peninsula (IP). It is well known that ocean temperature, particularly the tropical ocean, is one of the most convenient variables to be used as climate predictor. Oceans are considered as the principal heat storage of the planet due to the high heat capacity of water. When this energy is released, it alters the global atmospheric circulation regimes by teleconnection1 mechanisms. These changes in the general circulation of the atmosphere affect the regional temperature, precipitation, moisture, wind, etc., and those influence crop growth, development and yield. For the case of Europe, this implies that the atmospheric variability in a specific region is associated with the variability of others adjacent and/or remote regions as a consequence of Europe being affected by global circulations patterns which, in turn, are affected by oceanic patterns. The general objective of this Thesis is to analyze the variability of crop yields at climate time scales and its relation to the climate variability and teleconnections, as well as to evaluate their predictability. Moreover, this Thesis aims to establish a methodology to study the predictability of crop yield anomalies. The analysis focuses on wheat and maize as a reference crops for other field crops in the IP, for winter rainfed crops and summer irrigated crops respectively. Crop simulation experiments using a model chain methodology (climate + crop) are designed to evaluate the impacts of climate variability patterns on yield and its predictability. The present Thesis is structured in two parts. The first part is focused on the climate variability analyses, and the second part is an application of the quantitative crop forecasting for years that fulfill specific conditions identified in the first part. This Thesis is divided into 4 chapters, covering the specific objectives of the present research work. Part I. Climate variability analyses The first chapter shows an analysis of potential yield variability in one location, as a bioclimatic indicator of the El Niño teleconnections with Europe, putting forward its importance for improving predictability in both climate and agriculture. It also presents the chosen methodology to relate yield with atmospheric and oceanic variables. Crop yield is partially determined by atmospheric climate variability, which in turn depends on changes in the sea surface temperature (SST). El Niño is the leading mode of SST interannual variability, and its impacts extend worldwide. Nevertheless, the predictability of these impacts is controversial, especially those associated with European climate variability, which have been found to be non-stationary and non-linear. The study showed how potential2 crop yield obtained from reanalysis data and crop models serves as an alternative and more effective index of El Niño teleconnections because it integrates the nonlinearities between the climate variables in a unique time series. The relationships between El Niño and crop yield anomalies are more significant than the individual contributions of each of the atmospheric variables used as input in the crop model. Additionally, the non-stationarities between El Niño and European climate variability are more clearly detected when analyzing crop-yield variability. The understanding of this relationship allows for some predictability up to one year before the crop is harvested. This predictability is not constant, but depends on both high and low frequency modulation. The second chapter identifies the oceanic and atmospheric patterns of climate variability affecting summer cropping systems in the IP. Moreover, hypotheses about the eco-physiological mechanism behind crop response are presented. It is focused on an analysis of maize yield variability in IP for the whole twenty century, using a calibrated crop model at five contrasting Spanish locations and reanalyses climate datasets to obtain long time series of potential yield. The study tests the use of reanalysis data for obtaining only climate dependent time series of simulated crop yield for the whole region, and to use these yield to analyze the influences of oceanic and atmospheric patterns. The results show a good reliability of reanalysis data. The spatial distribution of the leading principal component of yield variability shows a similar behaviour over all the studied locations in the IP. The strong linear correlation between El Niño index and yield is remarkable, being this relation non-stationary on time, although the air temperature-yield relationship remains on time, being the highest influences during grain filling period. Regarding atmospheric patterns, the summer Scandinavian pattern has significant influence on yield in IP. The third chapter identifies the oceanic and atmospheric patterns of climate variability affecting winter cropping systems in the IP. Also, hypotheses about the eco-physiological mechanism behind crop response are presented. It is focused on an analysis of rainfed wheat yield variability in IP. Climate variability is the main driver of changes in crop growth, development and yield, especially for rainfed production systems. In IP, wheat yields are strongly dependent on seasonal rainfall amount and temporal distribution of rainfall during the growing season. The major source of precipitation interannual variability in IP is the North Atlantic Oscillation (NAO) which has been related in part with changes in the Tropical Pacific (El Niño) and Atlantic (TNA) sea surface temperature (SST). The existence of some predictability has encouraged us to analyze the possible predictability of the wheat yield in the IP using SSTs anomalies as predictor. For this purpose, a crop model with a site specific calibration for the Northeast of IP and reanalysis climate datasets have been used to obtain long time series of attainable wheat yield and relate their variability with SST anomalies. The results show that El Niño and TNA influence rainfed wheat development and yield in IP and these impacts depend on the concurrent state of the NAO. Although crop-SST relationships do not equally hold on during the whole analyzed period, they can be explained by an understood and stationary ecophysiological mechanism. During the second half of the twenty century, the positive (negative) TNA index is associated to a negative (positive) phase of NAO, which exerts a positive (negative) influence on minimum temperatures (Tmin) and precipitation (Prec) during winter and, thus, yield increases (decreases) in IP. In relation to El Niño, the highest correlation takes place in the period 1981-2001. For these decades, high (low) yields are associated with an El Niño to La Niña (La Niña to El Niño) transitions or to El Niño events finishing. For these events, the regional associated atmospheric pattern resembles the NAO, which also influences directly on the maximum temperatures (Tmax) and precipitation experienced by the crop during flowering and grain filling. The co-effects of the two teleconnection patterns help to increase (decrease) the rainfall and decrease (increase) Tmax in IP, thus on increase (decrease) wheat yield. Part II. Crop forecasting The last chapter analyses the potential benefits for wheat and maize yields prediction from using seasonal climate forecasts (precipitation), and explores methods to apply such a climate forecast to crop models. Seasonal climate prediction has significant potential to contribute to the efficiency of agricultural management, and to food and livelihood security. Climate forecasts come in different forms, but probabilistic. For this purpose, two methods were evaluated and applied for disaggregating seasonal climate forecast into daily weather realizations: 1) a conditioned stochastic weather generator (predictWTD) and 2) a simple forecast probability resampler (FResampler1). The two methods were evaluated in a case study where the impacts of three scenarios of seasonal rainfall forecasts on rainfed wheat yield, on irrigation requirements and yields of maize in IP were analyzed. In addition, we estimated the economic margins and production risks associated with extreme scenarios of seasonal rainfall forecasts (dry and wet). The predWTD and FResampler1 methods used for disaggregating seasonal rainfall forecast into daily data needed by the crop simulation models provided comparable predictability. Therefore both methods seem feasible options for linking seasonal forecasts with crop simulation models for establishing yield forecasts or irrigation water requirements. The analysis of the impact on gross margin of grain prices for both crops and maize irrigation costs suggests the combination of market prices expected and the seasonal climate forecast can be a good tool in farmer’s decision-making, especially on dry forecast and/or in locations with low annual precipitation. These methodologies would allow quantifying the benefits and risks of a seasonal weather forecast to farmers in IP. Therefore, we would be able to establish early warning systems and to design crop management adaptation strategies that take advantage of favorable conditions or reduce the effect of adverse conditions. The potential usefulness of this Thesis is to apply the relationships found to crop forecasting on the next cropping season, suggesting opportunity time windows for the prediction. The methodology can be used as well for the prediction of future trends of IP yield variability. Both public (improvement of agricultural planning) and private (decision support to farmers, insurance companies) sectors may benefit from such an improvement of crop forecasting.

Trans-Pacific migrations of the loggerhead turtle (Caretta caretta) demonstrated with mitochondrial DNA markers.

Juvenile loggerhead turtles (Caretta caretta) have recently been documented in the vicinity of Baja California and thousands of these animals have been captured in oceanic fisheries of the North Pacific. The presence of loggerhead turtles in the central and eastern North Pacific is a prominent enigma in marine turtle distribution because the nearest documented nesting concentrations for this species are in Australia and Japan, over 10,000 km from Baja California. To determine the origin of the Baja California feeding aggregate and North Pacific fishery mortalities, samples from nesting areas and pelagic feeding aggregates were compared with genetic markers derived from mtDNA control region sequences. Overall, 57 of 60 pelagic samples (95%) match haplotypes seen only in Japanese nesting areas, implicating Japan as the primary source of turtles in the North Pacific Current and around Baja California. Australian nesting colonies may contribute the remaining 5% of these pelagic feeding aggregates. Juvenile loggerhead turtles apparently traverse the entire Pacific Ocean, approximately one-third of the planet, in the course of developmental migrations, but mortality in high-seas fisheries raises concern over the future of this migratory population.

Variability of the ocean heat content during the last millennium - an assessment with the ECHO-g Model

Studies addressing climate variability during the last millennium generally focus on variables with a direct influence on climate variability, like the fast thermal response to varying radiative forcing, or the large-scale changes in atmospheric dynamics (e. g. North Atlantic Oscillation). The ocean responds to these variations by slowly integrating in depth the upper heat flux changes, thus producing a delayed influence on ocean heat content (OHC) that can later impact low frequency SST (sea surface temperature) variability through reemergence processes. In this study, both the externally and internally driven variations of the OHC during the last millennium are investigated using a set of fully coupled simulations with the ECHO-G (coupled climate model ECHAMA4 and ocean model HOPE-G) atmosphere-ocean general circulation model (AOGCM). When compared to observations for the last 55 yr, the model tends to overestimate the global trends and underestimate the decadal OHC variability. Extending the analysis back to the last one thousand years, the main impact of the radiative forcing is an OHC increase at high latitudes, explained to some extent by a reduction in cloud cover and the subsequent increase of short-wave radiation at the surface. This OHC response is dominated by the effect of volcanism in the preindustrial era, and by the fast increase of GHGs during the last 150 yr. Likewise, salient impacts from internal climate variability are observed at regional scales. For instance, upper temperature in the equatorial Pacific is controlled by ENSO (El Nino Southern Oscillation) variability from interannual to multidecadal timescales. Also, both the Pacific Decadal Oscillation (PDO) and the Atlantic Multidecadal Oscillation (AMO) modulate intermittently the interdecadal OHC variability in the North Pacific and Mid Atlantic, respectively. The NAO, through its influence on North Atlantic surface heat fluxes and convection, also plays an important role on the OHC at multiple timescales, leading first to a cooling in the Labrador and Irminger seas, and later on to a North Atlantic warming, associated with a delayed impact on the AMO.

Stable isotopes, Mg/Ca-ratios, and biogenic opal of five sediment cores from the subarctic northwest Pacific

Based on models and proxy data it has been proposed that salinity-driven stratification weakened in the subarctic North Pacific during the last deglaciation, which potentially contributed to the deglacial rise in atmospheric carbon dioxide. We present high-resolution subsurface temperature (TMg/Ca) and subsurface salinity-approximating (d18Oivc-sw) records across the last 20,000 years from the subarctic North Pacific and its marginal seas, derived from combined stable oxygen isotopes and Mg/Ca ratios of the planktonic foraminiferal species Neogloboquadrina pachyderma (sin.). Our results indicate regionally differing changes of subsurface conditions. During the Heinrich Stadial 1 and the Younger Dryas cold phases our sites were subject to reduced thermal stratification, brine rejection due to sea-ice formation, and increased advection of low-salinity water from the Alaskan Stream. In contrast, the Bølling-Allerød warm phase was characterized by strengthened thermal stratification, stronger sea-ice melting, and influence of surface waters that were less diluted by the Alaskan Stream. From direct comparison with alkenone-based sea surface temperature estimates (SSTUk'37), we suggest deglacial thermocline changes that were closely related to changes in seasonal contrasts and stratification of the mixed layer. The modern upper-ocean conditions seem to have developed only since the early Holocene.

Benthic foraminiferal assemblages and paleoproductivity reconstructions for the last 550 kyr of the Ontong Java Plateau, western equatorial Pacific

We analyzed foraminiferal and nannofossil assemblages and stable isotopes in samples from ODP Hole 807A on the Ontong Java Plateau in order to evaluate productivity and carbonate dissolution cycles over the last 550 kyr (kilo year) in the western equatorial Pacific. Our results indicate that productivity was generally higher in glacials than during interglacials, and gradually increased since MIS 13. Carbonate dissolution was weak in deglacial intervals, but often reached a maximum during interglacial to glacial transitions. Carbonate cycles in the western equatorial Pacific were mainly influenced by changes of deep-water properties rather than by local primary productivity. Fluctuations of the estimated thermocline depth were not related to glacial to interglacial alternations, but changed distinctly at ~280 kyr. Before that time the thermocline was relatively shallow and its depth fluctuated at a comparatively high amplitude and low frequency. After 280 kyr, the thermocline was deeper, and its fluctuations were at lower amplitude and higher frequency. These different patterns in productivity and thermocline variability suggest that thermocline dynamics probably were not a controlling factor of biological productivity in the western equatorial Pacific Ocean. In this region, upwelling, the influx of cool, nutrient-rich waters from the eastern equatorial Pacific or of fresh waters from rivers have probably never been important, and their influence on productivity has been negligible over the studied period. Variations in the inferred productivity in general are well correlated with fluctuations in the eolian flux as recorded in the northwestern Pacific, a proxy for the late Quaternary history of the central East Asian dust flux into the Pacific. Therefore, we suggest that the dust flux from the central East Asian continent may have been an important driver of productivity in the western Pacific.

Neogene radiolarian datum levels in the equatorial Indian and Pacific Oceans

Fifty radiolarian events of early Pleistocene and Neogene age were identified in an E-W transect of equatorial DSDP sites, extending from the Gulf of Panama to the western Pacific and eastern Indian Oceans. Our objective was to document the degree of synchroneity or time-transgressiveness of stratigraphically-useful datum levels from this geologic time interval. We restricted our study to low latitudes within which morphological variations of individual taxa are minimal, the total assemblage diversity remains high, and stratigraphic continuity is well-documented by an independent set of criteria. Each of the five sites chosen (503, 573, 289/586, 214) was calibrated to an "absolute" time scale, using a multiple of planktonic foraminiferal, nannofossil, and diatom datum levels which have been independently correlated to the paleomagnetic polarity time scale in piston core material. With these correlations we have assigned "absolute" ages to each radiolarian event, with a precision of 0.1-0.2 m.y. and an accuracy of 0.2-0.4 m.y. On this basis we have classified each of the events as either: (a) synchronous (range of ages <0.4 m.y.); (b) time-transgressive (i.e., range of ages >1.0 m.y.); and (c) not resolvable (range of ages 0.4-1.0 m.y.). Our results show that, among the synchronous datum levels, a large majority (15 out of 19) are last occurrences. Among those events which are clearly time-transgressive, most are first appearances (10 out of 13). In many instances taxa appear to evolve first in the Indian Ocean, and subsequently in the western and eastern Pacific Ocean. This pattern is particularly unexpected in view of the strong east-to-west zonal flow in equatorial latitudes. Three of the time-transgressive events have been used to define zonal boundaries: the first appearances of Spongaster pentas, Diartus hughesi, and D. petterssoni. Our results suggest that biostratigraphic non-synchroneity may be substantial (i.e., greater than 1 m.y.) within a given latitudinal zone; one would expect this effect to be even more pronounced across oceanographic and climatic gradients. We anticipate that the extent of diachroneity may be comparable for diatom, foraminiferal, and nannofossil datum levels as well. If this proves true, global "time scales" may need to be re-formulated on the basis of a smaller number of demonstrably synchronous events.

Strontium and boron geochemistry of ODP sites at Hydrate Ridge, eastern Pacific Ocean

ODP Leg 204, which drilled at Hydrate Ridge, provides unique insights into the fluid regime of an accretionary complex and delineates specific sub-seafloor pathways for fluid transport. Compaction and dewatering due to smectite-illite transition increase with distance from the toe of the accretionary prism and bring up fluids from deep within the accretionary complex to sampled depths (<= 600 mbsf). These fluids have a distinctly non-radiogenic strontium isotope signature indicating reaction with the oceanic basement. Boron isotopes are also consistent with a deep fluid source that has been modified by desorption of heavy boron as clay minerals change from smectite to illite. One of three major horizons serves as conduit for the transport of mainly fluid. Our results enable us to evaluate fluid migration pathways that play important roles on massive gas hydrate accumulations and seepage of methane-rich fluids on southern Hydrate Ridge.

Trace element abundance, and Sr and Nd isotope ratios of dust samples in the Pacific Ocean (Table 2)

Eolian dust preserved in deep-sea sediment cores provides a valuable indicator of past atmospheric circulation and continental paleoclimate. In order to identify the provenance of eolian dust, Nd and Sr isotopic compositions and Rb, Sr and rare earth element (REE) concentrations have been determined for the silicate fractions of deep-sea sediments from the north and central Pacific Ocean. Different regions of the Pacific Ocean are characterized by distinct air-borne inputs, producing a large range in epsolin-Nd (-10 to +1), 87Sr/86Sr (0.705-0.721), La/Yb (5-15), EuN/EuN* (0.6-1.0) and Sr/Nd (4-33). The average Nd isotopic composition of Pacific deep-sea sediments (epsilon-Nd = -6), is more radiogenic than the average from the Atlantic (epsilon-Nd = -8). In contrast, the average147Sm/144Nd ratio for Pacific sediments (0.114) is identical to that of Atlantic sediments and to that of global average riverine suspended material. The values of epsilon-Nd and147Sm/144Nd are positively correlated for the Pacific samples but negatively correlated for Atlantic samples, reflecting a fundamental difference between the dominant components in the end members with radiogenic Nd (island-arc components in the Pacific and LREE-enriched intraplate ocean island components in the Atlantic). Samples from the north central Pacific have distinctive unradiogenic epsilon-Nd values of -10, 87Sr/86Sr > 0.715, high La/Yb (> 12), and low EuN/EuN* (0.6) and Sr/Nd (3-6). These data are virtually identical to the values for loess from Asia and endorse the use of these sediments as indicators of Asian paleoclimate and paleowind directions. Island-arc contributions appear to dominate in the northwest Pacific, resulting in higher epsilon Nd (-1 to +1) and lower 87Sr/86Sr (~ 0.705) and La/Yb (~ 5). Sediments from the eastern Pacific tend to have intermediate Sr and Nd isotopic compositions but regionally variable Sr/Nd and REE patterns; they appear to be derived from the west margin of the North and South American continents, rather than from Asia. Our results confirm that dust provenance can be constrained by isotopic and geochemical analyses, which will facilitate reconstructions of past atmospheric circulation and continental paleoclimate.

Nd isotope ratios of surface sediments from the Pacific Ocean (Table 2)

The neodymium isotopic composition of the silicate fraction of Holocene pelagic sediments from the North Pacific define two provinces: a central North Pacific province characterized by unradiogenic and remarkably homogeneous end (-10.2 +/- 0.5) and a narrow circum-Pacific marginal province characterized by more radiogenic and variable end (-4.2 +/- 3.8). The silicate fraction in the central North Pacific is exclusively eolian; based on prevailing wind patterns, meteorological data, and neodymium isotopic data, the only significant sediment source is Chinese loess. Leaching experiments on Chinese loess confirm that leachable Nd is isotopically indistinguishable from bulk and residual silicate Nd. Silicates in the circum-North Pacific marginal province comprise eolian loess, volcanic ash, and hemipelagic sediments derived from volcanic arcs. A compilation of Pacific seawater and Mn nodule epsilon-Nd data shows no clear spatial variation except for a general decrease from surface to deep waters from -3 to -4 and slightly lower epsilon-Nd in bottom waters along the western North Pacific due to the incursion of Antarctic Bottom Water. The relative homogeneity of bottom water epsilon-Nd, which contrasts sharply with the distinctive variation in sediment epsilon-Nd, plus the large difference between the average end of bottom waters and the central North Pacific eolian silicates (-4 vs. -10), suggests that any contribution of REE to seawater from eolian materials is insignificant. Furthermore, leaching of REE from eolian particles as they sink though the water column must be insignificant because Nd in shallow waters is more radiogenic than Nd in deeper waters. That there is no contrast in the Nd isotopic composition of bottom waters that overlie the central and marginal sediment provinces suggests that the ash and hemipelagic sediments derived from Pacific rim volcanic arcs also contribute minimal REE to seawater. The elimination of eolian, ash, and hemipelagic sediments leaves only near-shore riverine particulates as a possibly significant particulate source of REE to seawater.