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.     

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

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

大西洋中脊热液硫化物的矿物学研究

利用１９８９年在大西洋热液活动区采集的块状硫化物样品，采用矿相显微镜、电子显微镜和电子探针等手段，进行了矿物共生组合、矿物形态与成分标型及其演化特征研究。结果表明，大西洋热液成矿作用可以划分为热液期与沉积期两个成矿期，和石英－黄铁矿阶段、黄铁矿阶段、多金属硫化物阶段、胶状黄铁矿阶段和非晶质ＳｉＯ２阶段等５个成矿阶段；不同成矿阶段黄铁矿的共生矿物和矿物特征不同，在晶体形态上具有从单形晶－聚形晶－胶状     

杰氏风暴潮预报方法在珠海及粤西海域的应用

本文通过杰氏风暴潮预报方法对历史上给珠海和粤西海域造成较大影响的台风增水事件进行后报试验,并通过统计分析方法对原有的预报公式进行订正,总结出适合珠海及粤西海域的风暴潮预报公式.利用5a的时间对珠海及粤西海域的台风增水进行试报,结果表明预报准确率分别达95.0%和87.6%,尤其对一些增水较大、影响较严重的台风增水预报相对更加准确,预报准确率均达90.0%以上.这说明订正后的公式对珠海及粤西海域的风暴潮预报能得到较好的预报效果.     

Estimating the geostrophic velocity obtained by HF radar observations in the upstream area of the Kuroshio

A method to extract geostrophic current in the daily mean HF radar data in the Kuroshio upstream region is established by comparison with geostrophic velocity determined from the along-track altimetry data. The estimated Ekman current in the HF velocity is 1.2% (1.5%) and 48° (38°)-clockwise rotated with respect to the daily mean wind in (outside) the Kuroshio. Furthermore, additional temporal smoothing is found necessary to remove residual ageostrophic currents such as the inertial oscillation. After removal of the ageostrophic components, the HF geostrophic velocity agrees well with that from the altimetry data with rms difference 0.14 (0.12) m/s in (outside) the Kuroshio.     

In situ and satellite observations of the eastward migration of the Western Alboran Sea Gyre

During October 2003 an intensive oceanographic survey (BIOMEGA) was carried out in the Alboran Sea, coinciding with a migration event of the Western Alboran Sea Gyre (WAG). The observations gathered during that cruise constitute the first field evidence of a migrated stage of the WAG. In this work we present the main differences between the 3D hydrodynamic fields observed during BIOMEGA and those corresponding to a WAG located at its usual position. The migration of the gyre was followed by satellite (altimetry and sea surface temperature) imagery. The causes of the gyre migration are explored in terms of the quasi-geostrophic tendency equation, in particular of the dynamics governing scales larger than the Rossby radius of deformation. It is shown that the steady state gyre must be almost equivalent barotropic and that the key process to break down the stationarity would be a density advection at gyre scale. The mechanisms to explain the migration of the WAG proposed by previous authors are discussed in light of the explanation proposed in this work.     

The POC / Th ratio of settling particles isolated using split flow-thin cell fractionation (SPLITT)

The common assumption that the ratio between particulate organic carbon (POC) and particulate ²³⁴Th obtained from shallow sediment traps and filterable particles are representative of the ratio in the total particle settling flux should be treated with caution in view of well-known biases associated with tethered shallow sediment traps and the decoupling between size and settling velocity of many natural particle regimes. To make progress toward reliably constraining the POC / ²³⁴Th ratio on truly settling particles, we have tested here a settling collection technique designed to remove any hydrodynamic bias; split flow-thin cell fractionation (SPLITT). These first results from a North Sea fjord and an open Baltic Sea time-series station indicates that the POC / ²³⁴Th ratio on the more complete particle-settling spectrum, isolated with SPLITT, was higher than the POC / ²³⁴Th ratio obtained simultaneously from tethered shallow sediment traps in seven out of seven parallel deployments with an average factor of 210%. The POC / ²³⁴Th ratio from the SPLITT was either in the same range or higher than that obtained on filtered “bulk” particles. To explain this novel data we hypothesize that the slowest settling fraction is organic-matter rich and does not strongly complex ²³⁴Th (i.e., high POC / ²³⁴Th). We suggest that this ultra-slow sinking fraction is better collected by SPLITT than with tethered sediment traps because of minimized hydrodynamic bias.This was tested using the ratio of POC / Al as a tracer of detrital mineral-ballast influenced settling velocity. The higher POC / Al ratios in SPLITT samples relative to in traps is consistent with the hypothesis that SPLITT is better suited for collecting also the slow-settling component of sinking particles. This important slow-settling component appears to here consist primarily of non-APS/TEP components of plankton exudates or other less-strongly ²³⁴Th-complexing organic matter. Further applications of the SPLITT technique are likely to return increasingly new insights on the composition (including “truly settling” POC / ²³⁴Th) of the total spectrum of particles settling out of the upper ocean.     

融冰季节北极破碎冰区热通量的初步研究

利用航空遥感数字影像的解析结果和实测气象，海洋和海冰资料，定量研究了夏季融冰期北极破碎冰区的热通量，计算了海洋对大气的热贡献，结果表明，在北极夏季海冰融化时，短波辐射远远大于感热和潜热通量，是表面热通量的决定因素，海洋对大气的热贡献主要由长波辐射决定，在观测期间，海洋对大气的热贡献为38～104Wm^-2，这部分热量的大小与海冰的密集度有关，当海冰密集度小于0.8时，海洋对大气的热贡献随海冰密度度的增大而减小，而当海冰密集度超过0.8以后，该热通量将随海冰密集度的增大而增大。