Present-day stress and pore pressure fields in the Cuu Long and Nam Con Son Basins, offshore Vietnam

Knowledge of the in situ, or contemporary stress field is vital for planning optimum orientations of deviated and horizontal wells, reservoir characterization and a better understanding of geodynamic processes and their effects on basin evolution.This study provides the first documented analysis of in situ stress and pore pressure fields in the sedimentary formations of the Cuu Long and Nam Con Son Basins, offshore Vietnam, based on data from petroleum exploration and production wells.In the Cuu Long Basin, the maximum horizontal stress is mainly oriented in NNW–SSE to N–S in the northern part and central high. In the Nam Con Son Basin, the maximum horizontal stress is mainly oriented in NE–SW in the northern part and to N–S in the central part of the basin.The magnitude of the vertical stress has a gradient of approximately 22.2 MPa/km at 3500 m depth. Minimum horizontal stress magnitude is approximately 61% of the vertical stress magnitude in normally pressured sequences.The effect of pore pressure change on horizontal stress magnitudes was estimated from pore pressure and fracture tests data in depleted zone caused by fluid production, and an average pore pressure–stress coupling ratio (ΔS_h/ΔP_p) obtained was 0.66. The minimum horizontal stress magnitude approaches the vertical stress magnitude in overpressured zones of the Nam Con Son Basin, suggesting that an isotropic or strike-slip faulting stress regime may exist in the deeper overpressured sequences.     

Petroleum composition in the Cuu Long Basin (Mekong Basin) offshore southern Vietnam

The Cuu Long Basin (Mekong Basin) is a rift basin off southern Vietnam, and the most important petroleum producing basin in the country. However, information on petroleum type and characteristics has hitherto been largely unavailable to the public. This paper presents petroleum geochemical data on nine oil samples from four different producing fields in the Cuu Long Basin: the Dragon (Rong), Black Lion (Sutu-Den), Sunrise (Rang ?ong) and White Tiger (Bach Ho) Fields. The oils are highly paraffinic with bimodal normal alkane distributions and show moderate pristane to phytane ratios and a conspicuous hyperbolic decrease in abundance with increasing carbon number of hopane homologues from C₃₀ to C₃₅. The TPP-index of Holba et al. (Holba, A.G., Dzou, L.I., Wood, G.D., Ellis, L., Adam, P., Schaeffer, P., Albrecht, P., Greene, T., Hughes, W.B., 2003. Application of tetracyclic polyprenoids as indicators of input from fresh–brackish water environments. Organic Geochemistry 34, 441–469) is equal to 1 in all samples which in combination with tricyclic triperpane T26/T25 ratios >1 and the n-alkane and hopane distributions mentioned above provide a strong indication of an origin from lacustrine source rocks. This is supported by the absence of marine C₃₀ desmethyl steranes (i.e. 24-n-propylcholestanes) and marine diatom-derived norcholestanes. Based on the overall biological marker distributions, the lakes probably belonged to the overfilled or balanced-fill types defined by Bohacs et al. (Bohacs, K.M., Carroll, A.R., Neal, J.E., Mankiewicz, P.J., 2000. Lake-basin type, source potential, and hydrocarbon character. An integrated sequence-stratigraphic–geochemical framework. AAPG Studies in Geology 46, 3–34). The oils were generated from source rocks at early- to mid-oil-window maturity, presumably Oligocene lacustrine shales that are present in the syn-rift succession. Oils from individual fields may, however, be distinguished by a combination of biological marker parameters, such as the oleanane index, the gammacerane index, the relative abundance of tricyclic terpanes, the proportions of diasteranes and 28-norspergulane, complemented by other parameters. The oils of the Cuu Long Basin show an overall similarity to the B-10 oil from the Song Hong Basin off northern Vietnam, but are markedly different from the seepage oils known from Dam Thi Nai on the coast of central Vietnam.     

Paleostress field reconstruction in the Oslo region

The geodynamic history of a region is archived in its geologic record which, in turn, may reflect deformation patterns that causally can be related to certain configurations of paleostresses. In the Oslo Region, the exposed geological record ranges from Precambrian high-grade metamorphic rocks through Cambro-Silurian sedimentary rocks to Permo-Carboniferous sedimentary and magmatic rocks, the latter being related to the development of the Oslo rift system. We investigate the kinematics of outcrop-scale faults to derive the diversity of paleostress states responsible for the observed strain. For this purpose, we combine different graphical and numerical approaches to separate heterogeneous fault-slip data sets and estimate the associated reduced stress tensors. A reduced stress tensor consists of the directions of the three principal stress axes with σ₁ ≥ σ₂ ≥ σ₃ and the ratio of principal stress differences, R = (σ₂ − σ₃)/(σ₁ − σ₃).     

Estimating horizontal stresses for mudrocks under one-dimensional compression

We examine horizontal stresses in mudrocks for the case of one-dimensional mechanical compression through laboratory measurements of the horizontal to vertical effective stress ratio, i.e. K₀. Previous empirical relationships which predict horizontal stresses as a function of depth or overburden stress are site specific and cannot easily be generalized. The K₀ values of various resedimented mudrocks from a diverse range of geologic origins have been investigated for effective stresses ranging from 0.1 to 100 MPa. The value of K₀ during normal compression (K_0NC) varies systematically as a function of both stress level and mudrock composition. K_0NC generally increases with stress level and values as high as 0.90 have been measured at effective stresses approaching 100 MPa. This finding implies that high horizontal stresses can develop at depth solely due to normal mechanical compression and that tectonic stresses, creep, or geologic unloading are not necessary. Correlations have been developed which allow K_0NC to be estimated from a mudrock's liquid limit, an easily and inexpensively measured property which reflects both the quantity and type of clay minerals present in a mudrock. High liquid limit, smectite-rich mudrocks display a more rapid increase in K_0NC with increasing stress compared to more silty, low liquid limit mudrocks. K₀ measurements made on laboratory resedimented specimens compare well with measurements made on intact core samples as well as in situ tests. Tests carried out to examine K₀ for overconsolidated materials show that K₀ during unloading can be well approximated by a power-law function originally proposed by Schmidt (1966). An original equation is proposed to describe K₀ during reloading. Our experimental results reveal sediment behavior at high stress levels that has not been previously observed and provide useful information on how horizontal stresses evolve within basins.     

The present-day state of tectonic stress in the Darling Basin,Australia: Implications for exploration and production

Knowledge of the full present-day stress tensor and pore pressure has significant applications in the exploration and production of conventional and unconventional hydrocarbon reservoirs. The Darling Basin of New South Wales, Australia, is an old sedimentary basin (Late Cambrian/Silurian to Early Carboniferous) in which there was limited information about the present-day stress field prior to this study. In this paper we evaluate the contemporary stress field of the Darling Basin using a dataset from recent exploration wells and perform a geomechanical risk assessment with respect to borehole stability, fracture/fault generation and reactivation. Our interpretations of borehole failures in borehole image logs reveal a prevailing east-west orientation of the maximum horizontal stress throughout the Darling Basin. The estimates of the magnitudes for the vertical, minimum and maximum horizontal stress in the studied wells indicate a transition between thrust and strike-slip faulting stress regime at 600–700 m depths, where the magnitude of vertical stress and minimum horizontal stress are close to each other. However, the presence of borehole breakouts and drilling-induced tensile fractures, that we observe in the image logs at greater depths (900–2100 m) indicate a transition into a strike-slip tectonic stress regime below a depth range of approximately 700–900 m. These findings are in agreement with overcoring stress measurements east and west of the investigated wells. Furthermore, there are several Neogene-to-Recent geological structures in the study area that indicate thrust faulting with an east-west oriented maximum horizontal stress orientation around this old sedimentary basin. The consistency between the orientation of maximum horizontal stress determined from wellbore data and neotectonic structures is significant, and implies that horizontal stress orientations derived from very recent geological features may be valuable inputs to geomechanical models in the absence of wellbore or other data. However, the recent surface geological structures suggest a thrust faulting stress regime that is in slight contrast to the transition between thrust and strike-slip stress regime (S_H > S_h ∼ S_v) indicated by petroleum data, and highlights a potential pitfall of using neotectonic structures in geomechanical models. In particular, careful attention and verification should be made when using neotectonic structures for input, calibration or confirmation of geomechanical models, especially in intraplate tectonic settings such as Australia.     

In-situ stress state in the Linxing region,eastern Ordos Basin,China: Implications for unconventional gas exploration and production

The Carboniferous and Permian sedimentary rocks (mainly the Shanxi and Taiyuan formations) in the Linxing region, eastern Ordos Basin, China, host a significant volume of unconventional gas resources (coalbed methane, shale gas and tight sandstone gas). Currently, the in-situ stress state is poorly understood but knowledge of this is extremely important for a range of applications, such as gas exploration and production, fracture stimulation and wellbore stability. The maximum horizontal stress (S_Hmax), minimum horizontal stress (S_hmin) and vertical stress (S_v) magnitudes, and the S_Hmax orientation in the Linxing region were systematically analyzed for the first time in the present study, which can provide a reference for subsequent numerical simulation and hydraulic fracturing design. Based on borehole breakouts and drilling-induced tensile fractures interpreted from borehole imaging logs, the S_Hmax orientation rotates from ∼NEE-SWW-trending in the southern part to ∼ NWW-SEE-trending in the northern part of the Linxing region. Both conventional logs and extended leak-off tests were used for stress magnitude determination. The results revealed three types of in-situ stress fields (S_v > S_Hmax > S_hmin, S_Hmax > S_v > S_hmin and S_Hmax > S_v ≈ S_hmin), and a dominant strike-slip stress regime (S_Hmax > S_v ≥ S_hmin) was found for the entire well section in the target Shanxi Formation and Taiyuan Formation in the Linxing region. In addition, differential stress increased with depth in the Linxing region, which indicates that wellbore instability might be a potentially significant problem when drilling wells that are vertical or ∼ N-S-trending.     

Adaptative responses of E. coli to marine environmental stresses: a modelling approach based on viability and dormancy concepts

In an attempt to synthesize conceptual and experimental information on the behaviour of enteric bacteria in seawater, a mathematical model has been developed. This model is based on changes in intracellular metabolisable components depending on physiological responses of the cell when subjected to nutrient starvation, salinity stress and solar radiation. Following a strategy which takes into account short- and mid-term physiological adaptation and reversal processes, the cells can develop in different states: culturable (B₁), viable and definitively nonculturable (B₂), and reversibly dormant (B₃) cells. Model parameters were deduced either from original microcosm experiments, from literature data or calibration procedures. Validation of the model was performed through specific data on the patterns of change in culturable cell abundances when Escherichia coli populations were subjected to separate or simultaneous experimental stresses. It was shown that model simulations fit these data satisfactorily. Simulation results clearly show the effectiveness of the antistress response and the importance of dormancy which confers higher resistance properties. The model also makes it possible to test the effect of the different stresses, and the role of pre-adaptation with regard to the dynamics of the different cellular states.     

Measurement and modeling of bed shear stress under solitary waves

Direct measurements of bed shear stresses (using a shear cell apparatus) generated by non-breaking solitary waves are presented. The measurements were carried out over a smooth bed in laminar and transitional flow regimes (~ 10⁴ < R_e < ~ 10⁵). Measurements were carried out where the wave height to water depth (h/d) ratio varied between 0.12 and 0.68; maximum near bed velocity varied between 0.16 m/s and 0.51 m/s and the maximum total shear stress (sum of skin shear stress and Froude–Krylov force) varied between 0.386 Pa and 2.06 Pa. The total stress is important in determining the stability of submarine sediment and in sheet flow regimes. Analytical modeling was carried out to predict total and skin shear stresses using convolution integration methods forced with the free stream velocity and incorporating a range of eddy viscosity models. Wave friction factors were estimated from skin shear stress at different instances over the wave (viz., time of maximum positive total shear stress, maximum skin shear stress and at the time of maximum velocity) using both the maximum velocity and the instantaneous velocity at that phase of the wave cycle. Similarly, force coefficients obtained from total stress were estimated at time of maximum positive and negative total stress and at maximum velocity. Maximum positive total shear stress was approximately 1.5 times larger than minimum negative total stress. Modeled and measured positive bed shear stresses are well correlated using the best convolution model, but the model underestimates the data by about 4%. Friction factors are dependent on the choice of normalizing using the maximum velocity, as is conventional, or the instantaneous velocity. These differ because the stress is not in phase with the velocity in general. Friction factors are consistent with previous data for monochromatic waves, and vary inversely with the square-root of the Reynolds number. The total shear stress leads the free stream fluid velocity by approximately 50°, whereas the skin friction shear stress leads by about 30°, which is similar to that reported by earlier researchers.     

Laboratory observation of boundary layer flow under spilling breakers in surf zone using particle image velocimetry

A boundary layer flow under spilling breakers in a laboratory surf zone with a smooth bottom is investigated using a high resolution particle image velocimetry (PIV) technique. By cross-correlating the images, oscillatory velocity profiles within a viscous boundary layer of O(1) mm in thickness are resolved over ten points. Using PIV measurements taken for an earlier study and the present study, flow properties in the wave bottom boundary layer (WBBL) over the laboratory surf zone are obtained, including the mean velocities, turbulence intensity, Reynolds stresses, and intermittency of coherent events. The data are then used to estimate the boundary layer thickness, phase variation, and bottom shear stress. It is found that while the time averaged mass transport inside the WBBL is onshore in the outer surf zone, it changes to offshore in the inner surf zone. The zero Eulerian mass transport occurs at h/h_b ≈ 0.92 in the outer surf zone. The maximum overshoot of the streamwise velocity and boundary layer thickness are not constant across the surf zone. The bottom shear stress is mainly contributed by the viscous stress through mean velocity gradient while the Reynolds stress is small and negligible. The turbulence level is higher in the inner surf zone than that in the outer surf zone, although only a slight increase of turbulent intensity is observed inside the WBBL from the outer surf zone to the inner surf zone. The variation of phase inside and outside the WBBL was examined through the spatial velocity distribution. It is found the phase lead is not constant and its value is significantly smaller than previous thought. By analyzing instantaneous velocity and vorticity fields, a remarkable number of intermittent turbulent eddies are observed to penetrate into the WBBL in the inner surf zone. The size of the observed large eddies is about 0.11 to 0.16 times the local water depth. Its energy spectra follow the − 5/3 slope in the inertial subrange and decay exponentially in the dissipation subrange.     

Distribution of bottom stress and tidal energy dissipation in a well-mixed estuary

Estimates of area-averaged tidal bottom stress are made for four channel segments of the Great Bay Estuary, N.H. Current and sealevel measurements are used to estimate acceleration and pressure gradient terms in the equation of motion, while the equation of motion itself is used to infer the remaining stress term. Dynamic terms, bottom stress values, friction coefficients and energy dissipation rates are estimated for each site. The analysis shows that while throughout the estuary the principal force balance is between the frictional stress and the pressure gradient forcing, RMS values of total bottom stress range from 2·67 to 10·38 Nm^?2 and friction coefficients vary from 0·015 to 0·054. Both stress and energy dissipation are largest in the seaward portion of the estuary with an order of magnitude decrease in dissipation at the most inland site.These distributions of stress and energy dissipation are consistent with cotidal charts of the principal semi-diurnal tidal constituent (M₂) which indicate that the estuary is composed of a highly dissipative more progressive tidal wave regime seaward and a less dissipative standing wave regime landward.     

Ring-Shaped Seismicity Structures in Southern California: Possible Preparation for Large Earthquake in the Los Angeles Basin

Some characteristics of seismicity in Southern California are studied. It is found that ring-shaped seismicity structures with threshold magnitudes M_th of 4.1, 4.1, and 3.8 formed prior to three large (M _w > 7.0) earthquakes in 1992, 1999, and 2010, respectively. The sizes of these structures are several times smaller than for intracontinental strike-slip events with similar magnitudes. Two ring-shaped structures are identified in areas east of the city of Los Angeles, where relatively large earthquakes have not occurred for at least 150 years. The magnitudes of large events which can occur in the areas of these structures are estimated on the basis of the previously obtained correlation dependence of ring sizes on magnitudes of the strike-slip earthquakes. Large events with magnitudes of M _w = 6.9 ± 0.2 and M _w = 8.6 ± 0.2 can occur in the area to the east of the city of Los Angeles and in the rupture zone of the 1857 great Fort Tejon earthquake, respectively. We believe that ring-structure formation, similarly to the other regions, is connected with deep-seated fluid migration.     

Effects of periodic saturation on stress–strain relationship of a sandstone mixture

While a crushed sandstone particle mixture, usually used as an engineering fill, is filled along bank of or in a large reservoir, it is subjected to periodic saturation induced by filling-drawdown cycles of reservoir water. The periodic saturation may induce post-construction settlement, and reduce stability of the filled body. In order to evaluate the effects of periodic saturation on stress–strain relationship of the mixture, several triaxial tests are carried out. According to experimental data, the periodic saturation may induce an increment of axial strain (Δε; <0.226%), a decrement of crest of deviator stress (Δ(σ₁?σ₃)_f; <0.192?MPa) and a decrement of angle of shearing resistance (Δφ; <3.70°). With increment of the number of periodic saturation cycles, the variations of the three parameters may be fitted by logarithmic curves. And with increase in the stress level for periodic saturation, the variations of the three parameters may be fitted by straight lines. Three fitting equations to predict the three parameters’ values, and an equation to calculate the settlement induced by the periodic saturation of a large-area filled foundation, are suggested.     

Effects of sea ice formation and diapycnal mixing on the Okhotsk Sea intermediate water clarified with oxygen isotopes

Processes relating to the formation of dense shelf water and intermediate water in the Okhotsk Sea were studied by examining oxygen isotope ratios (δ¹⁸O), salinity, and temperature. The salinity and δ¹⁸O of the cold dense shelf water on the northern continental shelf showed peculiar relationship. The relationship indicates that 3% of the mixed-layer water, having salinity of 32.6, froze and the remaining 97% became dense shelf water of salinities of more than 33.2 (σ_θ>26.7) during the sea ice formation. The salinity–δ¹⁸O relationship also shows that 20% of the Okhotsk Sea Intermediate Water at the σ_θ=26.8 level was derived from the dense shelf water. The remaining 80% came from the Western Subarctic Pacific water modified by diapycnal mixing of water affected by the surface cooling and freshening within the Okhotsk Sea. The mixing with dense shelf water contributes to only 26% of the temperature difference or 8% of the salinity difference between the original Pacific water and the Okhotsk Sea Intermediate Water at σ_θ=26.8. This result suggests that the cold and less saline properties of the Okhotsk Sea Intermediate Water are produced mainly by diapycnal mixing, rather than by mixing of the Pacific water with the dense shelf water.     

Deformation of a sandstone–mudstone particle mixture induced by periodic saturation

This study focuses on the evaluation on deformation induced by periodic saturation of a sandstone–mudstone particle mixture. Two types of triaxial tests, without and with periodic saturation, were performed. The strain–stress relationships from the two types of tests indicate that the periodic saturation may induce an increment of axial strain (Δε), and the Δε values are related to the ratio of confining to atmospheric pressure (σ₃/p_a), stress level for periodic saturation (L), and number of periodic saturation or cycles (N). The values of Δε are increasing along logarithmic curves with increment of N value from 1 to 20, and increase along straight lines with increasing L value from 0.18 to 0.82 or σ₃/p_a value from 1 to 4. Based on the analyses of experimental data, a logarithmic fitting equation, which is a function of N, L, and σ₃/p_a, is suggested to predict the Δε value. And based on the fitting equation and simple analyses on stress state, another equation, which may be used to estimate the settlement induced by periodic saturation of a large-area foundation filled using the sandstone–mudstone particle mixture, is also suggested.     

Current-induced scour beneath initially elevated subsea pipelines

When a subsea pipeline is laid on an uneven seabed, certain sections may have an initial elevation with respect to the far-field seabed, e_o, and thus potentially affecting the on-bottom stability of the pipeline. This paper focuses on quantifying the effects of the upstream dimensionless seabed shear stress, θ_∞, and Reynolds number, Re, on: (1) the maximum dimensionless seabed shear stress beneath the pipe, θ_max, to be compared to the critical shear stress in order to determine whether scour would occur and progress towards an equilibrium state; and, (2) the dimensionless equilibrium scour depth beneath the pipe, S_eq/D. Using a 2-D Reynolds averaged Navier-Stokes (RANS) approach along with the k-ω Shear Stress Transport (SST) turbulence model, a parametric study involving 243 computational fluid dynamics (CFD) simulations was conducted. The simulation results were used to develop a closed-form equation for the prediction of θ_max. Subsequently, experimental measurements of S_eq/D have been compiled from published literature, to develop a new closed-form equation for the prediction of S_eq/D with a high correlation to the experimental data. In summary, we present two closed-form equations for the prediction of θ_max and S_eq/D for pipelines with an initial e_o/D, which are applicable for both clear-water and live-bed conditions. The effects of θ_∞ and Re have been included, albeit Re having a small influence as compared to the other parameters.     

On the compression behavior of remolded cement-admixed soft clay

Although extensive research has been performed on the mechanical properties of cement-stabilized clays, quite a few attempts have been made on the compression behavior of remolded cement-admixed clays. The results from oedometer tests have been discussed to investigate the compressibility of remolded cement-admixed clays, taking into consideration cement amount and curing time. The findings show that the difference in shape and position of compression curves is attributed to cement amount and curing time. Most compression index (C_c) values of remolded cement-admixed clays are greater than those of untreated clay due to the presence of remolded yield stress σ′_yr that is closely related to initial water content and clay fabric. Based on the obtained test data, the relationships of C_c vs. e₀, C_c vs. w₀, C_c vs. e₁, C_c vs. e_yr, and σ′_yr vs. e_yr are preliminarily discussed and quantitatively established. Especially, an important divergence of void index I_v at effective stress σ′_v less than remolded yield stress σ′_yr can be observed at different cement amounts and curing durations. Being independent on cement amount, curing time, and initial state of soil, an excellent convergence occurs at stress σ′_v greater than yield stress σ′_yr. The normalized compression curves of I_v vs. σ′_v at σ′_v?>?σ′_y can be expressed by a unique line that agrees well with intrinsic compression line (ICL) and extended ICL.     

Interaction of dispersive water waves with weakly sheared currents of arbitrary profile

A set of depth-integrated equations describing combined wave–current flows is derived and validated. To account for the effect of turbulence induced by interactions between waves and currents with arbitrary horizontal vorticity, new additional stress terms are introduced. These stresses are functions of a parameter b that relates the relative importance of wave radiation stress and bottom friction stress to the wave–current interaction. To solve the equations, a fourth-order MUSCL-TVD scheme with an approximate Riemann solver is adopted. As a first-order check of the model, the Doppler shift effect and wave dispersion over linearly sheared currents are analytically shown to be retained appropriately in the equation set. The model results are then validated through comparisons with three experimental data sets. First, based on the experiments of Kemp and Simons (1982, 1983), a reasonable functional form of b is estimated. Second, simulations examining the propagation of a weakly dispersive wave over a depth-uniform or linearly sheared current are performed. Finally, the model is applied to a more complex configuration where bichromatic waves interact with spatially varying currents. Simulated results indicate that the model is capable of predicting nearshore interactions of waves with currents of arbitrary vertical structure. One of the unique properties of the developed model is its ability to assimilate an external current field from any source, be it from a circulation model or an observation, and predict the interaction of a nonlinear and dispersive wave field with that current.     

“Basin scale” versus “localized” pore pressure/stress coupling - Implications for trap integrity evaluation

In petroleum industry, the difference between pore pressure (P_p) and minimum horizontal stress S_h (termed the seal or retention capacity) is of major consideration because it is often assumed to represent how close a system is to hydraulic failure and thus the maximum hydrocarbon column height that can be maintained. While S_h and P_p are often considered to be independent parameters, several studies in the last decade have demonstrated that S_h and P_p are in fact coupled. However, the nature of this coupling relationship remains poorly understood. In this paper, we explore the influences of the spatial pore pressure distribution on S_h/P_p coupling and then on failure pressure predictions and trap integrity evaluation. With analytical models, we predict the fluid pressure sustainable within a reservoir before failure of its overpressured shale cover. We verify our analytical predictions with experiments involving analogue materials and fluids. We show that hydraulic fracturing and seal breach occur for fluid pressure greater than it would be expected from conventional retention capacity. This can be explained by the impact of the fluid overpressure field in the overburden and the pressure diffusion around the reservoir on the principal stresses. We calculate that supralithostatic pressure could locally be reached in overpressured covers. We also define the retention capacity of a cover (RC) surrounding a fluid source or reservoir as the difference between the failure pressure and the fluid overpressure prevailing in shale at the same depth. In response to a localized fluid pressure rise, we show that the retention capacity does not only depend on the pore fluid overpressure of the overburden but also on the tensile strength of the cover, its Poisson’s ratio, and the depth and width of the fluid source.     

Numerical study on spatially varying control parameters of a marine ecosystem dynamical model with adjoint method

Based on the simulation of a marine ecosystem dynamical model in the Bohai Sea,the Yellow Sea and the East China Sea,chlorophyll data are assimilated to study the spatially varying control parameters (CPs) by using the adjoint method.In this study,the CPs at some grid points are selected as the independent CPs,while the CPs at other grid points can be obtained through linear interpolation with the independent CPs.The independent CPs are uniformly selected from each 30 × 30 area,and we confirm that the optimal influence radius is 1.2° by a twin experiment.In the following experiments,when only the maximum growth rate of phytoplankton (V m) is estimated by two given types of spatially varying CPs,the mean relative errors of V m are 1.22% and 0.94% while the decrease rates of the mean error of chlorophyll in the surface are 94.6% and 95.8%,respectively.When the other four CPs are estimated respectively,the results are also satisfactory,which indicates that the adjoint method has a strong ability of optimizing the prescribed CP with spatial variations.However,when all these five most important CPs are estimated simultaneously,the collocation of the changing trend of each parameter influences the estimation results remarkably.Only when the collocation of the changing trend of each parameter is consistent with the ecological mechanisms which influence the growth of the phytoplankton in marine ecosystem,could the five most important CPs be estimated more accurately.     

Inertial and tidal shear variability above Reykjanes Ridge

Yearlong 75 kHz acoustic Doppler current profiler (ADCP) data were obtained well above Reykjanes Ridge (northern extension of the Mid-Atlantic Ridge (MAR)). The area is characterized by relatively large semidiurnal tidal (‘D₂’) currents that have (at lunar M₂) more than half a decade larger variance than inertial (f) currents. However, despite the relatively weak near-inertial kinetic energy, its vertical current shear shows larger magnitudes than at M₂ in an otherwise flat f−D₂ band limited between frequencies 0.74 and 1.35f, which equals the inertio-gravity wave bounds [σ_min, σ_max](N=f). N represents the buoyancy frequency. The shear in this band dominates all shear computed at 20 m effective vertical scale. As the kinetic energy spectrum peaks at M₂, but not (significantly) at S₂ and N₂, a difference in tidal (and inertial) scales and hence sources is observed. M₂-tides contribute mostly to large-scale coherent motions. The dominant incoherent f−D₂ shear is highly variable in time (∼2-day periodicity). Furthermore, inertial and tidal shear are more or less completely separated in space and time, each occurring in different layers in the vertical. The thin shear layers reflect the rapidly varying short vertical scale N_s profile, to within the ∼20 m limitation of ADCP data, rather than the large-scale smooth N_L. In each of large-N_s layers Ri≈1, probably. The yearlong smoothed shear magnitude follows N_L, but only as stable Ri≈5. The shear polarization is more circular than rectilinear, albeit varying with time, and highly symmetric around f. During transitions, e.g., between stratified and homogeneous layers and between waves from varying sources, near-circular motions can generate near-rectilinear shear in the direction of wave propagation (in the direction of the minor axis of the current ellipse). This contrasts with the possibility of near-rectilinear barotropic oscillatory motions generating near-circular shear under viscosity in shallow seas.