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.     

洋浦近岸海域污染控制规划方案优化——Ⅰ.排污口位置优选

在洋浦近岸海域计算潮流场基础上,计算预选排污口附近海水质点运动轨迹及预测污染物浓度分布,最后,从环保角度出发,推荐了洋浦地区各开发区的排污口位置,为洋浦近岸海域污染控制规划方案的制定及其优化提供了依据。     

环台湾岛海域半日潮波特征的三维模拟

用1997版POM海洋模式,首次应用于环台湾岛海域的潮波数值研究.得到该海域的半日潮波主要为23°N以南西太平洋传来的胁振潮.影响台湾海峡的半日潮波分别由海峡南北口传入的两支潮波,且北支强于南支.福建沿岸湄州湾-兴化湾为最强潮区,其M₂分潮最大振幅可达240cm.最强潮流区位于澎湖水道,M₂分潮最大潮流达196cm/s.环台湾岛海域潮波潮流水平结构上除海峡北部原有一个圆流点外,还发现另外存在4个新的圆流点.潮流垂直结构上主要为右偏,接近底层处为左偏.     

南海东北部台风潮的数值后报及极值增水估计

本研究对近40年来影响南海东北部陆架海区的28次台风引起的风暴潮进行了数值后报,其中8个过程的沿岸后报增水值与实测值进行了比较,表明后报值与实测值符合良好,90%以上的最大增水值偏差在30厘米以下。为了得出本海区多年一遇的台风增水极值,在后报台风路径密集处选择了9个不同水深点,对每点取出各次台风下的最大增水值,然后用Weibull分布进行拟合,得出了各点的极值分布。     

A Numerical Model for Natural Backfill of Pipeline Trenches Subjected to Unidirectional／Oscillatory Flows

A numerical model for the self-burial of a pipeline trench is developed. Morphological evolutions of a pipeline trench under steady-current or oscillatory-flow conditions are simulated with/without a pipeline inside the trench. The oscillatory flow in this study represents the action of waves. The two-dimensional Reynolds-averaged continuity and Navier-Stokes equations with the standard k-e turbulence closure, as well as the sediment transport equations, are solved with the finite difference method in a eurvilirrear coordinate system. Both bed and suspended loads of sediment transport are included in the morphological model. Because of the lack of experimental data on the backfilling of pipeline trenches, the numericalmodel is firstly verified against three closely-relevant experiments available in literature. A detailed measurement of the channel migration phenomenon under steady currents is employed for the assessment of the integral performance of themodel. The two experimental results from U-tube tests are used to validate the model‘s ability in predicting oscillatory flows. Different time-marching schemes are employed for the morphological computation under unidirectional and oscillatory conditions. It is found that vortex motions within the trench play an important role in the trench development.     

Effect of the Twin Jetties as A Regulating Structure at the Tao-er Estuary

The phase one project of the twin jetties at the Tao-er Estuary yields no effect on the flow, and there is a large area of deposition. A numerical model is employed to simulate the flow field, and the scouring and silting in the channel be-tween the two jetties are analyzed. The results show that the effective tidal flux between the two jetties is reduced due to the notch located at the heel of the east jetty, and the effect of contraction-flush is not obvious. As for the regulation of this kind of estuary under the tide action, the capability of tidal flux should be fully used.     

水平井含水率上升影响因素

利用塔里木油田塔中4油田(TZ4)底水油藏相关的地质、流体数据建立数值模型。在所建模型的基础上,应用数值模拟计算的累积产油、产水和产液量回归俞启泰水驱特征曲线,以反映水平井见水特征的参数b。以参数b为研究对象,采用正交试验的方法研究不同因素对b值的影响,筛选影响水平井见水特征的主要因素,认为原油粘度、油层厚度、非均质性及水平井水平段在油藏中的位置是影响水平井含水上升趋势的主要指标。最后建立主要因素与b值的关系式,结合俞启泰曲线b值图版,提出预测水平井见水规律的公式——图版法(F-b法)。应用F-b法对塔里木水平井含水率进行预测,并与其他相关方法及实际生产数据对比,认为F-b法可作为预测塔里木油田水平井含水率、估算可采储量的一种有效方法。     

On the multivariate total least-squares approach to empirical coordinate transformations. Three algorithms

The multivariate total least-squares (MTLS) approach aims at estimating a matrix of parameters, Ξ, from a linear model (Y−E _Y = (X−E _X ) · Ξ) that includes an observation matrix, Y, another observation matrix, X, and matrices of randomly distributed errors, E _Y and E _X . Two special cases of the MTLS approach include the standard multivariate least-squares approach where only the observation matrix, Y, is perturbed by random errors and, on the other hand, the data least-squares approach where only the coefficient matrix X is affected by random errors. In a previous contribution, the authors derived an iterative algorithm to solve the MTLS problem by using the nonlinear Euler–Lagrange conditions. In this contribution, new lemmas are developed to analyze the iterative algorithm, modify it, and compare it with a new ‘closed form’ solution that is based on the singular-value decomposition. For an application, the total least-squares approach is used to estimate the affine transformation parameters that convert cadastral data from the old to the new Israeli datum. Technical aspects of this approach, such as scaling the data and fixing the columns in the coefficient matrix are investigated. This case study illuminates the issue of “symmetry” in the treatment of two sets of coordinates for identical point fields, a topic that had already been emphasized by Teunissen (1989, Festschrift to Torben Krarup, Geodetic Institute Bull no. 58, Copenhagen, Denmark, pp 335–342). The differences between the standard least-squares and the TLS approach are analyzed in terms of the estimated variance component and a first-order approximation of the dispersion matrix of the estimated parameters.     

波流共存场中多向随机波浪传播变形数学模型

基于波作用量守恒方程建立了波流共存场中多向随机波浪传播变形数学模型,模型中考虑了波浪绕射的影响和水流引起的波浪弥散多普勒效应,应用包含水流和地形影响的激破波模式计算波浪破碎的能量耗散,采用一阶上迎风有限差分格式离散控制方程。分别计算了有无近岸流情况下单向和多向随机波浪的波高分布,考虑水流影响的数值计算结果与物理模型实验数据吻合良好,比较分析表明,所建立的数学模型能够复演由于离岸流引起的波高增大,可用于波流共存场多向随机波浪传播变形的模拟和预报。