(Table 1) Movement parameters of nine adult female snowy owls (Bubo scandiacus) tracked during the winter period 07/08 and 08/09 around Hudson Strait, Canada

Top predators of the arctic tundra are facing a long period of very low prey availability during winter and subsidies from other ecosystems such as the marine environment may help to support their populations. Satellite tracking of snowy owls, a top predator of the tundra, revealed that most adult females breeding in the Canadian Arctic overwinter at high latitudes in the eastern Arctic and spend several weeks (up to 101 d) on the sea-ice between December and April. Analysis of high-resolution satellite images of sea-ice indicated that owls were primarily gathering around open water patches in the ice, which are commonly used by wintering seabirds, a potential prey. Such extensive use of sea-ice by a tundra predator considered a small mammal specialist was unexpected, and suggests that marine resources subsidize snowy owl populations in winter. As sea-ice regimes in winter are expected to change over the next decades due to climate warming, this may affect the wintering strategy of this top predator and ultimately the functioning of the tundra ecosystem.

EOT11a - Global Empirical Ocean Tide model from multi-mission satellite altimetry, with links to model results

EOT11a is a global (E)mpirical (O)cean (T)ide model derived in 2011 by residual analysis of multi-mission satellite (a)ltimeter data. EOT11a includes amplitudes and phases of the main astronomical tides M2, S2, N2, K2, 2N2, O1, K1, P2, and Q1, the non-linear constituent M4, the long period tides Mm and Mf, and the radiational tide S1. Ocean tides as well as loading tides are provided. EOT11a was computed by means of residual tidal analysis of multi-mission altimeter data from TOPEX/Poseidon, ERS-2, ENVISAT, and Jason-1/2, as far as acquired between September 1992 and April 2010. The resolution of 7.5'x7.5' is identical with FES2004 which was used as reference model for the residual tide analysis. The development of EOT11a was funded by the Deutsche Forschungsgemeinschaft (DFG) under grant BO1228/6-2.

Benthic and planktic foraminiferal stable isotope records combined with Nd isotope ratios of ODP Hole 198-1210B

The Pacific Ocean is the largest water body on Earth, and circulation in the Pacific contributed significantly to climate evolution in the latest Cretaceous, the culmination of a period of long-term cooling. Here, we present new high-resolution late Campanian to Maastrichtian benthic and planktic foraminiferal stable isotope data and a neodymium (Nd) isotope record obtained from sedimentary ferromanganese oxide coatings of Ocean Drilling Program Hole 1210B from the tropical Pacific Ocean (Shatsky Rise). These new records resolve 13 million years in the latest Cretaceous, providing insights into changes in surface and bottom water temperatures and source regions of deep to intermediate waters covering the carbon isotope excursions of the Campanian-Maastrichtian Boundary Event (CMBE) and the Mid-Maastrichtian event (MME). Our new benthic foraminiferal d18O and Nd isotope records together with published Nd isotope data show markedly parallel trends across the studied interval over a broad range of bathyal to abyssal water depths interpreted to reflect changes in the intensity of deep-ocean circulation in the tropical Pacific. In particular, we observe a three-million-year-long period of cooler conditions in the early Maastrichtian (72.5 to 69.5 Ma) when a concomitant change toward less radiogenic seawater Nd isotope signatures probably marks a period of enhanced admixture and northward flow of deep waters with Southern Ocean provenance. We suggest this change to have been triggered by intensified formation and convection of deep waters in the high southern latitudes, a process that weakened during the MME (69.5 to 68.5 Ma). The early Maastrichtian cold interval is closely related to the negative and positive carbon isotope trends of the CMBE and MME. The millions-of-years long duration of these carbon cycle perturbations suggests a tectonic forcing of climatic cooling, possibly related to changes in ocean basin geometry and bathymetry.

Stable isotope record of benthic and planktonic foraminifera from ODP Site 175-1087 in the southern Cape Basin, Atlantic Ocean

High-resolution planktonic and benthic stable isotope records from Ocean Drilling Program Site 1087 off southeast Africa provide the basis for a detailed study of glacial-interglacial (G-IG) cycles during the last 500 k.y. This site is located in the Southern Cape Basin at the boundary of the coastal upwelling of Benguela and close to the gateway between the South Atlantic and the Indian Oceans. It therefore monitors variations of the hydrological fronts associated with the upwelling system and the Atlantic-Indian Ocean interconnections, in relation to global climate change. The coldest period of the last 500 k.y. corresponds to marine isotope Stage (MIS) 12, when surface water temperature was 4°C lower than during the last glacial maximum (LGM) as recorded by the surface-dwelling foraminifer Globigerinoides ruber. The warmest periods occurred during MISs 5 and 11, a situation slightly different to that observed at Site 704, which is close to the Polar Front Zone, where there is no significant difference between the interglacial stages for the past 450 k.y., except the long period of warmth during MIS 11. The planktonic and benthic carbon isotope records do not follow the G-IG cycles but show large oscillations related to major changes in the productivity regime. The largest positive 13C excursion between 260 and 425 ka coincides with the global mid-Brunhes event of carbonate productivity. The oxygen and carbon isotopic gradients between surface and deep waters display long-term changes superimposed on rapid and high-frequency fluctuations that do not follow the regular G-IG pattern; these gradients indicate modifications of the temperature, salinity, and productivity gradients due to changes in the thermocline depth, the position of the hydrological fronts, and the strength of the Benguela Current.

Neogene magneto- and biostratigraphy of ODP Leg 104 holes

Silicoflagellate assemblages of ODP Leg 104 Neogene sequences are the basis of an interpretation of changes in the Neogene paleoenvironment of the Norwegian Sea. Fluctuations in the percentages of temperature and nutrient-sensitive taxonomic groups document major changes in sea-surface conditions. A brief, but distinct, cooling event occurred at 18.0-17.5 Ma which resulted in the disappearance of Naviculopsis. Following this early Miocene cooling a long period of increasing surface-water temperature occurred, leading up to a thermal high in the early middle Miocene (14.0 Ma). The early late Miocene (10.0-9.0 Ma) was distinctly cooler than the middle Miocene, but warmer than the remainder of the Neogene. Conditions between 13.0 and 10.0 Ma are unrecorded because of a regional hiatus, which is the earliest evidence for an end to the more temperate and stable conditions of the early to middle middle Miocene. A major plunge in temperatures occurred between 8.5 and 7.4 Ma and during the remainder of the late Miocene and Pliocene; from 7.4 to 2.65 Ma subpolar conditions prevailed. Silicoflagellates disappeared, except for sporadic occurrences, at 2.64 Ma with the beginning of dominant glacial sedimentation. Biogenic opal is absent in sediments younger than 0.76 Ma, indicating the dominance of glacial conditions with extensive sea ice.

Age model and Plio-Pleistocene benthic stable oxygen isotope ratios of ODP Site 184-1143 in the Sourh China Sea

Based on benthic foraminiferal delta18O from ODP Site 1143, a 5-Myr astronomical timescale for the West Pacific Plio-Pleistocene was established using an automatic orbital tuning method. The tuned Brunhes/Matuyama paleomagnetic polarity reversal age agrees well with the previously published age of 0.78 Ma. The tuned ages for several planktonic foraminifer bio-events also agree well with published dates, and new ages for some other bio-events in the South China Sea were also estimated. The benthic delta18O from Site 1143 is highly coherent with the Earth's orbit (ETP) both at the obliquity and precession bands for the last 5 Myr, and at the eccentricity band for the last 2 Myr. In general, the 41-kyr cycle was dominant through the Plio-Pleistocene although the 23-kyr cycle was also very strong. The 100-kyr cycle became dominant only during the last 1 Myr. A comparison of the benthic delta18O between the Atlantic (ODP 659) and the East and West Pacific (846 and 1143) reveals that the Atlantic-Pacific benthic oxygen isotope difference ratio (Delta delta18OAtl-Pac) displays an increasing trend in three time intervals: 3.6-2.7 Ma, 2.7-2.1 Ma and 1.5-0.25 Ma. Each of the intervals begins with a rapid negative shift in Delta delta18OAtl-Pac, followed by a long period with an increasing trend, corresponding to the growth of the Northern Hemisphere ice sheet. This means that all three intervals of ice sheet growth in the Northern Hemisphere were accompanied at the beginning by a rapid relative warming of deep water in the Atlantic as compared to that of the Pacific, followed by its gradual relative cooling. This general trend, superimposed on the frequent fluctuations with glacial cycles, should yield insights into the processes leading to the boreal glaciation. Cross-spectral analyses of the Delta delta18OAtl-Pac with the Earth's orbit suggests that after the initiation of Northern Hemisphere glaciation at about 2.5 Ma, obliquity rather than precession had become the dominant force controlling the vertical structure or thermohaline circulation in the paleo-ocean.

Phenotypic evolution in a fossil gastropod species lineage: evidence for adaptive radiation?

Detecting speciation in the fossil record is a particularly challenging matter. Palaeontologists are usually confronted with poor preservation and limited knowledge on the palaeoenvironment. Even in the contrary case of adequate preservation and information, the linkage of pattern to process is often obscured by insufficient temporal resolution. Consequently, reliable documentations of speciation in fossils with discussions on underlying mechanisms are rare. Here we present a well-resolved pattern of morphological evolution in a fossil species lineage of the gastropod Melanopsis in the long-lived Lake Pannon. These developments are related to environmental changes, documented by isotope and stratigraphical data. After a long period of stasis, the ancestral species experiences a phenotypic change expressed as shift and expansion of the morphospace. The appearance of several new phenotypes along with changes in the environment is interpreted as adaptive radiation. Lake-level high stands affect distribution and availability of habitats and, as a result of varied functional demands on shell geometry, the distribution of phenotypes. The on-going divergence of the morphospace into two branches argues for increasing reproductive isolation, consistent with the model of ecological speciation. In the latest phase, however, progressively unstable conditions restrict availability of niches, allowing survival of one branch only.

Pollen and spores record, and paleotemperature reconstructions for mid-Holocene sediments of the Northwest Pacific

Sannai-Maruyama is one of the most famous and best-researched mid-Holocene (mid-Jomon) archaeological sites in Japan, because of a large community of people for a long period. Archaeological studies have shown that the Jomon people inhabited the Sannai-Maruyama site from 5.9-4.2 +/- 0.1 cal. kyr B.P. However, a continuous record of the terrestrial and marine environments around the site has not been available. Core KT05-7 PC-02, was recovered from Mutsu Bay, only 20 km from the site, for the reconstruction of high-resolution time series of environmental records, including sea surface temperature (SST). C37 alkenone SSTs showed clear fluctuations, with four periods of high (8.4-7.9, 7.0-5.9, 5.1-4.1, and 2.3-1.4 cal. kyr B.P.) and four of low (-8.4, 7.9-7.0, 5.9-5.1, and 4.1-2.3 cal. kyr B.P.) SST. Thus, each SST cycle lasted 1.0-2.0 kyr, and the amplitude of fluctuation was about 1.5-2.0 °C. Total organic carbon (TOC) and C37 alkenone contents, and the TOC/total nitrogen ratio indicate that marine biogenic production was low before 7.0 cal. kyr B.P., but was clearly increased between 5.9 and 4.0 cal. kyr B.P., because of stronger vertical mixing. During the period when the community at the site prospered (between 5.9 and 4.2 +/- 0.1 cal. kyr B.P.), the terrestrial climate was relatively warm. The high relative abundance of pollen of both Castanea and Quercus subgen. Cyclobalanopsis supports the interpretation that the local climate was optimal for human habitation. Between 5.9 and 5.1 cal. kyr B.P., in spite of warm terrestrial climates, the C37 alkenone SST was low; this apparent discrepancy may be attributed to the water column structure in the Tsugaru Strait, which differed from the modern condition. The evidence suggests that at about 5.9 cal. kyr B.P, high productivity of marine resources such as fish and shellfish and a warm terrestrial climate led to the establishment of a human community at the Sannai-Maruyama site. Then, at about 4.1 +/- 0.1 cal. kyr B.P., abrupt marine and terrestrial cooling, indicated by a decrease of about 2 °C in the C37 alkenone SST and an increase in pollen of taxa of cooler climates, led to a reduced terrestrial food supply, causing the people to abandon the site. The timing of the abandonment is consistent with the timing (around 4.0-4.3 cal. kyr B.P.) of the decline of civilizations in north Mesopotamia and along the Yangtze River. These findings suggest that a temperature rise of ~2 °C in this century as a result of global warming could have a great impact on the human community and especially on agriculture, despite the advances of contemporary society.

Geochemistry and minerals at DSDP LEg 58 Holes

Leg 58 successfully recovered basalt at Sites 442, 443, and 444, in the Shikoku Basin, and at Site 446 in the Daito Basin. Only at Site 442 did penetration reach unequivocal oceanic layer 2; at the other sites, only off-axis sills and flows were sampled. Petrographic observations indicate that back-arc basalts from the Shikoku Basin, with the exception of the kaersutite-bearing upper sill at Site 444, are mineralogically similar to basalts being erupted at normal mid-ocean ridges. However, the Shikoku Basin basalts are commonly very vesicular, indicating a high volatile content in the magmas. Site 446 in the Daito Basin penetrated a succession of 23 sills which include both kaersutite-bearing and kaersutite-free basalt varieties. A total of 187 samples from the four sites has been analyzed for major and trace elements using X-ray-fluorescence techniques. Chemically, the basalts from Sites 442 and 443 and the lower sill of Site 444 are subalkaline tholeiites and resemble N-type ocean-ridge basalts found along the East Pacific Rise and at 22° N on the Mid-Atlantic Ridge (MAR), although they are not quite as depleted in certain hygromagmatophile (HYG) elements. They do not show any chemical affinities with island-arc tholeiites. The basalts from Site 446 and from the upper sill at Site 444 show alkaline and tholeiitic tendencies, and are enriched in the more-HYG elements; they chemically resemble enriched or E-type basalts and their differentiates found along sections of the MAR (e.g., 45°N) and on ocean islands (e.g., Iceland and the Azores). Most of the intra-site variation may be attributed to crystal settling within individual massive flows and sills, to high-level fractional crystallization in sub-ridge magma chambers, or, where there is evidence of a long period of magmatic quiescence between units, to batch partial melting. However, the basalts from Sites 442 and 443 and from the lower sill at Site 444 cannot easily be related to those from Site 446 and the upper sill at Site 444, and it is possible that the different basalt types were derived from chemically distinct mantle sources. From comparison of the Leg 58 data with those already available for other intra-oceanic back-arc basins, it appears that the mantle sources giving rise to back-arc-basin basalts are chemically as diverse as those for mid-ocean ridges. In addition, the high vesicularity of the Shikoku Basin basalts supports previous observations that the mantle source of back-arc-basin basalts may be contaminated by a hydrous component from the adjacent subduction zone.

Water chemistry of lakes in Zackenbergdalen between 1997-2003, East Greenland

The ecology of arctic lakes is strongly influenced by climate-generated variations in snow coverage and by the duration of the ice-free period, which, in turn, affect the physical and chemical conditions of the lakes (Wrona et al., 2005, http://www.acia.uaf.edu/PDFs/ACIA_Science_Chapters_Final/ACIA_Ch08_Final.pdf). Most arctic lakes are characterised by a long period (8-10 months) of ice-cover, cold water and low algal biomass. The water temperature and nutrient concentrations, and most probably the nutrient input from the catchments, are closely related to the duration of snow- and ice-cover in the lakes. In years when the ice-out is late, - that is, in late July, - phytoplankton photosynthesis is limited by the lack of light and nutrients. Less food is then available to the next link in the food chain, such as copepods and daphnids, with implication on their growth rates.

Carbonate content and d13C record of a Jordan Cretaceous outcrop (Section GM3 (Ghawr Al Mazar (Ghor al Mazrar))

An evaluation of the global synchronicity and duration of "3rd-order" sea-level fluctuations during the Cretaceous greenhouse has been hampered by poor constraints on potential climatic and tectonic drivers, and limitations of geochronology and chronostratigraphic correlation. To provide insight into the nature of such sea-level fluctuations, here we present a new Late Cretaceous record from the Jordanian Levant Platform, comprising a detailed physical-, bio-, chemo- and sequence stratigraphy. Carbonate content of these strata reflects overall sequence stratigraphic development, and demonstrates a dramatic 3rd-order-scale cycle that is also apparent in the d°C record. Updated radioisotopic constraints and astrochronologic testing provide support for the inference of an ~1 million year long sea-level oscillation associated with this 3rd-order cycle, which likely reflects a long-period obliquity (1.2 Myr) control on eustasy and stratigraphic sequence development, linked to the global carbon cycle. The observation of cyclic sea-level fluctuations on this time scale suggests sustained global modulation of continental fresh-water-storage. The hypothesized link between astronomical forcing and sea-level forms a baseline approach in the global correlation of sequence boundaries.

Accumulation rates, grain-size fractions and calcium carbonate at DSDP Sites 89-586 and 90-591

Carbonate oozes recovered by hydraulic piston coring at DSDP Site 586 on Ontong-Java Plateau and Site 591 on Lord Howe Rise have carbonate contents that are consistently higher than 90% with only minor variations. Consequently, paleoceanographic signals were not recorded in detail in the carbonate contents. However, mass accumulation rates of carbonate increased in the late Miocene to mid-Pliocene, reflecting an increase in productivity, then abruptly decreased from mid-Pliocene to the present. Variations in relative abundances of coarse material (foraminifers) and fine material (mostly calcareous nannofossils) do reflect histories of current winnowing and biogenic productivity at the two sites. The late Miocene from 10.5 to 6.5 m.y. ago was a time of relatively constant, quiet, pelagic sedimentation with typical southwest Pacific sedimentation rates of 20-25 m/m.y. The average coarse-fraction abundances are always higher at Site 586 than at Site 591, which reflects winnowing at Site 586. These conditions were interrupted between 6.5 to 4.0 m.y. ago when increased upwelling at the Subtropical Divergence and the Equatorial Divergence produced greater productivity of calcareous planktonic organisms. The increased productivity is suggested by large increases in both fineand coarse-fraction material and constant ratios of foraminifers to nannofossils. The maximum of productivity was about 4.0 m.y. ago. This period of increased upwelling is coincident with the inferred development of the West Antarctic ice sheet. The high productivity was followed by an abrupt increase in winnowing about 2.5 m.y. ago at Site 591, but not until about 2.0 m.y. ago at Site 586. By 2.0 m.y. ago in the late Pliocene, quiet, pelagic sedimentation conditions prevailed, similar to those of the late Miocene. The last 0.7 m.y. has been a period of relatively intense winnowing on Lord Howe Rise but not on Ontong-Java Plateau. The coarse-fraction data have both long- and short-period fluctuations. Long-period fluctuations at Site 591 average about 850 *10**3 yr./cycle and those at Site 586 average 430*10**3 yr./cycle. The highest amplitudes are found in the Pliocene and Quaternary sections. The short-period fluctuations range from 100 to 48*10**3 yr./cycle at Site 586 and from 250 to 33 *10**3 yr./cycle at Site 591. The effects of local fluctuations of productivity and winnowing have modified the primary orbital forcing signals at these two sites to yield complex paleoceanographic records.

Foraminiferal oxygen and carbon isotopic measurements from 5 sites of southern Atlantic and Pacific

The stratigraphy and paleoceanography of the late Miocene and early Pliocene have been examined at six sites in the South Atlantic and southwest Pacific oceans: Deep Sea Drilling Project (DSDP) sites 284, 516A, 519, 588, and 590 and two piston cores from Chain cruise 115. A consistent stratigraphy was developed among sites using graphic correlation, which resulted in age models for all sites that are tied to the revised paleomagnetic time scale of Berggren et al. (1985). Applying these chronologies, we assessed latitudinal and interocean contrasts in the stratigraphic ranges of late Miocene-early Pliocene planktonic foraminiferal and nanno - fossil datums. Salient stratigraphic results include (1) The last appearance datum (LAD) of Globoquadrina dehiscens is a late Miocene (approx. 6.4 Ma) event in the subtropics and is not useful for the placement of the Miocene/Pliocene (M/P) boundary in this biogeographic province. (2) The first appearance datum (FAD) of Globorotalia crassaformis occurred at 5.1 Ma in the South Atlantic near the M/P boundary, suggesting that Gr. crassaformis may have first evolved in the South Atlantic and later migrated to other regions. (3) In the southwest Pacific, the FADs of Gr. margaritae (5.97 Ma), Gr. puncticulata (5.09 Ma), and Gr. crassaformis (4.87 Ma) are significantly time transgressive between temperate and warm subtropical regions. Time lags of 1.0 m.y. were required for these species to adapt to physical and/or biotic conditions peripheral to their endemic biogeographic provinces. (4) Between the subtropics of the South Atlantic and southwest Pacific, many planktonic foraminiferal datums (FAD of Dentogloboquadrina altispira, Gr. cibaoensis, Gr. conomiozea, Gr. margaritae, and Gq. dehiscens and LAD of Gr. cibaoensis) markedly depart from the correlation suggested by magnetostratigraphy, indicating that these datum levels are unreliable for correlation between these ocean basins. (5) In contrast, available calcareous nannofossil datum levels fall on or near the paleomagnetic correlation line, indicating synchroneity of events within the subtropics. (6) Biostratigraphic, magnetic, and 87Sr/86Sr correlation between sites 588 and 519 and the M/P neostratotype at Capo Rossello, Sicily, suggests that the base of the Zanclean stratotype occurs at 5.1-5.0 Ma in the lower reversed subchron of the Gilbert, about 2-3 * 10**5 years above the Gilbert/Chron 5 boundary. Oxygen isotopic results from DSDP sites 284, 519, and CH115 piston cores confirm a prolonged benthic d18O increase in the latest Miocene between 5.6 and 5.0 Ma, as originally proposed by Shackleton and Kennett (1975). At DSDP site 588, the benthic d18O record in the latest Miocene is marked by high-frequency fluctuations with amplitude variations of 0.5per mill, and a long-period wavelength component of 400,000 years. Maximum d18O values, however, occurred during the late Miocene (Kapitean Stage) between 5.5 and 5.1 Ma. The late Miocene d18O changes resulted from mid- and high-latitude cooling and pulses of ice sheet expansion and contraction. Glacial events were most intense during the latest Miocene (Kapitean Stage), and occurred at 5.50-5.35 Ma and at 5.10 Ma. Glacial events are estimated to have lowered sea level by 40 to 60 m and contributed to the isolation and desiccation of the Mediterranean Basin during the late Messinian. Interglacial conditions prevailed at 5.2 Ma and between 5.0 and 4.1 Ma in the early Pliocene. The beginning of the Pliocene was marked by changes in many proxy climatic indicators at all sites, suggesting a prolonged interval of warm, interglacial conditions between 5.0 and 4.1 Ma during the earliest Pliocene.

Statistical grain size parameter of surface samples from Kiel Bay

Geological observations, using "free-diving" techniques (Figure I) were made in September, 1960 and March 1961 along two continuous profiles in the outer Kiel Harbor, Germany and at several other spot locations in the Western Baltic Sea. A distinct terrace, cut in Pleistocene glacial till, was found that was covered with varying amounts and types of recent deposits. Hand samples were taken of the sea-floor sediments and grainsize distribution determined for both the sediment as a whole and for its heavy mineral fraction. From the Laboratory and Field observations it was possible to recognize two distinct types of sand; Type I, Sand resulting from transportation over a long period of time and distance and Type 11, Sand resulting from little transportation and found today near to xvhere it was formed. Several criterea related to the agent of movement could be used to classify the nature of the sediment; (1) undisturbed (the sediment Cover of the Pleistocene Terrace is essentially undisturbed), (2) mixed by organisms, (3) transported by water movements (sediment found with ripple marks, etc., and (4) "Scoured" (the movement of individual particles of sediment from around larger boulders causes a slow downward movement or "Creeping" which is due to both the force of gravity and bottom currents. These observations and laboratory studies are discussed concerning their relationship to the formation of residual sediments, the direction of sand transportation, and the intensive erosion on the outer edge of the wave-cut platform found in this part of the Baltic Sea.

Half-degree gridded nitrogen and phosphorus fertilizer use for global agriculture production during 1900-2013

In addition to enhance agricultural productivity, synthetic nitrogen (N) and phosphorous (P) fertilizer application in croplands dramatically altered global nutrient budget, water quality, greenhouse gas balance, and their feedbacks to the climate system. However, due to the lack of geospatial fertilizer input data, current Earth system/land surface modeling studies have to ignore or use over-simplified data (e.g., static, spatially uniform fertilizer use) to characterize agricultural N and P input over decadal or century-long period. We therefore develop a global time-series gridded data of annual synthetic N and P fertilizer use rate in croplands, matched with HYDE 3,2 historical land use maps, at a resolution of 0.5º latitude by longitude during 1900-2013. Our data indicate N and P fertilizer use rates increased by approximately 8 times and 3 times, respectively, since the year 1961, when IFA (International Fertilizer Industry Association) and FAO (Food and Agricultural Organization) survey of country-level fertilizer input were available. Considering cropland expansion, increase of total fertilizer consumption amount is even larger. Hotspots of agricultural N fertilizer use shifted from the U.S. and Western Europe in the 1960s to East Asia in the early 21st century. P fertilizer input show the similar pattern with additional hotspot in Brazil. We find a global increase of fertilizer N/P ratio by 0.8 g N/g P per decade (p< 0.05) during 1961-2013, which may have important global implication of human impacts on agroecosystem functions in the long run. Our data can serve as one of critical input drivers for regional and global assessment on agricultural productivity, crop yield, agriculture-derived greenhouse gas balance, global nutrient budget, land-to-aquatic nutrient loss, and ecosystem feedback to the climate system.