Numerical modeling of shallow water table behavior with Lisse effect

Air entrapment is an important consideration in environments with shallow water tables and sandy soil, like the condition of highly conductive sandy soils and flat topography in Florida, USA. It causes water table rises in soils, which are significantly faster and higher than those in soils without air entrapment. Two numerical models, Integrated Hydrologic Model (IHM) and HYDRUS-1D (a single-phase, one-dimensional Richards′ equation model) were tested at an area of west central Florida to help further understanding the shallow water table behavior during a long term air entrapment. This investigation employed field data with two modeling approaches to quantify the variation of air pressurization values. It was found that the air pressurization effect was responsible at time up to 40 cm of water table rise being recorded by the observation well for these two models. The values of air pressurization calculated from IHM and HYDRUS-1D match the previously published values. Results also indicated that the two numerical models did not consider air entrapment effect (as the predictive parameters remain uncertain) and thus results of depth to water table from these models did not compare to the observations for these selected periods. Incorporating air entrapment in prediction models is critical to reproduce shallow water table observations.     

Using otolith weight to predict the age of <Emphasis Type="Italic">Pennahia macrocephalus</Emphasis> in the mouth of the Beibu Gulf

The relation between otolith weight (OW) and the age of marine fish is studied. A total of 222 individuals of bighead white croaker, Pennahia macrocephalus were sampled seasonally in the mouth of the Beibu Gulf, the South China Sea, in 2007. Since there are no significant differences in sagittal OW between otolith in pairs (P≥0.05), the undamaged left sagittal otolith is used for age determination. The highest correlations among standard length, OW and fish ages are confirmed by linear, exponential and multinomial regression. Results show that sagittal OW overlaps only occasionally among age groups, and to individuals with similar standard length, the older and slower-growing fish has a heavier otolith because of the continued otolith material deposition. There are differences in sagittal OW among different age groups and significant positive linear relationship with age (P<0.05). The age readings can be verified by plotting the sagittal OW versus the standard length for age groups, and the individuals with similar standard length but in different ages can be separated by sagittal OW frequency analysis. Mostly, the predicted ages using the regression between sagittal OW and ages are closed to the observed ages by counting annulus on scale. It indicates that the sagittal OW analysis is a useful technique for validating the accuracy of age determination by annuli counts, especially for individuals of similar size. Furthermore, the technique is applied for Pennahia macrocephalus with discussion in this paper. Supported by the National Natural Science Foundation of China (No. 30771653) and Bureau of Fisheries, Ministry of Agriculture     

Factors acquisition and content estimation of farmland soil organic carbon based upon Internet of Things

Aiming at the shortage of sufficient continuous parameters for using models to estimate farmland soil organic carbon (SOC) content, an acquisition method of factors influencing farmland SOC and an estimation method of farmland SOC content with Internet of Things (IOT) are proposed in this paper. The IOT sensing device and transmission network were established in a wheat demonstration base in Yanzhou Distict of Jining City, Shandong Province, China to acquire data in real time. Using real-time data and statistics data, the dynamic changes of SOC content between October 2012 and June 2015 was simulated in the experimental area with SOC dynamic simulation model. In order to verify the estimation results, potassium dichromate external heating method was applied for measuring the SOC content. The results show that: 1) The estimated value matches the measured value in the lab very well. So the method is feasible in this paper. 2) There is a clear dynamic variation in the SOC content at 0.2 m soil depth in different growing periods of wheat. The content reached the highest level during the sowing period, and is lowest in the flowering period. 3) The SOC content at 0.2 m soil depth varies in accordance with the amount of returned straw. The larger the amount of returned straw is, the higher the SOC content.     

Capacities and failure modes of suction bucket foundation with internal bulkheads

Suction bucket foundations can be divided into four compartments by cruciform internal bulkheads, thereby yielding better capacity in certain conditions than those without internal bulkheads. As yet, no systematic study has been conducted regarding the effects of cruciform internal bulkheads on the capacities of suction bucket foundations. In this study, we established a large number of finite element models of suction bucket foundations with and without cruciform internal bulkheads and of solid embedded circular foundations. We found the uniaxial capacities and failure modes of suction bucket foundations with various depth ratios to remain basically unaffected by internal bulkheads in uniform clays. However, in inhomogeneous clay with high strength heterogeneity, we observed the uniaxial moment and horizontal capacities and corresponding failure modes of suction bucket foundations with a low depth ratio to be obviously affected by internal bulkheads. In this case, the uniaxial moment capacities, in particular, as well as the horizontal capacities of suction bucket foundations with cruciform internal bulkheads become obviously greater than those without internal bulkheads. Under combined loading, we found the failure envelopes of suction bucket foundations with and without cruciform internal bulkheads and of solid circular foundation to also be basically consistent in uniform clays. However, in inhomogeneous clay with high strength heterogeneity, cruciform internal bulkheads can obviously change the shapes of the failure envelopes of bucket foundations with a small depth ratio. We conclude that when the acting vertical load or foundation depth is relatively small, suction bucket foundations with cruciform internal bulkheads can be subjected to larger moment and horizontal loads in soft clays with high strength heterogeneity.     

Origin of Radio Enhancements in Type II Bursts in the Outer Corona

We study interplanetary (IP) solar radio type II bursts from 2011?–?2014 in order to determine the cause of the intense enhancements in their radio emission. Type II bursts are known to be due to propagating shocks that are often associated with fast halo-type coronal mass ejections (CMEs). We analysed the radio spectral data and the white-light coronagraph data from 16 selected events to obtain directions and heights for the propagating CMEs and the type II bursts. CMEs preceding the selected events were included in the analysis to verify whether CME interaction was possible. As a result, we were able to classify the events into five different groups. 1) Events where the heights of the CMEs and type II bursts are consistent, indicating that the shock is located at the leading front of the CME. The radio enhancements are superposed on the type II lanes, and they are probably formed when the shock meets remnant material from earlier CMEs, but the shock continues to propagate at the same speed. 2) Events where the type II heights agree with the CME leading front and an earlier CME is located at a height that suggests interaction. The radio enhancements and frequency jumps could be due to the merging process of the CMEs. 3) Events where the type II heights are significantly lower than the CME heights almost from the start. Interaction with close-by streamers is probably the cause for the enhanced radio emission, which is located at the CME flank region. 4) Events where the radio enhancements are located within wide-band type II bursts and the causes for the radio enhancements are not clear. 5) Events where the radio enhancements are associated with later-accelerated particles (electron beams, observed as type III bursts) that stop at the type II burst emission lane, and no other obvious reason for the enhancement can be identified.Most of the events (38%) were due to shock–streamer interaction, while one quarter of the events was due to possible CME–CME interaction. The drift rates, bandwidth characteristics, or cross-correlations of various characteristics did not reveal any clear association with particular category types. The chosen atmospheric density model causes the largest uncertainties in the derived radio heights, although in some cases, the emission bandwidths also lead to relatively large error margins.Our conclusion is that the enhanced radio emission associated with CMEs and propagating shocks can have different origins, depending on their overall configuration and the associated processes.     

Is use of oil-field brine as a dust-abating agent really benign? Tracing the source and flowpath of contamination by oil brine in a shallow phreatic aquifer

In Ohio, 1985 H.R. Bill 501 authorizes the local authorities to issue permits to use oil-field brine surface spreading as a dust and ice control agent. Such permits are usually given without any necessary hydrogeological expertise as to the potential impact on the shallow aquifers, particularly if the brine spreading occurs in a recharge area. One such case occurred recently in the southwestern suburb of the City of Wooster, OH, USA, where a group of home owners were seeking judicial relief when their water wells began yielding salty water as a result of brine spreading on a nearby large open storage area for oil and gas well drilling supplies. The defendant, owner of the storage yard, acted in accordance with the permit issued by the local authorities. Yet, decentralized decision making and an increased emphasis on local and citizen involvement have created a gap between science and society. The local authorities were not required by law to condition the issuance of the permit on an analysis of the potential environmental impact. The decision to issue that permit was made purely on non-scientific grounds. Therefore, the residents had no choice but to embark on a costly process of proving the damages, causation, and liability in court of law. During the protracted, 6-year long conflict, an extensive hydrogeological and hydrochemical data set (including stable isotopes’ analyses along with the complete chemical analyses of major and minor constituents) was amassed at a great cost to both sides. The following article presents the hydrogeological and hydrochemical interpretation of this data set. In addition, the case illustrates a glaring deficiency in the Ohio, 1985 H.R. Bill 501, and one pitfall in the “politically correct” philosophy of delegating decision-making process entirely to the non-professional local authorities.     

Behaviour of deep immersed tunnel under combined normal fault rupture deformation and subsequent seismic shaking

Immersed tunnels are particularly sensitive to tensile and compressive deformations such as those imposed by a normal seismogenic fault rupturing underneath, and those generated by the dynamic response due to seismic waves. The paper investigates the response of a future 70 m deep immersed tunnel to the consecutive action of a major normal fault rupturing in an earthquake occurring in the basement rock underneath the tunnel, and a subsequent strong excitation from a different large-magnitude seismic event that may occur years later. Non-linear finite elements model the quasi-static fault rupture propagation through the thick soil deposit overlying the bedrock and the ensuing interaction of the rupture with the immersed tunnel. It is shown that despite imposed bedrock offset of 2 m, net tension or excessive compression between tunnel segments could be avoided with a suitable design of the joint gaskets. Then, the already deformed (“injured”) structure is subjected to strong asynchronous seismic shaking. The thick-walled tunnel is modelled as a 3-D massive flexural beam connected to the soil through properly-calibrated nonlinear interaction springs and dashpots, the supports of which are subjected to the free-field acceleration time histories. The latter, obtained with 1-D wave propagation analysis, are then modified to account for wave passage effects. The joints between tunnel segments are modeled with special non-linear hyper-elastic elements, properly accounting for their 7-bar longitudinal hydrostatic pre-stressing. Sliding is captured with special gap elements. The effect of segment length and joint properties is explored parametrically. A fascinating conclusion emerges in all analysed cases for the joints between segments that were differentially deformed after the quasi-static fault rupture: upon subsequent very strong seismic shaking, overstressed joints de-compress and understressed joints re-compress—a “healing” process that leads to a more uniform deformation profile along the tunnel. This is particularly beneficial for the precariously de-compressed joint gaskets. Hence, the safety of the immersed tunnel improves with “subsequent” strong seismic shaking!     

Leak-off mechanism and pressure prediction for shallow sediments in deepwater drilling

Deepwater sediments are prone to loss circulation in drilling due to a low overburden gradient. How to predict the magnitude of leak-off pressure more accurately is an important issue in the protection of drilling safety and the reduction of drilling cost in deep water. Starting from the mechanical properties of a shallow formation and based on the basic theory of rock-soil mechanics, the stress distribution around a borehole was analyzed. It was found that the rock or soil on a borehole is in the plastic yield state before the effective tensile stress is generated, and the effective tangential and vertical stresses increase as the drilling fluid density increases; thus, tensile failure will not occur on the borehole wall. Based on the results of stress calculation, two mechanisms and leak-off pressure prediction models for shallow sediments in deepwater drilling were put forward, and the calculated values of these models were compared with the measured value of shallow leak-off pressure in actual drilling. The results show that the MHPS (minimum horizontal principle stress) model and the FIF (fracturing in formation) model can predict the lower and upper limits of leak-off pressure. The PLC (permeable lost circulation) model can comprehensively analyze the factors influencing permeable leakage and provide a theoretical basis for leak-off prevention and plugging in deepwater drilling.     

Long-Period Variations of the Eccentricity Vector Valid also for Near Circular Orbits around a Non-Spherical Body

This paper studies the long period variations of the eccentricity vector of the orbit of an artificial satellite, under the influence of the gravity field of a central body. We use modified orbital elements which are non-singular at zero eccentricity. We expand the long periodic part of the corresponding Lagrange equations as power series of the eccentricity. The coefficients characterizing the differential system depend on the zonal coefficients of the geopotential, and on initial semi-major axis, inclination, and eccentricity. The differential equations for the components of the eccentricity vector are then integrated analytically, with a definition of the period of the perigee based on the notion of “free eccentricity”, and which is also valid for circular orbits. The analytical solution is compared to a numerical integration. This study is a generalization of (Cook, Planet. Space Sci., 14, 1966): first, the coefficients involved in the differential equations depend on all zonal coefficients (and not only on the very first ones); second, our method applies to nearly circular orbits as well as to not too eccentric orbits. Except for the critical inclination, our solution is valid for all kinds of long period motions of the perigee, i.e., circulations or librations around an equilibrium point.     

Laboratory Validation of Ultra-Soft Soil Deformation Model

Ultra-soft soil with high moisture content will experience large strain deformation under one-dimensional compression with little or no gain in effective stress. Such deformation behaviour does not comply with Terzaghi’s effective stress gain theory. The e-log s_v^￠ \sigma_{v}^{\prime } relationship of ultra-soft soil is non-linear with large compression index in the first order of log cycle. This paper proposes three compression indices (C_c1^* C_{c1}^{*} , C_c2^* C_{c2}^{*} and C_c3^* C_{c3}^{*} ) for stresses covering three log cycles. Good prediction of settlement magnitude is possible with these newly proposed compression parameters for ultra-soft soil. In addition, implicit finite difference model applying the large strain theory is also proposed and validated with results from laboratory measurements. The time factor curves for ultra-soft soil with large strain compression are also proposed and validated.