A singular evolutive extended Kalman filter to assimilate ocean color data in a coupled physical–biochemical model of the North Atlantic ocean

Within the European DIADEM project, a data assimilation system for coupled ocean circulation and marine ecosystem models has been implemented for the North Atlantic and the Nordic Seas. One objective of this project is to demonstrate the relevance of sophisticated methods to assimilate satellite data such as altimetry, surface temperature and ocean color, into realistic ocean models. In this paper, the singular evolutive extended Kalman (SEEK) filter, which is an advanced assimilation scheme where three-dimensional, multivariate error statistics are taken into account, is used to assimilate ocean color data into the biological component of the coupled system. The marine ecosystem model, derived from the FDM model [J. Mar. Res. 48 (1990) 591], includes 11 nitrogen and carbon compartments and describes the synthesis of organic matter in the euphotic zone, its consumption by animals of upper trophic levels, and the recycling of detritic material in the deep ocean. The circulation model coupled to the ecosystem is the Miami isopycnic coordinate ocean model (MICOM), which covers the Atlantic and the Arctic Oceans with an enhanced resolution in the North Atlantic basin. The model is forced with realistic ECMWF ocean/atmosphere fluxes, which permits to resolve the seasonal variability of the circulation and mixed layer properties. In the twin assimimation experiments reported here, the predictions of the coupled model are corrected every 10 days using pseudo-measurements of surface phytoplankton as a substitute to chlorophyll concentrations measured from space. The diagnostics of these experiments indicate that the assimilation is feasible with a reduced-order Kalman filter of small rank (of order 10) as long as a sufficiently good identification of the error structure is available. In addition, the control of non-observed quantities such as zooplankton and nitrate concentrations is made possible, owing to the multivariate nature of the analysis scheme. However, a too severe truncation of the error sub-space downgrades the propagation of surface information below the mixed layer. The reduction of the actual state vector to the surface layers is therefore investigated to improve the estimation process in the perspective of sea-viewing wide field-of-view sensor (SeaWiFS) data assimilation experiments.     

Model simulations of carbon sequestration in the northwest Pacific by patch fertilization

Iron fertilization of nutrient-rich surface waters of the ocean is one possible way to help slow the rising levels of atmospheric CO₂ by sequestering it in the oceans via biological carbon export. Here, I use an ocean general circulation model to simulate a patch of nutrient depletion in the subpolar northwest Pacific under various scenarios. Model results confirm that surface fertilization is an inefficient way to sequester carbon from the atmosphere (Gnanadesikan et al., 2003), since only about 20% of the exported carbon comes initially from the atmosphere. Fertilization reduces future production and thus CO₂ uptake by utilizing nutrients that would otherwise be available later. Effectively, this can be considered as leakage when compared to a control run. This “effective” leakage and the actual leakage of sequestered CO₂ cause a significant, rapid decrease in carbon retention (only 30–45% retained after 10 years and less than 20% after 50 years). This contrasts markedly with the almost 100% retention efficiency for the same duration using the same model, when carbon is disposed directly into the northwest Pacific (Matsumoto and Mignone, 2005). As a consequence, the economic effectiveness of patch fertilization is poor in two limiting cases of the future price path of carbon. Sequestered carbon in patch fertilization is lost to the atmosphere at increasingly remote places as time passes, which would make monitoring exceedingly difficult. If all organic carbon from one-time fertilization reached the ocean bottom and remineralized there, acidification would be about −0.05 pH unit with O₂ depletion about −20 μmol kg⁻¹. These anomalies are probably too small to seriously threaten deep sea biota, but they are underestimated in the model because of its large grid size. The results from this study offer little to advocate purposeful surface fertilization as a serious means to address the anthropogenic carbon problem.     

Ocean Waves: Half-a-Century of Discovery

While the nature of most ocean waves has long been known and their basic physics understood since the nineteenth century, intense study of ocean waves during the second half of the twentieth century has taken the subject from the realm of mathematical exercises to that of practical engineering. Modern understanding of the generation, propagation and interactions of ocean waves with each other and with oceanic features has advanced to a quantitative level offering predictive capacity. This paper presents a brief qualitative review of advances in knowledge of sound waves, wind waves, tsunamis, tides, internal waves and long-period vorticity waves. The review is aimed at non-specialists who may benefit from an overview of the current state of the subject and access to a bibliography of general-interest references. This revised version was published online in August 2006 with corrections to the Cover Date.     

Topographic Effect of Izu Ridge on the Distribution of the North Pacific Intermediate Water South of Japan

The topographic effect of the Izu Ridge on the horizontal distribution of the North Pacific Intermediate Water (NPIW) south of Japan has been studied using observational data obtained by the Seisui-Maru of Mie University (Mie Univ. data) and those compiled by Japan Oceanographic Data Center (JODC data). Both data sets show that water of salinity less than 34.1 psu on potential density () surface of 26.8 is confined to the eastern side of the Izu Ridge, while water of salinity less than 34.2 psu is confined to the southern area over the Izu Ridge at a depth greater than 2000 m and to the southeastern area in the Shikoku Basin. It is also shown by T-S analysis of Mie Univ. data over the Izu Ridge that water of salinity less than 34.2 psu dominates south of 30°N, where the depth of the Izu Ridge is deeper than 2000 m and NPIW can intrude westward over the Izu Ridge. JODC data reveal that relatively large standard deviations of the salinity on surface of 26.7, 26.8 and 26.9 are detected along the mean current path of the Kuroshio and the Kuroshio Extension. Almost all of the standard deviations are less than 0.05 psu in other area with the NPIW, which shows that the time variation in the salinity can be neglected. This observational evidence shows that the topographic effect of the Izu Ridge on the horizontal distribution of the NPIW, which is formed east of 145°E by the mixing of the Kuroshio water and the Oyashio water, is prominent north of 30°N with a depth shallower than 2000 m.     

中太平洋CP25岩心的矿物、稀土元素及Sr、Nd同位素组成--晚新生代海底火山活动的证据

中太平洋ＣＰ２５岩心沉积物中钙十字沸石、蒙脱石／伊利石比值（Ｍ／Ｉ）、稀土元素含量、铈异常（δＣｅ）、铕异常（δＥｕ）及Ｓｒ、Ｎｄ同位素特征作为海底火山活动的证据,揭示了中太平洋海盆晚新生代以来火山活动和沉积的历史。其结果是,早中新世约２０～１７ＭａＢ．Ｐ．间,火山活动逐渐强烈,在１７ＭａＢ．Ｐ．最强烈,使得沉积物具有高Ｍ／Ｉ和∈Ｎｄ、低δＣｅ和８７Ｓｒ／８６Ｓｒ。１７～１ＭａＢ．Ｐ．发生沉积间断。１ＭａＢ．Ｐ．至现代,火山活动较早中新世减弱     

红河断裂构造带在东南亚的延伸特征及其大地构造意义

根据地质、地貌及地球物理资料分析,探讨展布于云南边陲哀牢山两侧的北西向断裂组成的红河断裂构造带,在东南亚的延伸特征。提出在该延伸带两侧沉积建造、构造活动,地球物理场和大地构造发展上均有显著差异。并进一步探讨该断裂构造带在大地构造上的意义。     

太平洋中部表层沉积物镁铝含量比的变化

本文根据中、日调查资料,研究了太平洋中部表层沉积物镁铝含量比的变化,并探讨其影响因素。     

Tectonics of a fast spreading center: A Deep-Tow and sea beam survey on the East Pacific rise at 19°30′ S

Sea Beam and Deep-Tow were used in a tectonic investigation of the fast-spreading (151 mm yr^-1) East Pacific Rise (EPR) at 19°30 S. Detailed surveys were conducted at the EPR axis and at the Brunhes/Matuyama magnetic reversal boundary, while four long traverses (the longest 96 km) surveyed the rise flanks. Faulting accounts for the vast majority of the relief. Both inward and outward facing fault scarps appear in almost equal numbers, and they form the horsts and grabens which compose the abyssal hills. This mechanism for abyssal hill formation differs from that observed at slow and intermediate spreading rates where abyssal hills are formed by back-tilted inward facing normal faults or by volcanic bow-forms. At 19°30 S, systematic back tilting of fault blocks is not observed, and volcanic constructional relief is a short wavelength signal (less than a few hundred meters) superimposed upon the dominant faulted structure (wavelength 2–8 km). Active faulting is confined to within approximately 5–8 km of the rise axis. In terms of frequency, more faulting occurs at fast spreading rates than at slow. The half extension rate due to faulting is 4.1 mm yr^-1 at 19°30 S versus 1.6 mm yr^-1 in the FAMOUS area on the Mid-Atlantic Ridge (MAR). Both spreading and horizontal extension are asymmetric at 19°30 S, and both are greater on the east flank of the rise axis. The fault density observed at 19°30 S is not constant, and zones with very high fault density follow zones with very little faulting. Three mechanisms are proposed which might account for these observations. In the first, faults are buried episodically by massive eruptions which flow more than 5–8 km from the spreading axis, beyond the outer boundary of the active fault zone. This is the least favored mechanism as there is no evidence that lavas which flow that far off axis are sufficiently thick to bury 50–150 m high fault scarps. In the second mechanism, the rate of faulting is reduced during major episodes of volcanism due to changes in the near axis thermal structure associated with swelling of the axial magma chamber. Thus the variation in fault spacing is caused by alternate episodes of faulting and volcanism. In the third mechanism, the rate of faulting may be constant (down to a time scale of decades), but the locus of faulting shifts relative to the axis. A master fault forms near the axis and takes up most of the strain release until the fault or fault set is transported into lithosphere which is sufficiently thick so that the faults become locked. At this point, the locus of faulting shifts to the thinnest, weakest lithosphere near the axis, and the cycle repeats.     

Instrumentation and Numerical Analysis of Cylindrical Diaphragm Wall Movement During Deep Excavation at Coastal Area

The lateral deflection of a cylindrical diaphragm wall and the associated ground movement induced by deep excavation are analyzed by performing site instrumentations and numerical analyses in the coastal area of Korea. Wall lateral deflection, rebar stress, and pore water pressure were measured and analyzed in eight directions. Variations of soil properties with the decrease of confining pressure are compared by performing various in situ tests before ad after excavation. To calculate the wall lateral deflection accurately, the effects of small strain nonlinearity, confining pressure, and the hysteresis loading/unloading loop developed during excavation are considered in the proposed numerical analysis. By comparing numerical results with measured ones, the importances of considering small strain nonlinearity and confining pressure reduction in the nonlinear (FEM) are emphasized. Also, the effects of wall stiffness on the performance of cylindrical diaphragm walls are studied for future similar excavation in the coastal area.     

The Aleutian Low, storm tracks, and winter climate variability in the Bering Sea

Previous studies have found inconsistent results regarding how wintertime conditions in the Bering Sea relate to variations in the North Pacific climate system. This problem is addressed through analysis of data from the NCEP/NCAR Reanalysis for the period 1950–2003. Composite patterns of sea-level pressure, 500 hPa geopotential heights, storm tracks and surface air temperature are presented for four situations: periods of strong Aleutian Low, weak Aleutian Low, warm Bering Sea air temperatures, and cold Bering Sea air temperatures. Winter temperatures in the Bering Sea are only marginally related to the strength of the Aleutian Low, and are much more sensitive to the position of the Aleutian Low and to variations in storm tracks. In particular, relatively warm temperatures are associated with either an enhanced storm track off the coast of Siberia, and hence anomalous southerly low-level flow, or an enhanced storm track entering the eastern Bering Sea from the southeast. These latter storms do not systematically affect the mean meridional winds, but rather serve to transport mild air of maritime origin over the Bering Sea. The leading indices for the North Pacific, such as the NP and PNA, are more representative of the patterns of tropospheric circulation and storm track anomalies associated with the strength of the Aleutian Low than patterns associated with warm and cold wintertime conditions in the Bering Sea.