Two-dimensional inversion of magnetotelluric data with consecutive use of conjugate gradient and least-squares solution with singular value decomposition algorithms

I investigated the two‐dimensional magnetotelluric data inversion algorithms in studying two significant aspects within a linearized inversion approach. The first one is the method of minimization and second one is the type of stabilizing functional used in parametric functionals. The results of two well‐known inversion algorithms, namely conjugate gradient and the least‐squares solution with singular value decomposition, were compared in terms of accuracy and CPU time. In addition, magnetotelluric data inversion with various stabilizers, such as L2‐norm, smoothing, minimum support, minimum gradient support and first‐order minimum entropy, were examined. A new inversion algorithm named least‐squares solution with singular value decomposition and conjugate gradient is suggested in seeing the outcomes of the comparisons carried out on least‐squares solutions with singular value decomposition and conjugate gradient algorithms subject to a variety of stabilizers. Inversion results of synthetic data showed that the newly suggested algorithm yields better results than those of the individual implementations of conjugate gradient and least‐squares solution with singular value decomposition algorithms. The suggested algorithm and the above‐mentioned algorithms inversion results for the field data collected along a line crossing the North Anatolian Fault zone were also compared each other and results are discussed.     

A case study for the effect of a flow improvement device (a partial wake equalizing duct) on ship powering characteristics

A case study for the energy saving in the powering characteristics of a river going general cargo ship has been carried out. Two different hull forms were generated from the original hull form of the vessel to optimise the stern flow of the vessel. A possible energy saving concept, such as partial wake equalizing duct was investigated in this manner. Resistance, self-propulsion and flow visualization measurements were performed with the hull models to explore the effect of partial wake equalizing ducts on the powering characteristics of the hull form. Analysis of the results indicates that the partial wake equalizing duct concept with an appropriate stern design affect not only the flow characteristics at aft-end, but also the propulsion characteristics. In order to identify effect of each component (partial wake equalizing duct and stern form) a further investigation is needed.     

Geoelectrical survey over perched aquifers in the northern part of Upper Sakarya River Basin,Türkiye

In this study, a groundwater exploration survey was conducted using the DC Resistivity (DCR) method in a hydrogeological setting containing a perched aquifer. DCR data were gathered and an electrical tomography section was recovered using conventional four-electrode instruments with a Schlumberger array and a two-dimensional (2D) inversion scheme. The proposed scheme was tested over a synthetic three-dimensional (3D) subsurface model before deploying it in a field situation. The proposed method indicated that gathering data with simple four-electrode instruments at stations along a line and 2D inversion of datasets at multiple stations can recover depth intervals of the studied aquifer in the hydrogeological setting even if it has a 3D structure. In this study, 2D inversion of parallel profiles formed a pseudo-3D volume of the subsurface resistivity structures and mapped out multiple resistive (>25 ohm·m) bodies at shallow (between 50–100 m) and deep sections (>150 m). In general, the proposed method is convenient to encounter geological units that have limited vertical and spatial extensions in any direction and presents resistivity contrast from groundwater-bearing geologic materials.     

Polycyclic aromatic hydrocarbons, polychlorinated biphenyls and organochlorine pesticides in urban air of Konya, Turkey

Polycyclic aromatic hydrocarbons (PAHs), polychlorinated biphenyls (PCBs) and organochlorine pesticides (OCPs) were determined in urban air samples of Konya, Turkey between August 2006 and May 2007. The concentrations of pollutants in both the gas and particulate phase were separately analysed. The average total (gas + particulate) concentrations of PAHs, PCBs and OCPs were determined as 206 ng m^− 3, 0.106 ng m^− 3, 4.78 ng m^− 3 respectively. All of the investigated target compounds were dominantly found in the gas phase except OCPs. Higher air concentrations of PAHs were found at winter season while the highest concentrations of PCBs were determined in September. The highest OCPs were detected in October and in March. In urban air of Konya, PCB 28 and PCB 52 congeners represent 46% and 35% of total PCBs while Phenanthrene, Fluoranthene, Pyrene accounted for 29%, 13%, 10% of total PAHs. HCH compounds (α + β + γ + δ-HCH), total DDTs (p,p′-DDE, p,p′-DDD, p,p′-DDT), Endosulfan compounds (Endosulfan I, Endosulfan II, Endosulfan sulfate) were dominantly determined as 30%, 21%, 20% of total OCPs respectively. Considering the relation between these compounds with temperature, there was no significant correlation observed. Despite banned/restricted use in Turkey, some OCPs were determined in urban air. These results demonstrated that they are either illegally being used in the course of agricultural activity and gardens in Konya or they are residues of past use in environment. According to these results, it can be suggested that Konya is an actively contributing region to persistent organic pollutants in Turkey.     

Deep resistivity structure of the Dikili-Bergama region, west Anatolia, revealed by two- dimensional inversion of vertical electrical sounding data

Within the framework of the National Marine Geological and Geophysical Program, we re‐examined deep vertical electrical sounding (VES) data. The data, measured in 1968 by the General Directorate of Mineral Research and Exploration (MTA) of Turkey with the aim of exploring the deep resistivity structure of the Dikili–Bergama region, focus on the geothermal potential. The geoelectrical resistivity survey was conducted using a Schlumberger array with a maximum electrode half‐spacing of 4.5 km. The two‐dimensional (2D) inversion was utilized to interpret the VES data that were collected along 15‐ to 30‐km profiles. The 2D resistivity–depth cross‐sections obtained show very low resistivity values near the Dikili and Kaynarca hot springs. The 2D inversion results also indicate the presence of fault zones striking nearly N–S and E–W, and fault‐bounded graben‐horst structures that show promising potential for geothermal field resources. The 2D gravity model, which is in good agreement with the density variation of the region, supports the resistivity structure revealed by 2D inversion. The lithology information obtained from the borehole near Kaynarca also confirms the results of the resistivity interpretation and the density model.     

Two dimensional combined inversion of short- and long-normal dc resistivity well log data

A method is presented for the two-dimensional combined inversion of short- and long-normal tool direct current resistivity data with symmetry. The forward problem is solved using the finite element method in the cylindrical coordinates system. The inverse problem is solved using a conjugate gradient technique with the partial derivatives obtained using reciprocity. The parameters were obtained by means of both conjugate gradient relaxation and conventional conjugate gradient method. The solution of this highly underdetermined inverse problem is stabilized using Tikhonov regularization and the scheme yields a blurred image of the subsurface. The scheme is tested using synthetic data and field data. Tests using synthetic data suggest that traces of the horizontal boundaries are delineated in the range of the exploration distance while the resolution of vertical boundaries depends upon the solution regularization. Application to field data shows that additional information is necessary for resolving the resistivity structure when there are low resistivity contrasts between formation units.     

Error control of iterative linear solvers for integrated groundwater models

An open problem that arises when using modern iterative linear solvers, such as the preconditioned conjugate gradient method or Generalized Minimum RESidual (GMRES) method, is how to choose the residual tolerance in the linear solver to be consistent with the tolerance on the solution error. This problem is especially acute for integrated groundwater models, which are implicitly coupled to another model, such as surface water models, and resolve both multiple scales of flow and temporal interaction terms, giving rise to linear systems with variable scaling. This article uses the theory of "forward error bound estimation" to explain the correspondence between the residual error in the preconditioned linear system and the solution error. Using examples of linear systems from models developed by the US Geological Survey and the California State Department of Water Resources, we observe that this error bound guides the choice of a practical measure for controlling the error in linear systems. We implemented a preconditioned GMRES algorithm and benchmarked it against the Successive Over-Relaxation (SOR) method, the most widely known iterative solver for nonsymmetric coefficient matrices. With forward error control, GMRES can easily replace the SOR method in legacy groundwater modeling packages, resulting in the overall simulation speedups as large as 7.74×. This research is expected to broadly impact groundwater modelers through the demonstration of a practical and general approach for setting the residual tolerance in line with the solution error tolerance and presentation of GMRES performance benchmarking results.     

Degradation of Malathion and Parathion by Ozonation,Photolytic Ozonation,and Heterogeneous Catalytic Ozonation Processes

The oxidation of organophosphorus pesticides (OPPs), such as malathion and parathion, in aqueous solution was studied using conventional ozonation (O₃), photolytic ozonation (O₃/UV, O₃/UV/H₂O₂), and heterogeneous catalytic ozonation (O₃/TiO₂/UV) processes. Experiments were performed in batch mode at laboratory scale and processes were compared in terms of disappearance kinetics. The best results of pesticide mineralization were obtained when TiO₂ particles in combination with ozone (O₃) and UV photolysis (λ = 254 nm) were applied. Decomposition of 99% of parent compounds were achieved in 10 min and oxon derivatives were completely removed in 30 min. The initial reaction rate increases linearly with increasing catalyst amount. Toxicity measurements of the treated solutions were carried out in order to evaluate the efficiency of the treatment methods. No detoxification was achieved for O₃ and O₃/UV applications. Heterogeneous photocatalytic ozonation was shown to be feasible for achieving complete decomposition of OPPs and their oxon intermediates.     

Constructing a 3D structural block diagram of the Central Basin in Marmara Sea by means of bathymetric and seismic data

In this study we made a comparative interpretation of multibeam bathymetric and seismic reflection data with different resolutions and penetration properties collected in the Central Basin of the Marmara Sea. Our main objectives were (i) to investigate and compare the active tectonic deformation observed on the sea bottom and within the uppermost sedimentary layers to that of the deep-seated deformation within the limits of resolution and penetration of the available geophysical data and (ii) to build a three-dimensional (3D) block diagram of the active tectonic and buried features by means of a sliced mapping technique. In this approach, we produced slice maps of the active and buried structural features at selected depths and then combined them to form a 3D structural block diagram. Motivation for our work was to produce a 3D structural diagram to derive a more detailed image of the structural features in the Central Basin where there is no available 3D seismic data. The observations from the bathymetry and seismic data and developed 3D diagram support the presence of a through-going strike-slip fault that forms a rotational depression zone against a right-stepping strike-slip faulting causing a pull-apart basin in the Central Depression zone.     

River suspended sediment modelling using a fuzzy logic approach

Correct estimation of sediment volume carried by a river is very important for many water resources projects. Conventional sediment rating curves, however, are not able to provide sufficiently accurate results. In this paper, a fuzzy logic approach is proposed to estimate suspended sediment concentration from streamflow. This study provides forecasting benchmarks for sediment concentration prediction in the form of a numerical and graphical comparison between fuzzy and rating‐curve models. Benchmarking was based on a 5‐year period of continuous streamflow and sediment concentration data of Quebrada Blanca Station operated by the United States Geological Survey. The benchmark results showed that the fuzzy model was able to produce much better results than rating‐curve models. The fuzzy model proposed in the study is site specific and does not simulate the hysteresis effects. Copyright © 2006 John Wiley & Sons, Ltd.