Particulate organic carbon fluxes to the ocean interior and factors controlling the biological pump: A synthesis of global sediment trap programs since 1983

Particulate organic carbon (POC) is vertically transported to the oceanic interior by aggregates and their ballasts, mainly CaCO₃ and biogenic opal, with a smaller role for lithogenic aerosols through the mesopelagic zone. Diel migrating zooplankton communities effect vertical transport and remineralization of POC in the upper layers of the ocean. Below 1.5 km, the presence of zooplankton is reduced and thus the aggregates travel mainly by gravitational transport. We normalized the fluxes of POC, CaCO₃, and biogenic opal from data published on samples collected at 134 globally distributed, bottom-tethered, time-series sediment trap (TS-trap) stations to annual mole fluxes at the mesopelagic/bathypelagic boundary (m/b) at 2 km and defined them as F_m/bC_org, F_m/bC_inorg, and F_m/bSi_bio. Using this global data set, we investigated (1) the geographic contrasts of POC export at m/b and (2) the supply rate of ∑CO₂ to the world mesopelagic water column. F_m/bC_org varies from 25 (Pacific Warm Pool) to 605 (divergent Arabian Sea) mmolC m⁻² yr⁻¹; F_m/bC_inorg varies from >8 (high latitude Polar Oceans) or 15 (Pacific Warm Pool) to 459 (divergent Arabian Sea) mmolC m⁻² yr⁻¹; and F_m/bSi_bio, the most spatially/temporally variable flux, ranges from 6 (North Atlantic Drift) to 1118 (Pacific Subarctic Gyre) mmolSi m⁻² yr⁻¹. The oceanic region exhibiting the highest POC flux over a significantly large region is the area of the North Pacific Boreal Gyres where the average F_m/bC_org = 213, F_m/bC_inorg = 126, and F_m/bSi_bio = 578 mmol m⁻² yr⁻¹. F_m/bC_org and F_m/bC_inorg are particularly high in large upwelling margins, including the divergent Arabian Sea and off Cape Verde. One of the data sets showing the lowest flux over a significant region/basin is F_m/bC_org = 39, F_m/bC_inorg = 69, and F_m/bSi_bio = 22 mmol m⁻² yr⁻¹ in the North Pacific subtropical/tropical gyres; Pan-Atlantic average fluxes are similar except F_m/bSi_bio fluxes are even lower. Where C_org/C_inorg and Si_bio/C_inorg are <1 defines the “Carbonate Ocean”, and where these ratios are ?1 defines the “Silica Ocean”. The Carbonate Ocean occupies about 80% of the present world pelagic ocean between the two major oceanographic fronts, the North Pacific Polar Front and the Antarctic Polar Front, and the Silica Ocean is found on the polar sides of these fronts. The total global annual fluxes of F_m/bC_org, F_m/bC_inorg, and F_m/bSi_bio at m/b calculated by parameterizations of the export flux data from 134 stations are surprisingly similar; 36.2, 33.8, and 34.6 teramol yr⁻¹ (120, 112, and 114 mmol m⁻² yr⁻¹), respectively, resulting in a near uniform binary ratio between the above three elements of about one. The global ternary % ratios estimated from 152 TS-trap samples of the three elements are 35:32:33. From our global F_m/bC_org and a published model estimate of the global export production, we estimate the regeneration rate of CO₂ through the mesopelagic zone by the biological pump is 441 teramolC yr⁻¹. Based on our global F_m/bC_inorg and recently estimated global primary production of PIC, 36-86 teramolC yr⁻¹ of PIC is assumed to be dissolved within the upper 2 km of the water column.     

Physical oceanographic conditions in the deepwater Gulf of Mexico in summer 2000–2002

The circulation and distribution of water properties in the water column of the Gulf of Mexico influence the flux of carbon to the benthic environment. The eddy field of the upper 1000 m creates environmental conditions that are favorable for biological productivity in an otherwise oligotrophic subtropical ocean. This eddy field results in the transport of nutrients and organic matter into the photic zone through cross-margin flow of shelf waters, upwelling in cyclones, and uplift from the interaction of anticyclones with bathymetry. These conditions then allow the productivity that becomes a possible source of carbon to the benthos.Data from four cruises during summers of 2000–2002 are used to describe the currents and water property distributions in the deepwater Gulf of Mexico, which consists of water depths greater than 400 m. Comparisons are made to historical data sets to provide an understanding of the persistence of the characteristics of the Gulf and the processes that occur there.The currents in the Gulf are surface intensified, have minimum in 800–1000 m depths, and also exhibit bottom intensification, especially near sloping topography. Historical time series records show current speeds near-bottom reach 50–100 cm s⁻¹. At basin scales, these currents tend to flow cyclonically (counter-clockwise) along the bathymetry. These near-bottom, episodic, high-speed currents provide a mechanism for the transport of organic material, in both large and small particle sizes, from one benthic area to another.The distributions of temperature, salinity, nutrients, and dissolved oxygen during the study appear to be unchanged from historical findings. The source waters for the deep Gulf are the water masses brought into the Gulf by the Loop Current system. The properties in the upper 100–200 m are the most variable of the water column, consistent with their proximity to wind mixing, river discharge mixing, and atmospheric influences. Below 1500 m, there are no major horizontal variations in these water properties.     

Formation of the Summertime Ozone Valley over the Tibetan Plateau: The Asian Summer Monsoon and Air Column Variations

The summertime ozone valley over the Tibetan Plateau is formed by two influences,the Asian summer monsoon(ASM) and air column variations.Total ozone over the Tibetan Plateau in summer was ～33 Dobson units(DU) lower than zonal mean values over the ocean at the same latitudes during the study period 2005-2009.Satellite observations of ozone profiles show that ozone concentrations over the ASM region have lower values in the upper troposphere and lower stratosphere(UTLS) than over the non-ASM region.This is caused by frequent convective transport of low-ozone air from the lower troposphere to the UTLS region combined with trapping by the South Asian High.This offset contributes to a ～20-DU deficit in the ozone column over the ASM region.In addition,along the same latitude,total ozone changes identically with variations of the terrain height,showing a high correlation with terrain heights over the ASM region,which includes both the Tibetan and Iranian plateaus.This is confirmed by the fact that the Tibetan and Iranian plateaus have very similar vertical distributions of ozone in the UTLS,but they have different terrain heights and different total-column ozone levels.These two factors(lower UTLS ozone and higher terrain height) imply 40 DU in the lower-ozone column,but the Tibetan Plateau ozone column is only ～33 DU lower than that over the non-ASM region.This fact suggests that the lower troposphere has higher ozone concentrations over the ASM region than elsewhere at the same latitude,contributing ～7 DU of total ozone,which is consistent with ozonesonde and satellite observations.     

Hydrothermal uranium deposits containing molybdenum and fluorite in the Marysvale volcanic field, west-central Utah

Uranium deposits containing molybdenum and fluorite occur in the Central Mining Area, near Marysvale, Utah, and formed in an epithermal vein system that is part of a volcanic/hypabyssal complex. They represent a known, but uncommon, type of deposit; relative to other commonly described volcanic-related uranium deposits, they are young, well-exposed and well-documented. Hydrothermal uranium-bearing quartz and fluorite veins are exposed over a 300 m vertical range in the mines. Molybdenum, as jordisite (amorphous MoS₂), together with fluorite and pyrite, increase with depth, and uranium decreases with depth. The veins cut 23-Ma quartz monzonite, 20-Ma granite, and 19-Ma rhyolite ash-flow tuff. The veins formed at 19-18 Ma in a 1 km² area, above a cupola of a composite, recurrent, magma chamber at least 24 × 5 km across that fed a sequence of 21- to 14-Ma hypabyssal granitic stocks, rhyolite lava flows, ash-flow tuffs, and volcanic domes. Formation of the Central Mining Area began when the intrusion of a rhyolite stock, and related molybdenite-bearing, uranium-rich, glassy rhyolite dikes, lifted the fractured roof above the stock. A breccia pipe formed and relieved magmatic pressures, and as blocks of the fractured roof began to settle back in place, flat-lying, concave-downward, “pull-apart” fractures were formed. Uranium-bearing, quartz and fluorite veins were deposited by a shallow hydrothermal system in the disarticulated carapace. The veins, which filled open spaces along the high-angle fault zones and flat-lying fractures, were deposited within 115 m of the ground surface above the concealed rhyolite stock. Hydrothermal fluids with temperatures near 200 °C, ¹⁸O_H2O∼−1.5, D_H2O∼−130, log f O₂ about −47 to −50, and pH about 6 to 7, permeated the fractured rocks; these fluids were rich in fluorine, molybdenum, potassium, and hydrogen sulfide, and contained uranium as fluoride complexes. The hydrothermal fluids reacted with the wallrock resulting in precipitation of uranium minerals. At the deepest exposed levels, wallrocks were altered to sericite; and uraninite, coffinite, jordisite, fluorite, molybdenite, quartz, and pyrite were deposited in the veins. The fluids were progressively oxidized and cooled at higher levels in the system by boiling and degassing; iron-bearing minerals in wall rocks were oxidized to hematite, and quartz, fluorite, minor siderite, and uraninite were deposited in the veins. Near the ground surface, the fluids were acidified by condensation of volatiles and oxidation of hydrogen sulfide in near-surface, steam-heated, ground waters; wall rocks were altered to kaolinite, and quartz, fluorite, and uraninite were deposited in veins. Secondary uranium minerals, hematite, and gypsum formed during supergene alteration later in the Cenozoic when the upper part of the mineralized system was exposed by erosion. Received: 23 June 1997 / Accepted: 15 October 1997     

Representativeness of wind measurements on a mesoscale grid with station separations of 312 m to 10 km

A field experiment was carried out in which wind speed and direction were measured over flat terrain at a height of 10 m using 13 identical instruments spaced logarithmically along two perpendicular 10 km lines. Station separations ranged from 312 m to 10 km. One-minute data from 11 sampling periods of duration 6 to 10 h were studied. p ]The statistics showed little dependence on whether the line of instruments was oriented along the wind or across the wind. The correlation coefficients between wind fluctuations at two stations separated by distance x were found to vary exponentially with x, with an integral distance scale on the order of 1 km. The integral time scale derived from the variation of the single station variances with averaging time was found to equal several minutes. At a station separation of 10 km, the correlation coefficients between the wind components at the two sites were calculated to be 0.24, 0.37, and 0.47 for averaging times of 1, 10, and 60 min, respectively. These values for the correlation coefficients correspond to root-mean-square differences in wind speed at the two stations of about 1.3, 1.0, and 0.7 m/s, respectively.Exponential formulas based on dimensional analysis are suggested for fitting these observations. It is found that the observations of spatial correlations are best fit if two independent integral distance scales are used — a boundary-layer distance scale of about 300 m that best applies to small station separations and a mesoscale distance scale of about 10 km that applies to larger station separations.     

The impact of the MJO on clusters of wintertime circulation anomalies over the North American region

Recent studies have shown that the Madden–Julian Oscillation (MJO) impacts the leading modes of intraseasonal variability in the northern hemisphere extratropics, providing a possible source of predictive skill over North America at intraseasonal timescales. We find that a k-means cluster analysis of mid-level geopotential height anomalies over the North American region identifies several wintertime cluster patterns whose probabilities are strongly modulated during and after MJO events, particularly during certain phases of the El Niño-Southern Oscillation (ENSO). We use a simple new optimization method for determining the number of clusters, k, and show that it results in a set of clusters which are robust to changes in the domain or time period examined. Several of the resulting cluster patterns resemble linear combinations of the Arctic Oscillation (AO) and the Pacific/North American (PNA) teleconnection pattern, but show even stronger responses to the MJO and ENSO than clusters based on the AO and PNA alone. A cluster resembling the positive (negative) PNA has elevated probabilities approximately 8–14 days following phase 6 (phase 3) of the MJO, while a negative AO-like cluster has elevated probabilities 10–20 days following phase 7 of the MJO. The observed relationships are relatively well reproduced in the 11-year daily reforecast dataset from the National Centers for Environmental Prediction (NCEP) Climate Forecast System version 2 (CFSv2). This study statistically links MJO activity in the tropics to common intraseasonal circulation anomalies over the North American sector, establishing a framework that may be useful for improving extended range forecasts over this region.     

Does increasing the spatial resolution of a regional climate model improve the simulated daily precipitation?

Three different resolution (50, 12, and 1.5 km) regional climate model simulations are compared in terms of their ability to simulate moderate and high daily precipitation events over the southern United Kingdom. The convection-permitting 1.5-km simulation is carried out without convective parametrisation. As in previous studies, increasing resolution (especially from 50 to 12 km) is found to improve the representation of orographic precipitation. The 50-km simulation underestimates mean precipitation over the mountainous region of Wales, and event intensity tends to be too weak; this bias is reduced in both the 12- and 1.5-km simulations for both summer and winter. In south–east England lowlands where summer extremes are mostly convective, increasing resolution does not necessary lead to an improvement in the simulation. For the 12-km simulation, simulated daily extreme events are overly intense. Even though the average intensity of summer daily extremes is improved in the 1.5-km simulation, this simulation has a poorer mean bias with too many events exceeding high thresholds. Spatial density and clustering of summer extremes in south–east England are poorly simulated in both the 12- and 1.5-km simulations. In general, we have not found any clear evidence to show that the 1.5-km simulation is superior to the 12-km simulation, or vice versa at the daily level.     

Examination of ‘global atmospheric temperature monitoring with satellite microwave measurements’: 2. Analysis of satellite data

Using Microwave Sounding Unit (MSU) channel 2 (Ch. 2, 53.74 GHz) data, Spencer and Christy (1992a) determined that the earth exhibits no temperature trend in the period 1979–90, while other authors find a temperature increase of roughly 0.1 K. Based on a theoretical analysis Prabhakara et al. (1995) showed that the information about the global atmospheric temperature deduced from MSU Ch. 2 observations has a small contamination, T ₂, as a result of the attenuation due to hydrometeors in the atmosphere. A method is developed in this study, that utilizes coincident measurements made by MSU in Ch. 1 (50.3 GHz), to estimate this T ₂ over the global oceans. The magnitude of T ₂ is found to be about 1 K over significant parts of the tropical oceanic rain belts and about 0.25 K over minor portions of the mid-latitude oceanic storm tracks. Due to events such as El Niôo, there is variability from year to year in the rain areas and rain intensity leading to significant change in the patterns of T ₂. The patterns of T ₂ derived for March 82 and March 83 reveal such a change. When averaged over the global oceans, from 50° N to 50° S, T ₂ has a value of 0.25 and 0.29 K for March 1982 and 1983, respectively. Due to these reasons the interannual temperature change derived by Spencer and Christy from MSU Ch. 2 will contain a residual hydrometeor effect. Thus in evaluating decadal trend of the global mean temperature of the order of 0.1 K from MSU Ch. 2 data one has to take into account completely the contamination due to hydrometeors.     

Zonal winds and southeast Australian rainfall in global and regional climate models

Southeast Australia is a region of high rainfall variability related to major climate drivers, with a long-term declining trend in cool-season rainfall. Projections of future rainfall trends are uncertain in this region, despite projected southward shifts in the subtropical ridge and mid-latitude westerlies. This appears to be related to a poor representation of the spatial relationships between rainfall variability and zonal wind patterns across southeast Australia in the latest Coupled Model Intercomparison Project ensemble, particularly in the areas where weather systems embedded in the mid-latitude westerlies are the main source of cool-season rainfall. Downscaling with regional climate models offers improvements in the mean rainfall climatology, and shows some ability to correct for poor modelled relationships between rainfall and zonal winds along the east coast of Australia. However, it provides only minor improvements to these relationships in southeast Australia, despite the improved representation of topographic features. These results suggest that both global and regional climate models may fail to translate projected circulation changes into their likely rainfall impacts in southeast Australia.     

Rossby wave propagation into the tropics in two GFDL general circulation models

During the northern winter the eastern Pacific is characterized by upper level westerly flow extending from the equator into the midlatitudes of both hemispheres. Theoretical and simple modeling studies suggest that such a region should favor the penetration of Rossby waves into the tropics from higher latitudes. Observational results by Kiladis and Weickmann using 200 mb data indicate that Rossby waves do indeed propagate freely into the tropical eastern Pacific during the northern winter from the Asian jet exit region. They also confirmed that cross-equatorial dispersion of energy from the Northern into the Southern Hemisphere occurs frequently. The present study examines these interactions in climatological runs of two GFDL GCMs. The northern wintertime mean states of these models are characterized by a rather realistically simulated upper level westerly regime in the tropical Pacific. Despite the relative weakness of the Asian jet and wave activity with respect to observations, propagation of Rossby waves into the tropics is present in both models, and these waves are strongly positively tilted as seen in the observations. A momentum budget of the zonal wind and E vector diagnostics over the tropical Pacific indicate that these transients are an important component of the momentum balance of the equatorial westerlies in both the observations and in the models.This paper was presented at the Second International Conference on Modelling of Global Climate Variability, held in Hamburg 7–11 September 1992 under the auspices of the Max Planck Institute for Meteorology. Guest Editor for these papers is L. Dümenil