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.

Time-series of seagrass and land cover maps from 1972 to 2010, in South East Queensland, Australia, with links to GeoTIFFs

Long term global archives of high-moderate spatial resolution, multi-spectral satellite imagery are now readily accessible, but are not being fully utilised by management agencies due to the lack of appropriate methods to consistently produce accurate and timely management ready information. This work developed an object-based remote sensing approach to map land cover and seagrass distribution in an Australian coastal environment for a 38 year Landsat image time-series archive (1972-2010). Landsat Multi-Spectral Scanner (MSS), Thematic Mapper (TM) and Enhanced Thematic Mapper (ETM+) imagery were used without in situ field data input (but still using field knowledge) to produce land and seagrass cover maps every year data were available, resulting in over 60 map products over the 38 year archive. Land cover was mapped annually using vegetation, bare ground, urban and agricultural classes. Seagrass distribution was also mapped annually, and in some years monthly, via horizontal projected foliage cover classes, sand and deep water. Land cover products were validated using aerial photography and seagrass maps were validated with field survey data, producing several measures of accuracy. An average overall accuracy of 65% and 80% was reported for seagrass and land cover products respectively, which is consistent with other studies in the area. This study is the first to show moderate spatial resolution, long term annual changes in land cover and seagrass in an Australian environment, created without the use of in situ data; and only one of a few similar studies globally. The land cover products identify several long term trends; such as significant increases in South East Queensland's urban density and extent, vegetation clearing in rural and rural-residential areas, and inter-annual variation in dry vegetation types in western South East Queensland. The seagrass cover products show that there has been a minimal overall change in seagrass extent, but that seagrass cover level distribution is extremely dynamic; evidenced by large scale migrations of higher seagrass cover levels and several sudden and significant changes in cover level. These mapping products will allow management agencies to build a baseline assessment of their resources, understand past changes and help inform implementation and planning of management policy to address potential future changes.

SHRIMP U-Pb zircon data from Gjelsvikfjella, Dronning Maud Land, East Antarctica

The Maud Belt in Dronning Maud Land (western East Antarctic Craton) preserves a high-grade polyphase tectono-thermal history with two orogenic episodes of Mesoproterozoic (1.2-1.0 Ga) and Neoproterozoic (0.6-0.5 Ga) age. New SHRIMP U-Pb zircon data from southern Gjelsvikfjella in the northeastern part of the belt make it possible to differentiate between a series of magmatic and metamorphic events. The oldest event recorded is the formation of an extensive 1140-1130 Ma volcanic arc. This was followed by 1104 ± 8 Ma granitoids that might represent, together with so far undated mafic dykes, part of a decompression melting-related bimodal suite that reflects the sub-continental Umkondo igneous event. The first high-grade metamorphism is constrained at 1070 Ma. The metamorphic age data are similar to those obtained from other parts of the Maud Belt, but also from the Namaqua-Natal Belt in South Africa, but the preceding arc formation was diachronous in the two belts. This indicates that the two belts did not form a continuous volcanic arc unit as suggested in previous models, but became connected only at the end of the Mesoproterozoic. Intense reworking during the Neoproterozoic, probably as a result of continent-continent collision between components of Gondwana, is indicated by ductile refliation, further high-grade metamorphic recrystallisation and metamorphic zircon overgrowths at approximately 530 Ma. This was followed by late- to post-tectonic magmatism, reflected by 500 Ma granite bodies and 490 Ma aplite dykes as well as a 480 Ma gabbro body.

Tab. 1: U-Pb analyses of zircons from rock samples in MT. Newton and Cumpston Massif, East Antarctica

We present new U-Pb zircon (SHRIMP) data on rocks from Mt Newton and Cumpston Massif in the southern Prince Charles Mountains. Our data demonstrate that Mt Newton was affected by a newly proposed Palaeoproterozoic "Newton" Orogeny at c. 2100-2200 Ma. Sedimentation, felsic volcanism (c. 2200 Ma), metamorphism and folding, followed by granite intrusion (c. 2100 Ma), suggest development of a trough or aulacogene in the area during the early Palaeoproterozoic. An orthogneiss from Cumpston Massif yielded an age of c. 3180 Ma for granitic protolith emplacement, which is in good agreement with many U-Pb zircon ages from similar rocks in the southern Mawson Escarpment. A syn- to late-tectonic muscovite-bearing pegmatite from Cumpston Massif yielded a c. 2500 Ma date of emplacement, which indicates early Palaeoproterozoic activity in this block, probably in response to a tectono-magmatic episode in the Lambert Terrane bordering the Ruker Terrane in the northeast. The correlation of tectono-magmatic events in both the Ruker and Lambert terranes of the southern Prince Charles Mountains provides evidence for their common evolution during the Proterozoic.

Maps (pdf) and raster images (tif) of predicted rural land cover suitability under current (2010) conditions and future climate scenarios

This study projects land cover probabilities under climate change for corn (maize), soybeans, spring and winter wheat, winter wheat-soybean double cropping, cotton, grassland and forest across 16 central U.S. states at a high spatial resolution, while also taking into account the influence of soil characteristics and topography. The scenarios span three oceanic-atmospheric global circulation models, three Representative Concentration Pathways, and three time periods (2040, 2070, 2100). As climate change intensifies, the suitable area for all six crops display large northward shifts. Total suitable area for spring wheat, followed by corn and soybeans, diminish. Suitable area for winter wheat and for winter wheat-soybean double-cropping expand northward, while cotton suitability migrates to new, more northerly, locations. Suitability for forest intensifies in the south while yielding to crops in the north; grassland intensifies in the western Great Plains as crop suitability diminishes. To maintain current broad geographic patterns of land use, large changes in the thermal response of crops such as corn would be required. A transition from corn-soybean to winter wheat-soybean doubling cropping is an alternative adaptation.

SHRIMP U-Pb dating of high-grade migmatites and related magmatites from Northwestern Oates Land, East Antarctica

High- to very-high-grade migmatitic basement rocks of the Wilson Hills area in northwestern Oates Land (Antarctica) form part of a low-pressure high-temperature belt located at the western inboard side of the Ross-orogenic Wilson Terrane. Zircon, and in part monazite, from four very-high grade migmatites (migmatitic gneisses to diatexites) and zircon from two undeformed granitic dykes from a central granulite-facies zone of the basement complex were dated by the SHRIMP U-Pb method in order to constrain the timing of metamorphic and related igneous processes and to identify possible age inheritance. Monazite from two migmatites yielded within error identical ages of 499 +/- 10 Ma and 493 +/- 9 Ma. Coexisting zircon gave ages of 500 +/- 4 Ma and 484 +/- 5 Ma for a metatexite (two age populations) and 475 +/- 4 Ma for a diatexite. Zircon populations from a migmatitic gneiss and a posttectonic granitic dyke yielded well-defined ages of 488 +/- 6 Ma and 482 +/- 4 Ma, respectively. There is only minor evidence of age inheritance in zircons of these four samples. Zircon from two other samples (metatexite, posttectonic granitic dyke) gave scattered 206Pb-238U ages. While there is a component similar in age and in low Th/U ratio to those of the other samples, inherited components with ages up to c. 3 Ga predominate. In the metatexite, a major detrital contribution from 545 - 680 Ma old source rocks can be identified. The new age data support the model that granulite- to high-amphibolite-facies metamorphism and related igneous processes in basement rocks of northwestern Oates Land were confined to a relatively short period of time of Late Cambrian to early Ordovican age. An age of approximately 500 Ma is estimated for the Ross-orogenic granulite-facies metamorphism from consistent ages of monazite from two migmatites and of the older zircon age population in one metatexite. The variably younger zircon ages are interpreted to reflect mineral formation in the course of the post-granulite-facies metamorphic evolution, which led to a widespread high-amphibolite-facies retrogression and in part late-stage formation of ms+bi assemblages in the basement rocks and which lasted until about 465 Ma. The presence of inherited zircon components of latest Neoproterozoic to Cambrian age indicates that the high- to very-grade migmatitic basement in northwestern Oates Land originated from clastic series of Cambrian age and, therefore, may well represent the deeper-crustal equivalent of lower-grade metasedimentary series of the Wilson Terrane.

Geological and geochemical characteristics of the Tersang gold deposit, Malaysia: implications for genesis and gold exploration

The Central gold belt of peninsular Malaysia comprises a number of gold deposits located in the east of the N-S striking Bentong-Raub Suture Zone. The Tersang gold deposit is one of the gold deposits in the gold belt and hosted in sandstone, rhyolite and breccia units. The deposit has an inferred resource of 528,000 ounces of gold. The geochronology of the Tersang deposit has been newly constrained by LA ICP-MS U-Pb zircon dating. The maximum depositional age of the host sedimentary rocks ranges from Early Carboniferous to Early Permian (261.5 ± 4.9 Ma to 333.5 ± 2.5 Ma) for the host sandstone and Late Triassic for the rhyolite intrusion (218.8 ± 1.7 Ma). Textural characteristics of pyrite have revealed five types including (1) Euhedral to subhedral pyrite with internal fracturing and porous cores located in the sandstone layers (pyrite 1); (2) Anhedral pyrite overgrowths on pyrite 1 and disseminated in stage 1 vein (pyrite 2); (3) Fracture-filled or vein pyrite located in stages 1 and 2 vein (pyrite 3); (4) Euhedral pyrite with internal fractures also located in stage 2 vein (pyrite 4); and (5) Subhedral clean pyrite located in the rhyolite intrusion (pyrite 5). Based on pyrite mapping and spot analyses, two main stages of gold enrichment are documented from the Tersang gold deposit. Gold in sandstone-hosted pyrite 1 (mean 4.3 ppm) shows best correlation with Bi and Pb (as evidenced on pyrite maps). In addition, gold in pyrite 3 (mean 8 ppm) located in stage 2 vein shows a good correlation with As, Ag, Sb, Cu, Tl, and Pb. In terms of gold exploration, we suggest that elements such as As, Ag, Sb, Cu, Tl, Bi, and Pb associated with Au may serve as vectoring tools in gold exploration. Our new geological, structural, geochemical and isotopic data together with mineral paragenesis, pyrite chemistry and ore fluid characteristics indicate that the Tersang gold deposit is comparable to a sediment-hosted gold deposit. Our new genetic model suggests deposition of the Permo-Carboniferous sediments followed by intrusion of rhyolitic magma in the Late Triassic. At a later stage, gold mineralisation overprinted the rhyolite intrusion and the sandstone.

Benthic foraminiferal stable isotope record from ODP Site 138-849 on the East Pacific Rise

Benthic foraminifer and delta13C data from Site 849, on the west flank of the East Pacific Rise (0°11 'N, 110°31'W; 3851 m), give relatively continuous records of deep Pacific Ocean stable isotope variations between 0 and 5 Ma. The mean sample spacing is 4 k.y. Most analyses are from Cibicides wuellerstorfi, but isotopic offsets relative to Uvigerina peregrina appear roughly constant. Because of its location west of the East Pacific Rise, Site 849 yields a suitable record of mean Pacific Ocean delta13C, which approximates a global oceanic signal. The ~100-k.y.-period climate cycle, which is prevalent in delta18O does not dominate the long-term delta13C record. For delta13C, variations in the ~400- and 41-k.y. periods are more important. Phase lags of delta13C relative to ice volume in the 41- and 23-k.y. bands are consistent with delta13C as a measure of organic biomass. A model-calculated exponential response time of 1-2 k.y. is appropriate for carbon stored in soils and shallow sediments responding to glacial-interglacial climate change. Oceanic delta13C leads ice volume slightly in the 100-k.y. band, and this suggests another process such as changes in continental weathering to modulate mean river delta13C at long periods. The delta13C record from Site 849 diverges from that of Site 677 in the Panama Basin mostly because of decay of 13C-depleted organic carbon in the relatively isolated Panama Basin. North Atlantic to Pacific delta13C differences calculated using published data from Sites 607 and 849 reveal variations in Pliocene deep water within the range of those of the late Quaternary. Maximum delta13C contrast between these sites, which presumably reflects maximum influx of high-delta13C northern source water into the deep North Atlantic Ocean, occurred between 1.3 and 2.1 Ma, well after the initiation of Northern Hemisphere glaciation. Export of high-delta13C North Atlantic Deep Water from the Atlantic to the circumpolar Antarctic, as recorded by published delta13C data from Subantarctic Site 704, appears unrelated to the North Atlantic-Pacific delta13C contrast. To account for this observation, we suggest that deep-water formation in the North Atlantic reflects northern source characteristics, whereas export of this water into the circumpolar Antarctic reflects Southern Hemisphere wind forcing. Neither process appears directly linked to ice-volume variations.

Great Barrier Reef coral element/Ca and stable isotope data and U-Th ages from IODP Expedition 325

Tropical south-western Pacific temperatures are of vital importance to the Great Barrier Reef (GBR), but the role of sea surface temperatures (SSTs) in the growth of the GBR since the Last Glacial Maximum remains largely unknown. Here we present records of Sr/Ca and d18O for Last Glacial Maximum and deglacial corals that show a considerably steeper meridional SST gradient than the present day in the central GBR. We find a 1-2 °C larger temperature decrease between 17° and 20°S about 20,000 to 13,000 years ago. The result is best explained by the northward expansion of cooler subtropical waters due to a weakening of the South Pacific gyre and East Australian Current. Our findings indicate that the GBR experienced substantial meridional temperature change during the last deglaciation, and serve to explain anomalous deglacial drying of northeastern Australia. Overall, the GBR developed through significant SST change and may be more resilient than previously thought.

(Table 1) Calcium isotope composition, Mg/Ca, Sr/Ca and U/Ca ratios of calcareous cysts of Thoracosphaera heimii

We investigated the influences of temperature, salinity and pH on the calcium isotope as well as trace and minor element (uranium, strontium, magnesium) to Ca ratios on calcium carbonate cysts of the calcareous dinoflagellate species Thoracosphaera heimii grown in laboratory cultures. The natural habitat of this species is the photic zone (preferentially at the chlorophyll maximum depth) of temperate to tropical oceans, and it is abundant in deep-sea sediments over the entire Cenozoic. In our experiments, temperatures ranged from 12 to 30 °C, salinity from 36.5 to 38.8 and pH from 7.9 to 8.4. The delta44/40Ca of T. heimii cysts resembles that of other marine calcifiers, including coccolithophores, foraminifers and corals. However, its temperature sensitivity is considerably smaller and statistically insignificant, and T. heimii might serve as a recorder of changes in seawater delta44/40Ca over geologic time. The Sr/Ca ratios of T. heimii cysts show a pronounced temperature sensitivity (0.016 mmol/mol °C**-1) and have the potential to serve as a palaeo-sea surface temperature proxy. No clear temperature- and pH-dependences were observed for Mg/Ca. U/Ca seems to be influenced by temperature and pH, but the correlations change sign at 23 °C and pH 8.2, respectively.

(Table 1) U and Th isotope analyses of pyrite and anhydrite of ODP Site 158-957

Hydrothermal circulation at oceanic spreading ridges causes sea water to penetrate to depths of 2 to 3 km in the oceanic crust where it is heated to ~400 °C before venting at spectacular 'black smokers'. These hydrothermal systems exert a strong influence on ocean chemistry (Edmond et al., 1979, doi:10.1016/0012-821X(79)90061-X), yet their structure, longevity and magnitude remain largely unresolved (Elderfield and Schultz., 1996, doi:10.1146/annurev.earth.24.1.191). The active Transatlantic Geotraverse (TAG) deposit, at 26° N on the Mid-Atlantic Ridge, is one of the largest, oldest and most intensively studied of the massive sulphide mounds that accumulate beneath black-smoker fields. Here we report ages of sulphides and anhydrites from the recently drilled (Humphris et al., 1995, doi:10.1038/377713a0) TAG substrate structures -determined from 234U-230Th systematics analysed by thermal ionization mass spectrometry. The new precise ages combined with existing data (Lalou et al., 1993, doi:10.1029/92JB01898; 1998, doi:10.2973/odp.proc.sr.158.214.1998) show that the oldest material (11,000 to 37,000 years old) forms a layer across the centre of the deposit with younger material (2,300-7,800 years old) both above and below. This stratigraphy confirms that much of the sulphide and anhydrite are precipitated within the mound by mixing of entrained sea water with hydrothermal fluid (James and Elderfield, 1996, doi:10.1130/0091-7613(1996)024<1147:COOFFA>2.3.CO;2). The age distribution is consistent with episodic activity of the hydrothermal system recurring at intervals of up to 2,000 years.

Net primary productivity, POC export, Th234 and U238 activities and abundance of autotrophs during U.S.C.G.C. Healy cruises HLY0802 and HLY0803

Seasonal patterns in the partitioning of phytoplankton carbon during receding sea ice conditions in the eastern Bering Sea water column are presented using rates of 14C net primary productivity (NPP), phototrophic plankton carbon content, and POC export fluxes from shelf and slope waters in the spring (March 30-May 6) and summer (July 3-30) of 2008. At ice-covered and marginal ice zone (MIZ) stations on the inner and middle shelf in spring, NPP averaged 76 ± 93 mmol C/m**2/d, and in ice-free waters on the outer shelf NPP averaged 102 ± 137 mmol C/m**2/d. In summer, rates of NPP were more uniform across the entire shelf and averaged 43 ± 23 mmol C/m**2/d over the entire shelf. A concomitant shift was observed in the phototrophic pico-, nano-, and microplankton community in the chlorophyll maximum, from a diatom dominated system (80 ± 12% autotrophic C) in ice covered and MIZ waters in spring, to a microflagellate dominated system (71 ± 31% autotrophic C) in summer. Sediment trap POC fluxes near the 1% PAR depth in ice-free slope waters increased by 70% from spring to summer, from 10 ± 7 mmol C/m**2/d to 17 ± 5 mmol C/m**2/d, respectively. Over the shelf, under-ice trap fluxes at 20 m were higher, averaging 43 ± 17 mmol C/m**2/d POC export over the shelf and slope estimated from 234Th deficits averaged 11 ± 5 mmol C/m**2/d in spring and 10 ± 2 mmol C/m**2/d in summer. Average e-ratios calculated on a station-by-station basis decreased by ~ 30% from spring to summer, from 0.46 ± 0.48 in ice-covered and MIZ waters, to 0.33 ± 0.26 in summer, though the high uncertainty prevents a statistical differentiation of these data.