Institute of Geophysics, University of TehranJournal of the Earth and Space Physics2538-371X48420230220Geophysical Surveys for Saltwater Intrusion Assessment Using Electrical Resistivity Tomography and Electromagnetic Induction MethodsGeophysical Surveys for Saltwater Intrusion Assessment Using Electrical Resistivity Tomography and Electromagnetic Induction Methods1208927510.22059/jesphys.2022.324755.1007328FADanialMansourianCorresponding Author, Department of Geology and Environment, Faculty of Science, Ghent University, Ghent, Belgium. E-mail: dmansourian@yahoo.co.uk0000-0003-1503-6199AsiehHamidiDepartment of Geology, Faculty of Science, Ferdowsi University of Mashhad, Mashhad, Iran. E-mail: asiehhamidi@yahoo.com0000-0001-6026-4348EsmaelMakarianDepartment of Petroleum Exploration, Faculty of Mining Engineering, Sahand University of Technology, Tabriz, Iran. E-mail: esmael.makarian@gmail.com0000-0001-5198-9398PedramNamazifardDepartment of Petroleum Exploration, Faculty of Mining Engineering, Sahand University of Technology, Tabriz, Iran. E-mail: pe_namazifard@sut.ac.ir0000-0002-4455-1926MaryamMirhashemiDepartment of Petroleum Exploration, Faculty of Mining Engineering, Sahand University of Technology, Tabriz, Iran. E-mail: mirhashemi.maryamm@yahoo.com0000-0002-2458-9337Journal Article20210714Saltwater intrusion is as an environmental hazard in coastal lines if not appropriately managed. The over-exploitation, over-population and climate change have invited and pushed the saltwater landwards and polluted the freshwater aquifers. This research studies the results of the implemented project at the coast of Saint Andre' located in Koksijde, Belgium, to study this phenomenon through near-surface geophysics. Two geophysical methods, including Electrical Resistivity Tomography (ERT) and Electromagnetic Induction (EMI) were used to identify the saltwater intrusion. The present study aimed to investigate the possibility of saltwater intrusion, its extension and assess the government reclamation attempts to push back the saltwater. In the inversions, the Depth of Investigation Index (DOI) and the topography effect were evaluated. The subsurface conductivity of both methods was compared. The reliability of both methods to identify the saltwater intrusion has been established; however, the ERT survey provided a more comprehensive visualization than the EMI. The saltwater intrusion was found in the first 80 m of the coastal line with resistivity values of 2 to 5 Ohm.m; however, the infiltration of freshwater and the reclamation operation have stopped the further progress salinity into the dunes. Local possibilities of brackish water or clay lenses were identified with 7 to 25 Ohm.m resistivity values. The freshwater body was observed at distances between 120 and 220 m of the ERT line with values between 46 and 136 Ohm.m. The results were correlated with other studies, proving the reliability of the models.Saltwater intrusion is as an environmental hazard in coastal lines if not appropriately managed. The over-exploitation, over-population and climate change have invited and pushed the saltwater landwards and polluted the freshwater aquifers. This research studies the results of the implemented project at the coast of Saint Andre' located in Koksijde, Belgium, to study this phenomenon through near-surface geophysics. Two geophysical methods, including Electrical Resistivity Tomography (ERT) and Electromagnetic Induction (EMI) were used to identify the saltwater intrusion. The present study aimed to investigate the possibility of saltwater intrusion, its extension and assess the government reclamation attempts to push back the saltwater. In the inversions, the Depth of Investigation Index (DOI) and the topography effect were evaluated. The subsurface conductivity of both methods was compared. The reliability of both methods to identify the saltwater intrusion has been established; however, the ERT survey provided a more comprehensive visualization than the EMI. The saltwater intrusion was found in the first 80 m of the coastal line with resistivity values of 2 to 5 Ohm.m; however, the infiltration of freshwater and the reclamation operation have stopped the further progress salinity into the dunes. Local possibilities of brackish water or clay lenses were identified with 7 to 25 Ohm.m resistivity values. The freshwater body was observed at distances between 120 and 220 m of the ERT line with values between 46 and 136 Ohm.m. The results were correlated with other studies, proving the reliability of the models.https://jesphys.ut.ac.ir/article_89275_a7dd1123e0d3aecbe3f1554dec705c45.pdfInstitute of Geophysics, University of TehranJournal of the Earth and Space Physics2538-371X48420230220Airborne Geophysical Data Interpretation of Nkalagu and Abakaliki Regions of the Lower Benue Trough, Nigeria: Implication for Mineral PotentialityAirborne Geophysical Data Interpretation of Nkalagu and Abakaliki Regions of the Lower Benue Trough, Nigeria: Implication for Mineral Potentiality21328691610.22059/jesphys.2022.331487.1007368FAEjike K.NnaemekaDepartment of Physics and Astronomy, Faculty of Physical Sciences, University of Nigeria, Nsukka, Enugu State, Nigeria. E-mail: ejikekings24@gmail.com0000-0002-8050-9047Johnson C.IbuotCorresponding Author, Department of Physics and Astronomy, Faculty of Physical Sciences, University of Nigeria, Nsukka, Enugu State, Nigeria. E-mail: johnson.ibuot@unn.edu.ng0000-0002-1655-6471Daniel N.ObioraDepartment of Physics and Astronomy, Faculty of Physical Sciences, University of Nigeria, Nsukka, Enugu State, Nigeria. E-mail: daniel.obiora@unn.edu.ng0000-0002-9893-0673SuleimanTaufiqDepartment of Science Education, Waziri Umaru Federal Polytechnic, Birnin Kebbi, Kebbi State, E-mail: suleimantaufiq@gmail.com0000-0003-2924-8678Journal Article20211012This study was carried out to delineate possible mineralized zones within Nkalagu and Abakaliki<strong> </strong>by mapping the structural and hydrothermal alteration zones deduced from the available aeromagnetic and aeroradiometric datasets. Magnetic enhancement techniques such as total magnetic intensity (TMI), reduction to the equator (RTE), analytical signal (AS) and center for exploration targeting (CET) were utilized. The Potassium, Thorium, Uranium, ternary image maps and the K/Th ratio map were produced to aid the interpretation process of alteration areas. The results obtained from the analysis of the airborne magnetic data identified and enhanced the magnetic signatures that reflect the structural features (lineament) of the study area, by revealing the structural trends at the ENE-WSW, NE-SW and WNW-ESE as major trends and NNE-SSW, NW-SE, NNW-SSE as minor trends. The analysis of radiometric data revealed the concentrations of Thorium (eTh), Uranium (eU), and Potassium (%K) concentrations, which were used for the classification of the rock types present in the area. These classifications identified intrusions of basically igneous rocks such as granite, gabbro, rhyolite, diabase, and metamorphic rocks such as quartzite and schist. The areas believed to be hydrothermally altered, aligned NW-SW, NNE-SSW, SE, central potions, and NE borders were identified based on the concentration of radioelements, using K/Th ratio and the ternary maps. Hence, the results obtained from the analysis of the two methods mapped geological structures, geological boundaries, and alteration areas that could be target areas of possible mineral deposits.This study was carried out to delineate possible mineralized zones within Nkalagu and Abakaliki<strong> </strong>by mapping the structural and hydrothermal alteration zones deduced from the available aeromagnetic and aeroradiometric datasets. Magnetic enhancement techniques such as total magnetic intensity (TMI), reduction to the equator (RTE), analytical signal (AS) and center for exploration targeting (CET) were utilized. The Potassium, Thorium, Uranium, ternary image maps and the K/Th ratio map were produced to aid the interpretation process of alteration areas. The results obtained from the analysis of the airborne magnetic data identified and enhanced the magnetic signatures that reflect the structural features (lineament) of the study area, by revealing the structural trends at the ENE-WSW, NE-SW and WNW-ESE as major trends and NNE-SSW, NW-SE, NNW-SSE as minor trends. The analysis of radiometric data revealed the concentrations of Thorium (eTh), Uranium (eU), and Potassium (%K) concentrations, which were used for the classification of the rock types present in the area. These classifications identified intrusions of basically igneous rocks such as granite, gabbro, rhyolite, diabase, and metamorphic rocks such as quartzite and schist. The areas believed to be hydrothermally altered, aligned NW-SW, NNE-SSW, SE, central potions, and NE borders were identified based on the concentration of radioelements, using K/Th ratio and the ternary maps. Hence, the results obtained from the analysis of the two methods mapped geological structures, geological boundaries, and alteration areas that could be target areas of possible mineral deposits.https://jesphys.ut.ac.ir/article_86916_aeae3421b985d1d3f7534bb9a2af5544.pdfInstitute of Geophysics, University of TehranJournal of the Earth and Space Physics2538-371X48420230220Spatial and Temporal Variation of Seismic Ambient Noise in Tehran Region for Frequency Range 1-30 HzSpatial and Temporal Variation of Seismic Ambient Noise in Tehran Region for Frequency Range 1-30 Hz33459078410.22059/jesphys.2023.332817.1007377FAMehrasaMasihDepartment of Seismology, Institute of Geophysics, University of Tehran, Tehran, Iran. E-mail: mehrasa.masih@ut.ac.ir0000-0002-8976-7501Zaher HosseinShomaliCorresponding Author, Department of Seismology, Institute of Geophysics, University of Tehran, Tehran, Iran.
E-mail: shomali@ut.ac.ir0000-0001-6254-7560EsmaeilBayramnejadDepartment of Seismology, Institute of Geophysics, University of Tehran, Tehran, Iran. E-mail: ebayram@ut.ac.ir0000-0002-9683-459XJournal Article20211106Because of the use of seismic ambient noise as a low-cost tool for researching subsurface structure and hazard assessments in recent decades, urban seismology has become an active research subject both with seismological<strong> </strong>objectives, as obtaining better microzonation maps in highly populated areas, and with engineering<strong> </strong>objectives, as the monitoring of traffic or the surveying of historical buildings. As a result, urban seismology has been used in the metropolis of Tehran, which is one of the world's most populated cities. The city is situated on the southern slopes of the Central Alborz Mountains on Quaternary alluvial deposits, and its southern section is situated on the northwest side of Central Iran's Great Desert, surrounded by active faults such as the Mosha Fault, North Tehran Fault and South Rey Fault.<br />By using cultural noises caused by human activities such as traffic, subway, concerts, and rituals in metropolitan settings, urban seismology explains underlying structures, enhances seismic hazard management, and zoning. As a result of the high level of noise created by human activities in Tehran, as well as the presence of subterranean structures in this metropolis, the necessity of urban seismology and seismic ambient noise approaches is clear. The data from Tehran's accelerometer networks that named Tehran Disaster Mitigation and Management Organization (TDMMO) and the Road, Housing and Urban Development Research Center (BHRC) which are equipped with a Guralp CMG-5T three component accelerometer and the power spectral density-probability density function (PSD-PDF) approach by PQLX software, which was made available to the seismologist community in June 2004 were utilized in this study to analyze the spatial and temporal fluctuations of seismic ambient noise in the frequency domain. For this purpose, the period of 10 December to 24 December 2020 was investigated. In addition, the effect of the Corona virus pandemic on the frequency and time domain level of seismic ambient noise was explored in period of 16 March to 26 March 2019 and 15 March to 25 March 2020. This period was chosen because the ancient Nowruz festival is in this period.<br />Finally, considering the energy level of cultural noise denoted by the large daytime/nighttime variation with<strong> </strong>large energy during working hours and much less during nighttime and weekends, most of the stations have been experiencing the lowest level of cultural noise between the hours of 3:30-4:30 AM across all three frequency ranges that include human movement, traffic and subway, and in terms of spatial variation, D011 station has been experiencing the lowest noise levels in all three frequency ranges due to its construction on a stone structure. Furthermore, investigations on the change of noise levels during the Corona epidemic revealed a drop of 1-4% and a rise of 1-15% for certain stations depending on their location in all three frequency bands.Because of the use of seismic ambient noise as a low-cost tool for researching subsurface structure and hazard assessments in recent decades, urban seismology has become an active research subject both with seismological<strong> </strong>objectives, as obtaining better microzonation maps in highly populated areas, and with engineering<strong> </strong>objectives, as the monitoring of traffic or the surveying of historical buildings. As a result, urban seismology has been used in the metropolis of Tehran, which is one of the world's most populated cities. The city is situated on the southern slopes of the Central Alborz Mountains on Quaternary alluvial deposits, and its southern section is situated on the northwest side of Central Iran's Great Desert, surrounded by active faults such as the Mosha Fault, North Tehran Fault and South Rey Fault.<br />By using cultural noises caused by human activities such as traffic, subway, concerts, and rituals in metropolitan settings, urban seismology explains underlying structures, enhances seismic hazard management, and zoning. As a result of the high level of noise created by human activities in Tehran, as well as the presence of subterranean structures in this metropolis, the necessity of urban seismology and seismic ambient noise approaches is clear. The data from Tehran's accelerometer networks that named Tehran Disaster Mitigation and Management Organization (TDMMO) and the Road, Housing and Urban Development Research Center (BHRC) which are equipped with a Guralp CMG-5T three component accelerometer and the power spectral density-probability density function (PSD-PDF) approach by PQLX software, which was made available to the seismologist community in June 2004 were utilized in this study to analyze the spatial and temporal fluctuations of seismic ambient noise in the frequency domain. For this purpose, the period of 10 December to 24 December 2020 was investigated. In addition, the effect of the Corona virus pandemic on the frequency and time domain level of seismic ambient noise was explored in period of 16 March to 26 March 2019 and 15 March to 25 March 2020. This period was chosen because the ancient Nowruz festival is in this period.<br />Finally, considering the energy level of cultural noise denoted by the large daytime/nighttime variation with<strong> </strong>large energy during working hours and much less during nighttime and weekends, most of the stations have been experiencing the lowest level of cultural noise between the hours of 3:30-4:30 AM across all three frequency ranges that include human movement, traffic and subway, and in terms of spatial variation, D011 station has been experiencing the lowest noise levels in all three frequency ranges due to its construction on a stone structure. Furthermore, investigations on the change of noise levels during the Corona epidemic revealed a drop of 1-4% and a rise of 1-15% for certain stations depending on their location in all three frequency bands.https://jesphys.ut.ac.ir/article_90784_7af01703b1c9f44cd8a54cb5e0347c16.pdfInstitute of Geophysics, University of TehranJournal of the Earth and Space Physics2538-371X48420230220The Magnitude Threshold for Detecting Recorded Earthquakes in Tehran's Accelerometer NetworksThe Magnitude Threshold for Detecting Recorded Earthquakes in Tehran's Accelerometer Networks47549062110.22059/jesphys.2023.332846.1007378FAMehrasaMasihDepartment of Seismology, Institute of Geophysics, University of Tehran, Tehran, Iran. E-mail: mehrasa.masih@ut.ac.ir0000-0002-8976-7501HosseinKianimehrIranian Seismological Center, Kerman, Iran. E-mail: kianimehrhossein@gmail.com0000-0002-5552-3584Zaher HosseinShomaliCorresponding Author, Department of Seismology, Institute of Geophysics, University of Tehran, Tehran, Iran. E-mail: shomali@ut.ac.ir0000-0001-6254-7560EsmaeilBayramnejadDepartment of Seismology, Institute of Geophysics, University of Tehran, Tehran, Iran. E-mail: ebayram@ut.ac.ir0000-0002-9683-459XJournal Article20211106Due to the population growth in metropolitan regions such as Tehran and the existence of the underground constructions, the importance of seismic investigation is evident to reduce damages caused by probable earthquakes. Accordingly, the precise detection of micro to medium earthquakes is effective tool for tracking fault dynamics in seismic cycles, as well as for earthquake prediction and seismic hazard assessment. In this study, the recorded ambient noise at Tehran Disaster Mitigation and Management Organization (TDMMO) as well as Road, Housing and Urban Development Research Center (BHRC) networks as an accelerometer network installed in Tehran city, have been used on the point of characterizing the noise spectrum for each station as a function of time for obtaining the detection threshold of these networks. Therefore, an indirect approach based on the signal-to-noise ratio (SNR) in the time domain, with parameterization in the frequency domain is applied. Based on SNR method, the source signature is simulated by a simple source model called a circular fault model. Thus, the signal is estimated via the Brune function as most common models for circular faults. While, to determine the noise, the real data of 13 accelerometer stations of the TDMMO and seven joint stations with the BHRC are used. In this respect, the Power Spectral Density (PSD) of noise is calculated using PQLX software in the frequency domain and then is transferred to the time domain by the Parsville theorem. Eventually, the SNR value is acquired for each station by dividing these two quantities. As a result, the minimum detectable magnitude in at least five stations with an SNR larger than 5 is 3.0 for S-waves and 3.3 for P-waves, which frequently occurs in the center of the network.<br />Another finding of these studies is to analyze the effect of spatial variations of the noise on the detection ability. For this, a constant noise is allotted to all stations, lowest observed noise level, as a result of which, the smallest magnitude detectable is 1.7 for S-waves and 2.2 for P-waves.<br />At last, the sensitivity of the detection capability to three fundamental parameters, including stress drop, focal depth and reduced time, which are assumed as constant values within the network, are investigated. In fact, these parameters are strongly affected by uncertainty and are not absolute values. Consequently, the impact of their changes was studied. In our case, it is implied that the variation in the stress drop has no effect on the detection threshold, but the focal depth and the reduced time are effectual. A 15 km variation in the focal depth, the detectable magnitude changes by 0.3 units, and by changing the reduced time from 0.015 s to 0.035 s, the detectable magnitude varies by 0.4 units in M<sub>w</sub>.Due to the population growth in metropolitan regions such as Tehran and the existence of the underground constructions, the importance of seismic investigation is evident to reduce damages caused by probable earthquakes. Accordingly, the precise detection of micro to medium earthquakes is effective tool for tracking fault dynamics in seismic cycles, as well as for earthquake prediction and seismic hazard assessment. In this study, the recorded ambient noise at Tehran Disaster Mitigation and Management Organization (TDMMO) as well as Road, Housing and Urban Development Research Center (BHRC) networks as an accelerometer network installed in Tehran city, have been used on the point of characterizing the noise spectrum for each station as a function of time for obtaining the detection threshold of these networks. Therefore, an indirect approach based on the signal-to-noise ratio (SNR) in the time domain, with parameterization in the frequency domain is applied. Based on SNR method, the source signature is simulated by a simple source model called a circular fault model. Thus, the signal is estimated via the Brune function as most common models for circular faults. While, to determine the noise, the real data of 13 accelerometer stations of the TDMMO and seven joint stations with the BHRC are used. In this respect, the Power Spectral Density (PSD) of noise is calculated using PQLX software in the frequency domain and then is transferred to the time domain by the Parsville theorem. Eventually, the SNR value is acquired for each station by dividing these two quantities. As a result, the minimum detectable magnitude in at least five stations with an SNR larger than 5 is 3.0 for S-waves and 3.3 for P-waves, which frequently occurs in the center of the network.<br />Another finding of these studies is to analyze the effect of spatial variations of the noise on the detection ability. For this, a constant noise is allotted to all stations, lowest observed noise level, as a result of which, the smallest magnitude detectable is 1.7 for S-waves and 2.2 for P-waves.<br />At last, the sensitivity of the detection capability to three fundamental parameters, including stress drop, focal depth and reduced time, which are assumed as constant values within the network, are investigated. In fact, these parameters are strongly affected by uncertainty and are not absolute values. Consequently, the impact of their changes was studied. In our case, it is implied that the variation in the stress drop has no effect on the detection threshold, but the focal depth and the reduced time are effectual. A 15 km variation in the focal depth, the detectable magnitude changes by 0.3 units, and by changing the reduced time from 0.015 s to 0.035 s, the detectable magnitude varies by 0.4 units in M<sub>w</sub>.https://jesphys.ut.ac.ir/article_90621_f80d4bad7ea1db52694b1857c5ce2312.pdfInstitute of Geophysics, University of TehranJournal of the Earth and Space Physics2538-371X48420230220Magnetotelluric Evidence for Electrical Anisotropy in the Sabalan Geothermal Reservoir, Northwestern IranMagnetotelluric Evidence for Electrical Anisotropy in the Sabalan Geothermal Reservoir, Northwestern Iran55658923810.22059/jesphys.2022.333571.1007381FAZahra SadatMirhadiDepartment of Earth Physics, Institute of Geophysics, University of Tehran, Tehran, Iran. E-mail: zahramirhadi@ut.ac.ir0000-0003-2359-7391BanafshehHabibian DehkordiCorresponding Author, Department of Earth Physics, Institute of Geophysics, University of Tehran, Tehran, Iran.
E-mail: bhabibian@ut.ac.ir0000-0002-6847-4291Journal Article20211116Ignoring electrical anisotropy, if present, results in inaccurate modeling of electromagnetic data and unreliable subsequent interpretation. Its identification through the data analysis procedure, therefore, can lead to the selection of the correct algorithm for modeling and inversion and ultimately trustworthy interpretation. In this study, a part of the magnetotelluric data acquired on Sabalan volcano, located in northwestern Iran, was examined in terms of the presence of electrical anisotropy. For this purpose, penetration depths, anisotropy coefficients, phase tensor, induction vectors, and distortion parameters were considered. The results confirm significant signatures for anisotropic features in the area in the form of different depths of two polarizations, high anisotropy magnitudes, and consistent deviation of the phase tensor main axes and anisotropy directions from the regional trend. This is consistent with the outflow direction towards the Moil valley, indicated by previous studies and the high density of fractures and faults related to the setting of hydrothermal reservoirs as the main path of fluid flow.Ignoring electrical anisotropy, if present, results in inaccurate modeling of electromagnetic data and unreliable subsequent interpretation. Its identification through the data analysis procedure, therefore, can lead to the selection of the correct algorithm for modeling and inversion and ultimately trustworthy interpretation. In this study, a part of the magnetotelluric data acquired on Sabalan volcano, located in northwestern Iran, was examined in terms of the presence of electrical anisotropy. For this purpose, penetration depths, anisotropy coefficients, phase tensor, induction vectors, and distortion parameters were considered. The results confirm significant signatures for anisotropic features in the area in the form of different depths of two polarizations, high anisotropy magnitudes, and consistent deviation of the phase tensor main axes and anisotropy directions from the regional trend. This is consistent with the outflow direction towards the Moil valley, indicated by previous studies and the high density of fractures and faults related to the setting of hydrothermal reservoirs as the main path of fluid flow.https://jesphys.ut.ac.ir/article_89238_75ced66cae14720eb4eb05b9613284e9.pdfInstitute of Geophysics, University of TehranJournal of the Earth and Space Physics2538-371X48420230220Long-Short Rupture Boundary of Coseismic Displacement Estimation Based on 30-Seconds GNSS ObservationLong-Short Rupture Boundary of Coseismic Displacement Estimation Based on 30-Seconds GNSS Observation67768927210.22059/jesphys.2022.334610.1007386FAAgan AulRizkiDepartment of Geodetic Engineering, Faculty of Engineering, University Gadjah Mada, Yogyakarta, Indonesia. E-mail: agan.aul.rizki@mail.ugm.ac.id0000-0001-8935-9121CecepPratamaCorresponding Author, Department of Geodetic Engineering, Faculty of Engineering, University Gadjah Mada, Yogyakarta, Indonesia. E-mail: cecep.pratama@ugm.ac.id0000-0002-6841-9970Leni SophiaHelianiDepartment of Geodetic Engineering, Faculty of Engineering, University Gadjah Mada, Yogyakarta, Indonesia. E-mail: lheliani@ugm.ac.id0000-0002-4242-2936AdiWibowoDepartment of Informatics, Faculty of Science and Mathematics, Diponegoro University, Semarang, Indonesia. E-mail: bowo.adi@live.undip.ac.id0000-0002-7966-1017David PrambudiSaharaGlobal Geophysics Research Group, Institut Teknologi Bandung, Bandung, Indonesia. E-mail: david.sahara@gf.itb.ac.id0000-0003-1785-8962Journal Article20211130Indonesia is located at the Pacific Ring of Fire and the meeting place of the world's four tectonic plates, which makes Indonesia to have a high tectonic activity and to be prone to earthquakes. Currently, early earthquake detection service in Indonesia is based on seismometers network. However, seismometer instruments that observe seismic waveforms might become saturated, and as such may lead to incorrect earthquake magnitude detection at an early stage. Therefore, a new approach is needed to detect earthquake coseismic information. Global Navigation Satellite System (GNSS) is a good instrument to measure the surface displacement due to an earthquake. However, previous studies in Indonesia still predominantly used daily solution data. To carry out early detection, it is not possible to use daily solutions. Therefore it is needed to use the data with a higher frequency solution than the daily solution. In this study, we used 30-second sampling rate data available from Indonesian Continuously Operating Reference Station (Ina-CORS) and Sumatran GPS Array (SuGAr). We will see how the 30-second GNSS data responds to earthquakes to estimate the value of coseismic displacement compared to daily solution data. The estimated value of this coseismic displacement can be used for earthquake early detection.Indonesia is located at the Pacific Ring of Fire and the meeting place of the world's four tectonic plates, which makes Indonesia to have a high tectonic activity and to be prone to earthquakes. Currently, early earthquake detection service in Indonesia is based on seismometers network. However, seismometer instruments that observe seismic waveforms might become saturated, and as such may lead to incorrect earthquake magnitude detection at an early stage. Therefore, a new approach is needed to detect earthquake coseismic information. Global Navigation Satellite System (GNSS) is a good instrument to measure the surface displacement due to an earthquake. However, previous studies in Indonesia still predominantly used daily solution data. To carry out early detection, it is not possible to use daily solutions. Therefore it is needed to use the data with a higher frequency solution than the daily solution. In this study, we used 30-second sampling rate data available from Indonesian Continuously Operating Reference Station (Ina-CORS) and Sumatran GPS Array (SuGAr). We will see how the 30-second GNSS data responds to earthquakes to estimate the value of coseismic displacement compared to daily solution data. The estimated value of this coseismic displacement can be used for earthquake early detection.https://jesphys.ut.ac.ir/article_89272_be6b3fd2dd9642bca7a6411fc237f91f.pdfInstitute of Geophysics, University of TehranJournal of the Earth and Space Physics2538-371X48420230220Prediction of Water Saturation by FSVM using Well Logs in a Gas FieldPrediction of Water Saturation by FSVM using Well Logs in a Gas Field77868690510.22059/jesphys.2022.334938.1007389FANastaranMoosaviDepartment of Earth Sciences, Science and Research Branch, Islamic Azad University, Tehran, Iran. E-mail: nastaran.moosavi@gmail.com0000-0002-0891-358XMajidBagheriCorresponding Author, Department of Earth Physics, Institute of Geophysics, University of Tehran, Tehran, Iran.
E-mail: majidbagheri@ut.ac.ir0000-0003-2059-0194MajidNabi-BidhendiDepartment of Earth Physics, Institute of Geophysics, University of Tehran, Tehran, Iran. E-mail: mnbhendi@ut.ac.ir0000-0002-9555-8327RezaHeidariDepartment of Earth Sciences, Science and Research Branch, Islamic Azad University, Tehran, Iran. E-mail: r.heidari61@gmail.com0000-0003-3927-5848Journal Article20211211Water saturation is one of the key petrophysical parameters that mainly affects the accuracy of initial oil estimation related to a hydrocarbon reservoir. Approximation of this parameter is inevitable since it has a high effect on economic development of hydrocarbon reservoirs. In this paper, we approximate a function, using two wells with two core data sets belonging to each well, to predict water saturation by means of Support Vector Machine (SVM) algorithm in one of the gas reservoirs in the Persian Gulf. Due to the inevitability of noise and outliers in the measured data, SVM is modified to Fuzzy SVM (FSVM). For this purpose, a membership function is applied on the points, so each data point receives a membership degree. In this case, each input point is able to contribute to the learning of decision function. In other words, FSVM is able to enhance SVM by devoting less value to noise and outliers, as a result, better models compared to SVM can be produced. In this study, application of SVM for regression purpose (Support Vector Regression) is carried out on eight logs of DT, GR, RHOB, NPHI, LLD, LLS, MSFL, PEF as input with relevant core data belonging to a gas zone. Then, we determine the coefficients based on the comparison between predicted water saturation (using both SVR and fuzzy SVR algorithm) and core data. Our results show that the predicted water saturation from fuzzy SVR and SVR are 95% and 71%, respectively (higher for fuzzy SVR than SVR).Water saturation is one of the key petrophysical parameters that mainly affects the accuracy of initial oil estimation related to a hydrocarbon reservoir. Approximation of this parameter is inevitable since it has a high effect on economic development of hydrocarbon reservoirs. In this paper, we approximate a function, using two wells with two core data sets belonging to each well, to predict water saturation by means of Support Vector Machine (SVM) algorithm in one of the gas reservoirs in the Persian Gulf. Due to the inevitability of noise and outliers in the measured data, SVM is modified to Fuzzy SVM (FSVM). For this purpose, a membership function is applied on the points, so each data point receives a membership degree. In this case, each input point is able to contribute to the learning of decision function. In other words, FSVM is able to enhance SVM by devoting less value to noise and outliers, as a result, better models compared to SVM can be produced. In this study, application of SVM for regression purpose (Support Vector Regression) is carried out on eight logs of DT, GR, RHOB, NPHI, LLD, LLS, MSFL, PEF as input with relevant core data belonging to a gas zone. Then, we determine the coefficients based on the comparison between predicted water saturation (using both SVR and fuzzy SVR algorithm) and core data. Our results show that the predicted water saturation from fuzzy SVR and SVR are 95% and 71%, respectively (higher for fuzzy SVR than SVR).https://jesphys.ut.ac.ir/article_86905_59612aa15982101481dd64d35a9975df.pdfInstitute of Geophysics, University of TehranJournal of the Earth and Space Physics2538-371X48420230220Source Fault Analyses from InSAR Data and Aftershocks for the Fin Doublet Earthquakes on 14 November 2021 in Hormozgan Province, South IranSource Fault Analyses from InSAR Data and Aftershocks for the Fin Doublet Earthquakes on 14 November 2021 in Hormozgan Province, South Iran87978927810.22059/jesphys.2022.337959.1007399FAMehdiRezapourCorresponding Author, Department of Seismology, Institute of Geophysics, University of Tehran, Tehran, Iran.
E-mail: rezapour@ut.ac.ir0000-0002-8112-7083Mohammad RezaJamalreyhaniDepartment of Seismology, Institute of Geophysics, University of Tehran, Tehran, Iran. E-mail: m.jamalreyhani@ut.ac.ir0000-0003-4181-7175Journal Article20220126The Fin doublet earthquakes with magnitudes <em>M<sub>N</sub></em> 6.2 and 6.3 occurred on November 14, 2021 in Hormozgan province, south of Iran, and were followed by many aftershocks with magnitude range between <em>M<sub>N</sub></em> 2.5 and <em>M<sub>N</sub></em> 5.0. In this study, first we conducted the InSAR analysis using the satellite data from European Space Agency (ESA) to the Fin doublet earthquakes. The displacement field obtained from interferometric observations in ascending and descending directions, exhibit ~E-W oriented pattern reverse fault. Then, we relocated aftershocks that were recorded between November 14, 2021 and January 14, 2022, using the double-difference relocation method (hypoDD). The general pattern of relocated aftershocks distribution shows a seismic zone covering an area of approximately 25×17 km<sup>2</sup>. The cross-sections through aftershock locations show that the dominant depth range of aftershocks is from ~4 to ~23 km. The results show that the focal depths decrease toward the northeastern part of the region and deeper-aftershocks are located in the southwestern part. The interferometric observations and distribution of relocated aftershocks suggest the source fault(s) of the Fin doublet earthquakes strikes along a nearly west–east line and dip toward south-to-southeast.The Fin doublet earthquakes with magnitudes <em>M<sub>N</sub></em> 6.2 and 6.3 occurred on November 14, 2021 in Hormozgan province, south of Iran, and were followed by many aftershocks with magnitude range between <em>M<sub>N</sub></em> 2.5 and <em>M<sub>N</sub></em> 5.0. In this study, first we conducted the InSAR analysis using the satellite data from European Space Agency (ESA) to the Fin doublet earthquakes. The displacement field obtained from interferometric observations in ascending and descending directions, exhibit ~E-W oriented pattern reverse fault. Then, we relocated aftershocks that were recorded between November 14, 2021 and January 14, 2022, using the double-difference relocation method (hypoDD). The general pattern of relocated aftershocks distribution shows a seismic zone covering an area of approximately 25×17 km<sup>2</sup>. The cross-sections through aftershock locations show that the dominant depth range of aftershocks is from ~4 to ~23 km. The results show that the focal depths decrease toward the northeastern part of the region and deeper-aftershocks are located in the southwestern part. The interferometric observations and distribution of relocated aftershocks suggest the source fault(s) of the Fin doublet earthquakes strikes along a nearly west–east line and dip toward south-to-southeast.https://jesphys.ut.ac.ir/article_89278_1f713138b0f2f151d3d15d812cac5cbb.pdfInstitute of Geophysics, University of TehranJournal of the Earth and Space Physics2538-371X48420230220Modelling Thermal Convection of Earth Mantle with Aspect CodeModelling Thermal Convection of Earth Mantle with Aspect Code991068924310.22059/jesphys.2022.338634.1007402FANazilaAsaadiCorresponding Author, Department of Physics, Faculty of Science, University of Zanjan, Zanjan, Iran. E-mail: nazilaasaadi@znu.ac.ir0000-0001-5297-1237SamanehNorouziDepartment of Physics, Faculty of Science, University of Zanjan, Zanjan, Iran. E-mail: samane.norouzi7@gmail.com0000-0002-6285-5739Journal Article20220212The study of mantle convection is one of the most important topics in geodynamics. Mantle convection causes the transfer of internal heat to the cold parts of the Earth, and the effects of this heat transfer are observed as the motion of tectonic plates on the Earth's surface. Earthquakes, volcanism, and mountain building at the plate margins result from the movement of tectonic plates. Although the mantle occupies a large volume of the Earth, there are many fundamental questions about mantle composition, rheology, dynamics, and history. Many of these questions remain unanswered due to our indirect observations of the mantle. A major tool to study mantle dynamics is numerically analyzing mantle convection equations. In this work, we used Aspect -short for Advanced Solver for Problems in Earth's Convection- code to simulate mantle convection. The geometric model used in the simulation is a box of 4200 km by 3000 km. Using this code, we investigated the effect of different Rayleigh numbers on controlling the mantle convection and creating mantle plumes. Results show that the number of mantle plumes increases with increasing Rayleigh number, and the rising mantle plumes become thinner with the Rayleigh number increasing. Finally, we studied the relationship between the Rayleigh number and the Nusselt number (surface heat flux). We conclude that there is a power-law relation between Rayleigh and Nusselt numbers.The study of mantle convection is one of the most important topics in geodynamics. Mantle convection causes the transfer of internal heat to the cold parts of the Earth, and the effects of this heat transfer are observed as the motion of tectonic plates on the Earth's surface. Earthquakes, volcanism, and mountain building at the plate margins result from the movement of tectonic plates. Although the mantle occupies a large volume of the Earth, there are many fundamental questions about mantle composition, rheology, dynamics, and history. Many of these questions remain unanswered due to our indirect observations of the mantle. A major tool to study mantle dynamics is numerically analyzing mantle convection equations. In this work, we used Aspect -short for Advanced Solver for Problems in Earth's Convection- code to simulate mantle convection. The geometric model used in the simulation is a box of 4200 km by 3000 km. Using this code, we investigated the effect of different Rayleigh numbers on controlling the mantle convection and creating mantle plumes. Results show that the number of mantle plumes increases with increasing Rayleigh number, and the rising mantle plumes become thinner with the Rayleigh number increasing. Finally, we studied the relationship between the Rayleigh number and the Nusselt number (surface heat flux). We conclude that there is a power-law relation between Rayleigh and Nusselt numbers.https://jesphys.ut.ac.ir/article_89243_274f2c491626b596bc6ffc7ea48c7331.pdfInstitute of Geophysics, University of TehranJournal of the Earth and Space Physics2538-371X48420230220Markov Chain Monte Carlo Non-linear Geophysical Inversion with an Improved Proposal Distribution: Application to Geo-electrical DataMarkov Chain Monte Carlo Non-linear Geophysical Inversion with an Improved Proposal Distribution: Application to Geo-electrical Data1071248924010.22059/jesphys.2022.339477.1007407FAZahraTafaghod KhabazDepartment of Earth Physics, Institute of Geophysics, University of Tehran, Tehran, Iran. E-mail: zahratafaghod73@gmail.com0000-0003-1969-1545RezaGhanatiCorresponding Author, Department of Earth Physics, Institute of Geophysics, University of Tehran, Tehran, Iran. E-mail: rghanati@ut.ac.ir0000-0003-1138-7442Journal Article20220306Geophysical inverse problems seek to provide quantitative information about geophysical characteristics of the Earth’s subsurface for indirectly related data and measurements. It is generally formulated as an ill-posed non-linear optimization problem commonly solved through deterministic gradient-based approaches. Using these methods, despite fast convergence properties, may lead to local minima as well as impend accurate uncertainty analysis. On the contrary, formulating a geophysical inverse problem in a probabilistic framework and solving it by constructing the multi-dimensional posterior probability density (PPD) allow for complete sampling of the parameter space and the uncertainty quantification. The PPD is numerically characterized using Markov Chain Monte Carlo (MCMC) approaches. However, the convergence of the MCMC algorithm (i.e. sampling efficiency) toward the target stationary distribution highly depends upon the choice of the proposal distribution. In this paper, we develop an efficient proposal distribution based on perturbing the model parameters through an eigenvalue decomposition of the model covariance matrix in a principal component space. The covariance matrix is retrieved from an initial burn-in sampling, which is itself initiated using a linearized covariance estimate. The proposed strategy is first illustrated for inversion of hydrogeological parameters and then applied to synthetic and real geo-electrical data sets. The numerical experiments demonstrate that the presented proposal distribution takes advantage of the benefits from an accelerated convergence and mixing rate compared to the conventional Gaussian proposal distribution.Geophysical inverse problems seek to provide quantitative information about geophysical characteristics of the Earth’s subsurface for indirectly related data and measurements. It is generally formulated as an ill-posed non-linear optimization problem commonly solved through deterministic gradient-based approaches. Using these methods, despite fast convergence properties, may lead to local minima as well as impend accurate uncertainty analysis. On the contrary, formulating a geophysical inverse problem in a probabilistic framework and solving it by constructing the multi-dimensional posterior probability density (PPD) allow for complete sampling of the parameter space and the uncertainty quantification. The PPD is numerically characterized using Markov Chain Monte Carlo (MCMC) approaches. However, the convergence of the MCMC algorithm (i.e. sampling efficiency) toward the target stationary distribution highly depends upon the choice of the proposal distribution. In this paper, we develop an efficient proposal distribution based on perturbing the model parameters through an eigenvalue decomposition of the model covariance matrix in a principal component space. The covariance matrix is retrieved from an initial burn-in sampling, which is itself initiated using a linearized covariance estimate. The proposed strategy is first illustrated for inversion of hydrogeological parameters and then applied to synthetic and real geo-electrical data sets. The numerical experiments demonstrate that the presented proposal distribution takes advantage of the benefits from an accelerated convergence and mixing rate compared to the conventional Gaussian proposal distribution.https://jesphys.ut.ac.ir/article_89240_079b4882d3c1064f5110c0d2bc9d1fdd.pdfInstitute of Geophysics, University of TehranJournal of the Earth and Space Physics2538-371X48420230220Calibration of IASPEI Standard Broad-band Magnitude mB for Iranian Plateau EarthquakesCalibration of IASPEI Standard Broad-band Magnitude mB for Iranian Plateau Earthquakes1251358926710.22059/jesphys.2022.345010.1007444FAMehranKianiDepartment of Seismology, Institute of Geophysics, University of Tehran, Tehran, Iran. E-mail: mehran.kiani@alumni.ut.ac.ir0000-0003-0084-226XAliMoradiCorresponding Author, Department of Seismology, Institute of Geophysics, University of Tehran, Tehran, Iran.
E-mail: asmoradi@ut.ac.ir0000-0002-0836-9027Journal Article20220626Located in Alp-Himalayan belt and an active tectonic plate, Iran is annually struck by major earthquakes. Since shallow earthquakes cause considerable loss of lives and property in this region, using any method to decrease the time of magnitude estimation of great earthquakes is very important for making a prompt decision about what to do. To achieve this aim, m<sub>B</sub> was computed as a rapid estimator for 38 earthquakes with magnitudes greater than 6 occurred in Iran and adjacent areas (24°-44°N, 42°-66°E) from 1990 to 2018. The magnitudes that estimated by using the calibration function by Saul & Bormann (2007) have a standard error of 0.49 from Mw (in this study). Therefore, m<sub>B</sub>’s calibration function was modified. As a result, the magnitudes obtained are approximately equal to those of reported Mw (a standard error of 0.18). The calibration function acquired in this study for Iran’s earthquakes is lower than the m<sub>B</sub>’s global calibration function obtained by Saul & Bormann (2007). Their difference is nearly one unit at short distances, which can be related to the earthquakes located in subduction zones and plates boundaries used by Saul & Bormann (2007) that systemically have lower stress drops than intraplate earthquakes considered here. Thus it is needed to develop improved region-specific calibration functions for m<sub>B</sub>. However, the difference became smaller at distances greater than 20°. Consequently, this method and new calibration function can be employed to estimate magnitudes as early as possible across Iran plateau.Located in Alp-Himalayan belt and an active tectonic plate, Iran is annually struck by major earthquakes. Since shallow earthquakes cause considerable loss of lives and property in this region, using any method to decrease the time of magnitude estimation of great earthquakes is very important for making a prompt decision about what to do. To achieve this aim, m<sub>B</sub> was computed as a rapid estimator for 38 earthquakes with magnitudes greater than 6 occurred in Iran and adjacent areas (24°-44°N, 42°-66°E) from 1990 to 2018. The magnitudes that estimated by using the calibration function by Saul & Bormann (2007) have a standard error of 0.49 from Mw (in this study). Therefore, m<sub>B</sub>’s calibration function was modified. As a result, the magnitudes obtained are approximately equal to those of reported Mw (a standard error of 0.18). The calibration function acquired in this study for Iran’s earthquakes is lower than the m<sub>B</sub>’s global calibration function obtained by Saul & Bormann (2007). Their difference is nearly one unit at short distances, which can be related to the earthquakes located in subduction zones and plates boundaries used by Saul & Bormann (2007) that systemically have lower stress drops than intraplate earthquakes considered here. Thus it is needed to develop improved region-specific calibration functions for m<sub>B</sub>. However, the difference became smaller at distances greater than 20°. Consequently, this method and new calibration function can be employed to estimate magnitudes as early as possible across Iran plateau.https://jesphys.ut.ac.ir/article_89267_e4f774fb589def729007dce668517973.pdfInstitute of Geophysics, University of TehranJournal of the Earth and Space Physics2538-371X484202302203D Surface Heat Flow, Low-Temperature Basins and Curie Point Depth of the Iranian Plateau: Hydrocarbon Reservoirs and Iron Deposits3D Surface Heat Flow, Low-Temperature Basins and Curie Point Depth of the Iranian Plateau: Hydrocarbon Reservoirs and Iron Deposits1371509103110.22059/jesphys.2023.348000.1007453FANaeimMousaviCorresponding Author, Department of Earth Physics, Institute of Geophysics, University of Tehran, Tehran, Iran.
E-mail: moosavi.naeim@gmail.com0000-0002-1863-4867VahidEbrahimzadeh ArdestaniDepartment of Earth Physics, Institute of Geophysics, University of Tehran, Tehran, Iran. E-mail: ebrahimz@ut.ac.ir0000-0003-3936-201XJournal Article20220912While Surface Heat Flow (SHF) is an important indicator of the hydrocarbon reservoirs and mineral potentials, the measurements over the Iranian plateau are very sparse. In light of accessing the crustal and lithospheric structure derived from a well-constrained geophysical-petrological model, this study provides a 3D SHF, Curie depth isotherm (580 ºC), Moho temperature and low-temperature sedimentary basins (<150 ºC) over the Iranian plateau and surrounding areas. We solve heat transfer equation using certain thermal boundary condition and user-defined thermophysical parameters for crust. Thermal conductivity of the lithosphere is calculated iteratively. The results indicate that the iron deposits (within the igneous provinces) are spatially correlated with highs in the 3D map of SHF (>60 mW/m<sup>2</sup>), the shallow Curie isotherm (<40 km) and warm Moho boundary (>800 ºC) where lithospheric thinning or crustal thickening occurs. SHF highs are observed in the northern part of the Zagros collision zone, Central Iran micro-continent and Kopet Dagh. The low-temperature sedimentary basins (<150 ºC) are illustrated by the lows in the 3D map of SHF (<60 mW/m<sup>2</sup>), deep Curie isotherm (>40 km), and cold Moho boundary (<800 ºC) where lithosphere thickening or crustal thinning is taken place. These basins are distributed in the Oman Sea, Persian Gulf, northern margin of the Arabian plate (Mesopotamian foreland basin), the Caspian Sea and Turan platform.While Surface Heat Flow (SHF) is an important indicator of the hydrocarbon reservoirs and mineral potentials, the measurements over the Iranian plateau are very sparse. In light of accessing the crustal and lithospheric structure derived from a well-constrained geophysical-petrological model, this study provides a 3D SHF, Curie depth isotherm (580 ºC), Moho temperature and low-temperature sedimentary basins (<150 ºC) over the Iranian plateau and surrounding areas. We solve heat transfer equation using certain thermal boundary condition and user-defined thermophysical parameters for crust. Thermal conductivity of the lithosphere is calculated iteratively. The results indicate that the iron deposits (within the igneous provinces) are spatially correlated with highs in the 3D map of SHF (>60 mW/m<sup>2</sup>), the shallow Curie isotherm (<40 km) and warm Moho boundary (>800 ºC) where lithospheric thinning or crustal thickening occurs. SHF highs are observed in the northern part of the Zagros collision zone, Central Iran micro-continent and Kopet Dagh. The low-temperature sedimentary basins (<150 ºC) are illustrated by the lows in the 3D map of SHF (<60 mW/m<sup>2</sup>), deep Curie isotherm (>40 km), and cold Moho boundary (<800 ºC) where lithosphere thickening or crustal thinning is taken place. These basins are distributed in the Oman Sea, Persian Gulf, northern margin of the Arabian plate (Mesopotamian foreland basin), the Caspian Sea and Turan platform.https://jesphys.ut.ac.ir/article_91031_ef69bd969772894f5aee6a95f0f6fd9b.pdfInstitute of Geophysics, University of TehranJournal of the Earth and Space Physics2538-371X48420230220Internet Application for Interactive Visualization of Geophysical and Space Data: Approach, Architecture, TechnologiesInternet Application for Interactive Visualization of Geophysical and Space Data: Approach, Architecture, Technologies1511609078510.22059/jesphys.2023.350281.1007467FAAndrei V.VorobevCorresponding Author, The Geophysical Center of the Russian Academy of Sciences, Moscow, Russia. E-mail: geomagnet@list.ru0000-0002-9680-5609Anatoliy A.SolovievThe Geophysical Center of the Russian Academy of Sciences, Moscow, Russia. E-mail: a.soloviev@gcras.ru0000-0002-6476-9471Vyacheslav A.PilipenkoSсhmidt Institute of Physics of the Earth of the Russian Academy of Sciences, Moscow, Russia. E-mail: pilipenko_va@mail.ru0000-0003-3056-7465Gulnara R.VorobevaDepartment of Geophysical Research Methods (Department of Geophysics), Faculty of Mining and Oil, Ufa State Petroleum Technological University, Ufa, Russia. E-mail: gulnara.vorobeva@gmail.com0000-0001-7878-9724Journal Article20221105The proposed software solution is a tool developed for the analysis, forecast and visualization of geophysical data, which is collected and provided by a set of spatially distributed heterogeneous data repositories via standard web protocols (HTTP, HTTPs, FTP, etc.). They include ground magnetic observatories and stations, satellites, as well as various numerical models based on geophysical standards and specifications. The technological stack is limited with the tool’s web-based implementation and represented by integrated client- and server-side technologies with specialized frameworks and APIs. Client-side implementation is represented by several markup, styling and interaction software technologies, which are HTML5/CSS3/JavaScript with geospatial ESRI ArcGIS API for JavaScript available as the RESTful resources. Django web framework based on the “Model – View – Controller” architectural model represents server-side implementation, where Python is the main programming language used for the application’s business logic. The complete Web-based GIS represents a web portal with a set of services providing a rich instrumentation for the appropriate geophysical data analysis, processing, and visualization. Each tool upon execution provides an interactive geospatial image, which is generated according to the user request parameters or by default date-time settings. The proposed web services are freely available at https://aurora-forecast.ru and https://geomagnetic.ru through the web browsers.The proposed software solution is a tool developed for the analysis, forecast and visualization of geophysical data, which is collected and provided by a set of spatially distributed heterogeneous data repositories via standard web protocols (HTTP, HTTPs, FTP, etc.). They include ground magnetic observatories and stations, satellites, as well as various numerical models based on geophysical standards and specifications. The technological stack is limited with the tool’s web-based implementation and represented by integrated client- and server-side technologies with specialized frameworks and APIs. Client-side implementation is represented by several markup, styling and interaction software technologies, which are HTML5/CSS3/JavaScript with geospatial ESRI ArcGIS API for JavaScript available as the RESTful resources. Django web framework based on the “Model – View – Controller” architectural model represents server-side implementation, where Python is the main programming language used for the application’s business logic. The complete Web-based GIS represents a web portal with a set of services providing a rich instrumentation for the appropriate geophysical data analysis, processing, and visualization. Each tool upon execution provides an interactive geospatial image, which is generated according to the user request parameters or by default date-time settings. The proposed web services are freely available at https://aurora-forecast.ru and https://geomagnetic.ru through the web browsers.https://jesphys.ut.ac.ir/article_90785_99baecd2aaa0b6e411293f89c90a19d1.pdfInstitute of Geophysics, University of TehranJournal of the Earth and Space Physics2538-371X48420230220The Comparison of MODIS Land Surface Temperature with Meteorological Stations Measurements in IranThe Comparison of MODIS Land Surface Temperature with Meteorological Stations Measurements in Iran1611729062010.22059/jesphys.2023.343651.1007433FAMasoudMoradiDepartment of Zrebar Lake Environmental Research, Kurdistan Studies Institute, University of Kurdistan, Sanandaj, Iran. E-mail: moradimasood@ymail.com0000-0001-8270-3057BromandSalahiCorresponding Author, Department of Physical Geography, Faculty of Social Science, University of Mohaghegh Ardabili, Ardabil, Iran. E-mail: bromand416@yahoo.com0000-0003-4826-6185Journal Article20220608Land surface temperature (LST) plays a key role in the transfer of heat to the atmosphere and to the subsurface layers of soil. This study aims at examining the determination coefficient of MODIS LST on air temperature and soil temperature at different depths of Iran. A new method was employed to create a time consistent LST from Terra and Aqua MODIS products, to eliminate the observation differences in local solar time. Preceding the production of time consistent MODIS LST for 12:30 PM, a comparison was carried out with temperature measurements of meteorological stations. The correlation of MODIS LST and Meteorological Station Measurement (hereafter MSM) demonstrate high values, especially for air temperature and 5cm-deep subsurface soil temperature (R<sup>2</sup>>0.95). The lowest value was obtained for 100cm-deep soil temperature (R<sup>2</sup>=0.83). The results of intra annual analysis revealed significant relationship between MODIS LST and MSM temperatures. In the comparison of MODIS LST with subsurface soil temperatures, the scatter plot changes from 1:1 to fusiform due to the delay in heat transfer from surface to the subsurface of soil layers. This result postulates that MODIS LST is consonant with MSM temperatures in arid and semiarid regions of Iran. Spatial variation of correlation is higher for 100cm-deep soil temperature (16%). On the contrary, for air temperature and 5cm-deep soil temperature showing the highest correlation, the spatial variation is negligible throughout Iran (6.2%). However, Root Mean Square Error (RMSE) analysis revealed LST differences from 2.43 to 24.88 ˚C throughout Iran rather than MSM temperatures.Land surface temperature (LST) plays a key role in the transfer of heat to the atmosphere and to the subsurface layers of soil. This study aims at examining the determination coefficient of MODIS LST on air temperature and soil temperature at different depths of Iran. A new method was employed to create a time consistent LST from Terra and Aqua MODIS products, to eliminate the observation differences in local solar time. Preceding the production of time consistent MODIS LST for 12:30 PM, a comparison was carried out with temperature measurements of meteorological stations. The correlation of MODIS LST and Meteorological Station Measurement (hereafter MSM) demonstrate high values, especially for air temperature and 5cm-deep subsurface soil temperature (R<sup>2</sup>>0.95). The lowest value was obtained for 100cm-deep soil temperature (R<sup>2</sup>=0.83). The results of intra annual analysis revealed significant relationship between MODIS LST and MSM temperatures. In the comparison of MODIS LST with subsurface soil temperatures, the scatter plot changes from 1:1 to fusiform due to the delay in heat transfer from surface to the subsurface of soil layers. This result postulates that MODIS LST is consonant with MSM temperatures in arid and semiarid regions of Iran. Spatial variation of correlation is higher for 100cm-deep soil temperature (16%). On the contrary, for air temperature and 5cm-deep soil temperature showing the highest correlation, the spatial variation is negligible throughout Iran (6.2%). However, Root Mean Square Error (RMSE) analysis revealed LST differences from 2.43 to 24.88 ˚C throughout Iran rather than MSM temperatures.https://jesphys.ut.ac.ir/article_90620_f3b57fb8f0f6e772f04b038a490cf5bc.pdfInstitute of Geophysics, University of TehranJournal of the Earth and Space Physics2538-371X48420230220Spatial Analysis of PM2.5 Concentration over Iraq during 2003-2020Spatial Analysis of PM2.5 Concentration over Iraq during 2003-20201731828922810.22059/jesphys.2022.344671.1007441FAM. WidSadarDepartment of Atmospheric Sciences, College of Science, Mustansiriyah University, Baghdad, Iraq. E-mail: wid9777@yahoo.com0000-0003-0493-505XH. MonimAl-JibooriCorresponding Author, Department of Atmospheric Sciences, College of Science, Mustansiriyah University, Baghdad, Iraq. E-mail: mhaljiboori@gmail.com0000-0002-0816-3918K. YaseenAl-TimimiDepartment of Atmospheric Sciences, College of Science, Mustansiriyah University, Baghdad, Iraq. E-mail: yaseen.k.abbas@gmail.com0000-0001-5820-0345Journal Article20220621The particulate matter with a diameter of 2.5 μm (PM<sub>2.5</sub>) concentration seriously impacts the environment, climate, and human health. PM<sub>2.5</sub> emissions are caused by anthropogenic or natural sources and are still a problem worldwide. In this study, monthly, seasonal, and annual spatial distributions of PM<sub>2.5</sub> concentrations were analyzed over Iraq for the period 2003 to 2020 by use of the geographic information system technique. The results indicated that the PM<sub>2.5</sub> concentration was higher in summer than in winter and autumn. The monthly mean maximum values of PM<sub>2.5</sub> occurred during June and July with 8.4 and 8.7 µg/m<sup>3</sup> over central and southern regions of Iraq, respectively, while monthly means of minimum values were observed during November over northern and western regions. These results conclude that the magnitude of outdoor PM<sub>2.5</sub> concentrations varies among seasons and regions. Also, the annual means of PM<sub>2.5</sub> were less than the standard maximum permissible limits, and their seasonal means were smaller than this limit at all seasons.The particulate matter with a diameter of 2.5 μm (PM<sub>2.5</sub>) concentration seriously impacts the environment, climate, and human health. PM<sub>2.5</sub> emissions are caused by anthropogenic or natural sources and are still a problem worldwide. In this study, monthly, seasonal, and annual spatial distributions of PM<sub>2.5</sub> concentrations were analyzed over Iraq for the period 2003 to 2020 by use of the geographic information system technique. The results indicated that the PM<sub>2.5</sub> concentration was higher in summer than in winter and autumn. The monthly mean maximum values of PM<sub>2.5</sub> occurred during June and July with 8.4 and 8.7 µg/m<sup>3</sup> over central and southern regions of Iraq, respectively, while monthly means of minimum values were observed during November over northern and western regions. These results conclude that the magnitude of outdoor PM<sub>2.5</sub> concentrations varies among seasons and regions. Also, the annual means of PM<sub>2.5</sub> were less than the standard maximum permissible limits, and their seasonal means were smaller than this limit at all seasons.https://jesphys.ut.ac.ir/article_89228_798ebdd186b3fc96d51b6e863575df27.pdfInstitute of Geophysics, University of TehranJournal of the Earth and Space Physics2538-371X48420230220The Effect of Dust Aerosols on Some Meteorological Parameters in Two Dry and Humid AreasThe Effect of Dust Aerosols on Some Meteorological Parameters in Two Dry and Humid Areas1831959096810.22059/jesphys.2023.349589.1007461FARoohalahSolgiDepartment of Space Physics, Institute of Geophysics, University of Tehran, Tehran, Iran. E-mail: roohalah.solgi@alumni.ut.ac.ir0000-0002-2635-577XMajidMazraeh FarahaniDepartment of Space Physics, Institute of Geophysics, University of Tehran, Tehran, Iran. E-mail: mazraeh@ut.ac.ir0000-0001-9449-2314MaryamGharaylouCorresponding Author, Department of Space Physics, Institute of Geophysics, University of Tehran, Tehran, Iran. E-mail: gharaylo@ut.ac.ir0000-0002-0014-9409Journal Article20221011The current study compares the effect of dust aerosols on two meteorological variables, temperature and relative humidity, in two different regions. For this purpose, AOD data from the Moderate Resolution Imaging Spectroradiometer (MODIS) were used for Kermanshah and Ahvaz from 2010 to 2015. In a subjective review, a day with the highest AOD value was highlighted. The effects of dust on temperature and relative humidity variations were investigated on the selected day and compared with a clean day. The effect of aerosols on the vertical profile of temperature shows that increasing aerosol concentrations in Kermanshah causes a rise in temperature at lower atmosphere during the day due to the absorption of solar radiation by dust aerosols and a decrease in temperature at night due to the longwave radiative cooling. Because of the high seasonal humidity in Ahvaz, the nature of the aerosols has resulted in the greenhouse effect, which has raised the temperature by absorbing radiation at night. The effect of aerosols on the vertical profile of relative humidity differs between Kermanshah and Ahvaz. The relative humidity has risen, particularly at lower levels in Ahvaz during the dusty days and nights. The increase in aerosols in both Ahvaz and Kermanshah regions had no effect on precipitation based on data from the Iran Meteorological Organization. The reason could be lack of precipitating systems in the two regions during the warm seasons.The current study compares the effect of dust aerosols on two meteorological variables, temperature and relative humidity, in two different regions. For this purpose, AOD data from the Moderate Resolution Imaging Spectroradiometer (MODIS) were used for Kermanshah and Ahvaz from 2010 to 2015. In a subjective review, a day with the highest AOD value was highlighted. The effects of dust on temperature and relative humidity variations were investigated on the selected day and compared with a clean day. The effect of aerosols on the vertical profile of temperature shows that increasing aerosol concentrations in Kermanshah causes a rise in temperature at lower atmosphere during the day due to the absorption of solar radiation by dust aerosols and a decrease in temperature at night due to the longwave radiative cooling. Because of the high seasonal humidity in Ahvaz, the nature of the aerosols has resulted in the greenhouse effect, which has raised the temperature by absorbing radiation at night. The effect of aerosols on the vertical profile of relative humidity differs between Kermanshah and Ahvaz. The relative humidity has risen, particularly at lower levels in Ahvaz during the dusty days and nights. The increase in aerosols in both Ahvaz and Kermanshah regions had no effect on precipitation based on data from the Iran Meteorological Organization. The reason could be lack of precipitating systems in the two regions during the warm seasons.https://jesphys.ut.ac.ir/article_90968_b0e6e7c3beb62791c2d7671a7127a8d5.pdfInstitute of Geophysics, University of TehranJournal of the Earth and Space Physics2538-371X48420230220The Effect of Variations of EEJ on the Ionospheric TEC at Different Longitudinal Sectors using Ground-based ObservationThe Effect of Variations of EEJ on the Ionospheric TEC at Different Longitudinal Sectors using Ground-based Observation1972209086910.22059/jesphys.2023.350185.1007464FAAlemayehu MengeshaCherkosCorresponding Author, Department of Physics, CNCS, Institute of Geophysics Space Science and Astronomy, Addis Ababa University, Addis Ababa, Ethiopia. E-mail: alexye9@yahoo.com0000-0002-1809-9059MelessewNigussieDepartment of Physics, Washera Geospace and Radar Science Research Laboratory, Bahir Dar University, Bahir Dar, Ethiopia. E-mail: melessewnigussie@yahoo.com0000-0001-7827-5117Journal Article20221114In this work, the longitudinal variations of equatorial electrojet (EEJ) and its effect on the diurnal behavior of the EIA during quiet days in the period of 2011- 2012 were investigated. EEJ has been estimated using a pair of ground-based magnetometers data from six longitudinal sectors, and the Global Positioning System (GPS) TEC have also been obtained at each longitudinal sector from three stations at Northern and Southern crests and trough regions. The statistical results show that the monthly mean variations of EIA crest are consistent with that of the strength of equatorial electrojet in most regions of the investigation. The mean EEJ and EIA crests are strongest around equinoctial months in the Peruvian and Southeast Asian sectors followed by the West African regions throughout the years investigated. The weakest EEJ peaks and TEC of EIA are observed over the Pacific sectors throughout the periods of investigation. The monthly mean characteristics of EEJ/counter electrojets (CEJ) and EIA are also presented. The results also show that the CEJ events occur more frequently in the Brazilian sectors followed by in the Peruvian and West African sectors. However, in most of the equinoctial months, the strongest equatorial EIA trough and weakest of EIA crests are observed in the Brazilian sector. The temporal extent of the well-developed EIA crest and its properties show a substantial dependence on the diurnal characteristics of the EEJ for each specific day.In this work, the longitudinal variations of equatorial electrojet (EEJ) and its effect on the diurnal behavior of the EIA during quiet days in the period of 2011- 2012 were investigated. EEJ has been estimated using a pair of ground-based magnetometers data from six longitudinal sectors, and the Global Positioning System (GPS) TEC have also been obtained at each longitudinal sector from three stations at Northern and Southern crests and trough regions. The statistical results show that the monthly mean variations of EIA crest are consistent with that of the strength of equatorial electrojet in most regions of the investigation. The mean EEJ and EIA crests are strongest around equinoctial months in the Peruvian and Southeast Asian sectors followed by the West African regions throughout the years investigated. The weakest EEJ peaks and TEC of EIA are observed over the Pacific sectors throughout the periods of investigation. The monthly mean characteristics of EEJ/counter electrojets (CEJ) and EIA are also presented. The results also show that the CEJ events occur more frequently in the Brazilian sectors followed by in the Peruvian and West African sectors. However, in most of the equinoctial months, the strongest equatorial EIA trough and weakest of EIA crests are observed in the Brazilian sector. The temporal extent of the well-developed EIA crest and its properties show a substantial dependence on the diurnal characteristics of the EEJ for each specific day.https://jesphys.ut.ac.ir/article_90869_8465260dfede6d2d51fb31d550ffa103.pdfInstitute of Geophysics, University of TehranJournal of the Earth and Space Physics2538-371X48420230220Effect of Nonextensive Perturbation on Ion Acoustic SolitonsEffect of Nonextensive Perturbation on Ion Acoustic Solitons2212299062310.22059/jesphys.2023.351257.1007470FAHamid RezaPakzad1. Corresponding Author, Department of Physics, Bojnourd Branch, Islamic Azad University, Bojnourd, Iran. E-mail: pakzad@bojnourdiau.ac.ir0000-0001-8514-3155Journal Article20221119The behavior of ion acoustic wave (IAW) is studied in an electron-ion plasma consisting of cold ions and nonextensive electrons. In this study, the reductive perturbation method is used and the new point is the non-uniformity of the nonextensive parameter in the media. We want to achieve more realistic results of ion acoustic wave behavior by better using the reductive method. In fact, the variation in the behavior of ion acoustic wave when it encounters the nonextensity perturbation region is examined. Perturbation area is a part of plasma where the nonextensivity changes slightly. Therefore, the presence of nonextensivity is introduced as the first order perturbation and the phase velocity is applied as a fixed parameter in the calculations. The modified KdV (mKdV) equation is derived to describe the behavior of the ion acoustic wave propagation in this model. The obtained equation clarifies the change of the soliton profile when moving in all through the perturbation area. Our numerical results show that part of ion acoustic waves propagates as oscillatory shock wave in the perturbed area. The results of this investigation can be helpful for understanding the behavior of ion acoustic waves in an astrophysical environment and space plasmas with varying nonextensivity (Qiu <em>et al.</em>, 2020; Silva et al., 1998; Lima <em>et al.</em>, 2000).The behavior of ion acoustic wave (IAW) is studied in an electron-ion plasma consisting of cold ions and nonextensive electrons. In this study, the reductive perturbation method is used and the new point is the non-uniformity of the nonextensive parameter in the media. We want to achieve more realistic results of ion acoustic wave behavior by better using the reductive method. In fact, the variation in the behavior of ion acoustic wave when it encounters the nonextensity perturbation region is examined. Perturbation area is a part of plasma where the nonextensivity changes slightly. Therefore, the presence of nonextensivity is introduced as the first order perturbation and the phase velocity is applied as a fixed parameter in the calculations. The modified KdV (mKdV) equation is derived to describe the behavior of the ion acoustic wave propagation in this model. The obtained equation clarifies the change of the soliton profile when moving in all through the perturbation area. Our numerical results show that part of ion acoustic waves propagates as oscillatory shock wave in the perturbed area. The results of this investigation can be helpful for understanding the behavior of ion acoustic waves in an astrophysical environment and space plasmas with varying nonextensivity (Qiu <em>et al.</em>, 2020; Silva et al., 1998; Lima <em>et al.</em>, 2000).https://jesphys.ut.ac.ir/article_90623_82a38d3501d0dc02bfacf7bcc2c1e986.pdf