ALBERT

Description: The seasonal variation of the intrusion of the Philippine Sea Water into the South China Sea was studied by analyzing the historical hydrographic station data in the northern South China Sea and the Philippine Sea. Water masses at 150, 200, and 250 m were classified by discriminant analysis according to their temperature-salinity characteristics. At each depth, most water in the study region was classified into two groups representing the Philippine Sea Water and the South China Sea Water, respectively. The geographic distribution of water masses in the South China Sea shows that the Philippine Sea Water was present along the continental margin south of China between October and January. A westward current in the northern South China Sea in winter was inferred from the distribution of the intrusion water.

Description: Multichannel seismic reflection data were used to determine the evolutionary history of the forearc region of the central Aleutian Ridge. Since at least late Miocene time this sector of the ridge has been obliquely underthrust 30° west of orthogonal convergence by the northwestward converging Pacific plate at a rate of 80–90 km/m.y. Our data indicate that prior to late Eocene time the forearc region was composed of rocks of the arc massif thinly mantled by slope deposits; the forearc region probably lacked both major depositional basins and a tectonically attached accretionary prism of offscraped oceanic deposits. Beginning in latest Miocene or earliest Pliocene time, a zone of outer-arc structural highs and a forearc basin began to form. Formation of these companion intraarc structures may be linked to the late Neogene growth of an accretionary wedge that formed as the result of the deposition of a thick turbidite wedge in the Aleutian Trench. Initial structures of the zone of outer-arc highs formed as the thickening wedge underran, compressively deformed, and uplifted the seaward edge of the arc massif above a landward dipping backstop thrust. Forearc basin strata ponded arcward of the elevating zone of outer-arc highs. However, most younger structures of the zone of outer-arc highs cannot be ascribed simply to the orthogonal effects of an underrunning wedge. Oblique convergence created a major right-lateral shear zone (the Hawley Ridge shear zone) that longitudinally disrupted the zone of outer-arc highs, truncating the seaward flank of the forearc basin and shearing the southern limb of Hawley Ridge, an exceptionally large antiformal outer-arc high structure. Slivers of forearc basement rocks and overlying strata have been transported along the shear zone that is flanked by differentially elevated structures attributed to localized transpressive and transtensional processes. Uplift of Hawley Ridge may be related to the thickening of the arc massif by westward directed basement duplexes. In addition, the forearc is disrupted by structures transverse to the margin that occur where unusually high-stress accumulations have resulted in the rupture of repeated great earthquakes. It is likely that many ancient active margins evolved in tectonic and depositional settings similar to those of the central Aleutian Ridge. Great structural complexity, including the close juxtaposition of coeval structures recording compression, extension, differential vertical movements, and strike-slip displacement, should be expected, even within areas of generally kindred tectonostratigraphic terranes.

Description: Several geochemical anomalies were observed before the Haichen, Longling, Tangshan, and Songpan earthquakes and their strong aftershocks. They included changes in groundwater radon levels; chemical composition of the groundwater (concentration of Ca++, Mg++, Cl−, SO4= and HCO3− ions); conductivity; and dissolved gases such as H2, CO2, etc. In addition, anomalous changes in water color and quality were observed before these large earthquakes. Before some events gases escaped from the surface, and there were reports of "ground odors" being smelled by local residents. The large amount of radon data can be grouped into long-term and short-term anomalies. The long-term anomalies have a radon emission build up time of from a few months to more than a year. The short-term anomalies have durations from a few hours or less to a few months.

Description: Until reliable procedures have been developed to preserve the phosphorus contained in particulate matter captured by in situ pumps and sediment traps and until these procedures are applied over a wide range of locations and depths in the sea, indirect methods will have to be used to determine the C/P ratio in marine detritus. We have taken two such approaches: (1) the use of C/N ratios for particulates captured in the upper thermocline in conjunction with 02/P and N/P ratios obtained from deconvolutions of ocean chemical data and (2) regression along isopycnals in the deep‐sea waters free of fossil fuel CO2. While neither approach yields a definitive answer, both suggest that a value of 127 carbon atoms per phosphorus atom would be a more appropriate interim value than that of 106 adopted long ago by A. C. Redfield and his associates.

Description: Axial volcanic ridges (AVRs) are found on most slow-spreading mid-ocean ridges and are thought to be the main locus of volcanism there. In this study we present high-resolution mapping of a typical, well-defined AVR on the Mid-Atlantic Ridge at 45°N. The AVR is characterized by “hummocky terrain,” composed typically of hummocks with pillowed or elongate pillowed flanks with pillowed or lobate lava flow summits, often with small haystacks sitting on their highest points. The AVR is surrounded by several areas of “flat seafloor,” composed of lobate and sheet lava flows. The spatial and morphological differences between these areas indicate different eruption processes operating on and off the AVR. Volcanic fissures are found all around and on the AVR, although those with the greatest horizontal displacement are found on the ridge crest and flat seafloor. Clusters of fissures may represent volcanic vents. Extremely detailed comparisons of sediment coverage and examination of contact relations around the AVR suggest that many of the areas of flat seafloor are of a similar age or younger than the hummocky terrain of the AVR. Additionally, all the lavas surveyed have similar degrees of sediment cover, suggesting that the AVR was either built or resurfaced in the same 50 ka time frame as the flat seafloor.

Description: The Norwegian Atlantic Current (NwAC) and its eddy field are examined using data from surface drifters. The data set used spans nearly 20 years, from June 1991 to December 2009. The results are largely consistent with previous estimates, which were based on data from the first decade only. With our new data set, statistical analysis of the mean fields can be calculated with larger confidence. The two branches of the NwAC, one over the continental slope and a second further offshore, are clearly captured. The Norwegian Coastal Current is also resolved. In addition, we observe a semipermanent anticylonic eddy in the Lofoten Basin, a feature seen previously in hydrography and in models. The eddy kinetic energy (EKE) is intensified along the path of the NwAC, with the largest values occurring in the Lofoten Basin. The strongest currents, exceeding 100 cm s−1, occur west of Lofoten. Lateral diffusivities were computed in five domains and ranged from 1–5 × 107 cm2 s−1. The Lagrangian integral time and space scales are 1–2 days and 7–23 km, respectively. The data set allows studies of seasonal and interannual variations as well. The strongest seasonal signal is in the NwAC itself, as the mean flow strengthens by approximately 20% in winter. The EKE and diffusivities on the other hand do not exhibit consistent seasonality in the sampled regions. There are no consistent indications of changes in either the mean or fluctuating surface velocities between the 1990s and 2000s.

Description: Samples collected at hourly intervals on May 18–19, 1980, at three sites 200 km downwind from Mount St. Helens, have made possible a detailed reconstruction of the conditions that contribute to the compositional heterogeneity of mineral and glass components observed in distal tephra layers. The air fall tephra deposited at the sites during the first 7 hours of the May 18 eruption is mostly coarse grained, microlite-rich, nonjuvenile glass and feldspar. Grain-size maxima in this initial tephra can be related to the cataclysmic blast at 0832 and a subsequent pulse of the eruption at 1200. Juvenile, microlite-free glass increases in relative abundance at the sampling sites beginning at about 1900. Such a change between nonjuvenile and juvenile tephra can be related to a 5-km increase in column height associated with the last major pulse of the eruption which occurred at 1700 at the volcano. Electron microprobe study of both microlite-rich and microlite-free pumice in the time series samples reveals significant compositional differences. Interstitial glass in nonjuvenile pumice deposited during the first few hours at the sampling sites is enriched in SiO2 and K2O and depleted in TiO2, FeO*, and MgO relative to juvenile glass. By comparison, major element composition of the least evolved juvenile glass sampled during the last several hours of the eruption displays a slight trend toward less evolved composition. Least squares calculations suggest that the more evolved character of the nonjuvenile glass can be explained by greater fractional crystallization brought about by enhanced cooling in a cryptodome prior to eruption, whereas the temporal changes observed in juvenile glass composition during the last several hours of the eruption suggest the presence of a small, slightly zoned magma chamber at depth. Electron microprobe study of glass-coated ilmenites, magnetites, and plagioclases provides the following estimates of the physical conditions in this reservoir: 865°±50°C, PH2O = 2.2 kbar and -log ƒO2 = 11.7. Analyses of bulk pumice, glass and selected mineral phases from May 25, June 12, July 22, and October 16–18 pumices erupted from Mount St. Helens indicate that the bulk pumice (magma) compositions have become slightly more andesitic with time, while mineral and co-existing glass compositions have changed significantly in post-May 18 eruptions with both being more highly evolved than those associated with the May 18 eruption. An application of the magnetite-ilmenite geothermometer to June 12 and July 22 samples indicates temperatures of 919°±30°C and 930°±50°C, respectively. Least squares calculations suggest that such evolved post-May 18 glass and mineral phases can be derived by fractional crystallization of a magma composition like bulk May 18 pumice into approximately 50% crystals and 50% residual liquid. Such partitioning between crystals and residual liquid appears to have occurred on the scale of centimeters and is interpreted as a consequence of accelerated crystallization under reduced water pressure.

Description: The well-known caldera of Thira (Santorini), Greece, was not formed during a single eruption but is composed of two overlapping calderas superimposed upon a complex volcanic field that developed along a NE trending line of vents. Before the Minoan eruption of 1400 B.C., Thira consisted of three Java shields in the northern half of the island and a flooded depression surrounded by tuff deposits in the southern half. Andesitic lavas formed the overlapping shields of the north and were contemporaneous with and, in many places, interbedded with the southern tuff deposits. Although there appears to be little difference between the composition of magmas erupted, differences in eruption style indicate that most of the activity in the northern half of the volcanic field was subaerial, producing lava flows, whereas in the south, eruptions within a flooded depression produced a sequence of mostly phreatomagmatic tuffs. Many of these tuffs are plastered onto the walls of what appears to have been an older caldera, most probably associated with an eruption of rhyodacitic tephra 100,000 years ago. The Minoan eruption of about 1400 B.C. had four distinct phases, each reflecting a different vent geometry and eruption mechanism. The Minoan activity was preceded by minor eruptions of fine ash. (1) The eruption began with a Plinian phase, from subaerial vent(s) located on the easternmost of the lava shields. (2) Vent(s) grew toward the SW into the flooded depression. Subsequent activity deposited large-scale base surge deposits during vent widening by phreatomagmatic activity. (3) The third eruptive phase was also phreatomagmatic and produced 60% of the volume of the Minoan Tuff. This activity was nearly continuous and formed a large featureless tuff ring with poorly defined bedding. This deposit contains 5–40% lithic fragments that are typical of the westernmost lava shield and appears to have been erupted when caldera collapse began. (4) The last phase consisted of eruption of ignimbrites from vent(s) on the eastern shield, not yet involved in collapse. Collapse continued after eruption of the ignimbrites with foundering of the eastern half of the caldera. Total volume of the collapse was about 19 km3, overlapping the older caldera to form the caldera complex visible today. Intracaldera eruptions have formed the Kameni Islands along linear vents concomitant with vents that may have been sources for the Minoan Tuff.

Description: Seismic techniques provide the highest-resolution measurements of the structure of the crust and have been conducted on a worldwide basis. We summarize the structure of the continental crust based on the results of seismic refraction profiles and infer crustal composition as a function of depth by comparing these results with high-pressure laboratory measurements of seismic velocity for a wide range of rocks that are commonly found in the crust. The thickness and velocity structure of the crust are well correlated with tectonic province, with extended crust showing an average thickness of 30.5 km and orogens an average of 46.3 km. Shields and platforms have an average crustal thickness nearly equal to the global average. We have corrected for the nonuniform geographical distribution of seismic refraction profiles by estimating the global area of each major crustal type. The weighted average crustal thickness based on these values is 41.1 km. This value is 10% to 20% greater than previous estimates which underrepresented shields, platforms, and orogens. The average compressional wave velocity of the crust is 6.45 km/s, and the average velocity of the uppermost mantle (Pn velocity) is 8.09 km/s. We summarize the velocity structure of the crust at 5-km depth intervals, both in the form of histograms and as an average velocity-depth curve, and compare these determinations with new measurements of compressional wave velocities and densities of over 3000 igneous and metamorphic rock cores made to confining pressures of 1 GPa. On the basis of petrographic studies and chemical analyses, the rocks have been classified into 29 groups. Average velocities, densities, and standard deviations are presented for each group at 5-km depth intervals to crustal depths of 50 km along three different geotherms. This allows us to develop a model for the composition of the continental crust. Velocities in the upper continental crust are matched by velocities of a large number of lithologies, including many low-grade metamorphic rocks and relatively silicic gneisses of amphibolite facies grade. In midcrustal regions, velocity gradients appear to originate from an increase in metamorphic grade, as well as a decrease in silica content. Tonalitic gneiss, granitic gneiss, and amphibolite are abundant midcrustal lithologies. Anisotropy due to preferred mineral orientation is likely to be significant in upper and midcrustal regions. The bulk of the lower continental crust is chemically equivalent to gabbro, with velocities in agreement with laboratory measurements of mafic granulite. Garnet becomes increasingly abundant with depth, and mafic garnet granulite is the dominant rock type immediately above the Mohorovicic discontinuity. Average compressional wave velocities of common crustal rock types show excellent correlations with density. The mean crustal density calculated from our model is 2830 kg/m3, and the average SiO2 content is 61.8%.

Description: We present a conceptual model of fluid circulation in a ridge flank hydrothermal system, the Mariana Mounds. The model is based on chemical data from pore waters extracted from piston cores and from push cores collected by deep-sea research vessel Alvin in small, meter-sized mounds situated on a local topographic high. These mounds are located within a region of heat flow exceeding that calculated from a conductive model and are zones of strong pore water upflow. We have interpreted the chemical data with time-dependent transport-reaction models to estimate pore water velocities. In the mounds themselves pore water velocities reach several meters per year to kilometers per year. Within about 100 m from these zones of focused upflow velocities decrease to several centimeters per year up to tens of centimeters per year. A larger area of low heat flow surrounds these heat flow and topographic highs, with upwelling pore water velocities less than 2 cm/yr. In some nearby cores, downwelling of bottom seawater is evident but at speeds less than 2 cm/yr. Downwelling through the sediments appears to be a minor source of seawater recharge to the basaltic basement. We conclude that the principal source of seawater recharge to basement is where basement outcrops exist, most likely a scarp about 2–4 km to the east and southeast of the study area.

Description: During a 10-year study more than 2,000 phytoplankton samples were collected from the entire coast of Peru and analyzed. In general, diatoms were the most abundant group of organisms in all seasons. Predominant species were Rhizosolenia delicatula, Skeletonema costatum Thalassiosira subtilis, Thalassionema nitzschioides and several species of the genus Chaetooeros. Dinoflagellates and flagellates were observed frequently during summer. The mean distribution of the phytoplankton concentration during the 10 years shows the existence of several centers with higher cell densities along the coast, coinciding with the areas of more intense and persistent upwelling. Four major centers have been identified: Pimentel (˜6°S), Chimbote (˜9°S), Callao (˜12°S), and Tambo de Mora-Pisco (˜15°S); and two minor centers, Talara (˜4°S) and Ilo (˜17°S). The relative importance of each center seems to change according to the season. The highest phytoplankton concentration tended to be in the northern part of the coast during fall and winter and in the south through spring and summer.

Description: Methane is a potent greenhouse gas and large-scale rapid release of methane from hydrate may have contributed to past abrupt climate change inferred from the geological record. The discovery in 2008 of over 250 plumes of methane gas escaping from the seabed of the West Svalbard continental margin at ~400 m water depth (mwd) suggests that hydrate is dissociating in the present-day Arctic. Here we model the dynamic response of hydrate-bearing sediments over a period of 2300 years and investigate ocean warming as a possible cause for present-day and likely future dissociation of hydrate, within 350–800 mwd, west of Svalbard. Future temperatures are given by two climate models, HadGEM2 and CCSM4, and scenarios, Representative Concentration Pathways (RCPs) 8.5 and 2.6. Our results suggest that over the next three centuries 5.3–29 Gg yr−1 of methane may be released to the Arctic Ocean on the West Svalbard margin.

Description: Aftershocks of shallow earthquakes larger than magnitude 7 in the Aleutians, southern Alaska, southeast Alaska, and offshore British Columbia from 1920 to 1970 were relocated by computer in an attempt to delineate the rupture zones of large earthquakes. Plate tectonic theory indicates that gaps in activity for large earthquakes for the past 10's to 100's of years are likely sites of future large earthquakes. Three prominent gaps of this type are delineated: one in southeast Alaska; another in southern Alaska near the epicenters of the great earthquakes of 1899 and 1900; and one in the far western Aleutians. These gaps deserve high priority for study and instrumentation. Large earthquakes appear to be much more regular than smaller shocks in their distributions with respect to space, time, and size. Aftershock zones of events since 1930 that are larger than magnitude 7.8 are longer than 250 km and those less than 7.5 are shorter than 125 km. The rupture zones of events that occurred before 1930 could not be delineated from aftershock locations. Aftershock zones of large earthquakes tend to abut without significant overlap even for rupture zones as long as 1200 km. Nearly the entire Alaska-Aleutian zone from 145°W to 171°E has broken since 1938 in a series of large earthquakes. The rupture zones of five large events appear to form a space-time sequence that progressed from 155°W in 1938 to 171°E in 1965. This sequence is much like the well-known westward progression of activity since 1939 along the North Anatolian fault. Shocks with long rupture zones tend to occur along those parts of the Alaska-Aleutian zone that are relatively simple tectonically. The ends of many aftershock zones of large earthquakes are located at the intersection of major transverse features with the Aleutian arc. Large earthquakes rarely, if ever, reoccur along the same part of a fault zone in less than several tens of years, i.e. within a time less than that for substantial strain accumulation. Events of comparable magnitude that occur soon after some great earthquakes usually involve rupture in a region adjacent to but different from that of the main shock. The March 30, 1965, earthquake of magnitude 7.5, which involved normal faulting in the Aleutian trench, appears to have been triggered by thrust faulting along the adjacent inner margin of the trench in the magnitude 7.9 earthquake of February 4, 1965. Large events of the thrust type are commonly followed within ten years by events involving normal faulting in the adjacent part of the trench. Estimates of average displacements and of the repeat times of great earthquakes from measurements of 20-sec surface waves are systematically too small and do not agree with the meager historic record of great shocks. Other estimates of repeat times vary from 30 to 850 years, but neither of these extremes appears to be typical. The aftershock zone of the April 1, 1946, Aleutian earthquake, which generated one of the largest and most widespread seismic sea waves in the Pacific during this century, was very small. A large displacement of the ocean floor may be responsible for the generation of the large sea wave. An average displacement of 2.4 to 4.1 meters was calculated from amplitudes of 100-sec waves.

Description: We measured the respiratory isotope effect ϵresp for seven representative unicellular marine organisms. The bacterium Pseudomonas halodurans, the diatom Phaeodactylum tricornutum, the phytoflagellates Cryptomonas baltica and Dunaliella tertiolecta, the heterotrophic flagellates Paraphysomonas imperforata and Bodo sp., and the ciliate Uronema sp. exhibit ϵresp values in the range 14-26‰. We also measured ϵresp for three metazoans. The ϵresp for the copepod Acartia tonsa ranged from 17 to 25‰, while two larger organisms, the mollusk Mercenaria mercenaria and the salmon Salmo salmar, respire with a smaller ϵresp of 5-10‰. The average respiratory isotope effect of the dominant marine respirers (the bacteria, microalgae and zooplankton) is about 20 ± 3‰. An ϵresp of this magnitude supports the hypothesis that the photosynthesis-respiration cycle is responsible for the 23.5‰ enrichment in the δ18O ratio of atmospheric O2 relative to seawater (the Dole effect). The large value and high variability in the average ϵresp limits the usefulness of a proposed method using the δ18O of naturally fractionated dissolved O2 in seawater as a tracer of primary production in the oligotrophic ocean.

Description: We review the current understanding of the Dole effect (the observed difference between the δ18O of atmospheric O2 and that of seawater) and its causes, extend the record of variations in the Dole effect back to 130 kyr before present using data on the δ18O of O2 obtained from studying the Vostok ice core (Sowers et al., 1993), and discuss the significance of temporal variations. The Dole effect reflects oxygen isotope fractionation during photosynthesis, respiration, and hydrologic processes (evaporation, precipitation, and evapotranspiration). Our best prediction of the present-day Dole effect, +20.8‰, is considerably lower than the observed value, +23.5‰, and we discuss possible causes of this discrepancy. During the past 130 kyr, the Dole effect has been 0.05‰ lower than the present value, on average. The standard deviation of the Dole effect from the mean has been only ±0.2‰, and the Dole effect is nearly unchanged between glacial maxima and interglacial periods. The small variability in the Dole effect suggests that relative rates of primary production in the land and marine realms have been relatively constant. Most periodic variability in the Dole effect is in the precession band, suggesting that changes in this global biogeochemical term reflects variations in low-latitude land hydrology and productivity or possibly variability in low-latitude oceanic productivity.

Description: Much attention has been paid, in recent years, to the potential application of the Ce anomaly, measured in various marine phases, as a paleoceanographic indicator of widespread marine anoxia. In this paper we present and discuss results from recent studies of present‐day rare earth element (REE) distributions (and hence Ce anomaly distributions) in the marine environment which are particularly pertinent to paleoceanography. Subsequently, we review and discuss the validity of the recent literature in which Ce anomalies have been employed as paleoredox indicators.

Description: On Syros, high‐pressure metamorphism affects a lithological pile that is composed of, from base to top: (1) the Komito‐Vari granitic basement, (2) a margin sedimentary sequence that is predominantly made of marbles and schists (the Pyrgos and Kastri units), and (3) the Kambos metaophiolitic mélange. The tectonic history occurred in three main stages. During the first stage, in the mid‐Eocene, the Kambos oceanic unit was thrust southward on top of the sedimentary pile. Top‐to‐the‐south‐southwest ductile senses of shear are synchronous with prograde high‐pressure metamorphism and associated with this thrusting event. The second stage corresponds to a top‐to‐the‐northeast ductile shear that affects the whole metamorphic pile and is synchronous with the metamorphic retrogression from eclogite to greenschist facies. However, the Kambos oceanic unit remained partly undeformed, as shown by significant volumes containing undeformed lawsonite pseudomorphs. No major extensional detachment related to this exhumation event outcrops on the island. The localized semibrittle to brittle deformation of the third stage is associated with the postmetamorphic development of (1) a ramp‐flat extensional system at the island scale, whose southward minimum displacement is estimated at approximately 7 km, and (2) two sets of steeply dipping strike‐slip faults with a normal component, trending either east–west or around north–south, indicating that the mean stretching and shortening directions are trending NNE–SSW and ESE–WNW, respectively. This sequence of major tectonic events and their relationship to metamorphism are interpreted within the framework of the subduction of the Pindos Ocean and then of the Adria continental passive margin.

Description: Abundances of 12 species of planktonic foraminifera collected in two plankton tows from the east tropical Atlantic are compared to the chlorophyll content and the temperature of the sea water from which they were collected. As expected from previous work in the tropics, all dominant tropical species occur in greatest abundance within the photic zone. Many species occur in greatest abundances in the seasonal thermocline in association with the maximum chlorophyll concentration, while a few algal symbiont-bearing species occur in greatest abundance in the mixed layer. The δ18O measurements of planktonic foraminifera shells from core top sediment samples confirm the vertical stratification within the photic zone that is suggested by the relationship between hydrography and abundances found in the plankton tows and found in the statistical study by Ravelo et al. [1990]. Comparison between the measured δ18O values of planktonic foraminifera with the predicted δ18O profiles of the overlying water column at three core locations indicate that species abundances in the sediment record the seasonally integrated conditions of the photic zone and suggests that the abundance of a species in the sediment depends on whether the preferred ecological conditions of that species may be found within the photic zone of the overlying water column sometime during the year. Species which calcify below the photic zone have only trace relative abundances. Finally, it appears that the total range of δ18O values of the dominant species approximates the predicted annual δ18O of calcite range in the upper 80 m of the water column.

Description: Pore water has been analyzed from sediment cores taken from three areas on the eastern flank of the Juan de Fuca Ridge as part of FlankFlux 90, a study of hydrothermal circulation through mid-ocean ridge flanks. Seismic reflection and heat flow surveys (Davis et al., 1992a) indicate that the three areas differ in sediment thickness, basement topography, abundance of outcrops, basement temperature, and fraction of heat lost by advection versus conduction. Area 1 is on 0.6 Ma crust with nearly continuous basement outcrop, area 2 is on 1.3 Ma crust over the first buried ridge parallel to the present ridge axis, and area 3 is on 3.5–3.8 Ma crust over two axis-parallel buried ridges that penetrate the sediment cover in three locations. Each area includes a hydrothermal system in which seawater flows into basement, reacts with crustal basalt, and then exits basement either through the sediment or directly into the overlying water column. As constrained by concentrations of sulfate and lithium in the pore waters, at least some seawater enters basement in all three areas without reacting fully with the overlying sediment, even where no outcrops are known nearby. Speeds of up welling of pore water through the sediment have been estimated by fitting profiles of dissolved magnesium and chlorinity, which behave conservatively in these areas, to numerical time-dependent transport models. The estimated velocities range from 〈0.1 to 7.4 cm/yr; faster flows probably occur but were not sampled. Upwelling speed correlates positively with heat flow and basement highs and negatively with sediment thickness. The correlation with heat flow differs from area 2 to area 3 along with differences in physical properties of the turbidite sediment. We have documented pore water upwelling through sediment up to 100 m thick. We estimate that upwelling continues at decreasing speeds through sediment up to 160 m thick, corresponding to a heat flow of 0.44 W/m2 in area 2 and 0.3 W/m2 in area 3. Concentrations of magnesium and chlorinity in the altered seawater upwelling from basement are uniform within each area but differ from one area to the next. Both species remain at the bottom seawater concentration in area 1, where basement is cooled to 〈10°C at the base of the sediments mainly by advection. The concentration of magnesium decreases with increasing basement temperature in areas 2 and 3 to a minimum of 2.5 mmol/kg at about 90°C in area 3. The transition from largely advective to largely conductive heat loss occurs over only 20 km between areas 1 and 2 and corresponds to a dramatic change in the composition of fluid circulating through basement, as the uppermost basement is heated from 〈10° to 40–50°C. Chlorinity of the basement fluid increases above the present-day bottom seawater concentration in areas 2 and 3 and in nearly all other mid-ocean ridge flanks studied to date, as a result of rock hydration and the higher chlorinity of bottom seawater during the last glacial period. While chlorinity generally correlates positively with uppermost basement temperature in various ridge flank hydrothermal systems, it reaches a maximum in area 2 at only 40°C, probably because alteration there occurs at a lower water/rock ratio than elsewhere. For all mid-ocean ridge flanks studied to date, the temperature at the basement interface correlates better with the fraction of heat lost by advection versus conduction and with the average thickness of the sediment cover than with crustal age.

Description: An examination of 311 intraplate earthquakes in the Australian plate portion of the Pacific Ocean basin reported from 1918 to 1990 reveals that only 113 events are reliably intraplate, with most of the rest relocating to active trenches and transforms. The non-random distribution of the reliably intraplate events gives insight into the tectonic stresses present. The central Tasman Sea is mostly aseismic except for a swarm of activity at the predicted site of the Tasmantid hot spot. To the north, the broad regions of the Coral Sea, South Fiji Basin and Lord Howe Rise show very little intraplate seismicity, yet the narrow Norfolk Ridge and Three Kings Rise, caught between the double convergence of the New Hebrides and Tonga subduction zones, support many more earthquakes. High levels of intraplate seismicity in the southern Tasman Sea adjacent to the Macquarie Ridge Complex (MRC) indicate that this region may be undergoing internal deformation due to the unusual nature of the Australia-Pacific plate boundary. Additional support exists in the form of intraplate focal mechanisms similar to those at the plate boundary and a set of parallel gravity rolls which are observed in recent geoid maps. Some aftershocks of the Mw = 8.2 Macquarie Ridge earthquake of 1989 occurred in a fracture zone west of the Macquarie Ridge Complex [Das, 1992], but we have found several earthquakes from as early as 1924 which relocate to this feature, suggesting that its reactivation may be more significant than previously thought. This reactivation of a fossil fracture zone may be the result of the increasing amount of oblique convergence between the Australia and Pacific plates at the Macquarie Ridge Complex, formerly a spreading center, and the stresses associated with subducting recently formed Australian ocean crust beneath the older Pacific plate.

Description: The eastern Mediterranean transient (EMT) was caused by a combination of high‐salinity waters intruding into the Aegean Sea and the two particularly strong winters of 1991–1992 and 1992–1993. The approach in this chapter is to search for specific signatures in the historic hydrographic observations, which date back to 1910. To deal with the problem that up into the 1950s the data not only are of limited precision but also have gaps of about 20 years, it is advantageous to consider the fact that the evolution of the actual EMT is rather well documented over a similar time span. The chapter begins by outlining the characteristics of the current EMT. Thereafter, a selection of suitable hydrographic observations among the available historic data is provided to compare these with signatures expected from the evolution of the actual EMT.

Description: The Mauritanian coastal area is one of the most biologically productive upwelling regions in the world ocean. Shipboard observations carried out during maximum upwelling season and short-term moored observations are used to investigate diapycnal mixing processes and to quantify diapycnal fluxes of nutrients. The observations indicate strong tide-topography interactions that are favored by near-critical angles occurring on large parts of the continental slope. Moored velocity observations reveal the existence of highly nonlinear internal waves and bores and levels of internal wave spectra are strongly elevated near the buoyancy frequency. Dissipation rates of turbulent kinetic energy at the slope and shelf determined from microstructure measurements in the upper 200 m averages to ɛ = 5 × 10−8 W kg−1. Particularly elevated dissipation rates were found at the continental slope close to the shelf break, being enhanced by a factor of 100 to 1000 compared to dissipation rates farther offshore. Vertically integrated dissipation rates per unit volume are strongest at the upper continental slope reaching values of up to 30 mW m−2. A comparison of fine-scale parameterizations of turbulent dissipation rates for shelf regions and the open ocean to the measured dissipation rates indicates deficiencies in reproducing the observations. Diapycnal nitrate fluxes above the continental slope at the base of the mixed layer yielding a mean value of 12 × 10−2 μmol m−2 s−1 are amongst the largest published to date. However, they seem to only represent a minor contribution (10% to 25%) to the net community production in the upwelling region.

Description: Water transported by slabs into the mantle at subduction zones plays key roles in tectonics, magmatism, fluid and volatiles fluxes, and most likely in the chemical evolution of the Earth's oceans and mantle. Yet, incorporation of water into oceanic plates before subduction is a poorly understood process. Several studies suggest that plates may acquire most water at subduction trenches because the ocean crust and uppermost mantle there are intensely faulted caused by bending and/or slab pull, and display anomalously low seismic velocities. The low velocities are interpreted to arise from a combination of fluid-filled fractures associated to normal faulting and mineral transformation by hydration. Mantle hydration by transformation of nominally dry peridotite to water-rich serpentinite could potentially create the largest fluid reservoir in slabs and is therefore the most relevant for the transport of water in the deep mantle. The depth of fracturing by normal-fault earthquakes is usually not well constrained, but could potentially create deep percolation paths for water that might hydrate up to tens of kilometers into the mantle, restrained only by serpentine stability. Yet, interpretation of deep intraplate mineral alteration remains speculative because active-source seismic experiments have sampled only the uppermost few kilometers of mantle, leaving the depth-extent of anomalous velocities and their relation to faulting unconstrained. Here we use a joint inversion of active-source seismic data, and both local and regional earthquakes to map the three dimensional distribution of anomalous velocities under a seismic network deployed at the trench seafloor. We found that anomalous velocities are restrained to the depth of normal-fault micro-earthquake activity recorded in the network, and are considerably shallower than either the rupture depth of teleseismic, normal-fault earthquakes, or the limit of serpentine stability. Extensional micro-earthquakes indicate that each fault in the region slips every 2–3 months which may facilitate regular water percolation. Deeper, teleseismic earthquakes are comparatively infrequent, and possibly do not cause significant fracturing that remains open long enough to promote alteration detectable with our seismic study. Our results show that the stability field of serpentine does not constrain the depth of potential mantle hydration.

Description: Understanding the pre-anthropogenic Pb cycle of central North Pacific deep water has attracted a lot of attention in recent years, partly because of its unique geographical location in that it is a remote gyre system characterized by high dust fluxes and sluggish overturning circulation. However, the factors controlling Pb isotope evolution in this area over the Cenozoic are still controversial and various mechanisms have been proposed in previous studies. Here we report new Pb and Nd isotope time series of four ferromanganese crusts (two from the western Pacific near the Mariana arc and the other two from the central Pacific). Together with previously published records, we discuss for the first time the significance of a persistent and systematic Pb isotopic provinciality recorded by central North Pacific crusts over the Cenozoic. We propose that globally well mixed stratosphere volcanic aerosols could contribute Pb but have not been the major factors controlling the Pb isotope distribution in the central North Pacific over time. Island arc input (and probably enhanced hydrothermal input between about 45 and 20 Ma) likely controlled the Pb isotope provinciality and evolution prior to ~20 Ma, when coeval Pb isotope records in different crusts showed large differences and atmospheric silicate dust flux was extremely low. After the Eocene, in particular after 20 Ma, Asian dust input has become an isotopically resolvable source, while island arc-derived Pb has remained important to balance the dust input and to produce the observed Pb isotope distribution in the central North Pacific during this period.

Description: A realistic primitive-equation model of the Southern Ocean at eddying spatial resolution is used to examine the effect of ocean-surface-velocity dependence of the wind stress on the strength of near-inertial oscillations. Accounting for the ocean-surface-velocity dependence of the wind stress leads to a large reduction of wind-induced near-inertial energy of approximately 40 percent and of wind power input into the near-inertial frequency band of approximately 20 percent. A large part of this reduction can be explained by the leading-order modification to the wind stress if the ocean-surface velocity is included. The strength of the reduction is shown to be modulated by the inverse of the ocean-surface-mixed-layer depth. We conclude that the effect of surface-velocity dependence of the wind stress should be taken into account when estimating the wind-power input into the near-inertial frequency band and when estimating near-inertial energy levels in the ocean due to wind forcing.

Description: This study presents a new estimate of the oceanic anthropogenic CO2(Cant) sink over the industrial era (1780 to present), from assimilation of potential temperature, salinity, radiocarbon, and CFC-11 observations in a global steady state ocean circulation inverse model (OCIM). This study differs from previous data-based estimates of the oceanic Cant sink in that dynamical constraints on ocean circulation are accounted for, and the ocean circulation is explicitly modeled, allowing the calculation of oceanic Cant storage, air-sea fluxes, and transports in a consistent manner. The resulting uncertainty of the OCIM-estimated Cant uptake, transport, and storage is significantly smaller than the comparable uncertainty from purely data-based or model-based estimates. The OCIM-estimated oceanic Cant storage is 160–166 PgC in 2012, and the oceanic Cant uptake rate averaged over the period 2000–2010 is 2.6 PgC yr−1 or about 30% of current anthropogenic CO2 emissions. This result implies a residual (primarily terrestrial) Cant sink of about 1.6 PgC yr−1 for the same period. The Southern Ocean is the primary conduit for Cant entering the ocean, taking up about 1.1 PgC yr−1 in 2012, which represents about 40% of the contemporary oceanic Cant uptake. It is suggested that the most significant source of remaining uncertainty in the oceanic Cant sink is due to potential variability in the ocean circulation over the industrial era.

Description: Sedimentation rates (corrected for compaction) from along the passive continental margin of Africa between the Equatorial Fracture Zone and Somalia are used to compare the rates of subsidence of the continental crust since early Mesozoic time. Three distinctive subsidence histories can be identified which correspond with basinal areas that have different structural styles: rifted (west coast), sheared (Equatorial and Agulhas fracture zones) and sunk (zones of vertical tectonics in eastern Africa). A comparison of subsidence rates with other tensional margins (NE USA and the North Sea) and a consideration of the plate tectonic history of the African margins leads to the proposal of a geo and thermodynamic model that takes cognizance of the worldwide mid-Cretaceous rheological discontinuity between taphrogenic and epeirogenic basin formation recognized by Kent, and the more generally accepted, purely plate tectonic driven model of margin subsidence. The new suggestion involves a lower Mesozoic worldwide rise in the geothermal gradient in the lithosphere which produces metamorphism of the base of the continental crust and initiates taphrogenesis along lineaments throughout Gondwanaland. A lowering of the geothermal gradient in the lower Cretaceous produces a switch to epeirogenic subsidence, driven solely by sediment loading and thermal contraction, by Aptian/Albian times. The thermal event facilitated continental separation, and sea floor spreading commenced locally at various times along the active taphrogenic belts. Local thermal and tectonic aberrations associated with this phenomenon over print onto the worldwide pattern of marginal basin subsidence. A further rise in the geothermal gradient may have been responsible for renewed taphrogenesis in eastern Africa in Tertiary times.

Description: Stratospheric ozone depletion and emission of greenhouse gases lead to a trend of the southern annular mode (SAM) toward its high-index polarity. The positive phase of the SAM is characterized by stronger than usual westerly winds that induce changes in the physical carbon transport. Changes in the natural carbon budget of the upper 100 m of the Southern Ocean in response to a positive SAM phase are explored with a coupled ecosystem-general circulation model and regression analysis. Previously overlooked processes that are important for the upper ocean carbon budget during a positive SAM period are identified, namely, export production and downward transport of carbon north of the polar front (PF) as large as the upwelling in the south. The limiting micronutrient iron is brought into the surface layer by upwelling and stimulates phytoplankton growth and export production but only in summer. This leads to a drawdown of carbon and less summertime outgassing (or more uptake) of natural CO2. In winter, biological mechanisms are inactive, and the surface ocean equilibrates with the atmosphere by releasing CO2. In the annual mean, the upper ocean region south of the PF loses more carbon by additional export production than by the release of CO2 into the atmosphere, highlighting the role of the biological carbon pump in response to a positive SAM event.

Description: We combined data sets of measured sedimentary calcium carbonate (CaCO3) and satellite-derived pelagic primary production to parameterize the relation between CaCO3 content on the Antarctic shelves and primary production in the overlying water column. CaCO3 content predicted in this way was in good agreement with the measured data. The parameterization was then used to chart CaCO3 content on the Antarctic shelves all around the Antarctic, using the satellite-derived primary production. The total inventory of CaCO3 in the bioturbated layer of Antarctic shelf sediments was estimated to be 0.5 Pg C. This quantity is comparable to the total CO2 uptake by the Southern Ocean in only one to a few years (dependent on the uptake estimate and area considered), indicating that the dissolution of these carbonates will neither delay ocean acidification in this area nor augment the Southern Ocean CO2 uptake capacity.

Description: The evolution of ridge-hotspot systems is not well understood. In this investigation, satellite-derived marine gravity data are used in conjunction with underway bathymetric and magnetic anomaly profiles to investigate the nature of ridge-hotspot interaction at four sparsely explored systems in the Southern Ocean. These systems illustrate three different stages of ridge-hotspot interaction in which a migrating spreading center approaches a hotspot (Pacific-Antarctic/Louisville), passes over or is captured by the hotspot (Mid-Atlantic/Shona-Discovery), and ultimately migrates away from the hotspot (Southeast Indian/Kerguelen). All of these systems show some evidence of discrete ridge jumps in the direction of the hotspot as the spreading center attempts to relocate toward the hotspot by asymmetric spreading. Interestingly, these ridge jumps show no evidence of propagating offsets as have been seen on many other ridge-hotspot systems. A simple model predicts that typical plume excess temperatures can weaken the lithosphere sufficiently to promote asymmetric spreading and possibly allow a discrete ridge jump. The presence of previously uncharted, obliquely oriented aseismic ridges and gravity lineations between the ridge and the hotspot supports the notion of asthenospheric flux from the plume to the spreading center both before and after the time when the hotspot is ridge centered. The azimuths of the aseismic ridges cannot be explained by plate kinematics alone; they consistently extend from the ends toward the centers of the adjacent spreading segments suggesting some interaction between plume derived asthenospheric flux and local lithospheric structure. The features discussed here also indicate that the transfer of asthenospheric material from the plume to the spreading center is influenced by the local plate boundary configuration and interaction with transform offsets.

Description: We have compiled both laboratory and worldwide field data on electrical conductivity to help understand the physical implications of deep crustal electrical profiles. Regional heat flow was used to assign temperatures to each layer in regional electrical conductivity models; we avoided those data where purely conductive heat flow suggested temperatures more than about 1000°C, substantially higher than solidus temperatures and outside the range of validity of heat flow models. The resulting plots of log conductivity σ versus 1/T demonstrate that even low-conductivity layers (LCL) have conductivities several orders of magnitude higher than dry laboratory samples and that the data can be represented by straight line fits. In addition, technically active regions show systematically higher conductivities than do shield areas. Because volatiles are usually lost in laboratory measurements and their absence is a principal difference between laboratory and field conditions, these materials probably account for the relatively higher conductivities of rocks in situ in the crust; free water in amounts of 0.01–0.1% in fracture porosity could explain crustal conductivities. Other possibilities are graphite, hydrated minerals in rare instances, and sulfur in combination with other volatiles. As most of the temperatures are less than 700°C, partial melting seems likely only in regions of highest heat flow where the conductive temperature profiles are inappropriate. Another result is that at a given temperature, crustal high-conductivity layers (HCL) are more conductive by another order of magnitude and show more scatter than do LCL's. Because the differences between HCL's and LCL's are independent of temperature, we must invoke more than temperature increases as a cause for large conductivity increases; increased fluid concentration in situ seems a probable cause for enhanced conductivities in HCL's. From the point of view of these observations, it does not matter whether the fluids are in communication with the surface or trapped at lithostatic pressures.

Description: The combination of the Sunda megathrust and the (strike-slip) Sumatran Fault (SF) represents a type example of slip-partitioning. However, superimposed on the SF are geometrical irregularities that disrupt the local strain field. The largest such feature is in central Sumatra where the SF splits into two fault strands up to 35 km apart. A dense local network was installed along a 350 km section around this bifurcation, registering 1016 crustal events between April 2008 and February 2009. 528 of these events, with magnitudes between 1.1 and 6.0, were located using the double-difference relative location method. These relative hypocentre locations reveal several new features about the crustal structure of the SF. Northwest and southeast of the bifurcation, where the SF has only one fault strand, seismicity is strongly focused below the surface trace, indicating a vertical fault that is seismogenic to ∼15 km depth. By contrast intense seismicity is observed within the bifurcation, displaying streaks in plan and cross-section that indicate a complex system of faults bisecting the bifurcation. In combination with analysis of topography and focal mechanisms, we propose that the bifurcation is a strike-slip duplex system with complex faulting between the two main fault branches.

Description: Correlations between particulate organic carbon (POC) and mineral fluxes in the deep ocean have inspired the inclusion of “ballast effect” parameterizations in carbon cycle models. A recent study demonstrated regional variability in the effect of ballast minerals on the flux of POC in the deep ocean. We have undertaken a similar analysis of shallow export data from the Arctic, Atlantic, and Southern Oceans. Mineral ballasting is of greatest importance in the high-latitude North Atlantic, where 60% of the POC flux is associated with ballast minerals. This fraction drops to around 40% in the Southern Ocean. The remainder of the export flux is not associated with minerals, and this unballasted fraction thus often dominates the export flux. The proportion of mineral-associated POC flux often scales with regional variation in export efficiency (the proportion of primary production that is exported). However, local discrepancies suggest that regional differences in ecology also impact the magnitude of surface export. We propose that POC export will not respond equally across all high-latitude regions to possible future changes in ballast availability.

Description: The role of biominerals in driving carbon export from the surface ocean is unclear. We compiled surface particulate organic carbon (POC), and mineral ballast export fluxes from 55 different locations in the Atlantic and Southern Oceans. Substantial surface POC export accompanied by negligible mineral export was recorded implying that association with mineral phases is not a precondition for organic export to occur. The proportion of non-mineral associated sinking POC ranged from 0 to 80% and was highest in areas previously shown to be dominated by diatoms. This is consistent with previous estimates showing that transfer efficiency in such regions is low. However we propose that, rather than the low transfer efficiency arising from diatom blooms being inherently characterized by poorly packaged aggregates which are efficiently exported but which disintegrate readily in mid water, it is due to such environments having very high levels of unballasted organic C export.

Description: A numerical algorithm based on Fermat's Principle was developed to simulate the propagation of Global Positioning System (GPS) radio signals in the refractivity field of a numerical weather model. The unique in the proposed algorithm is that the ray-trajectory automatically involves the location of the ground-based receiver and the satellite, i.e. the posed two-point boundary value problem is solved by an implicit finite difference scheme. This feature of the algorithm allows the fast and accurate computation of the signal travel-time delay, referred to as Slant Total Delay (STD), between a satellite and a ground-based receiver. We provide a technical description of the algorithm and estimate the uncertainty of STDs due to simplifying assumptions in the algorithm and due to the uncertainty of the refractivity field. In a first application, we compare STDs retrieved from GPS phase-observations at the German Research Centre for Geosciences Potsdam (GFZ STDs) with STDs derived from the European Center for Medium-Range Weather Forecasts analyses (ECMWF STDs). The statistical comparison for one month (August 2007) for a large and continuously operating network of ground-based receivers in Germany indicates good agreement between GFZ STDs and ECMWF STDs; the standard deviation is 0.5% and the mean deviation is 0.1%.

Description: In the vast Low Nutrient Low-Chlorophyll (LNLC) Ocean, the vertical nutrient supply from the subsurface to the sunlit surface waters is low, and atmospheric contribution of nutrients may be one order of magnitude greater over short timescales. The short turnover time of atmospheric Fe and N supply (〈1 month for nitrate) further supports deposition being an important source of nutrients in LNLC regions. Yet, the extent to which atmospheric inputs are impacting biological activity and modifying the carbon balance in oligotrophic environments has not been constrained. Here, we quantify and compare the biogeochemical impacts of atmospheric deposition in LNLC regions using both a compilation of experimental data and model outputs. A metadata-analysis of recently conducted field and laboratory bioassay experiments reveals complex responses, and the overall impact is not a simple “fertilization effect of increasing phytoplankton biomass” as observed in HNLC regions. Although phytoplankton growth may be enhanced, increases in bacterial activity and respiration result in weakening of biological carbon sequestration. The application of models using climatological or time-averaged non-synoptic deposition rates produced responses that were generally much lower than observed in the bioassay experiments. We demonstrate that experimental data and model outputs show better agreement on short timescale (days to weeks) when strong synoptic pulse of aerosols deposition, similar in magnitude to those observed in the field and introduced in bioassay experiments, is superimposed over the mean atmospheric deposition fields. These results suggest that atmospheric impacts in LNLC regions have been underestimated by models, at least at daily to weekly timescales, as they typically overlook large synoptic variations in atmospheric deposition and associated nutrient and particle inputs. Inclusion of the large synoptic variability of atmospheric input, and improved representation and parameterization of key processes that respond to atmospheric deposition, is required to better constrain impacts in ocean biogeochemical models. This is critical for understanding and prediction of current and future functioning of LNLC regions and their contribution to the global carbon cycle.

Description: The importance of the Gulf Stream Extension region in climate and seasonal prediction research is being increasingly recognised. Here we use satellite-derived eddy momentum fluxes to drive a shallow water model for the North Atlantic Ocean that includes the realistic ocean bottom topography. The results show that the eddy momentum fluxes can drive significant transport, sufficient to explain the observed increase in transport of the Gulf Stream following its separation from the coast at Cape Hatteras, as well as the observed recirculation gyres. The model also captures recirculating gyres seen in the mean sea surface height field within the North Atlantic Current system east of the Grand Banks of Newfoundland, including a representation of the Mann Eddy.

Description: A decade of weak convection in the Labrador Sea associated with decreasing water mass transformation, in combination with advective and eddy fluxes into the convection area, caused significant warming of the deep waters in both the central Labrador Sea and boundary current system along the Labrador shelf break. The connection to the export of Deep Water was studied based on moored current meter stations between 1997 and 2009 at the exit of the Labrador Sea, near the shelf break at 5˚3N. More than 100 year -long current meter records spanning the full water column have been analyzed with respect to high frequency variability, decaying from the surface to the bottom layer, and for the annual mean flow, showing intra- to interannual variability but no detectable decadal trend in the strength of the deep and near-bottom flow out of the Labrador Sea.

Description: Recent work shows that multichannel seismic (MCS) systems provide detailed information on the oceans' finestructure. The aim of this paper is to analyze if high order numerical algorithms are suitable to accurately model the extremely weak wavefield scattered by the oceans' finestructures. For this purpose, we generate synthetic shot records along a coincident seismic and oceanographic profile acquired across a Mediterranean salt lens in the Gulf of Cadiz. We apply a 2D finite-difference time-domain propagation model, together with second-order Complex Frequency Shifted Perfectly Matched Layers at the numerical boundaries, using as reference a realistic sound speed map with the lateral resolution of the seismic data. We show that our numerical propagator creates an acoustical image of the ocean finestructures including the salt lens that reproduces with outstanding detail the real acquired one

Description: Recently seismic reflection methods have been successfully applied to oceanographic issues. Here, we present a new approach, combining XBT and CTD surveys with seismic observations, to visualize long sections with a resolution down to a few meters. The challenge to a full investigation of mixing processes has been the tremendous span of spatial scales ranging from hundreds of kilometers to centimeters. Traditional hydrographic observations could only resolve the large scale effects by measuring temperature and salinity profiles at discrete locations typically several kilometers apart, whereas dedicated localized measurements allowed investigation of the ocean fine structure at the other end of the spatial spectrum. The intermediate scales have in contrast been difficult to observe systematically. Here we present temperature and salinity data inverted from seismic observations that cover the intermediate scales and provide a new approach to image mesoscale processes and allow the investigation of their dynamics at unprecedented resolution.

Description: Nitrogen (N) fixation by specialized microorganisms (diazotrophs) influences global plankton productivity because it provides the ocean with most of its bio-available N. However, its global rate and large-scale spatial distribution is still regarded with considerable uncertainty. Here we use a global ocean nitrogen isotope model, in comparison with δ15NO3− observations, to constrain the pattern of N2 fixation across the Pacific Ocean. N2 fixation introduces isotopically light atmospheric N2 from to the ocean (δ15N = 0‰) relative to the oceanic average near 5‰, which makes nitrogen isotopes suitable to infer patterns of N2 fixation. Including atmospheric iron limitation of diazotrophy in the model shifts the pattern of simulated N2 fixation from the South Pacific to the North Pacific and from the East Pacific westward. These changes considerably improve the agreement with meridional transects of available δ15NO3− observations, as well as excess P (PO43− − NO3−/16), suggesting that atmospheric iron deposition is indeed important for N fixation in the Pacific Ocean. This study highlights the potential for using δ15N observations and model simulations to constrain patterns and rates of N fixation in the ocean.

Description: The Atlantic‐Mediterranean exchange of water at Gibraltar represents a significant heat and freshwater sink for the North Atlantic and is a major control on the heat, salt and freshwater budgets of the Mediterranean Sea. Consequently, an understanding of the response of the exchange system to external changes is vital to a full comprehension of the hydrographic responses in both ocean basins. Here, we use a synthesis of empirical (oxygen isotope, planktonic foraminiferal assemblage) and modeling (analytical and general circulation) approaches to investigate the response of the Gibraltar Exchange system to Atlantic freshening during Heinrich Stadials (HSs). HSs display relatively flat W–E surface hydrographic gradients more comparable to the Late Holocene than the Last Glacial Maximum. This is significant, as it implies a similar state of surface circulation during these periods and a different state during the Last Glacial Maximum. During HS1, the gradient may have collapsed altogether, implying very strong water column stratification and a single thermal and d18Owater condition in surface water extending from southern Portugal to the eastern Alboran Sea. Together, these observations imply that inflow of Atlantic water into the Mediterranean was significantly increased during HS periods compared to background glacial conditions. Modeling efforts confirm that this is a predictable consequence of freshening North Atlantic surface water with iceberg meltwater and indicate that the enhanced exchange condition would last until the cessation of anomalous freshwater supply into to the northern North Atlantic. The close coupling of dynamics at Gibraltar Exchange with the Atlantic freshwater system provides an explanation for observations of increased Mediterranean Outflow activity during HS periods and also during the last deglaciation. This coupling is also significant to global ocean dynamics, as it causes density enhancement of the Atlantic water column via the Gibraltar Exchange to be inversely related to North Atlantic surface salinity. Consequently, Mediterranean enhancement of the Atlantic Meridional Overturning Circulation will be greatest when the overturning itself is at its weakest, a potentially critical negative feedback to Atlantic buoyancy change during times of ice sheet collapse.

Description: The stratospheric climate and variability from simulations of sixteen chemistryclimate models is evaluated. On average the polar night jet is well reproduced though its variability is less well reproduced with a large spread between models. Polar temperature biases are less than 5 K except in the Southern Hemisphere (SH) lower stratosphere in spring. The accumulated area of low temperatures responsible for polar stratospheric cloud formation is accurately reproduced for the Antarctic but underestimated for the Arctic. The shape and position of the polar vortex is well simulated, as is the tropical upwelling in the lower stratosphere. There is a wide model spread in the frequency of major sudden stratospheric warnings (SSWs), late biases in the breakup of the SH vortex, and a weak annual cycle in the zonal wind in the tropical upper stratosphere. Quantitatively, “metrics” indicate a wide spread in model performance for most diagnostics with systematic biases in many, and poorer performance in the SH than in the Northern Hemisphere (NH). Correlations were found in the SH between errors in the final warming, polar temperatures, the leading mode of variability, and jet strength, and in the NH between errors in polar temperatures, frequency of major SSWs, and jet strength. Models with a stronger QBO have stronger tropical upwelling and a colder NH vortex. Both the qualitative and quantitative analysis indicate a number of common and long‐standing model problems, particularly related to the simulation of the SH and stratospheric variability.

Description: The phosphorus budget of the pre-human modern ocean is constrained applying the most recent estimates of the natural riverine, eolian, and ice-rafted input fluxes, the phosphorus burial in marine sediments, and the hydrothermal removal of dissolved phosphate from the deep ocean. This review of current flux estimates indicates that the phosphorus budget of the ocean is unbalanced since the accumulation of phosphorus in marine sediments and altered oceanic crust exceeds the continental input of particulate and dissolved phosphorus. The phosphorus mass balance is further tested considering the dissolved phosphate distribution in the deep water column, the marine export production of particulate organic matter, rain rates of phosphorus to the seafloor, benthic dissolved phosphate fluxes, and the organic carbon to phosphorus ratios in marine particles. These independent data confirm that the phosphate and phosphorus budgets were not at steadystate in the pre-human global ocean. The ocean is losing dissolved phosphate at a rate of ≥ 11.6 x 1010 mol yr-1 corresponding to a decline in the phosphate inventory of ≥ 4.5 % kyr-1. Benthic data show that phosphate is preferentially retained in pelagic deep-sea sediments where extended oxygen exposure times favor the degradation of particulate organic matter and the up-take of phosphate in manganese and iron oxides and hydroxides. Enhanced C : P regeneration ratios observed in the deep water column (〉400 m water depth) probably reflect the preferential burial of phosphorus in pelagic sediments. Excess phosphate is released from continental margin sediments deposited in low-oxygen environments. The dissolved oxygen threshold value for the enhanced release of dissolved phosphate is ~20 μM. Benthic phosphate fluxes increase drastically when oxygen concentrations fall below this value.

Description: The decay kinetics of superoxide (O2−) reacting with organic matter was examined in oligotrophic waters at, and nearby, the TENATSO ocean observatory adjacent to the Cape Verde archipelago. Superoxide is the short-lived primary photochemical product of colored dissolved organic matter (CDOM) photolysis and also reacts with CDOM or trace metals (Cu, Fe) to form H2O2. In the present work we focused our investigations on reactions between CDOM and superoxide. O2− decay kinetics experiments were performed by adding KO2 to diethylenetriaminepentaacetic acid (DTPA) amended seawater and utilizing an established chemiluminescence technique for the detection of O2− at nM levels. In Cape Verdean waters we found a significant reactivity of superoxide with CDOM with maximal rates adjacent to the chlorophyll maximum, presumably from production of new CDOM from bacteria/phytoplankton. This work highlights a poorly understood process which impacts on the biogeochemical cycling of CDOM and trace metals in the open ocean.

Description: Geophysical datasets sensitive to different physical parameters can be used to improve resolution of Earth's internal structure. Herein, we jointly invert long-period magnetotelluric (MT) data and surface-wave dispersion curves. Our approach is based on a joint inversion using a genetic algorithm for a one-dimensional (1-D) isotropic structure, which we extend to 1-D anisotropic media. We apply our new anisotropic joint inversion to datasets from Central Germany demonstrating the capacity of our joint inversion algorithm to establish a 1-D anisotropic model that fits MT and seismic datasets simultaneously and providing new information regarding the deep structure in Central Germany. The lithosphere/asthenosphere boundary is found at approx. 84 km depth and two main anisotropic layers with coincident most conductive/seismic fast-axis direction are resolved at lower crustal and asthenospheric depths. We also quantify the amount of seismic and electrical anisotropy in the asthenosphere showing an emerging agreement between the two anisotropic coefficients.

Description: A natural carbon dioxide (CO2) seep was discovered during an expedition to the southern German North Sea (October 2008). Elevated CO2 levels of ∼10–20 times above background were detected in seawater above a natural salt dome ∼30 km north of the East-Frisian Island Juist. A single elevated value 53 times higher than background was measured, indicating a possible CO2 point source from the seafloor. Measured pH values of around 6.8 support modeled pH values for the observed high CO2 concentration. These results are presented in the context of CO2 seepage detection, in light of proposed subsurface CO2 sequestering and growing concern of ocean acidification. We explore the boundary conditions of CO2 bubble and plume seepage and potential flux paths to the atmosphere. Shallow bubble release experiments conducted in a lake combined with discrete-bubble modeling suggest that shallow CO2 outgassing will be difficult to detect as bubbles dissolve very rapidly (within meters). Bubble-plume modeling further shows that a CO2 plume will lose buoyancy quickly because of rapid bubble dissolution while the newly CO2-enriched water tends to sink toward the seabed. Results suggest that released CO2 will tend to stay near the bottom in shallow systems (〈200 m) and will vent to the atmosphere only during deep water convection (water column turnover). While isotope signatures point to a biogenic source, the exact origin is inconclusive because of dilution. This site could serve as a natural laboratory to further study the effects of carbon sequestration below the seafloor.

Description: Reconstructions of the spatial pattern of recent multi-decadal sea level trends in the Indian Ocean (IO) indicate a zonally-extended band in the southern tropics where sea level has substantially fallen between the 1960s and 1990s; the decline is consistent with the observed subsurface cooling associated with a shoaling thermocline in this region. Here the origin and spatio-temporal characteristics of these trends are elucidated by a sequence of ocean model simulations. Whereas interannual variability in the southwestern tropical IO appears mainly governed by IO atmospheric forcing, longer term changes in the south tropical IO involve a strong contribution from the western Pacific via wave transmission of thermocline anomalies through the Indonesian Archipelago, and their subsequent westward propagation by baroclinic Rossby waves. The late 20th-century IO subsurface cooling trend reversed in the 1990s, reflecting the major regime shift in the tropical Pacific easterlies associated with the Pacific Decadal Oscillation.

Description: We determined methane (CH4) emissions in a field enclosure experiment in a littoral freshwater marsh under the influence of experimentally simulated warming and enhanced nitrogen deposition. Methane emissions by ebullition from the marsh composed of Phragmites australis were measured with funnel traps deployed in a series of enclosures for two 3 week periods. Diffusive fluxes were estimated on the basis of measured CH4 concentrations and application of Fick's law. Neither diffusive nor ebullitive fluxes of methane were significantly affected by warming or nitrate enrichment, possibly because variability both within and among replicate experimental enclosures was high. Average emission rates resulted primarily from ebullition (0.2–30.3 mmol CH4 m−2 d−1), which were 4 orders of magnitude higher than estimated diffusive fluxes and were of similar importance as the coarsely estimated advective methane transport through plants. Significant correlations between dissolved oxygen and dissolved methane and ebullition flux suggest that methane release from the sediment might feed back positively on methane production by reducing dissolved oxygen in the water column and oxygen flux into the sediment. Nitrate may have a similar effect. Extrapolation of our limited data indicates that total methane fluxes from vegetated littoral zones of temperate lakes may contribute 0.5%–7% of the global natural CH4 emissions. These results emphasize the importance of freshwater marshes as sources of methane emissions to the atmosphere, even when they occupy only relatively small littoral areas. More detailed investigations are clearly needed to assess whether global warming and nitrogen deposition can have climate feedbacks by altering methane fluxes from these wetlands.

Description: The sensitivity of the hydrological cycle to changes in orbital forcing and atmospheric greenhouse gas (GHG) concentrations is assessed using a fully coupled atmosphere-ocean-sea ice general circulation model (Kiel Climate Model). An orbitally-induced intensification of the summer monsoon circulation during the Holocene and Eemian drives enhanced water vapor advection into the Northern Hemisphere, thereby enhancing the rate of water vapor changes by about 30% relative to the rate given by the Clausius-Clapeyron Equation, assuming constant relative humidity. Orbitally-induced changes in hemispheric-mean precipitation are fully attributed to inter-hemispheric water vapor exchange in contrast to a GHG forced warming, where enhanced precipitation is caused by increased both the moisture advection and evaporation. When considering the future climate on millennial time scales, both forcings combined are expected to exert a strong effect.

Description: We present results from a seismic refraction and wide-angle experiment surveying an oceanic core complex on the Mid-Atlantic Ridge at 22°19′N. Oceanic core complexes are settings where petrological sampling found exposed lower crustal and upper mantle rocks, exhumed by asymmetric crustal accretion involving detachment faulting at magmatically starved ridge sections. Tomographic inversion of our seismic data yielded lateral variations of P wave velocity within the upper 3 to 4 km of the lithosphere across the median valley. A joint modeling procedure of seismic P wave travel times and marine gravity field data was used to constrain crustal thickness variations and the structure of the uppermost mantle. A gradual increase of seismic velocities from the median valley to the east is connected to aging of the oceanic crust, while a rapid change of seismic velocities at the western ridge flank indicates profound differences in lithology between conjugated ridge flanks, caused by un-roofing lower crust rocks. Under the core complex crust is approximately 40% thinner than in the median valley and under the conjugated eastern flank. Clear PmP reflections turning under the western ridge flank suggest the creation of a Moho boundary and hence continuous magmatic accretion during core complex formation.

Description: The Mid-Atlantic Ridge south of the equator is a key region for many aspects of spreading axis studies, from biogeography to ridge-hotspot interaction. Despite this, the ridge axis had, until 2004, seen little systematic study. Repeated trips to the area since then have mapped and explored some 900 km of ridge length, from 2° to 14°S. The result is complete bathymetric and side-scan coverage of the axial region and the discovery and characterization of the first hydrothermal vents south of the equator. Such multisegment detailed and interdisciplinary coverage allows us to formulate a general model for the interplay between volcanism, tectonics, and hydrothermalism on a slow spreading ridge. The model defines three basic types of ridge morphology with specific hydrothermal characteristics: (a) a deep, tectonically dominated rift valley where hydrothermalism is seldom associated with volcanism and much more likely confined to long-lived bounding faults; (b) a shallower, segment-center bulge where a combination of repeated magmatic activity and tectonism results in repeated, possibly temporally overlapping periods of hydrothermal activity on the ridge axis; and (c) a very shallow axis beneath which temperatures in all but the uppermost crust are so high that deformation is ductile, inhibiting the formation of high-porosity deep fractures and severely depressing hydrothermal circulation. This model is used together with satellitederived predicted bathymetry to provide forecasts of the best places to look for hydrothermal sites in the remaining unexplored regions of the South Atlantic.

Description: The Sahara Slide is a giant submarine landslide on the northwest African continental margin. The landslide is located on the open continental slope offshore arid Western Sahara, with a headwall at a water depth of ∼2000 m. High primary productivity in surface waters drives accumulation of thick fine-grained pelagic/hemipelagic sediment sequences in the slide source area. Rare but large-scale slope failures, such as the Sahara Slide that remobilized approximately 600 km3 of sediment, are characteristic of this sedimentological setting. Seismic profiles collected from the slide scar reveal a stepped profile with two 100 m high headwalls, suggesting that the slide occurred retrogressively as a slab-type failure. Sediment cores recovered from the slide deposit provide new insights into the process by which the slide eroded and entrained a volcaniclastic sand layer. When this layer was entrained at the base of the slide it became fluidized and resulted in low apparent friction, facilitating the exceptionally long runout of ∼900 km. The slide location appears to be controlled by the buried headwall of an older slope failure, and we suggest that the cause of the slide relates to differential sedimentation rates and compaction across these scarps, leading to local increases of pore pressure. Sediment cores yield a date of 50–60 ka for the main slide event, a period of global sea level rise which may have contributed to pore pressure buildup. The link with sea level rising is consistent with other submarine landslides on this margin, drawing attention to this potential hazard during global warming.

Description: The supply of oxygen-rich water to the oxygen minimum zones (OMZs) of the eastern North and South Pacific via zonal tropical currents is investigated using shipboard acoustic Doppler current profiler and hydrographic section data. Near the equator, the Equatorial Undercurrent (EUC), Northern and Southern Subsurface Countercurrents (SCCs), and the Northern and Southern Intermediate Countercurrents (ICCs) all carry water that is oxygen richer than adjacent westward flows, thereby providing a net oxygen supply to the eastern Pacific OMZs. The synoptic velocity-weighted oxygen concentration difference between eastward and westward flows is typically 10–50 μmol kg−1. Subthermocline zonal oxygen fluxes reflect decreasing oxygen concentrations of the EUC, the SCCs, and the ICCs as they flow eastward. Approximately 30 year time series in well-sampled regions of the equatorial Pacific show oxygen content decreasing as rapidly as −0.55 μmol kg−1 yr−1 in the major oxygen supply paths of the OMZs for a 200–700 m layer and similar trends for a density layer spanning roughly these depths. This finding is in gross agreement with climate models, which generally predict expanding OMZs.

Description: Near the western boundary of the tropical North Atlantic, where the North Brazil Current (NBC) retroflects into the North Equatorial Countercurrent, large anticyclonic rings are shed. After separating from the retroflection region, the so-called NBC rings travel northwestward along the Brazilian coast, until they reach the island chain of the Lesser Antilles and disintegrate. These rings contribute substantially to the upper limb return flow of the Atlantic Meridional Overturning Circulation by carrying South Atlantic Water into the northern subtropical gyre. Their relevance for the northward transport of South Atlantic Water depends on the frequency of their generation as well as on their horizontal and vertical structure. The ring shedding and propagation and the complex interaction of the rings with the Lesser Antilles are investigated in the inline equation Family of Linked Atlantic Model Experiments (FLAME) model. The ring properties simulated in FLAME reach the upper limit of the observed rings in diameter and agree with recent observations on seasonal variability, which indicates a maximum shedding during the first half of the year. When the rings reach the shallow topography of the Lesser Antilles, they are trapped by the island triangle of St. Lucia, Barbados and Tobago and interact with the island chain. The model provides a resolution that is capable of resolving the complex topographic conditions at the islands and illuminates various possible fates for the water contained in the rings. It also reproduces laboratory experiments that indicate that both cyclones and anticyclones are formed after a ring passes through a topographic gap. Trajectories of artificial floats, which were inserted into the modeled velocity field, are used to investigate the pathways of the ring cores and their fate after they encounter the Lesser Antilles. The majority of the floats entered the Caribbean, while the northward Atlantic pathway was found to be of minor importance. No prominent pathway was found east of Barbados, where a ring could avoid the interaction with the islands and migrate toward the northern Lesser Antilles undisturbed.

Description: The Southern Hemisphere winter stratosphere exhibits prominent traveling planetary-scale Rossby waves, which generally are not able to induce Stratospheric Sudden Warmings. A series of runs of a simplified general circulation model is presented, aimed at better understanding the generation of these waves. While the generation of planetary-scale traveling waves through the interaction of synoptic-scale waves is observed in a control run, when the model is truncated to permit only waves with zonal wave number 1 or 2, the long waves are found to increase in strength, leading to a considerably more active stratosphere including Sudden Warmings comparable in strength to Northern Hemisphere winter. This finding suggests that the role of tropospheric synoptic eddies is two-fold: while generating a weak planetary-scale wave flux into the stratosphere, their main effect is to suppress baroclinic instability of planetary-scale waves by stabilizing the tropospheric mean state.

Description: The Mw = 9.0 earthquake of 11 March 2011 at the Japan Trench and its devastating tsunami underscore the importance of understanding seismogenic behavior of subduction faults and realistically estimating the potential size of future earthquakes and tsunamis. For the Cascadia subduction zone (Figure 1a), a critical knowledge gap is the level of microseismicity offshore, especially near the megathrust, needed to better understand the state of the locked zone. In 2010 the first detailed seafloor earthquake monitoring campaign along the northern Cascadia subduction zone recorded nearby earthquakes in the local magnitude (ML) range from possibly around zero to 3.8 (Figures 1b and 1c) and larger earthquakes from outside this region. Preliminary analyses indicate that the network appears to have yielded a fairly complete catalog for events with ML 〉 1.2. Only a few tens of these events occurred beneath the continental shelf and slope (Figure 1a). The majority of the earthquakes were located along the margin-perpendicular Nootka fault zone. The relatively low seismicity away from the Nootka fault is consistent with a fully locked megathrust. Land-based GPS measurements cannot resolve the question of whether the offshore part of the megathrust seismogenic zone is narrow and fully locked or wider and only partially locked (slowly creeping). If it were only partially locked, the seafloor seismometer data should show many more small earthquakes along the interface than were actually detected.

Description: The accurate dating of meteorite impact structures on Earth has proven to be challenging. Melt sheets are amenable to high-precision dating by the U-Pb and 40Ar/39Ar methods, but many impact events do not produce them, or they are not preserved. In cases where high-temperature shock metamorphism of the target materials has occurred without widespread melting, these isotopic chronometers may be partially reset and yield dates that are difficult to interpret unambiguously as the age of impact. However, the (U-Th)/He chronometer is sensitive to thermal resetting and can provide a powerful new tool for dating impactites. We report (U-Th)/He dates for accessory minerals from the Manicouagan impact structure in Quebec, Canada. Nine zircons from a melt sheet sample yield a weighted mean age of 213.2 ± 5.4 Ma (2SE), indistinguishable from the published 214 ± 1 Ma (2σ) U-Pb zircon age for the impact. In contrast, five apatites from this sample yield dates between 205.9 ± 6.5 and 162.0 ± 5.3 Ma (2σ), indicating variable postimpact helium loss due to low-temperature thermal disturbance. Preimpact titanite crystals from a shocked meta-anorthosite sample yield two dates consistent with the impact age, at 212 ± 27 and 214 ± 13 Ma (2σ), and two younger dates of 189.6 ± 6.9 and 192.2 ± 9.8 Ma (2σ), suggestive of postimpact helium loss. These results indicate that (U-Th)/He chronometry is a suitable method for dating impact events, although interpretation of the results requires recognition of possible 4He loss related to reheating subsequent to impact.

Description: The tropospheric response to sudden stratospheric warmings (SSWs) is analyzed in an idealized model setup regarding the respective roles of planetary-scale and synoptic-scale waves. The control model run includes a full interactive wave spectrum, while a second run includes interactive planetary-scale waves but only the time-mean synoptic-scale wave forcing from the control run. In both runs, the tropospheric response is characterized by the negative phase of the respective tropospheric annular mode. But given their different latitudinal structure, the control run shows the expected response, i.e., an equatorward shift of the tropospheric jet, whereas the response in the absence of interactive synoptic eddies is characterized by a poleward jet shift. This opposite jet shift is associated with a different planetary wave variability that couples with the zonal flow between the stratosphere and the surface. These results indicate that the synoptic eddy feedback is necessary for the observed tropospheric response to SSWs.

Description: At 8:04 P.M. Pacific daylight time (PDT) on 27 October 2012 (03:04 universal time (UT), 28 October), Canada's second largest instrumentally recorded earthquake rocked Haida Gwaii (formerly Queen Charlotte Islands) and the mainland coast of British Columbia. The M 7.7 event off the west coast of Moresby Island caused a tsunami with local runup of more than 7 meters and amplitudes up to 0.8 meter on tide gauges 4000 kilometers away in Hawaii. Shaking was felt as far away as the Yukon, Alberta, Washington, and Montana, up to 1500 kilometers away. Little damage was caused, as the immediate region is an uninhabited National Park Reserve. The closest point of the rupture zone, as defined by aftershocks (Figures 1a and 1c), was 50 kilometers from the nearest community, Queen Charlotte, where damage was confined to a few chimneys and slumped roads.

Description: To reconstruct the history of water mass exchange between the NE Atlantic and the Nordic seas, sediment cores from ∼2 km water depth were studied across Termination II (TII) and through the last interglaciation (MIS5e). During early TII the sudden appearance of the low-latitude planktonic foraminifera Beella megastoma is noted in both regions along with a steep decrease in benthic foraminiferal δ18O. Since other proxies indicate that surface waters were cold and stratified because of meltwater, conditions which prevented near-surface thermohaline circulation and vertical convection in the Nordic seas, water mass exchange between the two areas occurred at the subsurface. During later TII, surface conditions changed, and this subsurface circulation style was eventually replaced by vertical convection. In the Nordic seas, B. megastoma vanished from the record together with ice-rafted debris (IRD) at the end of TII, while subpolar foraminiferal abundance rose. Peak interglacial conditions with intensive vertical convection now fully developed, generating a bottom water temperature gradient of ∼4°C between the two areas. However, surface water temperatures deteriorated in the Nordic seas already notably before IRD recurred, and δ18O increased at the end of MIS5e.

Description: This chapter contains sections titled: Introduction Temperature Data Temperature Observations Stratospheric Water Vapor Summary and Discussion Appendix A: Linear Regression Analysis

Description: [1] To assess the temporal variability in the methane fluxes from marine sediments that overly gas hydrate bearing sediments and the factors that might control its rate, in situ methane concentrations were measured near Bullseye Vent on the Northern Cascadia continental margin. A long-term sampling device collected overlying water and pore-fluid samples from 25 cm above seafloor, at the sediment-water interface (SWI), and 7 cmbsf (centimeters below seafloor) over a 9 month period (August 2009–May 2010). These samples provide a record at ∼4 day resolution of in situ methane, ethane, propane, sulfate, and chloride concentrations, as well as stable carbon isotope ratios of methane (δ13C-CH4) and dissolved inorganic carbon (δ13C-DIC). We show that pore fluids near the SWI are saturated or supersaturated with respect to methane (∼80 mM) and the methane flux from the seabed is variable over time. We hypothesized that regional seismic activity controlled this variable CH4 flux in the Northern Cascadia continental margin setting. However, we found no direct correlation between earthquakes and CH4 flux. We also posited alternative controls on CH4 flux variability, such as storms, regional oceanography and microbial activity. Again, no direct correlation was seen. This study takes first steps toward exploring which physical factors play a role in methane flux from hydrate-bearing sediments.

Description: The recent discovery of large ionospheric disturbances associated with sudden stratospheric warmings (SSW) has challenged the current understanding of mechanisms coupling the stratosphere and ionosphere. Non-linear interaction of planetary waves and tides has been invoked as a primary mechanism for such coupling. Here we show that planetary waves may play a more complex role than previously thought. Planetary wave forcing induces a global circulation that leads to the build-up of ozone density in the tropics at 30–50 km altitude, the primary region responsible for the generation of the migrating semidiurnal tide. The increase in the ozone density reaches 25% and lasts for ∼35 days following the SSW, long after the collapse of the planetary waves. Ozone enhancements are not only associated with SSW but are also observed after other amplifications in planetary waves. In addition, the longitudinal distribution of the ozone becomes strongly asymmetric, potentially leading to the generation of non-migrating semidiurnal tides. We report a persistent increase in the variability of ionospheric total electron content that coincides with the increase in stratospheric ozone and we suggest that the ozone fluctuations affect the ionosphere through the modified tidal forcing.

Description: There is an association between the polar stratospheric temperature in the northern winter and the solar cycle in the winters when the 50-mb equatorial winds are westerly: The lower the sunspot number in such winters, the lower is the temperature. No major mid-winter warmings occurred in these winters when the sunspot number was below about 100. There is no such relationship in the easterly phase of the QBO. In that phase the temperatures are generally higher than in the westerly phase, and major mid-winter warmings occur regardless of the state of the solar cycle.

Description: Gridded datasets have been constructed for the physical parameters controlling the formation of gas hydrate. Simultaneous solutions to the gas hydrate phase and pressure-temperature equations were obtained for each node of these grids. The results show sub-bottom depths of potential methane hydrate BSRs. These are presented as colour contoured grids for the thickness of the hydrate stability zone for the European margins. The chart proposes that if gas hydrates exist over these areas, then this is the potential depth to the BSR. This depth is generally greater along the continental margins and also increases with increasing age of the margin. Shallowing of the BSR can be seen over areas of thinned continental crust and plateaus. The geological factors controlling gas hydrate formation determine areas of likely occurrence so the apparent paucity of identified hydrate layer BSRs off the European margins, particularly in the Mediterranean, is most notable.

Description: We present data on the burial, displacement and exhumation history of the Himalayan fold-thrust belt in eastern Bhutan. These data document the magnitude and timing of displacement of large, discrete structures and highlight temporal variability in shortening rates. Eight new40Ar/39Ar ages from white mica, 32 new zircon (U-Th)/He ages, 7 new apatite fission track ages, and 1 new U-Pb zircon (LA-MC-ICP-MS) metamorphic rim growth age are combined with published cooling ages and deformation temperatures, and incremental shortening magnitudes from restorations of two published balanced cross sections, to illustrate the kinematic and temporal development of the Bhutan thrust belt. Integrating these data from ∼23 Ma to the present illustrates rapid horizontal shortening rates (28–35 mm/yr) between 23–20 Ma and 15–10 Ma, separated by more moderate rates (10–23 mm/yr). Shortening rates decrease significantly to 7–10 mm/yr (and possibly as low as 3–4 mm/yr) from 10 to 0 Ma. This decrease is interpreted to represent the onset of strain partitioning in the eastern part of the Himalayan-Tibetan orogenic system, between shortening in the Bhutan thrust belt, uplift of the Shillong Plateau, and deformation and outward growth of the northern and eastern Tibetan Plateau. Within estimated error, horizontal shortening rates during emplacement of the Main Central thrust sheet and during construction of the upper Lesser Himalayan duplex approached India-Asia tectonic velocities. Thus, for periods of time between ∼23–20 Ma and ∼15–10 Ma, the Bhutan thrust belt may have absorbed nearly all India-Asian convergence at this longitude.

Description: Archean cratons belong to the most remarkable features of our planet since they represent continental crust that has avoided reworking for several billions of years. Even more, it has become evident from both geophysical and petrological studies that cratons exhibit deep lithospheric keels which equally remained stable ever since the formation of the cratons in the Archean. Dating of inclusions in diamonds from kimberlite pipes gives Archean ages, suggesting that the Archean lithosphere must have been cold soon after its formation in the Archean (in order to allow for the existence of diamonds) and must have stayed in that state ever since. Yet, although strong evidence for the thermal stability of Archean cratonic lithosphere for billions of years is provided by diamond dating, the long-term thermal stability of cratonic keels was questioned on the basis of numerical modeling results. We devised a viscoelastic mantle convection model for exploring cratonic stability in the stagnant lid regime. Our modeling results indicate that within the limitations of the stagnant lid approach, the application of a sufficiently high temperature-dependent viscosity ratio can provide for thermal craton stability for billions of years. The comparison between simulations with viscous and viscoelastic rheology indicates no significant influence of elasticity on craton stability. Yet, a viscoelastic rheology provides a physical transition from viscously to elastically dominated regimes within the keel, thus rendering introduction of arbitrary viscosity cutoffs, as employed in viscous models, unnecessary.

Description: There is an incomplete description of the mid-depth circulation and its link to the oxygen minimum zone (OMZ) in the eastern tropical South Pacific. Subsurface currents of the OMZ in the eastern tropical South Pacific are investigated with a focus at 400 m depth, close to the core of the OMZ, using several Acoustic Doppler Current Profiler sections recorded in January and February 2009. Five profiling floats with oxygen sensors were deployed along 85°50’W in February 2009 with a drift depth at 400 m. Their spreading paths are compared with the model flow field from a 1/10° Tropical Pacific model (TROPAC01) and the Simple Ocean Data Assimilation (SODA) model. Overall the mean currents in the eastern tropical South Pacific are weak so that eddy variability influences the flow and ultimately feed oxygen-poor water to the OMZ. The center of the OMZ is a stagnant area so that floats stay much longer in this region and can even reverse direction. In one case of one float deployed at 8°S returned to the same location after 15 month. On the northern side of the OMZ in the equatorial current system, floats move rapidly to the west. Most current bands reported for the near surface layer exist also in the depth range of the OMZ. A schematic circulation flow field for the OMZ core depth is derived which shows the northern part of the South Pacific subtropical gyre south of the OMZ and the complicated zonal equatorial flow field north of the OMZ.

Description: The ongoing warming of bottom water in the Arctic region is anticipated to destabilize some of the gas hydrate present in shallow seafloor sediment, potentially causing the release of methane from dissociating hydrate into the ocean and the atmosphere. Ocean-bottom seismometer (OBS) experiments were conducted along the continental margin of western Svalbard to quantify the amount of methane present as hydrate or gas beneath the seabed. P- and S-wave velocities were modeled for five sites along the continental margin, using ray-trace forward modeling. Two southern sites were located in the vicinity of a 30 km long zone where methane gas bubbles escaping from the seafloor were observed during the cruise. The three remaining sites were located along an E-W orientated line in the north of the margin. At the deepest northern site, Vp anomalies indicate the presence of hydrate in the sediment immediately overlying a zone containing free gas up to 100-m thick. The acoustic impedance contrast between the two zones forms a bottom-simulating reflector (BSR) at approximately 195 m below the seabed. The two other sites within the gas hydrate stability zone (GHSZ) do not show the clear presence of a BSR or of gas hydrate. However, anomalously low Vp, indicating the presence of free gas, was modeled for both sites. The hydrate content was estimated from Vp and Vs, using effective-medium theory. At the deepest northern site, modeling suggests a pore-space hydrate concentration of 7–12%, if hydrate forms as part of a connected framework, and about 22% if it is pore-filling. At the two other northern sites, located between the deepest site and the landward limit of the GHSZ, we suggest that hydrate is present in the sediment as inclusions. Hydrate may be present in small quantities at these two sites (4–5%) of the pore space. The variation in lithology for the three sites indicated by high-resolution seismic profiles may control the distribution, concentration and formation of hydrate and free gas.

Description: A 1.6 km riser borehole was drilled at site C0009 of the NanTroSEIZE, in the center of the Kumano forearc basin, as a landward extension of previous drilling in the southwest Japan Nankai subduction zone. We determined principal horizontal stress orientations from analyses of borehole breakouts and drilling-induced tensile fractures by using wireline logging formation microresistivity images and caliper data. The maximum horizontal stress orientation at C0009 is approximately parallel to the convergence vector between the Philippine Sea plate and Japan, showing a slight difference with the stress orientation which is perpendicular to the plate boundary at previous NanTroSEIZE sites C0001, C0004 and C0006 but orthogonal to the stress orientation at site C0002, which is also in the Kumano forearc basin. These data show that horizontal stress orientations are not uniform in the forearc basin within the surveyed depth range and suggest that oblique plate motion is being partitioned into strike-slip and thrusting. In addition, the stress orientations at site C0009 rotate clockwise from basin sediments into the underlying accretionary prism.

Description: [1] The Integrated Ocean Drilling Program (IODP) Expedition 310 recovered drill cores from the drowned reefs around the island of Tahiti (17°40′S, 149°30′W), many of which contained samples of massive corals from the genus Porites. Herein we report on one well-preserved fossil coral sample: a 13.6 cm long Porites sp. dated by uranium series techniques at 9523 ± 33 years. Monthly δ18O and Sr/Ca determinations reveal nine clear and robust annual cycles. Coral δ18O and Sr/Ca determinations estimate a mean temperature of ∼24.3°C (∼3.2°C colder than modern) for Tahiti at 9.5 ka; however, this estimate is viewed with caution since potential sources of cold bias in coral geochemistry remain to be resolved. The interannual variability in coral δ18O is similar between the 9.5 ka coral record and a modern record from nearby Moorea. The seasonal cycle in coral Sr/Ca is approximately the same or greater in the 9.5 ka coral record than in modern coral records from Tahiti. Paired analysis of coral δ18O and Sr/Ca indicates cold/wet (warm/dry) interannual anomalies, opposite from those observed in the modern instrumental record.

Description: Arctic Ocean freshening can exert a controlling influence on global climate, triggering strong feedbacks on ocean‐atmospheric processes and affecting the global cycling of the world’s oceans. Glacier‐fed ocean currents such as the Alaska Coastal Current are important sources of freshwater for the Bering Sea shelf, and may also influence the Arctic Ocean freshwater budget. Instrumental data indicate a multiyear freshening episode of the Alaska Coastal Current in the early 21st century. It is uncertain whether this freshening is part of natural multidecadal climate variability or a unique feature of anthropogenically induced warming. In order to answer this, a better understanding of past variations in the Alaska Coastal Current is needed. However, continuous long‐term high‐resolution observations of the Alaska Coastal Current have only been available for the last 2 decades. In this study, specimens of the long‐lived crustose coralline alga Clathromorphum nereostratum were collected within the pathway of the Alaska Coastal Current and utilized as archives of past temperature and salinity. Results indicate that coralline algal Mg/Ca ratios provide a 60 year record of sea surface temperatures and track changes of the Pacific Decadal Oscillation, a pattern of decadal‐to‐multidecadal ocean‐atmosphere climate variability centered over the North Pacific. Algal Ba/Ca ratios (used as indicators of coastal freshwater runoff) are inversely correlated to instrumentally measured Alaska Coastal Current salinity and record the period of freshening from 2001 to 2006. Similar multiyear freshening events are not evident in the earlier portion of the 60 year Ba/Ca record. This suggests that the 21st century freshening of the Alaska Coastal Current is a unique feature related to increasing glacial melt and precipitation on mainland Alaska.

Description: The Indian Ocean water that ends up in the Atlantic Ocean detaches from the Agulhas Current retroflection predominantly in the form of Agulhas rings and cyclones. Using numerical Lagrangian float trajectories in a high-resolution numerical ocean model, the fate of coherent structures near the Agulhas Current retroflection is investigated. It is shown that within the Agulhas Current, upstream of the retroflection, the spatial distributions of floats ending in the Atlantic Ocean and floats ending in the Indian Ocean are to a large extent similar. This indicates that Agulhas leakage occurs mostly through the detachment of Agulhas rings. After the floats detach from the Agulhas Current, the ambient water quickly looses its relative vorticity. The Agulhas rings thus seem to decay and loose much of their water in the Cape Basin. A cluster analysis reveals that most water in the Agulhas Current is within clusters of 180 km in diameter. Halfway in the Cape Basin there is an increase in the number of larger clusters with low relative vorticity, which carry the bulk of the Agulhas leakage transport through the Cape Basin. This upward cascade with respect to the length scales of the leakage, in combination with a power law decay of the magnitude of relative vorticity, might be an indication that the decay of Agulhas rings is somewhat comparable to the decay of two-dimensional turbulence.

Description: Interannual anomalies in vertical profiles of stratospheric ozone, in both equatorial and extratropical regions, have been shown to exhibit a strong seasonal persistence, namely, extended temporal autocorrelations during certain times of the calendar year. Here we investigate the relationship between this seasonal persistence of equatorial and extratropical ozone anomalies using the SAGE-corrected SBUV data set, which provides a long-term ozone profile time series. For the regions of the stratosphere where ozone is under purely dynamical or purely photochemical control, the seasonal persistence of equatorial and extratropical ozone anomalies arises from distinct mechanisms but preserves an anticorrelation between tropical and extratropical anomalies established during the winter period. In the 16–10 hPa layer, where ozone is controlled by both dynamical and photochemical processes, equatorial ozone anomalies exhibit a completely different behavior compared to ozone anomalies above and below in terms of variability, seasonal persistence, and especially the relationship between equatorial and extratropical ozone. Cross-latitude-time correlations show that for the 16–10 hPa layer, Northern Hemisphere (NH) extratropical ozone anomalies show the same variability as equatorial ozone anomalies but lagged by 3–6 months. High correlation coefficients are observed during the time frame of seasonal persistence of ozone anomalies, which is June–December for equatorial ozone and shifts by approximately 3–6 months when going from the equatorial region to NH extratropics. Thus in the transition zone between dynamical and photochemical control, equatorial ozone anomalies established in boreal summer/autumn are mirrored by NH extratropical ozone anomalies with a time lag similar to transport time scales. Equatorial ozone anomalies established in boreal winter/spring are likewise correlated with ozone anomalies in the Southern Hemisphere extratropics with a time lag comparable to transport time scales, similar to what is seen in the NH. However, the correlations between equatorial and SH extratropical ozone in the 10–16 hPa layer are weak.

Description: The linear relation between the strength of the Agulhas Current at nominal latitude 34°S and the gradient in sea level height anomaly across the current is investigated in a 1/10° resolution regional numerical ocean model. Our results show that the strength of the current can be estimated with reasonable accuracy using altimeter data, once it has been calibrated using in-situ transport measurements. Three years of transport measurements provide a calibration with worst-case correlation R = 0.78. In that case the errors in proxy transport have a standard deviation of 9.8 Sv, compared to a 20.2 Sv standard deviation of the transport time series itself. From these results we conclude that the design of the Agulhas Current Timeseries (ACT) experiment, a three-year deployment of moorings across the Agulhas Current and along a TOPEX/Jason altimeter ground track, will likely produce a good quality multi-decadal time series of Agulhas Current strength.

Description: The Ka band radar imaging mechanism of the submerged wreck/sand ribbon of the motor vessel (M/V) Birkenfels in the southern North Sea is investigated by applying the quasi‐specular scattering theory and considering the capillary as well as the gravity wave ranges of the wave energy density spectrum. For the imaging of wrecks and other oceanographic and meteorological phenomena at the sea surface it is assumed that quasi‐specular scattering becomes dominant at higher radar frequencies like Ka and X band and wind speeds ≥ 7–8 m s−1. Multibeam echo sounder images of the Birkenfels wreck and associated sand ribbons as well as other available environmental in situ data have been analyzed. The formation of sand ribbons at the sea bed and the manifestation of its radar signatures at the water surface are caused by an elliptical vortex or helical flow cell triggered by unidirectional tidal current flow interacting with the wreck. The difference between simulated and measured normalized radar cross section (NRCS) modulation as a function of the space variable is less than 31.6%. Results are presented for NRCS simulations dependent on position for different effective incidence angles, unidirectional current speeds, wind speeds, and relaxation rates. The calculated current gradient or strain rate of the imaging theory has the same order of magnitude as those obtained for marine sand waves. This implies that the responsible hydrodynamic interaction mechanism is able to produce radar signatures of submerged wrecks/sand ribbons and make them visible at the sea surface.

Description: There is an evidence of the increasing intensity as well as occurrence frequency of the so-called central Pacific (CP) El Niño events since the 1990s. We examine whether such an increase in the frequency of CP El Niño may be a manifestation of natural climate variability. A control simulation of the Kiel Climate Model, run for 4200 years with the present values of greenhouse gases, exhibit large variations of the occurrence frequency of the CP El Niño versus the eastern Pacific (EP) El Niño. A model simulates to some extent changes in the occurrence ratio of CP and EP El Niño in comparison with the observations. Therefore, we can not exclude the possibility that an increasing of occurrence frequency of CP El Niño during recent decades in the observation could be a part of natural variability in the tropical climate system.

Description: Using a global ocean model with regionally focused high resolution (1/10°) in the East China Sea (ECS), we studied the oceanic heat budget in the ECS. The modeled sea surface height variability and eddy kinetic energy are consistent with those derived from satellite altimetry. Significant levels of eddy kinetic energy are found east of the Ryukyu Islands and east of Taiwan, where the short-term variability is spawned by active mesoscale eddies coalescing with the circulation. Furthermore, the simulated vertical cross-stream structure of the Kuroshio (along the Pollution Nagasaki line) and the volume transport through each channel in the ECS are in good agreement with the observational estimates. The time-averaged temperature fluxes across the Taiwan Strait (TWS), Tsushima Strait (TSS), and the 200 m isobath between Taiwan and Japan are 0.20 PW, 0.21 PW, and 0.05 PW, respectively. The residual heat flux of 0.04 PW into the ECS is balanced by the surface heat loss. The eddy temperature flux across the 200 m isobath is 0.005 PW, which accounts for 11.2% of the total temperature flux. The Kuroshio onshore temperature flux has two major sources: the Kuroshio intrusion northeast of Taiwan and southwest of Kyushu. The Ekman temperature flux induced by the wind stress in the ECS shows the same seasonal cycle and amplitude as the onshore temperature flux, with a maximum in autumn and a minimum in summer. We conclude that the Ekman temperature flux dominates the seasonal cycle of Kuroshio onshore flux.

Description: At the Costa Rica margin along the Middle America Trench along‐strike variations in heat flow are well mapped. These variations can be understood in terms of either ventilated fluid flow, where exposed basement allows fluids to freely advect heat between the crustal aquifer and ocean, or insulated fluid flow where continuous sediment cover restricts heat advection to within the crustal aquifer. We model fluid flow within the subducting aquifer using Nusselt number approximations coupled with finite element models of subduction and explore its effect on temperatures along the subduction thrust. The sensitivity of these models to the initial thermal state of the plate and styles of fluid flow, either ventilated or insulated, is explored. Heat flow measurements on cool crust accreted at the East Pacific Rise are consistent with ventilated hydrothermal cooling that continues with subduction. These models yield much cooler temperatures than predicted from simulations initialized with conductive predictions and without hydrothermal circulation. Heat flow transects on warm crust accreted at the Cocos‐Nazca spreading center are consistent with models of insulated hydrothermal circulation that advects heat updip within the subducting crustal aquifer. Near the trench these models are warmer than conductive predictions and cooler than conductive predictions downdip of the trench. Comparisons between microseismicity and modeled isotherms suggest that the updip limit of microseismicity occurs at temperatures warmer than 100°C and that the downdip extent of microseismicity is bounded by the intersection of the subduction thrust with the base of the overriding crust.

Description: The thermal structure of convergent margins provides information related to the tectonics, geodynamics, metamorphism, and fluid flow of active plate boundaries. We report 176 heat flow measurements made with a violin bow style probe across the Costa Rican margin at the Middle America Trench. The probe measurements are collocated with seismic reflection lines. These seismic reflection lines show widespread distribution of bottom‐simulating reflectors (BSRs). To extend the spatial coverage of heat flow measurements we estimate heat flow from the depth of BSRs. Comparisons between probe measurements and BSR‐derived estimates of heat flow are generally within 10% and improve with distance landward of the deformation front. Together, these determinations provide new information on the thermal regime of this margin. Consistent with previous studies, the margin associated with the northern Nicoya Peninsula is remarkably cool. We define better the southern boundary of the cool region. The northern extent of the cool region remains poorly determined. A regional trend of decreasing heat flow landward of the deformation front is apparent, consistent with the downward advection of heat by the subducting Cocos Plate. High wave number variability at a scale of 5–10 km is significantly greater than the measurement uncertainty and is greater south of the northern Nicoya Peninsula. These heat flow anomalies vary between approximately 20 and 60 mW m−2 and are most likely due to localized fluid flow through mounds and faults on the margin. Simple one‐dimensional models show that these anomalies are consistent with flow rates of 7–15 mm yr−1. Across the margin toe variability is significant and likely due to fluid flow through deformation structures associated with the frontal sedimentary prism.

Description: A system for the shipboard measurement of air-sea CO2 fluxes by eddy covariance was developed and tested. The system was designed to reduce two major sources of experimental uncertainty previously reported. First, the correction for in situ water vapor fluctuations (the “Webb” correction) was reduced by 97% by drying the air sample stream. Second, motion sensitivity of the gas analyzer was reduced by using an open-path type sensor that was converted to a closed-path configuration to facilitate drying of the air stream. High-quality CO2 fluxes were obtained during 429 14 min flux intervals during two cruises in the North Atlantic. The results suggest that the gas analyzer resolved atmospheric CO2 fluctuations well below its RMS noise level. This noise was uncorrelated with the vertical wind and therefore filtered out by the flux calculation. Using climatological data, we estimate that the techniques reported here could enable high-quality measurements of air-sea CO2 flux over much of the world oceans.

Description: Blue mussel individuals (Mytilus edulis) were cultured at four different salinities (17, 20, 29, and 34). During the course of the experiment, temperature was gradually increased from 6°C to 14°C. Mg/Ca and Sr/Ca ratios of the shell calcite portions produced during the 9 weeks of experimental treatment as well parts that were precipitated before the treatment phase were measured by laser ablation–multicollector–inductively coupled plasma–mass spectrometry. Mg/Ca ratios show a positive correlation with temperature for individuals cultured at salinity 29 and 34 (Mg/Ca (mmol/mol) ∼ (0.2–0.3)*T (°C)), while for individuals cultured at low salinities (17, 20) no trend was observed. Sr/Ca ratios were not affected by temperature but strongly by salinity. The data show very strong biological influence (“individual differences” and “physiological variability”) on elemental ratios (79% on Mg/Ca and 41% on Sr/Ca) in M. edulis calcite. The results challenge the use of blue mussel shell data as environmental proxies.

Description: New high-resolution multichannel seismic data GWADASEIS-2009 and JC45/46-2010 cruises; 72 and 60 channels, respectively) combined with previous data(AGUADOMAR-1999 and CARAVAL-2002; 6 and 24 channels, respectively) allow a detailed investigation of mass-wasting processes around the volcanic island of Montserrat in the Lesser Antilles. Seven submarine deposits have sources on the flanks of Montserrat, while three are related to the nearby Kahouanne submarine volcanoes. The most voluminous deposit (∼20 km3) within the Bouillante-Montserrat half-graben has not been described previously and is probably related to a flank instability of the Centre Hills Volcano on Montserrat, while other events are related to the younger South Soufrière Hills-Soufrière Hills volcanic complex. All deposits are located to the south or southeast of the island in an area delimited by faults of the Bouillante-Montserrat half-graben. They cover a large part of the southeast quarter of the surrounding seafloor (∼520 km2), with a total volume of ∼40 km3. Our observations suggest that the Bouillante-Montserrat half-graben exerts a control on the extent and propagation of the most voluminous deposits. We propose an interpretation for mass-wasting processes around Montserrat similar to what has happened for the southern islands of the Lesser Antilles.

Description: A global Earth System model is employed to investigate the role of direct temperature effects in the response of marine ecosystems to climate change. While model configurations with and without consideration of explicit temperature effects can reproduce observed current biogeochemical tracer distributions and estimated carbon export about equally well, carbon flow through the model ecosystem reveals strong temperature sensitivities. Depending on whether biological processes are assumed temperature sensitive or not, simulated marine net primary production (NPP) increases or decreases under projected climate change driven by a business-as-usual CO2 emission scenario for the 21st century. This suggests that indirect temperature effects such as changes in the supply of nutrients and light are not the only relevant factors to be considered when modeling the response of marine ecosystems to climate change. A better understanding of direct temperature effects on marine ecosystems is required before even the direction of change in NPP can be reliably predicted.

Description: Mud volcanoes are seafloor expressions of focused fluid flow that are common in compressional tectonic settings. New high-resolution 3-D seismic data from the Mercator mud volcano (MMV) and an adjacent buried mud volcano (BMV) image the internal structure of the top 800 m of sediment at both mud volcanoes, revealing that both are linked and have been active episodically. The total volumes of extruded mud range between 0.15 and 0.35 km3 and 0.02–0.05 km3 for the MMV and the BMV, respectively. The pore water composition of surface sediment samples suggests that halokinesis has played an important role in the evolution of the mud volcanoes. We propose that erosion of the top of the Vernadsky Ridge that underlies the mud volcanoes activated salt movement, triggering deep migration of fluids, dissolution of salt, and sediment liquefaction and mobilization since the end of the Pliocene. Since beginning of mud volcanism in this area, the mud volcanoes erupted four times while there was only one reactivation of salt tectonics. This implies that there are other mechanisms that trigger mud eruptions. The stratigraphic relationship of mudflows from the MMV and BMV indicates that the BMV was triggered by the MMV eruptions. This may either be caused by loading-induced hydrofracturing within the BMV or due to a common feeder system for both mud volcanoes. This study shows that the mud volcanoes in the El Arraiche mud volcano field are long-lived features that erupt with intervals of several tens of thousands of years.

Description: Geological and hydrographic records contain evidence of substantial past variations in the oxygenation of the global ocean. Numerical models predicts a future decrease of marine oxygen levels under global warming. Using a global biogeochemical-climate model in which diapycnal mixing is parametrised as the sum of the regionally heterogenous tidal and homogenous background vertical mixing, we here show that simulated total oceanic oxygen content and the extent of marine suboxia are both sensitive to the strength of background vertical mixing. Eight otherwise identical configurations of the model were spun up under pre-industrial conditions for different vertical diffusivities ranging from background values of 0.01 cm 2/s to 0.5 cm 2/s. This range corresponds to various observational estimates and to values currently used in numerical ocean circulation models. Whereas the simulated total oceanic oxygen content is larger for larger mixing intensities, the simulated suboxic volume displays a maximum at intermediate diffusivities of about 0.2 cm 2/s. The intensity of vertical mixing also determines the evolution of suboxic areas under projected 21st century CO 2 emissions: while all model configurations predict a decline in total oceanic oxygen, the simulated extent of marine suboxia shows a 21st century expansion only for mixing rates higher than 0.2 cm 2/s, whereas the suboxic volume declines for lower mixing rates despite an overall loss of marine oxygen. Differences in the poorly constrained mixing parameterisation can thus lead to qualitatively different estimates about the future evolution of marine suboxia under projected climate change.

Description: Transient tracer measurements can constrain the rates and pathways of ocean ventilation and act as proxies for biogeochemically relevant gases such as CO2 and oxygen. Various techniques have deduced changes in ocean ventilation over decadal timescales using transient tracer measurements made on repeat sections, but these require a priori assumptions about mixing in the ocean interior. Here, we introduce a simple, direct observational method that takes advantage of the similar atmospheric increase rates of chlorofluorocarbon-12 and sulfur hexafluoride, but with a time lag (offset) of 1415 years. Such repeat measurements can be directly compared without prior assumptions about mixing. A difference larger than similar to 2 years between modern sulfur hexafluoride and historical chlorofluorocarbon-12 tracer ages implies a change in ventilation, although lack of difference does not necessarily imply no change. Several tracer data sets are presented, which suggest changes in ventilation in the South Pacific and North Atlantic Oceans.

Description: By comparison of the methane mixing ratio and the carbon isotope ratio (δ13CCH4) in Arctic air with regional background, the incremental input of CH4 in an air parcel and the source δ13CCH4 signature can be determined. Using this technique the bulk Arctic CH4 source signature of air arriving at Spitsbergen in late summer 2008 and 2009 was found to be −68‰, indicative of the dominance of a biogenic CH4 source. This is close to the source signature of CH4 emissions from boreal wetlands. In spring, when wetland was frozen, the CH4 source signature was more enriched in 13C at −53 ± 6‰ with air mass back trajectories indicating a large influence from gas field emissions in the Ob River region. Emissions of CH4 to the water column from the seabed on the Spitsbergen continental slope are occurring but none has yet been detected reaching the atmosphere. The measurements illustrate the significance of wetland emissions. Potentially, these may respond quickly and powerfully to meteorological variations and to sustained climate warming.

Description: Results from an interannually forced, 0.08 degrees eddy-resolving simulation based on the Hybrid Coordinate Ocean Model, in conjunction with a small but well-determined transport database, are used to investigate the currents and transports associated with the Atlantic meridional overturning circulation (AMOC) in the subpolar North Atlantic (SPNA). The model results yield a consistent warming in the western SPNA since the early 1990s, along with mean transports similar to those observed for the trans-basin AMOC across the World Ocean Circulation Experiment hydrographic section AR19 (16.4 Sv) and boundary currents at the exit of the Labrador Sea near 53 degrees N (39.0 Sv) and east of the Grand Banks near 43 degrees N (15.9 Sv). Over a 34 year integration, the model-determined AMOC across the AR19 section and the western boundary current near 53 degrees N both exhibit no systematic trend but some long-term (interannual and longer) variabilities, including a decadal transport variation of 3-4 Sv from relatively high in the 1990s to low in the 2000s. The decadal variability of the model boundary current transport near 53 degrees N lags the observed winter time North Atlantic Oscillation index by about 2 years and leads the model AMOC across the AR19 section by about 1 year. The model results also show that the long-term variabilities are low compared to those on shorter time scales. Thus, rapid sampling of the current over long time intervals is required to filter out high-frequency variabilities in order to determine the lower frequency variabilities of interest

Description: To investigate the mechanisms of the offset of seawater radiogenic Nd-Hf isotope compositions from those of the upper continental crust rocks, combined Nd-Hf isotope compositions of desert and loess samples from northern China (which integrate a wide range of lithologies and ages of continental rocks) are presented in this study. The results show significant and systematic fractionation of Hf isotopes between finegrained detritals/leachates (〈5 mm) and coarser fractions (〈75 mm) of the same samples. A small but systematic difference of Nd isotope compositions between leachates and detrital silicates is also revealed. Overall, the leaching data either plot along or slightly above the Nd-Hf seawater array, providing strong direct support that the seawater Nd-Hf isotope relationship is predominantly generated by weathering of upper continental crust. Our study supports the application of dissolved Hf isotopes as a proxy for different modes of weathering regimes rather than for continental source provenances.

Description: Large internal wave breaking is observed exceeding a vertical array of 61 high-resolution temperature sensors at 1 m intervals between 7 and 67 m above the bottom. The array was moored for 5 days at 969 m of Opouawe Bank, New Zealand, a known methane seep area. As breaking internal waves dominate sediment resuspension above sloping topography in other ocean areas, they are expected to also influence methane transport. Despite being visible in single beam echosounder data, indications for turbulence due to rising gas bubbles are not found in the present 1 Hz sampled temperature records. Likely, the mooring was too far away from the very localized bubble release spot. Instead, the temperature sensors show detailed internal wave-turbulence transitions. Every tidal cycle, a solibore (a frontal turbulent bore with a train of trailing solitary waves) changes shape and intensity. These solibores are highly turbulent and they restratify the bottom boundary layer, thereby maintaining efficient mixing. Details of different turbulent bore developments are discussed. Averaged over a few tidal cycles and over the sensors range, mean vertical eddy diffusivity amounts to 3 ± 1 × 10−3 m2 s−1 and mean turbulent kinetic energy dissipation to 1.6 ± 0.7 × 10−7 W kg−1, with variations over 4 orders of magnitude. Such turbulence will affect the distribution of dissolved methane and other geochemical species in the lower 100–150 m above the bottom and their release from the bottom. The above mean values are remarkably similar to those found at various other sites in the NE Atlantic Ocean.