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.

Sediment age, abundance of planktonic foraminifera shells, degree of their dissolution, percentage of benthic foraminifera in total foram fauna, Core GC-11, Sea of Okhotsk

Taxonomic composition and distribution of planktonic foraminifera are studied in section of Core GC-11 penetrated through Upper Quaternary sediments of the Bowers Ridge western slope, south Bering Sea. It is shown that structure of foraminiferal assemblage and productivity varied substantially during the last 32000 calendar years in response to changes in surface water temperatures and water mass circulation in the North Pacific including the Bering Sea. Productivity was maximal during the deglaciation epoch, being notably lower in Holocene and minimal at glaciation time.

(Table 2) Petrophysical, petrological, and geochemical characteristics of ODP Hole 163-990A sediments

Ocean Drilling Program Hole 990A penetrated 131 m of subaerially emplaced Paleocene flood basalts on the Southeast Greenland margin with a recovery of 74%. Shipboard P-wave velocity (Vp), density, and magnetic susceptibility were measured with 2- to 15-cm intervals on the core. Individual flow units were divided into four zones based on the observed petrophysical characteristics. From the top, these are Zone I (<7 m thick with a Vp of ~2.5 km/s), Zone II (3-5 m thick with a strongly increasing Vp from 2.5 to 5.5 km/s), Zone III (up to 20 m thick with a Vp of ~5.5-6.0 km/s), and Zone IV (<2 m thick with a strongly decreasing Vp from 6.0 to 2.5 km/s). Eighteen samples were selected from three of the fourteen penetrated basalt units for geochemical, petrological, and petrophysical studies focusing on the altered, low-velocity upper lava Zones I and II. Zone I is strongly altered to >50% clay minerals (smectite) and iron hydroxides, and the petrophysical properties are primarily determined by the clay properties. Zone II is intermediately altered with 5%-20% clay minerals, where the petrophysical properties are a function of both the degree of alteration and porosity variations. Shipboard and shore-based measurements of the same samples show that storage permanently lowers the elastic moduli of basalt from Zones I to III. This is related to the presence of even small quantities of swelling clays. The data show that alteration processes are important in determining the overall seismic properties of flood basalt constructions. The degree and depth of alteration is dependent on the primary lava flow emplacement structures and environment. Thus, the interplay of primary emplacement and secondary alteration structures determine the elastic properties of basalt piles. Rock property theories for sand-clay systems are further used to model the physical property variations in these altered crystalline rocks.