Estimation of Pullout Capacity of Vertical Plate Anchors in Cohesionless Soil Using MARS

Vertical plate anchors provide an economical solution to safely resist the large horizontal forces experienced by the foundation of different structures such as bulkheads, sheet piles, retaining walls and so forth. This paper develops a multivariate adaptive regression spline (MARS) model-based approach for the determination of horizontal pullout capacity (P _u ) of vertical plate anchors buried in cohesionless soil by utilizing experimental results reported by different researchers. Based on the collection of forty different pullout experimental test results reported in the literature for anchors buried in loose to dense cohesionless soil with an embedment ratio ranges from 1 to 5, a predictive approach for P _u of vertical plate anchors has been developed in terms of non-dimensional pullout coefficient (M _γq ). The capability of the proposed MARS model for estimating the values of M _γq is examined by comparing the results obtained in the present study with those methods available in the literature. Using different statistical error measure criteria, this study indicates that the present approach is efficient in estimating the horizontal pullout capacity of vertical plate anchors as compared to other methods. The sensitivity analysis indicates that the embedment ratio (H/h, where H = embedment depth of anchor, and h = height of anchor) and internal friction angle (?) of soil mass are the two most important parameters for the evaluation of non-dimensional pullout coefficient (M _γq ) using the proposed MARS model.     

Seismic Pullout Capacity of Inclined Anchor Plates in Sand

By using the lower bound finite elements limit analysis, the pullout capacity of an inclined strip anchor plate embedded in a cohesionless soil medium has been computed with an inclusion of pseudo-static horizontal earthquake body forces. The variation of the pullout capacity factor (F _γ ) with changes in horizontal earthquake acceleration co-efficient (α _h ) has been computed by varying the inclination angle (β) of the anchor plate between 0° and 90°. The results clearly reveal that the pullout capacity factor (F _γ ) decreases significantly with an increase in the value of α _h . The reduction in the pullout resistance due to seismic forces (1) becomes much more extensive for a vertical anchor plate as compared to the horizontal anchor, (2) decreases generally with increases in the soil friction angle (?) and (3) increases with an increase in friction angle between soil and anchor plate (δ). The developments of the failure zone around the anchor plate were also examined by varying α _h and β. The results obtained from the analysis compare well with the solutions reported in literature.     

Effect of 2-D Random Field Discretization on Failure Probability and Failure Mechanism in Probabilistic Slope Stability

This paper investigates, using the random field theory and Monte Carlo simulation, the effects of random field discretization on failure probability, p _f, and failure mechanism of cohesive soil slope stability. The spatial sizes of the discretized elements in random field Δx, Δy in horizontal and vertical directions, respectively, are assigned a series of combinational values in order to model the discretization accuracy. The p _f of deterministic critical slip surface (DCSS) and that of the slope system both are analyzed. The numerical simulation results have demonstrated that both the ratios of Δy/λ _y (λ _y  = scale of fluctuation in vertical direction) and Δx/λ _x (λ _x  = scale of fluctuation in horizontal direction) contribute in a similar manner to the accuracy of p _f of DCSS. The effect of random field discretization on the p _f can be negligible if both the ratios of Δx/λ _x and Δy/λ _y are no greater than 0.1. The normalized discrepancy tends to increase at a linear rate with Δy/λ _y when Δx/λ _x is larger than 0.1, and vice versa for p _f of DCSS. The random field discretization tends to have more considerable influence on the p _f of DCSS than on that of the slope system. The variation of p _f versus λ _x and λ _y may exhibit opposite trends for the cases where the limit state functions of slope failure are defined on DCSS and on the slope system as well. Finally, the p _f of slope system converges in a more rapid manner to that of DCSS than the failure mechanism does to DCSS as the spatial variability of soil property grows from significant to negligible.     

Enhancement of random finite element method in reliability analysis and risk assessment of soil slopes using Subset Simulation

Random finite element method (RFEM) provides a rigorous tool to incorporate spatial variability of soil properties into reliability analysis and risk assessment of slope stability. However, it suffers from a common criticism of requiring extensive computational efforts and a lack of efficiency, particularly at small probability levels (e.g., slope failure probability P _f ?<?0.001). To address this problem, this study integrates RFEM with an advanced Monte Carlo Simulation (MCS) method called “Subset Simulation (SS)” to develop an efficient RFEM (i.e., SS-based RFEM) for reliability analysis and risk assessment of soil slopes. The proposed SS-based RFEM expresses the overall risk of slope failure as a weighed aggregation of slope failure risk at different probability levels and quantifies the relative contributions of slope failure risk at different probability levels to the overall risk of slope failure. Equations are derived for integrating SS with RFEM to evaluate the probability (P _f ) and risk (R) of slope failure. These equations are illustrated using a soil slope example. It is shown that the P _f and R are evaluated properly using the proposed approach. Compared with the original RFEM with direct MCS, the SS-based RFEM improves, significantly, the computational efficiency of evaluating P _f and R. This enhances the applications of RFEM in the reliability analysis and risk assessment of slope stability. With the aid of improved computational efficiency, a sensitivity study is also performed to explore effects of vertical spatial variability of soil properties on R. It is found that the vertical spatial variability affects the slope failure risk significantly.     

Simulation of suspended sediment based on gamma test,heuristic, and regression-based techniques

In the present study, four different heuristic techniques viz. multi-layer perceptron (MLP), radial basis function (RBF), self-organizing maps (SOM), and co-active neuro-fuzzy inference system (CANFIS) with hyperbolic tangent and sigmoid transfer functions and two regression-based techniques, i.e., multiple linear regression (MLR) and sediment-rating curve (SRC), were used for suspended sediment modeling. Gamma test (GT), correlation function (CF), M test, and trail–error procedure were applied for estimation of appropriate input variables as well as training data length. The results of the GT and CF suggested the five input variables (Q_t, Q_t?1, Q_t?2, S_t?1, and S_t?2, where Q_t?1 and S_t?1 indicate the discharge and sediment values of one previous day) as the best combination. The optimal training data length (75% of total data) was estimated by M test and trail–error procedure for development of the applied models. The MLP with sigmoid transfer function (M-2) performed better than the all other models. The results of sensitivity analysis indicated that the present-day discharge (Q_t), 1-day lag discharge (Q_t?1) and 1-day lag suspended sediment (S_t?1) are the most influenced parameters in modeling current day suspended sediment (S_t).     

Numerical analysis under seismic loads of soils improvement with floating stone columns

Geotechnical Engineering has developed many methods for soil improvement so far. One of these methods is the stone column method. The structure of a stone column generally refers to partial change of suitable subsurface ground through a vertical column, poor stone layers which are completely pressed. In general terms, to improve bearing capacity of problematic soft and loose soil is implemented for the resolution of many problems such as consolidation and grounding problems, to ensure filling and splitting slope stability and liquefaction that results from a dynamic load such as earthquake. In this study, stone columns method is preferred as an improvement method, and especially load transfer mechanisms and bearing capacity of floating stone column are focused. The soil model, 32 m in width and 8 m in depth, used in this study is made through Plaxis 2D finite element program. The clay having 5° internal friction angle with different cohesion coefficients (c 10, c 15, c 20 kN/m²) are used in models. In addition, stone columns used for soil improvement are modeled at different internal friction angles (? 35°, ? 40°, ? 45°) and in different s/D ranges (s/D 2, s/D 3), stone column depths (B, 2B, 3B) and diameters (D 600 mm, D 800 mm, D 1000 mm). In the study, maximum acceleration (a _max = 1.785 m/s²) was used in order to determine the seismic coefficient used. In these soil models, as maximum acceleration, maximum east–west directional acceleration value of Van Muradiye earthquake that took place in October 23, 2011 was used. As a result, it was determined that the stone column increased the bearing capacity of the soil. In addition, it is observed that the bearing capacity of soft clay soil which has been improved through stone column with both static and earthquake load effect increases as a result of increase in the diameter and depth of the stone column and decreases as a result of the increase in the ranges of stone column. In the conducted study, the bearing capacity of the soil models, which were improved with stone column without earthquake force effect, was calculated as 1.01–3.5 times more on the average, compared to the bearing capacity of the soil models without stone column. On the other hand, the bearing capacity of the soil models with stone columns, which are under the effect of earthquake force, was calculated as 1.02–3.7 times more compared to the bearing capacity of the soil models without stone column.     

Failure modes and bearing capacity of strip footings on soft ground reinforced by floating stone columns

This study evaluates the failure modes and the bearing capacity of soft ground reinforced by a group of floating stone columns. A finite difference method was adopted to analyze the performance of reinforced ground under strip footings subjected to a vertical load. The investigation was carried out by varying the aspect ratio of the reinforced zone, the area replacement ratio, and the surface surcharge. General shear failure of the reinforced ground was investigated numerically without the surcharge. The results show the existence of an effective length of the columns for the bearing capacity factors N _c and N _γ. When certain surcharge was applied, the failure mode of the reinforced ground changed from the general shear failure to the block failure. The aspect ratio of the reinforced zone and the area replacement ratio also contributed to this failure mode transition. A counterintuitive trend of the bearing capacity factor N _q can be justified with a shift in the critical failure mode. An upper-bound limit method based on the general shear failure mode was presented, and the results agree well with those of the previous studies of reinforced ground. Equivalent properties based on the area-weighted average of the stone columns and clay parameters were used to convert the individual column model to an equivalent area model. The numerical model produced reasonable equivalent properties. Finally, a theoretical method based on the comparison of the analytical equations for different failure modes was developed for engineering design. Good agreement was found between the theoretical and numerical results for the critical failure mode and its corresponding bearing capacity factors.     

Radial pulsations of helium stars and a model for BX CIR

The results of hydrodynamical calculations of radially pulsating helium stars with masses 0.5M_⊙≤M≤0.9M_⊙, bolometric luminosities 600L_⊙≤5×10³L_⊙, and effective temperatures 1.5×10⁴ K≤T_eff≤3.5×10⁴ K are presented. The pulsation instability of these stars is due to the effects of ionization of iron-group elements in layers with temperatures T～2×10⁵ K. The calculations were carried out using opacities for the relative mass abundances of hydrogen and heavy elements X=0 and Z=0.01, 0.015, and 0.02. Approximate formulas for the pulsation constant Q over the entire range of pulsation instability of the hot helium stars in terms of the mass M, radius R, effective temperature T_eff, and heavy-element abundance Z are derived. The instability of BX Cir to radial pulsations with the observed period Π=0.1066 d occurs only for a mass M≥0.55M_⊙, effective temperature T_eff≥23000 K, and heavy-element abundance Z≥0.015. The allowed mass of BX Cir is in the range 0.55M_⊙≤M≤0.8M_⊙, which corresponds to luminosities 800L_⊙≤M≤1400L_⊙ and mean radii 1.7R_⊙?R?2.1R_⊙.     

Seismic hazard scenario and attenuation model of the Garhwal Himalaya using near-field synthesis from weak motion seismometry

In this paper, we present a seismic hazard scenario for the Garhwal region of the north-western Himalayan range, in terms of the horizontal Peak Ground Acceleration. The scenario earthquake of moment magnitude M _w 8.5 has a 10% exceedance probability over the next 50 years. These estimates, the first for the region, were calculated through a stepwise process based on:

• An estimation of the Maximum Credible Earthquake from the seismicity of the region and Global Seismic Hazard Assessment Program considerations, and
• four seismotectonic parameters abstracted from near field weak-motion data recorded at five stations installed in Chamoli District of the Garhwal region in the aftermath of the 1999 Chamoli earthquake. The latter include
• The frequency dependent power law for the shear wave quality factor, Q _S
• the site amplification at each station using horizontal-to-vertical-spectral ratio and generalized inversion technique
• source parameters of various events recorded by the array and application of the resulting relations between the scalar seismic moment M ₀ (dyne-cm) and moment magnitude M _w and the corner frequency, ? _c (Hz) and moment magnitude M _w to simulate spectral acceleration due to higher magnitude events corresponding to the estimated Maximum Credible Earthquake, and
• regional and site specific local spectral attenuation relations at different geometrically central frequencies in the low, moderate and high frequency bands.

     

Shear wave velocity as function of cone penetration resistance and grain size for Holocene-age uncemented soils: a new perspective

For feasibility studies and preliminary design estimates, field measurements of shear wave velocity, V _s, may not be economically adequate and empirical correlations between V _s and more available penetration measurements such as cone penetration test, CPT, data turn out to be potentially valuable at least for initial evaluation of the small-strain stiffness of soils. These types of correlations between geophysical (Vs) and geotechnical (N-SPT, q _c-CPT) measurements are also of utmost importance where a great precision in the calculation of the deposit response is required such as in liquefaction evaluation or earthquake ground response analyses. In this study, the stress-normalized shear wave velocity V _s1 (in m/s) is defined as statistical functions of the normalized dimensionless resistance, Q _tn-CPT, and the mean effective diameter, D ₅₀ (in mm), using a data set of different uncemented soils of Holocene age accumulated at various sites in North America, Europe, and Asia. The V _s1–Q _tn data exhibit different trends with respect to grain sizes. For soils with mean grain size (D ₅₀) < 0.2 mm, the V _s1/Q _tn ^0.25 ratio undergoes a significant reduction with the increase in D ₅₀ of the soil. This trend is completely reversed with further increase in D ₅₀ (D ₅₀ > 0.2 mm). These results corroborate earlier results that stressed the use of different CPT-based correlations with different soil types, and those emphasized the need to impose particle-size limits on the validity of the majority of available correlations.     

Towards a better understanding of the evolution of the flood risk in Mediterranean urban areas: the case of Barcelona

Dam failure constitutes a grave threat to human life. However, there is still a lack of systematic and comprehensive research on the loss of life (L) caused by dam break in China. From the perspective of protecting human life, a new calculation method for L occurred in dam break floods is put forward. Fourteen dam failure cases in China are selected as the basic data by three-dimensional stratified sampling, balancing spatial, vertical elevation and temporal representations, as well as considering various conditions of the dam collapse. The method includes three progressive steps: Firstly, some impact factors of loss of life (IFL) are selected by literature survey, i.e., severity of dam break flood (S _F), population at risk (P _R), understanding of dam break (U _B), warning time (T _W) and evacuation condition (E _C). And the other IFL of weather during dam break (W _B), dam break mode (M _B), water storage (S _W), building vulnerability (V _B), dam break time (T _B) and average distance from affected area to dam (D _D) are also taken into account to get a more comprehensive consideration. According to disaster system and disaster risk, these eleven IFL are divided into four categories. Through the improved entropy method, eight key IFL are further selected out of the eleven. Secondly, four L modules are built based on four categories, which are L-causing factor module (M ₁), L-prone environment module (M ₂), affected body module (M ₃) and rescue condition module (M ₄). Eventually, by using two methods of multivariate nonlinear regression and leave-one-out cross-validation in combination with coupled four modules, the calculation method for L is established. Compared with the results of Graham method and D&M method, the result of the proposed one is much closer to the actual value and performs better in fitting effect and regional applicability. In the application, L calculation and consequence assessment are carried out in the example of Hengjiang reservoir that has already broken down. At the same time, L calculation and risk prediction are used in the analysis of Yunshan reservoir, which is under planning. The proposed method can not only be applied to estimate L and its rate (f _L ) under various types of dam break conditions in China, but also provide a reliable consequence assessment and prediction approach to reduce the risk of L.     

Photosynthetic capability and Fe,Mn, Cu,and Zn contents in two Moraceae species under different phosphorus levels

The strong adaptability of Broussonetia papyrifera (L.) Vent. to low phosphorus (P) conditions can be attributed to the large amount of root-exuded organic acids and the high efficiency of P extraction. However, microelement contents are influenced by low-P stress, and their effects on the photosynthetic capability of B. papyrifera remain unknown. In this study, we investigated the effects of low-P treatment on net photosynthetic rate (P _N); chlorophyll a fluorescence (ChlF) characteristics; and Fe, Mn, Cu, and Zn contents of B. papyrifera and Morus alba L. seedlings. Results show that B. papyrifera exhibited better photosynthetic capability under moderate P deficiency (0.125, 0.063, and 0.031 mmol/L P treatments), whereas the photosynthetic capability of M. alba decreased under moderate and severe P deficiency (0.016 and 0 mmol/L P treatments). Under moderate P deficiency, the decrease in Cu and Zn contents in B. papyrifera was lower than that in M. alba. Under severe P deficiency, a considerable decrease of photosynthetic capability in B. papyrifera and M. alba was associated with low Cu and Zn contents. The P _N of the two Moraceae species exhibited a better correlation with Cu and Zn contents than with Fe or Mn content. P deficiency could not only decrease cyclic photophorylation and photosynthetic efficiency, but could also affect the stability of thylakoid membrane structure and electron transport efficiency by influencing the contents of Cu or Zn, thereby affecting photosynthesis.     

Masses of the visual components and black holes in X-ray novae: Effects of proximity of the components

It is shown that the approximation of the complex, tidally distorted shape of a star as a circular disc with local line profiles and a linear limb-darkening law, which is usually applied when deriving equatorial stellar rotation velocities from line profiles, leads to overestimation of the equatorial velocity V _rot sin i and underestimation of the component mass ratio q = M _x /M _v . A formula enabling correction of the effect of these simplifying assumptions on the shape of a star is used to re-determine the mass ratios q and the masses of the black holes M _x and visual components M _v in low-mass X-ray binary systems containing black holes. Taking into account the tidal–rotational distortion of the stellar shape can significantly increase the mass ratios q = M _x /M _v , reducing M _v , while M _x changes only slightly. The resulting distribution of M _v attains its maximum near M _v ? 0.35M _⊙, in disagreement with the results of population synthesis computations realizing standard models for Galactic X-ray novae with black holes. Possible ways to overcome this inconsistency are discussed. The derived distribution of M _x also differs strongly from the mass distribution for massive stars in the Galaxy.     

Dynamic response characteristics of a rock slope with discontinuous joints under the combined action of earthquakes and rapid water drawdown

In order to study the dynamic response characteristics of a rock slope with discontinuities under the combined action of earthquakes and rapid water drawdown, a large-scale shaking table test was performed on a rock slope with discontinuous joints. Wenchuan earthquake (WE) seismic records were performed to investigate the horizontal and vertical acceleration response and displacement response. In particular, three-dimensional optical measurement techniques was used to obtain the slope surface displacements. A comparison was made on the seismic response according to the analysis of PGD (peak ground displacement) and M _PGA (acceleration amplification coefficient) of the modeled slope. The results show that the experimental slope mainly underwent settlement and horizontal deformation when the WE records were applied in the z and x directions, respectively. The slope was first shaken by the P wave, which caused the differential settlement to occur at the surface slope; then, the slope was shaken more severely by the S wave, which led to a greater horizontal deformation. Moreover, analysis of the ΔPGD (increment of PGD) and ΔM _PGA (increment of M _PGA) under rapid drawdown suggests that the rapid water drawdown mainly impacts the deformation of the surface slope, particularly between the high and low water levels. The water infiltration through the cracks softened the material of the surface slope, and the rapid drawdown also enhanced the slope deformation. In addition, the damage evolution process of the slope can be identified, mainly including three stages: an elastic stage (<?0.168 g), a plastic stage (0.168–0.336 g), and a failure stage (>?0.336 g).     

Development of Analytical Models to Estimate Pile Setup in Cohesive Soils Based on FE Numerical Analyses

This paper presents the numerical simulation of pile installation and the subsequent increase in the pile capacity over time (or setup) after installation that was performed using the finite element software Abaqus. In the first part, pile installation and the following load tests were simulated numerically using the volumetric cavity expansion concept. The anisotropic modified Cam-Clay and Dracker–Prager models were adopted in the FE model to describe the behavior of the clayey and sandy soils, respectively. The proposed FE model proposed was successfully validated through simulating two full-scale instrumented driven pile case studies. In the second part, over 100 different actual properties of individual soil layers distracted from literature were used in the finite element analysis to conduct parametric study and to evaluate the effect of different soil properties on the pile setup behavior. The setup factor A was targeted here to describe the pile setup as a function of time after the end of driving. The selected soil properties in this study to evaluate the setup factor A include: soil plasticity index (PI), undrained shear strength (S _u ), vertical coefficient of consolidation (C _v ), sensitivity ratio (S _r ), and over-consolidation ratio (OCR). The predicted setup factor showed direct proportion with the PI and S _r and inverse relation with S _u , C _v and OCR. These soil properties were selected as independent variables, and nonlinear multivariable regression analysis was performed using Gauss–Newton algorithm to develop appropriate regression models for A. Best models were selected among all based on level of errors of prediction, which were validated with additional nineteen different site information available in the literature. The results indicated that the developed model is able to predict the setup behavior for individual cohesive soil layers, especially for values of setup factor greater than 0.10, which is the most expectable case in nature.     

Bearing capacity of strip footings on <Emphasis Type="Italic">c</Emphasis>–<Emphasis Type="Italic">φ</Emphasis> soils with square voids

The presence of underground voids has an adverse influence on the performance of shallow foundations. In this study, the bearing capacity and failure mechanism of footings placed on cohesive-frictional soils with voids are evaluated using discontinuity layout optimization. By introducing a reduction coefficient, a set of design charts that can be directly applied to the classical bearing capacity formulation is presented. The results indicate that the undrained bearing capacity with voids is sensitive to soil weight and cohesion, as both the bearing capacity and stability issues exist in the problem. The failure mechanism is directly related to a variety of soil properties, the locations of single voids, and the horizontal distance between two voids. The presence of voids has a more dominant effect on c–φ soils compared to that on undrained soil. An interpretation of the critical and adverse locations for single-void and dual-void cases with various soil strengths is presented.     

Coda <Emphasis Type="Italic">Q</Emphasis> in Qeshm Island,south of Iran,using aftershocks of the Qeshm earthquake of November 27, 2005

Coda wave attenuation is estimated for Qeshm Island which is located in the southeastern part of Zagros. For this purpose, the aftershocks of Qeshm earthquake in November 27, 2005, recorded within an epicentral distance less than 100 km, have been used. More than 829 earthquakes were recorded by a local temporary network consisting of 16 short period stations installed after a week after the main shock for ～10 weeks. The coda quality factor, Q _c, was estimated using the single-backscattering model in frequency bands of 0.5–24 Hz. In this research, lateral and vertical variations of coda Q in Qeshm Island are explored. In Qeshm Island, absence of significant lateral variation of coda Q is observed. To investigate the attenuation variation with depth, the coda Q value was calculated for coda time windows with different lengths (5, 10, 15, 20, 25, and 30 s). It is observed that coda Q increases with depth. However, in our study area, the rate of increase of coda Q with depth is not uniform. Beneath Qeshm Island, the rate of increase of coda Q is greater at depths less than ～40 km compared with those of larger depths. This is indicating the existence of a low attenuation anomalous structure under the ～40-km depth which may be correlated with the Moho depth in this region. The average frequency relation for this region is Q _c = 36 ± 1.2f ^{0.94 ± 0.039} at a 5 s-lapse time window length and Q _c = 110 ± 1.8f ^{0.88 ± 0.09} at a 30-s lapse time window length.     

Vertical stability and the Brunt-Väisäla frequency of deep natural waters by the example of Lake Baikal,Lake Tanganyika,and the World Ocean

Theoretical analysis, calculations, and comparison with the results of observations in Lake Baikal, Lake Tanganyika, and the World Ocean are performed for the vertical stability E and the Brunt-Väisäla frequency N in the form of N ² with regard to all components (at the constant temperature T and the salinity S, the common adiabatic form at T, S Const). The adiabatic stability E _ad and the Väisäla frequency N in the form of N _ad ² are always positive; at a change from the inverse to the direct temperature stratification, they have deep minimums reaching 10^?16 m^?1 and 10^?15 s^?2 and less; the minimums have the form of a special point, a reversal point of the first kind called a “cusp.” The reality of these reversal points is confirmed by the analysis of the investigation procedure, comparison with the results of previous theoretical (Sherstyankin, et al., 2007), and experimental (observations in Baikal, Shimaraev et al., 1994) works. The features of vertical profiles of E _ad , E andN _ad ² , N ², as well as the layers where the Brunt-Väisäla frequency is less than the inertial frequency, are studied. The analysis with regard to all components of the stability E _ad and the Brunt-Väisäla frequency N makes a great contribution to understanding of mixing processes in theoretical and experimental investigations; it is valid in all reservoirs of the Earth with inverse and direct temperature stratification, including Lake Baikal, Lake Tanganyika, and the World Ocean.     

Quantifying the Effect of Salinity Stratification on Phytoplankton Density Patterns in Estuaries

To quantify the effect of salinity stratification on phytoplankton density (denoted as P) patterns, experiments were conducted with an idealised model that couples physical and biological processes. Results show that the idealised model is capable of capturing the main features of observed P patterns in the Columbia River estuary during the spring season: during weak stratification, P is almost vertically uniform with values decreasing towards the estuary mouth, whereas during strong stratification, high values of P extend further seawards but are confined to the upper layer. Sensitivity studies reveal that the strong vertical gradients of P can only occur if the intensity of turbulence (measured by depth-averaged values of vertical eddy viscosity and eddy diffusivity) is weak. The advection of P by subtidal currents is important in obtaining a smaller along-estuary gradient of P during weak stratification and in obtaining a smaller horizontal gradient and a larger vertical gradient of P during strong stratification. Accounting for stratification controlled vertical distribution of vertical eddy viscosity and eddy diffusivity is necessary for obtaining realistic P patterns if stratification is strong, but not if stratification is weak. A higher osmotic stress, which leads to faster loss of phytoplankton in salt water, results in a larger along-estuary gradient of P if stratification is weak and in a larger vertical gradient of P if stratification is strong.     

Supermassive black holes: Relation to dark halos

Estimates of the masses of supermassive black holes (M _bh ) in the nuclei of disk galaxies with known rotation curves are compared with estimates of the rotational velocities V _m and the “indicative” masses of the galaxies M _i . Although there is a correlation between M _bh and V _m or M _i , it is appreciably weaker than the correlation with the central velocity dispersion. The values of M _bh for early-type galaxies (S0-Sab), which have more massive bulges, are, on average, higher than the values for late-type galaxies with the same rotational velocities. We conclude that the black-hole masses are determined primarily by the properties of the bulge and not the rotational velocity or the mass of the galaxy.