He, Ne and Ar isotopic composition of Fe-Mn crusts from the western and central Pacific Ocean and implications for their genesis

The noble gas nuclide abundances and isotopic ratios of the upmost layer of Fe-Mn crusts from the western and central Pacific Ocean have been determined. The results indicate that the He and Ar nuclide abundances and isotopic ratios can be classified into two types: low He-3/He-4 type and high He-3/He-4 type. The low He-3/He-4 type is characterized by high He-4 abundances of 191x10(-9) cm(3.)STP(.)g(-1) on average, with variable He-4, Ne-20 and Ar-40 abundances in the range (42.8-421)x10(-9) cm(3.)STP(.)g(-1), (5.40-141)x10(-9)cm(3.)STP(.)g(-1), and (773-10976)x10(-9) cm(3.)STP(.)g(-1), respectively. The high He-3/He-4 samples are characterized by low He-4 abundances of 11.7x10(-9) cm(3.)STP(.)g(-1) on average, with He-4, Ne-20 and Ar-40 abundances in the range of (7.57-17.4)x10(-9) cm(3.)STP(.)g(-1), (110.4-25.5)x10(-9) cm(3.)STP(.)g(-1) and (5354-9050)x10(-9) cm(3.)STP(.)g(-1), respectively. The low He-3/He-4 samples have He-3/He-4 ratios (with RIRA ratios of 2.04-2.92) which are lower than those of MORB (R/R-A=8 +/- 1) and Ar-40/Ar-36 ratios (447-543) which are higher than those of air (295.5). The high He-3/He-4 samples have He-3/He-4 ratios (with R/R-A ratios of 10.4-12.0) slightly higher than those of MORB (R/R-A=8 +/- 1) and Ar-40/Ar-36 ratios (293-299) very similar to those of air (295.5). The Ne isotopic ratios (Ne-20/Ne-22 and Ne-21/Ne-22 ratios of 10.3-10.9 and 0.02774-0.03039, respectively) and the Ar-38/Ar-36 ratios (0.1886-0.1963) have narrow ranges which are very similar to those of air (the Ne-20/Ne-22, Ne-21/Ne-22, Ar-38/Ar-36 ratios of 9.80, 0.029 and 0.187, respectively), and cannot be differentiated into different groups. The noble gas nuclide abundances and isotopic ratios, together with their regional variability, suggest that the noble gases in the Fe-Mn crusts originate primarily from the lower mantle. The low He-3/He-4 type and high He-3/He-4 type samples have noble gas characteristics similar to those of HIMU (High U/Pb Mantle)- and EM (Enriched Mantle)-type mantle material, respectively. The low He-3/He-4 type samples with HIMU-type noble gas isotopic ratios occur in the Magellan Seamounts, Marcus-Wake Seamounts, Marshall Island Chain and the Mid-Pacific Seamounts whereas the high He-3/He-4 type samples with EM-type noble gas isotopic ratios occur in the Line Island Chain. This difference in noble gas characteristics of these crust types implies that the Magellan Seamounts, Marcus-Wake Seamounts, Marshall Island Chain, and the Mid-Pacific Seamounts originated from HIMU-type lower mantle material whereas the Line Island Chain originated from EM-type lower mantle material. This finding is consistent with variations in the Pb-isotope and trace element signatures in the seamount lavas. Differences in the mantle surce may therefore be responsible for variations in the noble gas abundances and isotopic ratios in the Fe-Mn crusts. Mantle degassing appears to be the principal factor controlling noble gas isotopic abundances in Fe-Mn crusts. Decay of radioactive isotopes has a negligible influence on the nuclide abundances and isotopic ratios of noble gases in these crusts on the timescale of their formation.

Genesis of Th-230 excess in basalts from mid-ocean ridges and ocean islands: Constraints from the global U-series isotope database and major and rare earth element geochemistry

Based on Th-230-U-238 disequilibrium and major element data from mid-ocean ridge basalts (MORBs) and ocean island basalts (OIBs), this study calculates mantle melting parameters, and thereby investigates the origin of Th-230 excess. (Th-230/U-238) in global MORBs shows a positive correlation with Fe-8, P (o), Na-8, and F-melt (Fe-8 and Na-8 are FeO and Na2O contents respectively after correction for crustal fractionation relative to MgO = 8 wt%, P (o)=pressure of initial melting and F (melt)=degree of melt), while Th-230 excess in OIBs has no obvious correlation with either initial mantle melting depth or the average degree of mantle melting. Furthermore, compared with the MORBs, higher (Th-230/U-238) in OIBs actually corresponds to a lower melting degree. This suggests that the Th-230 excess in MORBs is controlled by mantle melting conditions, while the Th-230 excess in OIBs is more likely related to the deep garnet control. The vast majority of calculated initial melting pressures of MORBs with excess Th-230 are between 1.0 and 2.5 GPa, which is consistent with the conclusion from experiments in recent years that D (U)> D (Th) for Al-clinopyroxene at pressures of > 1.0 GPa. The initial melting pressure of OIBs is 2.2-3.5 GPa (around the spinel-garnet transition zone), with their low excess Ra-226 compared to MORBs also suggesting a deeper mantle source. Accordingly, excess Th-230 in MORBs and OIBs may be formed respectively in the spinel and garnet stability field. In addition, there is no obvious correlation of K2O/TiO2 with (Th-230/U-238) and initial melting pressure (P (o)) of MORBs, so it is proposed that the melting depth producing excess Th-230 does not tap the spinel-garnet transition zone. OIBs and MORBs in both (Th-230/U-238) vs. K2O/TiO2 and (Th-230/U-238) vs. P (o) plots fall in two distinct areas, indicating that the mineral phases which dominate their excess Th-230 are different. Ce/Yb-Ce curves of fast and slow ridge MORBs are similar, while, in comparison, the Ce/Yb-Ce curve for OIBs shows more influence from garnet. The mechanisms generating excess Th-230 in MORBs and OIBs are significantly different, with formation of excess Th-230 in the garnet zone only being suitable for OIBs.

The variation of upper ocean structure and paleoproductivity in the Kuroshio source region during the last 200 kyr

A core from the source region of the Kuroshio warm current (east of the Luzon Island) was analyzed using several proxies in order to study the variability of the Western Pacific Warm Pool (WPWP) during the last two glacial-interglacial cycles. Primary productivity (PP) variations were deduced from variations in the coccolith flora. Primary productivity was higher during glacial periods (the end of Marine Isotope Stage [MIS] 3, some periods in MIS 2 and 6), and decreased during interglacial periods (MIS 7, MIS Se and probably MIS 5c-5d), with the lowest PP in MIS 5e. variations in the delta C-13 difference in benthic and bulk carbonate, thus in the vertical gradient of delta C-13 in dissolved inorganic carbon (Delta delta C-13(c). (wuellerstorfi-N. dutertrei) and Delta delta C-13(c.) (wuellerstorfi-coccolith)) Coincided With the PP Changes, showing that export productivity was low during interglacial periods (MIS 7, MIS 5e and Holocene) and high during glacial periods (MIS 6, probably MIS 5c-5d, late MIS 4 and late MIS 3). Comparison of foraminiferal carbonate dissolution indicators and PP changes reveals that nannofossil assemblage in core Ph05-5 is not sensitive to carbonate dissolution intensity. The depth of the thermocline (DOT) was estimated from planktonic forminiferal assemblages, and was relatively greater during interglacial periods (MIS 7, MIS 5e, probably MIS 5c and Holocene) than during glacials (middle MIS 6, probably MIS 5b and 5d, some periods in MIS 4, MIS 3 and MIS 2). Good coherence between the paleoproductivity records and the DOT suggests that the DOT changes could be the primary control factor in changes of paleoproductivity, and the glacial high productivity in the Kuroshio source region could be associated with a global increase of nutrient concentration in the intermediate waters that upwelled into the photic zone. The low CO2 values derived for intervals of high productivity and a relatively shallow DOT suggest that the changes in biological productivity and DOT in the equatorial Pacific could have modified atmospheric CO2 concentrations. High Sea Surface Temperatures (SSTs) during the warm MIS 5e in combination with intensified monsoonal rain fall could have resulted in a more intense stratification of the upper waters, resulting in low nutrient supply to the surface waters and a resulting decrease in productivity. (C) 2010 Elsevier B.V. All rights reserved.

Living planktonic foraminifera from the northwestern Pacific Ocean: SSU rDNA sequences and their paleoceanographic significance

Living planktonic foraminifera (PF) samples from the Okinawa Trough of the northwestern Pacific Ocean were taken for DNA analysis. The SSU rDNA sequences of two PF species, Globigerina sp. and Pulleniatina obliquiloculata collected at Station WP01, were obtained and compared with those from the southwestern Pacific Ocean. Only small differences (< 0.7%-1.2% for P. obliquiloculata, and 0.3% for Globigerina sp.) were found between samples from the north- and south-western Pacific Ocean areas and this molecular evidence supported that these micropaleontological species are the same species, which implies that the West Pacific Ocean circulation system influences the planktonic foraminiferal gene communication.

Effect of low-frequency Rossby wave on thermal structure of the upper southwestern tropical Indian Ocean

We investigate the influence of low-frequency Rossby waves on the thermal structure of the upper southwestern tropical Indian Ocean (SWTIO) using Argo profiles, satellite altimetric data, sea surface temperature, wind field data and the theory of linear vertical normal mode decomposition. Our results show that the SWTIO is generally dominated by the first baroclinic mode motion. As strong downwelling Rossby waves reach the SWTIO, the contribution of the second baroclinic mode motion in this region can be increased mainly because of the reduction in the vertical stratification of the upper layer above thermocline, and the enhancement in the vertical stratification of the lower layer under thermocline also contributes to it. The vertical displacement of each isothermal is enlarged and the thermal structure of the upper level is modulated, which is indicative of strong vertical mixing. However, the cold Rossby waves increase the vertical stratification of the upper level, restricting the variability related to the second baroclinic mode. On the other hand, during decaying phase of warm Rossby waves, Ekman upwelling and advection processes associated with the surface cyclonic wind circulation can restrain the downwelling processes, carrying the relatively colder water to the near-surface, which results in an out-of-phase phenomenon between sea surface temperature anomaly (SSTA) and sea surface height anomaly (SSHA) in the SWTIO.

Two modes of dipole events in tropical Indian Ocean

By analyzing the distributions of subsurface temperature and the surface wind stress anomalies in the tropical Pacific and Indian Oceans during the Indian Ocean Dipole (IOD) events, two major modes of the IOD and their formation mechanisms are revealed. (1) The subsurface temperature anomaly (STA) in the tropical Indian Ocean during the IOD events can be described as a "<" -shaped and west-east-oriented dipole pattern; in the east side of the "<" pattern, a notable tongue-like STA extends westward along the equator in the tropical eastern Indian Ocean; while in the west side of the "<" pattern, the STA has opposite sign with two centers (the southern one is stronger than the northern one in intensity) being of rough symmetry about the equator in the tropical mid-western Indian Ocean. (2) The IOD events are composed of two modes, which have similar spatial pattern but different temporal variabilities due to the large scale air-sea interactions within two independent systems. The first mode of the IOD event originates from the air-sea interaction on a scale of the tropical Pacific-Indian Ocean and coexists with ENSO. The second mode originates from the air-sea interaction on a scale of the tropical Indian Ocean and is closely associated with changes in the position and intensity of the Mascarene high pressure. The strong IOD event occurs when the two modes are in phase, and the IOD event weakens or disappears when the two modes are out of phase. Besides, the IOD events are normally strong when either of the two modes is strong. (3) The IOD event is caused by the abnormal wind stress forcing over the tropical Indian Ocean, which results in vertical transports, leading to the upwelling and pileup of seawater. This is the main dynamic processes resulting in the STA. When the anomalous easterly exists over the equatorial Indian Ocean, the cold waters upwell in the tropical eastern Indian Ocean while the warm waters pileup in the tropical western Indian Ocean, hence the thermocline in the tropical Indian Ocean is shallowed in the east and deepened in the west. The off-equator component due to the Coriolis force in the equatorial area causes the upwelling of cold waters and the shallowing of the equatorial India Ocean thermocline. On the other hand, the anomalous anticyclonic circulations and their curl fields located on both sides of the equator, cause the pileup of warm waters in the central area of their curl fields and the deepening of the equatorial Indian Ocean thermocline off the equator. The above three factors lead to the occurrence of positive phase IOD events. When anomalous westerly dominates over the tropical Indian Ocean, the dynamic processes are reversed, and the negative-phase IOD event occurs.

Neural network retrieval of ocean surface parameters from SSM/I data

A new algorithm based on the multiparameter neural network is proposed to retrieve wind speed (WS), sea surface temperature (SST), sea surface air temperature, and relative humidity ( RH) simultaneously over the global oceans from Special Sensor Microwave Imager (SSM/I) observations. The retrieved geophysical parameters are used to estimate the surface latent heat flux and sensible heat flux using a bulk method over the global oceans. The neural network is trained and validated with the matchups of SSM/I overpasses and National Data Buoy Center buoys under both clear and cloudy weather conditions. In addition, the data acquired by the 85.5-GHz channels of SSM/I are used as the input variables of the neural network to improve its performance. The root-mean-square (rms) errors between the estimated WS, SST, sea surface air temperature, and RH from SSM/I observations and the buoy measurements are 1.48 m s(-1), 1.54 degrees C, 1.47 degrees C, and 7.85, respectively. The rms errors between the estimated latent and sensible heat fluxes from SSM/I observations and the Xisha Island ( in the South China Sea) measurements are 3.21 and 30.54 W m(-2), whereas those between the SSM/ I estimates and the buoy data are 4.9 and 37.85 W m(-2), respectively. Both of these errors ( those for WS, SST, and sea surface air temperature, in particular) are smaller than those by previous retrieval algorithms of SSM/ I observations over the global oceans. Unlike previous methods, the present algorithm is capable of producing near-real-time estimates of surface latent and sensible heat fluxes for the global oceans from SSM/I data.

Long-Wave Dynamics of Sea Level Variations during Indian Ocean Dipole Events

Long-wave dynamics of the interannual variations of the equatorial Indian Ocean circulation are studied using an ocean general circulation model forced by the assimilated surface winds and heat flux of the European Centre for Medium-Range Weather Forecasts. The simulation has reproduced the sea level anomalies of the Ocean Topography Experiment (TOPEX)/Poseidon altimeter observations well. The equatorial Kelvin and Rossby waves decomposed from the model simulation show that western boundary reflections provide important negative feedbacks to the evolution of the upwelling currents off the Java coast during Indian Ocean dipole (IOD) events. Two downwelling Kelvin wave pulses are generated at the western boundary during IOD events: the first is reflected from the equatorial Rossby waves and the second from the off-equatorial Rossby waves in the southern Indian Ocean. The upwelling in the eastern basin during the 1997-98 IOD event is weakened by the first Kelvin wave pulse and terminated by the second. In comparison, the upwelling during the 1994 IOD event is terminated by the first Kelvin wave pulse because the southeasterly winds off the Java coast are weak at the end of 1994. The atmospheric intraseasonal forcing, which plays an important role in inducing Java upwelling during the early stage of an IOD event, is found to play a minor role in terminating the upwelling off the Java coast because the intraseasonal winds are either weak or absent during the IOD mature phase. The equatorial wave analyses suggest that the upwelling off the Java coast during IOD events is terminated primarily by western boundary reflections.

Remote sensing of ocean waves by along-track interferometric synthetic aperture radar

A new wave retrieval method for the Along-Track Interferometric Synthetic Aperture Radar (AT-InSAR) phase image is presented. The new algorithm, named parametric retrieval algorithm (PRA), uses the full nonlinear mapping relations. It differs from previous retrieval algorithms in that it does not require a priori information about the sea state or the wind vector from scatterometer data. Instead, it combines the observed AT-InSAR phase spectrum and assumed wind vector to estimate the wind sea spectrum. The method has been validated using several C-band and X-band HH-polarized AT-InSAR observations collocated with spectral buoy measurements. In this paper, X-band and C-band HH-polarized AT-InSAR phase images of ocean waves are first used to study AT-InSAR wave imaging fidelity. The resulting phase spectra are quantitatively compared with forward-mapped in situ directional wave spectra collocated with the AT-InSAR observations. Subsequently, we combine the parametric retrieval algorithm (PRA) with X-band and C-band HH-polarized AT-InSAR phase images to retrieve ocean wave spectra. The results show that the ocean wavelengths, wave directions, and significant wave heights estimated from the retrieved ocean wave spectra are in agreement with the buoy measurements.

Validation of RADARSAT-2 fully polarimetric SAR measurements of ocean surface waves

C band RADARSAT-2 fully polarimetric (fine quad-polarization mode, HH+VV+HV+VH) synthetic aperture radar (SAR) images are used to validate ocean surface waves measurements using the polarimetric SAR wave retrieval algorithm, without estimating the complex hydrodynamic modulation transfer function, even under large radar incidence angles. The linearly polarized radar backscatter cross sections (RBCS) are first calculated with the copolarization (HH, VV) and cross-polarization (HV, VH) RBCS and the polarization orientation angle. Subsequently, in the azimuth direction, the vertically and linearly polarized RBCS are used to measure the wave slopes. In the range direction, we combine horizontally and vertically polarized RBCS to estimate wave slopes. Taken together, wave slope spectra can be derived using estimated wave slopes in azimuth and range directions. Wave parameters extracted from the resultant wave slope spectra are validated with colocated National Data Buoy Center (NDBC) buoy measurements (wave periods, wavelengths, wave directions, and significant wave heights) and are shown to be in good agreement.

Distribution of the nonlinear random ocean wave period

Because of the intrinsic difficulty in determining distributions for wave periods, previous studies on wave period distribution models have not taken nonlinearity into account and have not performed well in terms of describing and statistically analyzing the probability density distribution of ocean waves. In this study, a statistical model of random waves is developed using Stokes wave theory of water wave dynamics. In addition, a new nonlinear probability distribution function for the wave period is presented with the parameters of spectral density width and nonlinear wave steepness, which is more reasonable as a physical mechanism. The magnitude of wave steepness determines the intensity of the nonlinear effect, while the spectral width only changes the energy distribution. The wave steepness is found to be an important parameter in terms of not only dynamics but also statistics. The value of wave steepness reflects the degree that the wave period distribution skews from the Cauchy distribution, and it also describes the variation in the distribution function, which resembles that of the wave surface elevation distribution and wave height distribution. We found that the distribution curves skew leftward and upward as the wave steepness increases. The wave period observations for the SZFII-1 buoy, made off the coast of Weihai (37A degrees 27.6' N, 122A degrees 15.1' E), China, are used to verify the new distribution. The coefficient of the correlation between the new distribution and the buoy data at different spectral widths (nu=0.3-0.5) is within the range of 0.968 6 to 0.991 7. In addition, the Longuet-Higgins (1975) and Sun (1988) distributions and the new distribution presented in this work are compared. The validations and comparisons indicate that the new nonlinear probability density distribution fits the buoy measurements better than the Longuet-Higgins and Sun distributions do. We believe that adoption of the new wave period distribution would improve traditional statistical wave theory.

Modeling vertical heat transport in the wave affected surface layer of the ocean

In considering the vertical heat transport problems in the upper ocean, the flat upper boundary approximation for the free surface and the horizontal homogenous hypothesis are usually applied. However, due to the existence of the wave motion, the application of this approximation may result in some errors to the solar irradiation since it decays quickly in respect to the actual thickness of the water layer below the surface; on the other hand, due to the fluctuation of the water layer depth, it is improper to neglect the effects of the horizontal advection and turbulent diffusion since they also contribute to the vertical heat transport. A new model is constructed in this study to reflect these effects. The corresponding numerical simulations show that the wave motion may remarkably accelerate the vertical heat transferring process and the variation of the temperature in the wave affected layer appears in an oscillating manner.

Self-organized Criticality Model for Ocean Internal Waves

In this paper, we present a simple spring-block model for ocean internal waves based on the self-organized criticality (SOC). The oscillations of the water blocks in the model display power-law behavior with an exponent of -2 in the frequency domain, which is similar to the current and sea water temperature spectra in the actual ocean and the universal Garrett and Munk deep ocean internal wave model [Geophysical Fluid Dynamics 2(1972) 225; J. Geophys. REs. 80 (1975) 291]. The influence of the ratio of the driving force to the spring coefficient to SOC behaviors in the model is also discussed.

Interannual horizontal heat advection in the surface mixed layer over the equatorial Pacific Ocean: assimilation versus TAO analysis

The interannual anomalies of horizontal heat advection in the surface mixed layer over the equatorial Pacific Ocean in an assimilation experiment are studied and compared with existing observational analyses. The assimilation builds upon a hindcast study that has produced a good simulation of the observed equatorial currents and optimizes the simulation of the Reynolds sea surface temperature (SST) data. The comparison suggests that the assimilation has improved the simulation of the interannual horizontal heat advection of the surface mixed layer significantly. During periods of interrupted current measurements, the assimilation is shown to produce more meaningful anomalies of the heat advection than the interpolation of the observational data does. The assimilation also shows that the eddy heat flux due to the correlation between high-frequency current and SST variations, which is largely overlooked by the existing observational analyses, is important for the interannual SST balance over the equatorial Pacific. The interannual horizontal heat advection anomalies are found to be sensitive to SST errors where oceanic currents are strong, which is a challenge for ENSO prediction. The study further suggests that the observational analyses of the tropical SST balance based on the TAO and the Reynolds SST data contain significant errors due to the large gradient errors in the Reynolds SST data, which are amplified into the advection anomalies by the large equatorial currents.

Variations in the eastern Indian Ocean warm pool and its relation to the dipole in the tropical Indian Ocean

Based on the monthly average SST and 850 hPa monthly average wind data, the seasonal, interannual and long-term variations in the eastern Indian Ocean warm pool (EIWP) and its relationship to the Indian Ocean Dipole (IOD), and its response to the wind over the Indian Ocean are analyzed in this study. The results show that the distribution range, boundary and area of the EIWP exhibited obviously seasonal and interannual variations associated with the ENSO cycles. Further analysis suggests that the EIWP had obvious long-term trend in its bound edge and area, which indicated the EIWP migrated westwards by about 14 longitudes for its west edge, southwards by about 5 latitudes for its south edge and increased by 3.52x10(6) km(2) for its area, respectively, from 1950 to 2002. The correlation and composite analyses show that the anomalous westward and northward displacements of the EIWP caused by the easterly wind anomaly and the southerly wind anomaly over the eastern equatorial Indian Ocean played an important and direct role in the formation of the IOD.