A new ellipsoidal gravimetric, satellite altimetry and astronomic boundary value problem, a case study: The geoid of Iran

A new gravimetric, satellite altimetry, astronomical ellipsoidal boundary value problem for geoid computations has been developed and successfully tested. This boundary value problem has been constructed for gravity observables of the type (i) gravity potential, (ii) gravity intensity (i.e. modulus of gravity acceleration), (iii) astronomical longitude, (iv) astronomical latitude and (v) satellite altimetry observations. The ellipsoidal coordinates of the observation points have been considered as known quantities in the set-up of the problem in the light of availability of GPS coordinates. The developed boundary value problem is ellipsoidal by nature and as such takes advantage of high precision GPS observations in the set-up. The algorithmic steps of the solution of the boundary value problem are as follows:

- Application of the ellipsoidal harmonic expansion complete up to degree and order 360 and of the ellipsoidal centrifugal field for the removal of the effect of global gravity and the isostasy field from the gravity intensity and the astronomical observations at the surface of the Earth.

- Application of the ellipsoidal Newton integral on the multi-cylindrical equal-area map projection surface for the removal from the gravity intensity and the astronomical observations at the surface of the Earth the effect of the residual masses at the radius of up to 55 km from the computational point.

- Application of the ellipsoidal harmonic expansion complete up to degree and order 360 and ellipsoidal centrifugal field for the removal from the geoidal undulations derived from satellite altimetry the effect of the global gravity and isostasy on the geoidal undulations.

- Application of the ellipsoidal Newton integral on the multi-cylindrical equal-area map projection surface for the removal from the geoidal undulations derived from satellite altimetry the effect of the water masses outside the reference ellipsoid within a radius of 55 km around the computational point.

- Least squares solution of the observation equations of the incremental quantities derived from aforementioned steps in order to obtain the incremental gravity potential at the surface of the reference ellipsoid.

- The removed effects at the application points are restored on the surface of reference ellipsoid.

- Application of the ellipsoidal Bruns’ formula for converting the potential values on the surface of the reference ellipsoid into the geoidal heights with respect to the reference ellipsoid.

- Computation of the geoid of Iran has successfully tested this new methodology.

Performance of the horizontal drains in upstream shell of earth dams on the upstream slope stability during rapid drawdown conditions

Slope stability analysis during rapid drawdown is an important consideration in the design of embankment dams. During rapid drawdown, the stabilizing effect of the water on the upstream face is lost, but the pore water pressures within the embankment may remain high. As a result, the stability of the upstream face of the dam can be much reduced. Installing horizontal drains is a very efficient and cost-effective method for reducing the pore water pressure and increasing the stability of the upstream slope. The theory of horizontal drains in the upstream shell of earth dams is well established, but there seems to be limited resources available for the design of this type of horizontal drains. Hence, this study is focused on the performance of horizontal drains in the upstream shell of the slope of earth dams on the upstream slope stability during rapid drawdown conditions. The parametric study has been conducted on the variation of horizontal drain parameters such as the number of drains, their length, and their location. In this study, ten scenarios were analyzed based on different drainage configurations and the performance of each scenario is investigated on the seepage and the upstream slope stability during rapid drawdown conditions using finite element and limit equilibrium methods. The results demonstrated that the stability of the upstream slope during rapid drawdown conditions increases by increasing the number of drains. The length of drains extending further from its intersection with the critical failure surface does not provide any significant change in the factor of safety. Finally, the study also found that installing drains in the lower region of the upstream shell of earth dams gives more stability than those installed in higher elevations.     

Erosion Rate Prediction Model for Levee-Floodwall Overtopping Applications in Fine-Grained Soils

Characterizing soil erosion and predicting levee erosion rates for various levee soils and storm conditions during floodwall overtopping events is necessary in designing levee-floodwall systems. In this study, a series of laboratory scaled levee-floodwall erosion tests were conducted to determine erosion characteristics of fine grained soils subject to overtopping from different floodwall heights with variable flow-rates. A decreasing rate of erosion was observed as a pool of water was generated in the created scour hole at the crest of the levee model. The erosion rates were also assessed using jet erosion test (JET) and erosion function apparatus (EFA) tests. The results of levee-floodwall overtopping along with soil geotechnical characteristics such as plasticity index, compaction level, and saturation level of the levee soils as well as hydraulic parameters such as water overtopping velocity were used to develop a levee-floodwall erosion rate prediction model. Then, the results of JET and EFA were integrated to develop another prediction model for levee-floodwall erosion rate estimation. Consequently, the prediction models were evaluated by conducting additional tests and comparing the prediction results with the actual measured erosion rates.     

Contribution of serpentinized ultramafics to marine magnetic anomalies at slow and intermediate spreading centres: insights from the shape of the anomalies

Detailed characteristics of marine magnetic anomalies 33r and 20r suggest that the magnetization of the deeper magnetic layers, including the lower crust and possibly the uppermost mantle, is horizontally displaced with respect to that of the upper crust. We examine the possibility that serpentinization of ultramafics in the lower crust and possibly the uppermost mantle delays the acquisition of magnetization and introduces a shift between the upper- and lower-crustal magnetization patterns. Thermal evolution models and the resulting magnetization patterns of the oceanic lithosphere are calculated for a wide range of physical parameters such as the Nusselt number and the depth of hydrothermal circulation in the crust, and the temperature range of serpentinization. The models with moderate hydrothermal cooling of the whole crust and serpentinization temperatures ranging between 200 and 300 C successfully explain the anomalous skewness and the 'hook shape' of observed sea-level magnetic anomalies created at slow and intermediate spreading rates.     

IRG2018: A regional geoid model in Iran using Least Squares Collocation

Studia Geophysica et Geodaetica - The aim of this study is to determine an accurate geoid model for Iran based on the Least Squares Collocation method in the framework of the Remove — Compute...     

Partitioning of metals in sediments of the Haraz River (Southern Caspian Sea basin)

The Haraz River flows northwards through the Iranian Alborz mountains in the central region of Mazandaran province and empties into the Caspian Sea. This area has been a rich source of minerals from times immemorial. About 45 mines (coal, limestone, sand and gravel, etc.) have been operational for the last eight decades. Towards the estuary, the river receives a discharge of industrial, agricultural and urban wastes. Eight sediment samples from the Haraz River and its main tributaries were collected and analyzed for base metals as well as Sr and Fe. The chemical partitioning of metals (Cu, Zn, As, Cd, Pb, Fe, Ni, Cr, Co and Sr) in each sample was determined in four fractions (acid-soluble, reducible, oxidizable and residual). The total content of each metal was also determined. The results showed relatively higher concentrations of Cd, As, Sr and Pb in comparison to that of shale. However, based on the chemical partitioning of metals, it is found that Pb, Co, Cd and Sr are the most mobile metals. In spite of the high As concentrations in sediments, it is not likely that this element is a major hazard for the aquatic environment, since it is found mainly in the residual fraction. Cadmium is the metal that showed the highest percentages in the acid-soluble fraction (the most labile) and the lowest in the residual fraction. However, Fe, Cr and Ni are present in the greatest percentages in the residual fraction, which implies that these metals are strongly linked to the sediments.     

Three‐dimensional inverse analyses of a deep excavation in Chicago clays

Numerical models are commonly used to estimate excavation‐induced ground movements. Two‐dimensional (2D) plain strain assumption is typically used for the simulation of deep excavations which might not be suitable for excavations where three‐dimensional (3D) effects dominate the ground response. This paper adapts an inverse analysis algorithm to learn soil behavior from field measurements using a 3D model representation of an excavation. The paper describes numerical issues related to this development including the generation of the 3D model mesh from laser scan images of the excavation. The inverse analysis to extract the soil behavior in 3D is presented. The model captures the measured wall deflections. Although settlements were not sufficiently measured, the predicted settlements around the excavation site reflected strong 3D effects and were consistent with empirical correlations. Copyright © 2010 John Wiley & Sons, Ltd.     

Comparison of two inverse analysis techniques for learning deep excavation response

Performance observation is a necessary part of the design and construction process in geotechnical engineering. For deep urban excavations, empirical and numerical methods are used to predict potential deformations and their impacts on surrounding structures. Two inverse analysis approaches are described and compared for an excavation project in downtown Chicago. The first approach is a parameter optimization approach based on genetic algorithm (GA). GA is a stochastic global search technique for optimizing an objective function with linear or non-linear constraints. The second approach, self-learning simulations (SelfSim), is an inverse analysis technique that combines finite element method, continuously evolving material models, and field measurements. The optimization based on genetic algorithm approach identifies material properties of an existing soil model, and SelfSim approach extracts the underlying soil behavior unconstrained by a specific assumption on soil constitutive behavior. The two inverse analysis approaches capture well lateral wall deflections and maximum surface settlements. The GA optimization approach tends to overpredict surface settlements at some distance from the excavation as it is constrained by a specific form of the material constitutive model (i.e. hardening soil model); while the surface settlements computed using SelfSim approach match the observed ones due to its ability to learn small strain non-linearity of soil implied in the measured settlements.     

Tectonic implication of geomorphometric analyses along the Saravan Fault: evidence of a difference in tectonic movements between the Sistan Suture Zone and Makran Mountain Belt

In this paper, remote sensing techniques,as well as field studies, have been used to investigate the geomorphological processes and landscape evolution along the Saravan Fault, SE Iran to highlight how topographic features were influenced by active tectonics. Quantitative geomorphic analysis was carried out using mountain-front sinuosity(Smf),valley floor width-valley height ratio(Vf), drainage basin asymmetry factor(Af), Hypsometric integral(Hi), drainage basin shape index(Bs), mean axial slope of channel(MASC), standard deviation of topography(STD) and index of active tectonic(Iat).Remote sensing techniques, as well as field studies revealed that the Saravan Fault have three parts trending N-S, NW-SE, and E-W. Obtained results show that basins with high Iat index are located at where the strike of the Saravan Faults changes and where several strike-slip faults are crossed the Saravan fault.