موسسه ژئوفیزیک دانشگاه تهران فیزیک زمین و فضا 2538-371X 44 4 2018 12 22 Investigation of Geostrophic and Ekman Surface Current Using Satellite Altimetry Observations and Surface Wind in Persian Gulf and Oman Sea Investigation of Geostrophic and Ekman Surface Current Using Satellite Altimetry Observations and Surface Wind in Persian Gulf and Oman Sea 1 18 64857 10.22059/jesphys.2018.244979.1006938 FA Saeed Farzaneh Assistant Professor, Department of Surveying and Geomatics Engineering, Faculty of Engineering, University of Tehran, Iran Kamal Parvazi Ph.D. Student, Department of Surveying and Geomatics Engineering, Faculty of Engineering, University of Tehran, Iran Tayebe Noroozi M.Sc. Student, Department of Surveying Engineering, faculty of Engineering, University of Zanjan, Zanjan, Iran Journal Article 2017 11 11 The rise of satellite altimetry is a revolution in the ocean sciences. Due to its global coverage and its high resolution, altimetry classically outperforms in situ water level measurement. Ekman and geostrophic currents are large parts of the ocean’s current, playing a vital role in global climate variations. According to the classic oceanography, Ekman and geostrophic currents can be calculated through the pressure gradient force as well as the friction force assuming that the water’s density is constant. Investigation of Ekman and geostrophic currents existence along with the determination of their velocities can profoundly affect the various events of oceanography and different interactive processes between the atmosphere and the ocean. Additionally, the measurement of sea currents can be useful in determination of contamination transport, seawater exchange, fisheries, oil transfer, immigration of aquatic animals and several marine activities (e.g. military, telecommunication, fishing and research activities) and also has different effects on the regional climate. In the current study, local and climatic conditions, Ekman and geostrophic currents and their velocities have been investigated based on the solution of Ekman and geostrophic equilibrium equations in the region of the Persian Gulf and the Oman Sea. To this end, using data of Saral and Jeason-2 altimetry satellites and surface wind data measured by ASCAT satellite, velocities values of v and u as well as the value and the direction of Ekman and geostrophic currents were extracted in forms of monthly data. The results were compared with obtained measurements by AVISO and NOAA for the region of the Persian Gulf and the Oman Sea, and based on the obtained results of this study, the difference in the value of these currents is about 1 cm/s. The rise of satellite altimetry is a revolution in the ocean sciences. Due to its global coverage and its high resolution, altimetry classically outperforms in situ water level measurement. Ekman and geostrophic currents are large parts of the ocean’s current, playing a vital role in global climate variations. According to the classic oceanography, Ekman and geostrophic currents can be calculated through the pressure gradient force as well as the friction force assuming that the water’s density is constant. Investigation of Ekman and geostrophic currents existence along with the determination of their velocities can profoundly affect the various events of oceanography and different interactive processes between the atmosphere and the ocean. Additionally, the measurement of sea currents can be useful in determination of contamination transport, seawater exchange, fisheries, oil transfer, immigration of aquatic animals and several marine activities (e.g. military, telecommunication, fishing and research activities) and also has different effects on the regional climate. In the current study, local and climatic conditions, Ekman and geostrophic currents and their velocities have been investigated based on the solution of Ekman and geostrophic equilibrium equations in the region of the Persian Gulf and the Oman Sea. To this end, using data of Saral and Jeason-2 altimetry satellites and surface wind data measured by ASCAT satellite, velocities values of v and u as well as the value and the direction of Ekman and geostrophic currents were extracted in forms of monthly data. The results were compared with obtained measurements by AVISO and NOAA for the region of the Persian Gulf and the Oman Sea, and based on the obtained results of this study, the difference in the value of these currents is about 1 cm/s. https://jesphys.ut.ac.ir/article_64857_ab61f43564a0f42884d7588f142826a3.pdf

موسسه ژئوفیزیک دانشگاه تهران فیزیک زمین و فضا 2538-371X 44 4 2018 12 22 Transforming Geocentric Cartesian Coordinates to Geodetic Coordinates by a New Initial Value Calculation Paradigm Transforming Geocentric Cartesian Coordinates to Geodetic Coordinates by a New Initial Value Calculation Paradigm 19 28 65883 10.22059/jesphys.2018.246251.1006946 FA Norollah Tatar Ph.D. Student, Department of Surveying and Geomatics Engineering, Faculty of Engineering, University of Tehran, Iran Saeed Farzaneh Assistant Professor, Department of Surveying and Geomatics Engineering, Faculty of Engineering, University of Tehran, Iran Journal Article 2017 11 25 Transforming geocentric Cartesian coordinates (X, Y, Z) to geodetic curvilinear coordinates (φ, λ, h) on a biaxial ellipsoid is one of the problems used in satellite positioning, coordinates conversion between reference systems, astronomy and geodetic calculations. For this purpose, various methods including Closed-form, Vector method and Fixed-point method have been developed. In this paper, a new paradigm for calculation of initial values is presented. According to the new initial values, two state of the art iterative methods are modified to calculate the geodetic height and the geodetic latitude accurately and without iteration. The results show that for those points with height values between -10 to 1,000,000 km (30-fold more than the altitude of GPS satellites), the maximum error of the calculated height and geodetic latitude is less than 1.5×10^-8 m and 1×10^-14 rad (error lower than 0.001 mm in horizontal), respectively. Transforming geocentric Cartesian coordinates (X, Y, Z) to geodetic curvilinear coordinates (φ, λ, h) on a biaxial ellipsoid is one of the problems used in satellite positioning, coordinates conversion between reference systems, astronomy and geodetic calculations. For this purpose, various methods including Closed-form, Vector method and Fixed-point method have been developed. In this paper, a new paradigm for calculation of initial values is presented. According to the new initial values, two state of the art iterative methods are modified to calculate the geodetic height and the geodetic latitude accurately and without iteration. The results show that for those points with height values between -10 to 1,000,000 km (30-fold more than the altitude of GPS satellites), the maximum error of the calculated height and geodetic latitude is less than 1.5×10^-8 m and 1×10^-14 rad (error lower than 0.001 mm in horizontal), respectively. https://jesphys.ut.ac.ir/article_65883_777e6f29c32d6b5da94c599565a81494.pdf

موسسه ژئوفیزیک دانشگاه تهران فیزیک زمین و فضا 2538-371X 44 4 2018 12 22 Large-scale Inversion of Magnetic Data Using Golub-Kahan Bidiagonalization with Truncated Generalized Cross Validation for Regularization Parameter Estimation Large-scale Inversion of Magnetic Data Using Golub-Kahan Bidiagonalization with Truncated Generalized Cross Validation for Regularization Parameter Estimation 29 39 65881 10.22059/jesphys.2018.247879.1006954 FA Saeed Vatankhah Assistant Professor, Department of Earth Physics, Institute of Geophysics, University of Tehran, Iran Journal Article 2017 12 20 In this paper a fast method for large-scale sparse inversion of magnetic data is considered. The L₁-norm stabilizer is used to generate models with sharp and distinct interfaces. To deal with the non-linearity introduced by the L₁-norm, a model-space iteratively reweighted least squares algorithm is used. The original model matrix is factorized using the Golub-Kahan bidiagonalization that projects the problem onto a Krylov subspace with a significantly reduced dimension. The model matrix of the projected system inherits the ill-conditioning of the original matrix, but the spectrum of the projected system accurately captures only a portion of the full spectrum. Equipped with the singular value decomposition of the projected system matrix, the solution of the projected problem is expressed using a filtered singular value expansion. This expansion depends on a regularization parameter which is determined using the method of Generalized Cross Validation (GCV), but here it is used for the truncated spectrum. This new technique, Truncated GCV (TGCV), is more effective compared with the standard GCV method. Numerical results using a synthetic example and real data demonstrate the efficiency of the presented algorithm. In this paper a fast method for large-scale sparse inversion of magnetic data is considered. The L₁-norm stabilizer is used to generate models with sharp and distinct interfaces. To deal with the non-linearity introduced by the L₁-norm, a model-space iteratively reweighted least squares algorithm is used. The original model matrix is factorized using the Golub-Kahan bidiagonalization that projects the problem onto a Krylov subspace with a significantly reduced dimension. The model matrix of the projected system inherits the ill-conditioning of the original matrix, but the spectrum of the projected system accurately captures only a portion of the full spectrum. Equipped with the singular value decomposition of the projected system matrix, the solution of the projected problem is expressed using a filtered singular value expansion. This expansion depends on a regularization parameter which is determined using the method of Generalized Cross Validation (GCV), but here it is used for the truncated spectrum. This new technique, Truncated GCV (TGCV), is more effective compared with the standard GCV method. Numerical results using a synthetic example and real data demonstrate the efficiency of the presented algorithm. https://jesphys.ut.ac.ir/article_65881_6379712cc0856924f302337e3f0a59a9.pdf

موسسه ژئوفیزیک دانشگاه تهران فیزیک زمین و فضا 2538-371X 44 4 2018 12 22 Application of Single-Frequency Time-Space Filtering Technique for Seismic Ground Roll and Random Noise Attenuation Application of Single-Frequency Time-Space Filtering Technique for Seismic Ground Roll and Random Noise Attenuation 41 51 65887 10.22059/jesphys.2018.249021.1006959 FA Mohammad Radad Assistant Professor, Department of Petroleum and Geophysics, Faculty of Mining, Petroleum and Geophysics Engineering, Shahrood University of Technology, Shahrood, Iran Journal Article 2017 12 31 Time-frequency filtering is an acceptable technique for attenuating noise in 2-D (time-space) and 3-D (time-space-space) reflection seismic data. The common approach for this purpose is transforming each seismic signal from 1-D time domain to a 2-D time-frequency domain and then denoising the signal by a designed filter and finally transforming back the filtered signal to original time domain. The technique is efficient for ground roll and also random noise attenuation. However, if we deal with a large data set and a great number of contaminated signals with ground roll noise, a much move consuming time will be required. In this paper, time-frequency filtering is formulated and carried out by a different approach. The data is transformed from original time-space domain into several single-frequency time-space domains, and the filters to reduce noise is designed in the new domains. The transform is easily and completely invertible. The employed time frequency analysis method is a high-resolution version of S-transform. Application to synthetic and real shot gathers confirms the good performance and efficiency of the method for attenuating ground roll noise and random noise. Time-frequency filtering is an acceptable technique for attenuating noise in 2-D (time-space) and 3-D (time-space-space) reflection seismic data. The common approach for this purpose is transforming each seismic signal from 1-D time domain to a 2-D time-frequency domain and then denoising the signal by a designed filter and finally transforming back the filtered signal to original time domain. The technique is efficient for ground roll and also random noise attenuation. However, if we deal with a large data set and a great number of contaminated signals with ground roll noise, a much move consuming time will be required. In this paper, time-frequency filtering is formulated and carried out by a different approach. The data is transformed from original time-space domain into several single-frequency time-space domains, and the filters to reduce noise is designed in the new domains. The transform is easily and completely invertible. The employed time frequency analysis method is a high-resolution version of S-transform. Application to synthetic and real shot gathers confirms the good performance and efficiency of the method for attenuating ground roll noise and random noise. https://jesphys.ut.ac.ir/article_65887_cd59b4dd8167eb112faaebf31f5c4ad8.pdf

موسسه ژئوفیزیک دانشگاه تهران فیزیک زمین و فضا 2538-371X 44 4 2018 12 22 Sensitivity analysis of time lapse gravity for monitoring fluid saturation changes in a giant multi-phase gas reservoir located in south of Iran Sensitivity analysis of time lapse gravity for monitoring fluid saturation changes in a giant multi-phase gas reservoir located in south of Iran 53 61 65878 10.22059/jesphys.2018.249232.1006961 FA Abdolmalek Khosravi M.Sc. Graduated, Department of Earth Physics, Institute of Geophysics, University of Tehran, Iran Seyed Hani Motavalli-Anbaran Assistant Professor, Department of Earth Physics, Institute of Geophysics, University of Tehran, Iran Sajad Sarallah- Zabihi Expert, National Iranian Oil Company, Exploration Directorate, Tehran, Iran Mohammad Emami Niri Assistant Professor, Institute Of Petroleum Engineering, College of Engineering, University of Tehran, Iran Journal Article 2018 01 02 The time lapse gravity method is a widely used technique to monitor the subsurface density changes in time and space. In hydrocarbon reservoirs, the density variations are due to different factors, such as: substitution of fluids with high density contrast, water influx, gas injection, and the variation in reservoir geomechanical behavior. Considering the monitoring of saturation changes in the reservoir that cannot be inferred directly by seismic survey, a forward modelling followed by a sensitivity study is performed to examine that in what conditions the saturation changes are detectable by means of 4D gravity method in the understudy reservoir. Then static and dynamic models of a giant multi-phase gas reservoir are constructed. Then, synthetic gravity data are generated after variation of production time intervals and the number of production and injection wells. In addition to detecting the gravity signal for shallower reservoirs with similar characteristics to our reservoir, a sensitivity analysis was conducted for variation in depth of the reservoir. As either the depth of the reservoir decreases or the number of the production wells and production time periods increases, the produced gravity signal is more prone to be detectable by means of modern offshore gravimeters. The gravity signal could be detected with the maximum magnitude range of 9 -  in different scenarios as a consequence of gas-water substitution, which is consistent with water drive support from surrounding aquifers. Therefore, this method is applicable for providing complementary and even independent source of information about the saturation front changes in the under-study reservoir. The time lapse gravity method is a widely used technique to monitor the subsurface density changes in time and space. In hydrocarbon reservoirs, the density variations are due to different factors, such as: substitution of fluids with high density contrast, water influx, gas injection, and the variation in reservoir geomechanical behavior. Considering the monitoring of saturation changes in the reservoir that cannot be inferred directly by seismic survey, a forward modelling followed by a sensitivity study is performed to examine that in what conditions the saturation changes are detectable by means of 4D gravity method in the understudy reservoir. Then static and dynamic models of a giant multi-phase gas reservoir are constructed. Then, synthetic gravity data are generated after variation of production time intervals and the number of production and injection wells. In addition to detecting the gravity signal for shallower reservoirs with similar characteristics to our reservoir, a sensitivity analysis was conducted for variation in depth of the reservoir. As either the depth of the reservoir decreases or the number of the production wells and production time periods increases, the produced gravity signal is more prone to be detectable by means of modern offshore gravimeters. The gravity signal could be detected with the maximum magnitude range of 9 -  in different scenarios as a consequence of gas-water substitution, which is consistent with water drive support from surrounding aquifers. Therefore, this method is applicable for providing complementary and even independent source of information about the saturation front changes in the under-study reservoir. https://jesphys.ut.ac.ir/article_65878_3298014f4b11e9f75d155c18a3778ea7.pdf

موسسه ژئوفیزیک دانشگاه تهران فیزیک زمین و فضا 2538-371X 44 4 2018 12 22 Location and dimensionality estimation of geological bodies using eigenvectors of "Computed Gravity Gradient Tensor" Location and dimensionality estimation of geological bodies using eigenvectors of "Computed Gravity Gradient Tensor" 63 71 65888 10.22059/jesphys.2018.253742.1006984 FA Korosh Karimi M.Sc. Student, Department of Physics, faculty of Science, Razi University, Kermanshah, Iran Mohsen Oveisy Moakhar Assistant Professor, Department of Physics, faculty of Science, Razi University, Kermanshah, Iran Farzad Shirzaditabar Assistant Professor, Department of Physics, faculty of Science, Razi University, Kermanshah, Iran Journal Article 2018 03 11 One of the methodologies employed in gravimetry exploration is eigenvector analysis of Gravity Gradient Tensor (GGT) which yields a solution including an estimation of a causative body’s Center of Mass (COM), dimensionality and strike direction. The eigenvectors of GGT give very rewarding clues about COM and strike direction. Additionally, the relationships between its components provide a quantity (<em>I</em>), representative of a geologic body dimensions. Although this procedure directly measures derivative components of gravity vector, it is costly and demands modern gradiometers. This study intends to obtain GGT from an ordinary gravity field measurement (g_z). This Tensor is called Computed GGT (CGGT). In this procedure, some information about a geologic mass COM, strike and rough geometry, just after an ordinary gravimetry survey, is gained. Because of derivative calculations, the impacts of noise existing in the main measured gravity field (g_z) could be destructive in CGGT solutions. Accordingly, to adjust them, a “moving twenty-five point averaging” method, and “upward continuation” are applied. The methodology is tested on various complex isolated and binary models in noisy conditions. It is also tested on real geologic example from a salt dome, USA, and all the results are highly acceptable. One of the methodologies employed in gravimetry exploration is eigenvector analysis of Gravity Gradient Tensor (GGT) which yields a solution including an estimation of a causative body’s Center of Mass (COM), dimensionality and strike direction. The eigenvectors of GGT give very rewarding clues about COM and strike direction. Additionally, the relationships between its components provide a quantity (<em>I</em>), representative of a geologic body dimensions. Although this procedure directly measures derivative components of gravity vector, it is costly and demands modern gradiometers. This study intends to obtain GGT from an ordinary gravity field measurement (g_z). This Tensor is called Computed GGT (CGGT). In this procedure, some information about a geologic mass COM, strike and rough geometry, just after an ordinary gravimetry survey, is gained. Because of derivative calculations, the impacts of noise existing in the main measured gravity field (g_z) could be destructive in CGGT solutions. Accordingly, to adjust them, a “moving twenty-five point averaging” method, and “upward continuation” are applied. The methodology is tested on various complex isolated and binary models in noisy conditions. It is also tested on real geologic example from a salt dome, USA, and all the results are highly acceptable. https://jesphys.ut.ac.ir/article_65888_67bd0e81abfc725964e7de62cba1dab8.pdf

موسسه ژئوفیزیک دانشگاه تهران فیزیک زمین و فضا 2538-371X 44 4 2018 12 22 Attenuation of spatial aliasing in CMP domain by non-linear interpolation of seismic data along local slopes Attenuation of spatial aliasing in CMP domain by non-linear interpolation of seismic data along local slopes 73 85 67749 10.22059/jesphys.2018.257443.1007005 FA Mohammad Javad Khoshnavaz Post-Doc, Department of Earth Physics, Institute of Geophysics, University of Tehran, Iran 0000-0001-7821-5923 Hamid Reza Siahkoohi Professor, Department of Earth Physics, Institute of Geophysics, University of Tehran, Iran Andrej Bóna Professor, Department of Exploration Geophysics, Curtin University of Technology, Perth, Western Australia Journal Article 2018 05 14 Spatial aliasing is an unwanted side effect that produces artifacts during seismic data processing, imaging and interpolation. It is often caused by insufficient spatial sampling of seismic data and often happens in CMP (Common Mid-Point) gather. To tackle this artifact, several techniques have been developed in time-space domain as well as frequency domain such as frequency-wavenumber, frequency-space, and frequency-time. The main advantages of seismic interpolation in time-space domain over frequency domain are: a) frequency components of the initial signals are preserved, and b) the prior knowledge that a seismic event consists of many plane wave segments, can be used. Using the later advantage, a seismic event can be predicted by pursuing the continuity of seismic events in a trace-by-trace manner. This process, which has become popular in seismic data reconstruction and imaging within the past few years, is known as predictive painting. We use predictive painting to predict the wavefronts and two-way-travel time curves in regularly sampled CMP gathers followed by increasing the number of traces by cubic interpolation. Then, the amplitude of the interpolated trace is obtained by averaging the amplitudes of the neighbouring traces. Performance of the proposed method is demonstrated on several synthetic seismic data examples as well as a field data set. Spatial aliasing is an unwanted side effect that produces artifacts during seismic data processing, imaging and interpolation. It is often caused by insufficient spatial sampling of seismic data and often happens in CMP (Common Mid-Point) gather. To tackle this artifact, several techniques have been developed in time-space domain as well as frequency domain such as frequency-wavenumber, frequency-space, and frequency-time. The main advantages of seismic interpolation in time-space domain over frequency domain are: a) frequency components of the initial signals are preserved, and b) the prior knowledge that a seismic event consists of many plane wave segments, can be used. Using the later advantage, a seismic event can be predicted by pursuing the continuity of seismic events in a trace-by-trace manner. This process, which has become popular in seismic data reconstruction and imaging within the past few years, is known as predictive painting. We use predictive painting to predict the wavefronts and two-way-travel time curves in regularly sampled CMP gathers followed by increasing the number of traces by cubic interpolation. Then, the amplitude of the interpolated trace is obtained by averaging the amplitudes of the neighbouring traces. Performance of the proposed method is demonstrated on several synthetic seismic data examples as well as a field data set. https://jesphys.ut.ac.ir/article_67749_dc6c6dbbc2c8aa9d2734cc830934f465.pdf

موسسه ژئوفیزیک دانشگاه تهران فیزیک زمین و فضا 2538-371X 44 4 2018 12 22 Kalman filter and Neural Network methods for detecting irregular variations of TEC around the time of powerful Mexico (Mw=8.2) earthquake of September 08, 2017 Kalman filter and Neural Network methods for detecting irregular variations of TEC around the time of powerful Mexico (Mw=8.2) earthquake of September 08, 2017 87 97 67738 10.22059/jesphys.2018.258251.1007007 FA Mehdi Akhoondzadeh Hanzaei Assistant Professor, Department of Remote Sensing, School of Surveying and Geospatial Engineering, College of Engineering, University of Tehran, Iran Journal Article 2018 05 28 In 98 km SW of Tres Picos in Mexico (15.022^°N, 93.899^°W, 47.40 km depth) a powerful earthquake of M_w=8.2 took place at 04:49:19 UTC (LT=UTC-05:00) on September 8, 2017. In this study, using three standard, classical and intelligent methods including median, Kalman filter, and Neural Network, respectively, the GPS Total Electron Content (TEC) measurements of three months were surveyed to detect the potential unusual variations around the time and location of Mexico earthquake. Every three implemented methods indicated a striking irregular variation of TEC at the earthquake time. However, on the earthquake day, the geomagnetic indices D_st and A_p have exceeded the allowed ranges and even reached maximum values during the studied time period. Besides, the solar index of F10.7 showed high activity around the earthquake day. Therefore, it is difficult to acknowledge the seismicity nature of the detected TEC unusual variations on earthquake day. Therefore, in this case, we encounter a mixed and complex behavior of ionosphere. In 98 km SW of Tres Picos in Mexico (15.022^°N, 93.899^°W, 47.40 km depth) a powerful earthquake of M_w=8.2 took place at 04:49:19 UTC (LT=UTC-05:00) on September 8, 2017. In this study, using three standard, classical and intelligent methods including median, Kalman filter, and Neural Network, respectively, the GPS Total Electron Content (TEC) measurements of three months were surveyed to detect the potential unusual variations around the time and location of Mexico earthquake. Every three implemented methods indicated a striking irregular variation of TEC at the earthquake time. However, on the earthquake day, the geomagnetic indices D_st and A_p have exceeded the allowed ranges and even reached maximum values during the studied time period. Besides, the solar index of F10.7 showed high activity around the earthquake day. Therefore, it is difficult to acknowledge the seismicity nature of the detected TEC unusual variations on earthquake day. Therefore, in this case, we encounter a mixed and complex behavior of ionosphere. https://jesphys.ut.ac.ir/article_67738_ca53ea2073ad26fad3988d1e06142b16.pdf

موسسه ژئوفیزیک دانشگاه تهران فیزیک زمین و فضا 2538-371X 44 4 2018 12 22 Application of Wavelet Neural Networks for Improving of Ionospheric Tomography Reconstruction over Iran Application of Wavelet Neural Networks for Improving of Ionospheric Tomography Reconstruction over Iran 99 114 67736 10.22059/jesphys.2018.245567.1006940 FA Mir Reza Ghaffari Razin Assistant Professor, Department of Surveying Engineering, Arak University of Technology, Arak, Iran 0000-0002-5579-5889 Behzad Voosoghi Associate Professor, Department of Geodesy, Faculty of Geodesy and Geomatics Engineering, K. N. Toosi Univ. of Technology, Tehran, Iran Journal Article 2017 11 20 In this paper, a new method of ionospheric tomography is developed and evaluated based on the neural networks (NN). This new method is named ITNN. In this method, wavelet neural network (WNN) with particle swarm optimization (PSO) training algorithm is used to solve some of the ionospheric tomography problems. The results of ITNN method are compared with the residual minimization training neural network (RMTNN) and modified RMTNN (MRMTNN). In all three methods, empirical orthogonal functions (EOFs) are used as a vertical objective function. To apply the methods for constructing a 3D-image of the electron density, GPS measurements of the Iranian permanent GPS network (in three days in 2007) are used. Besides, two GPS stations from international GNSS service (IGS) are used as test stations. The ionosonde data in Tehran (φ=35.73820, λ=51.38510) has been used for validating the reliability of the proposed methods. The minimum RMSE for RMTNN, MRMTNN, ITNN are 0.5312, 0.4743, 0.3465 (10¹¹ele./m³) and the minimum bias are 0.4682, 0.3890, and 0.3368 (10¹¹ele./m³) respectively. The results indicate the superiority of ITNN method over the other two methods.   In this paper, a new method of ionospheric tomography is developed and evaluated based on the neural networks (NN). This new method is named ITNN. In this method, wavelet neural network (WNN) with particle swarm optimization (PSO) training algorithm is used to solve some of the ionospheric tomography problems. The results of ITNN method are compared with the residual minimization training neural network (RMTNN) and modified RMTNN (MRMTNN). In all three methods, empirical orthogonal functions (EOFs) are used as a vertical objective function. To apply the methods for constructing a 3D-image of the electron density, GPS measurements of the Iranian permanent GPS network (in three days in 2007) are used. Besides, two GPS stations from international GNSS service (IGS) are used as test stations. The ionosonde data in Tehran (φ=35.73820, λ=51.38510) has been used for validating the reliability of the proposed methods. The minimum RMSE for RMTNN, MRMTNN, ITNN are 0.5312, 0.4743, 0.3465 (10¹¹ele./m³) and the minimum bias are 0.4682, 0.3890, and 0.3368 (10¹¹ele./m³) respectively. The results indicate the superiority of ITNN method over the other two methods.   https://jesphys.ut.ac.ir/article_67736_f39b2acfb22bf1add208fdba10fb90db.pdf

موسسه ژئوفیزیک دانشگاه تهران فیزیک زمین و فضا 2538-371X 44 4 2018 12 22 Prediction of monthly rainfall using artificial neural network mixture approach, Case Study: Torbat-e Heydariyeh Prediction of monthly rainfall using artificial neural network mixture approach, Case Study: Torbat-e Heydariyeh 115 126 67737 10.22059/jesphys.2018.244511.1006941 FA Iman Zabbah Lecturer, Department of Computer, Torbat-e Heydariyeh branch, Islamic Azad University, Torbat-e Heydariyeh, Iran Ali Reza Roshani Assistant Professor, Department of Water Engineering, Torbat-e Heydariyeh branch, Islamic Azad University, Torbat-e Heydariyeh, Iran Amin Khafage M.Sc. Graduated, Department of Computer, Torbat-e Heydariyeh branch, Islamic Azad University, Torbat-e Heydariyeh, Iran Journal Article 2017 11 12 Rainfall is one of the most important elements of water cycle used in evaluating climate conditions of each region. Long-term forecast of rainfall for arid and semi-arid regions is very important for managing and planning of water resources. To forecast appropriately, accurate data regarding humidity, temperature, pressure, wind speed etc. is required.This article is analytical and its database includes 7336 records situated in 11 features from daily brainstorm data within a twenty year period. The samples were selected based on a case study in Torbat-e Heydariyeh. 70% were chosen for learning and 30% were chosen for taking tests. From 7181 available data, 75% and 25% were used for training and evaluating, respectively. This research studied the performance of different neural networks in order to predict precipitation and then presented an algorithm for combining neural networks with linear and nonlinear methods. After modeling and comparing their results using neural networks, the root mean square error was recorded for each method. In the first modeling, the artificial neural network error was 0.05, in the second modeling, linear combination of neural networks error was 0.07, and in the third model, nonlinear combination neural networks error was 0.001. Reducing the error of forecasting precipitation has always been one of the goals of the researchers. This study, with the forecast of precipitation by neural networks, suggested that the use of a more robust method called a nonlinear combination neural network can lead to improve men is in for cast diagnostic accuracy. Rainfall is one of the most important elements of water cycle used in evaluating climate conditions of each region. Long-term forecast of rainfall for arid and semi-arid regions is very important for managing and planning of water resources. To forecast appropriately, accurate data regarding humidity, temperature, pressure, wind speed etc. is required.This article is analytical and its database includes 7336 records situated in 11 features from daily brainstorm data within a twenty year period. The samples were selected based on a case study in Torbat-e Heydariyeh. 70% were chosen for learning and 30% were chosen for taking tests. From 7181 available data, 75% and 25% were used for training and evaluating, respectively. This research studied the performance of different neural networks in order to predict precipitation and then presented an algorithm for combining neural networks with linear and nonlinear methods. After modeling and comparing their results using neural networks, the root mean square error was recorded for each method. In the first modeling, the artificial neural network error was 0.05, in the second modeling, linear combination of neural networks error was 0.07, and in the third model, nonlinear combination neural networks error was 0.001. Reducing the error of forecasting precipitation has always been one of the goals of the researchers. This study, with the forecast of precipitation by neural networks, suggested that the use of a more robust method called a nonlinear combination neural network can lead to improve men is in for cast diagnostic accuracy. https://jesphys.ut.ac.ir/article_67737_05268b9a274487cf4c088a21757bff5c.pdf

موسسه ژئوفیزیک دانشگاه تهران فیزیک زمین و فضا 2538-371X 44 4 2018 12 22 Quantification and assessment of effective of global warming on the occurrence of heat and cold waves in some selected stations in Iran Quantification and assessment of effective of global warming on the occurrence of heat and cold waves in some selected stations in Iran 127 144 64870 10.22059/jesphys.2018.245853.1006943 FA Gholam Reza Roshan Assistant Professor, Department of Geography, Golestan University, Gorgan, Iran Abdolazim Ghanghermeh Assistant Professor, Department of Geography, Golestan University, Gorgan, Iran Elham Neyazmand M.Sc. Student, Department of Geography, Golestan University, Gorgan, Iran Journal Article 2017 11 20 One of the atmospheric hazards that seriously affect human life and health is the occurrence of thermal tensions and stress in the form of heat and cold waves. Iran is one of the areas of the planet that has climate variability due to its geographical characteristics; therefore, consequently, its different regions are not immune to heat and cold waves. On the other hand, Iran's climate variability is the factor causing the difference between thresholds of heat and cold wave occurrence for its different regions. Therefore, in this study, based on three different thresholds, spatial analysis of the frequency of occurrence of heat and cold waves has been done. Thus in this work, using average daily data from 1960 to 2014, PET (Physiological Equivalent Temperature) was used to monitor heat and cold waves of four stations in Iran. Results of this study showed that in the context of global warming, although significant differences in the frequency of cold waves cannot be seen, these changes are significant and increasing for the frequency of occurrence of heat waves of selected station. One of the atmospheric hazards that seriously affect human life and health is the occurrence of thermal tensions and stress in the form of heat and cold waves. Iran is one of the areas of the planet that has climate variability due to its geographical characteristics; therefore, consequently, its different regions are not immune to heat and cold waves. On the other hand, Iran's climate variability is the factor causing the difference between thresholds of heat and cold wave occurrence for its different regions. Therefore, in this study, based on three different thresholds, spatial analysis of the frequency of occurrence of heat and cold waves has been done. Thus in this work, using average daily data from 1960 to 2014, PET (Physiological Equivalent Temperature) was used to monitor heat and cold waves of four stations in Iran. Results of this study showed that in the context of global warming, although significant differences in the frequency of cold waves cannot be seen, these changes are significant and increasing for the frequency of occurrence of heat waves of selected station. https://jesphys.ut.ac.ir/article_64870_1e9cd5a2394bfdd0a9ef377b63b96017.pdf

موسسه ژئوفیزیک دانشگاه تهران فیزیک زمین و فضا 2538-371X 44 4 2018 12 22 Spatiotemporal Variations of Total Cloud Cover and Cloud Optical Thickness in Iran Spatiotemporal Variations of Total Cloud Cover and Cloud Optical Thickness in Iran 145 164 67759 10.22059/jesphys.2018.248041.1006956 FA Mahmoud Ahmadi Associate Professor, Department of physical geography, Faculty of Earth Sciences, Shahid Beheshti University, Tehran, Iran 0000-0001-8546-0361 Abbasali Dadashiroudbari Ph.D. Student, Department of physical geography, Faculty of Earth Sciences, Shahid Beheshti University, Tehran, Iran 0000-0002-9308-1019 Hamzeh Ahmadi Ph.D. Student, Department of physical geography, Faculty of Geography and Environmental Sciences, Hakim Sabzevari University, Iran Journal Article 2018 01 06 A knowledge of cloud properties and spatiotemporal variations of clouds is especially crucial to understand the radiative forcing of climate. This research aims to study cloudiness in Iran using the most recent satellite data, powerful databases, and regional and seasonal analyses. In this study, three data series were used for the spatiotemporal variations of cloudiness in the country: A) Cloudiness data of 42 synoptic stations in the country during the statistical period from 1970 to 2005, B) Cloud Optical Thickness (COT) of Terra and Aqua MODIS sensors for 2003-2015, and C) Total Cloud Cover (TCC) of ECMWF Database, ERA-Interim version, for 1979-2015. The values obtained in the country were located via the kriging geostatistical method by RMSE. The results showed that the highest TCC occurs during the winter months. At this time of the year, the cloud cover is reduced from North to South and from West to East. Besides, COT showed that in the cold months of the year, the highest COT is observed in January and the lowest in March. The west and central Zagros highlands have the highest COT. Incorporating COT and TCC results showed that the two factors of height and approximation and access to moisture sources contribute significantly to the regional differences of cloudiness in Iran. A knowledge of cloud properties and spatiotemporal variations of clouds is especially crucial to understand the radiative forcing of climate. This research aims to study cloudiness in Iran using the most recent satellite data, powerful databases, and regional and seasonal analyses. In this study, three data series were used for the spatiotemporal variations of cloudiness in the country: A) Cloudiness data of 42 synoptic stations in the country during the statistical period from 1970 to 2005, B) Cloud Optical Thickness (COT) of Terra and Aqua MODIS sensors for 2003-2015, and C) Total Cloud Cover (TCC) of ECMWF Database, ERA-Interim version, for 1979-2015. The values obtained in the country were located via the kriging geostatistical method by RMSE. The results showed that the highest TCC occurs during the winter months. At this time of the year, the cloud cover is reduced from North to South and from West to East. Besides, COT showed that in the cold months of the year, the highest COT is observed in January and the lowest in March. The west and central Zagros highlands have the highest COT. Incorporating COT and TCC results showed that the two factors of height and approximation and access to moisture sources contribute significantly to the regional differences of cloudiness in Iran. https://jesphys.ut.ac.ir/article_67759_a2006a1aaf922ab17b6c856a48ce34c9.pdf

موسسه ژئوفیزیک دانشگاه تهران فیزیک زمین و فضا 2538-371X 44 4 2018 12 22 Periodicity of Downward Longwave Radiation at an Equatorial Location Periodicity of Downward Longwave Radiation at an Equatorial Location 165 177 67740 10.22059/jesphys.2018.253506.1006981 FA Nsikan Obot Ph.D. Graduated, Department of Physics, University of Lagos, Akoka, Lagos State, Nigeria Ibifubara Humphrey Ph.D. Graduated, Department of Physics, University of Lagos, Akoka, Lagos State, Nigeria Michael Chendo Professor, Department of Physics, University of Lagos, Akoka, Lagos State, Nigeria Elijah Oyeyemi Professor, Department of Physics, University of Lagos, Akoka, Lagos State, Nigeria Sunday Udo Professor, Department of Physics, University of Calabar, Calabar, Cross River State, Nigeria Journal Article 2018 03 03 A good understanding of the diverse mechanisms in the atmosphere is required in modelling the climate. In this study, the diurnal and seasonal patterns of all-sky downward longwave radiation (DLR) at Ilorin (8^o 32^l N, 4^o 34^l E), Nigeria were investigated alongside relative humidity (RH) and temperature. The average diurnal pattern of DLR gives an arc that begins by increasing gradually with some inherent fluctuations from 01:00 hour to a maximum at 14:00 hour local time, and then gradually decreases to a minimum at 00:00 hour. However, sometimes erratic and double peak arc diurnal DLR patterns occur around the mid and the end of the year periods respectively. The seasonal, diurnal peak of temperature and the minimum of relative humidity (RH) occur approximately two hours after the peak of DLR. Besides, the seasonal trends of both DLR and RH match each other, except sometimes in June, which could be due to the midyear DLR erratic diurnal effect. Possibly, the mechanisms of the inter-tropical discontinuity (ITD) influence the particular diurnal patterns of DLR at the midyear and year-end periods.  Moreover, monthly dispersion of DLR is high during known months of high atmospheric aerosols, and annual maximum temperature occurs after the Harmattan season. Hence, variations in DLR are influenced mostly by mineral dust in the atmosphere, mechanisms of ITD and changes in the sun-earth distance, which reflects on the different seasons in Ilorin. A good understanding of the diverse mechanisms in the atmosphere is required in modelling the climate. In this study, the diurnal and seasonal patterns of all-sky downward longwave radiation (DLR) at Ilorin (8^o 32^l N, 4^o 34^l E), Nigeria were investigated alongside relative humidity (RH) and temperature. The average diurnal pattern of DLR gives an arc that begins by increasing gradually with some inherent fluctuations from 01:00 hour to a maximum at 14:00 hour local time, and then gradually decreases to a minimum at 00:00 hour. However, sometimes erratic and double peak arc diurnal DLR patterns occur around the mid and the end of the year periods respectively. The seasonal, diurnal peak of temperature and the minimum of relative humidity (RH) occur approximately two hours after the peak of DLR. Besides, the seasonal trends of both DLR and RH match each other, except sometimes in June, which could be due to the midyear DLR erratic diurnal effect. Possibly, the mechanisms of the inter-tropical discontinuity (ITD) influence the particular diurnal patterns of DLR at the midyear and year-end periods.  Moreover, monthly dispersion of DLR is high during known months of high atmospheric aerosols, and annual maximum temperature occurs after the Harmattan season. Hence, variations in DLR are influenced mostly by mineral dust in the atmosphere, mechanisms of ITD and changes in the sun-earth distance, which reflects on the different seasons in Ilorin. https://jesphys.ut.ac.ir/article_67740_ebc2472021c9ad6640f77874afb8c8e9.pdf