Institute of Geophysics, University of TehranJournal of the Earth and Space Physics2538-371X45220190723Complexity investigation of seismicity temporal distribution in the Zagros region by using fractal analysisComplexity investigation of seismicity temporal distribution in the Zagros region by using fractal analysis2472677098410.22059/jesphys.2019.254583.1006988FASaherehGolrizM.Sc. Graduated, School of Earth Sciences, Damghan University, Damghan, IranAmir PiroozKolahi-AzarAssistant Professor, School of Earth Sciences, Damghan University, Damghan, Iran0000-0001-8005-8680Journal Article20180418In this research it is tried to examine the fractal complexity of seismicity temporal dispersion in the Zagros Mountain range. The Wavelet Transform Modoulos Maxima (WTMM) as an innovative strang attractor formalism has been utilized for the multifractal investigation. Earthquakes that occurred from December 2003 to May 2016 have been collected from the master catalog of the International Institute of Earthquake Engineering and Seismology (IIEES). As all events in the master list are reported based on the local magnitude (M<sub>L</sub>), the achieved catalog is already homogeneous. M<sub>L</sub> is saturated for the earthquakes with magnitude greater than 5.5, so they are converted to the moment magnitude (M<sub>W</sub>) using some empirical relations. For a reliable and comprehensive seismicity examination, the Gutenberg–Richter analysis is performed over the cumulative distribution of events, and the minimum magnitude of completeness (M<sub>C</sub>) has been obtained. For M<sub>C</sub> calculation, the maximum curvature method is used and an overall Mc=3.1 is computed for the attained earthquake catalog. To complete the catalog, all events with M<sub>W</sub><M<sub>C</sub> have been removed from the earthquakes list. As the occurrence time is the most reliable seismicity parameter, the time-series are prepared as interevent times between the consecutive earthquakes for the different subzones of the Zagros region. The WTMM technique has been applied to each of the time-series and their fractal characteristics are gaind from the attributes of the related scaling and singularity spectrums. The obtained results revealed that the seismicity is scale invariant; however, its multifractal nature is not constant. There are some differences among the fractal aspects of seismicity temporal changes in the different portions of the belt. Chronological distribution of earthquakes in the simply-folded belt and Dezful embayment are remarkably more complex than the other portions of the Zagros Mountain range. Dezful embayment as an indenter plays an important role on deformation style in the Zagros Mountain. It causes crust materials to escape from the frontal regions toward the Fars-Arc and Lorestan side-salients. Our findings indicate a relatively complex and heterogeneous temporal variation of earthquakes in the salients and Dezful indenter with respect to those in high-Zagros and Izeh frontal subzones. Abadan plain is the quietest subzone seismically and it shows the least amuont of temporal complexity. From the dependency point of view, the seismicity of high-Zagros, Izeh, and Abadan plain has an anticorrelate sharing. On the contrary, Fars-Arc and Lorestan salients have correlated seismic activities and in Dezful embayment the seismicity behaves in a random (stochastic) manner. These findings reveal that the seismicity offers relatively inconsistent configuration in regions with a high-stress concentration and in contrary, earthquakes work dependably in other calm areas. Generally, in the Zagros region independent (scattered) earthquakes are more heterogeneous with respect to the dependent (clustered) seismicity. In other words, the Zagros tectonic setting is such that the independent earthquakes have more intricate temporal spreading with respect to the affiliated temblors. The results of this study are in agreement with Kolahi-Azar and Golriz (2018) examination. In the mentioned work topography complexity has been measured for the different subzones of the Zagros region. Assuming the topography is affected by the superficial tectonic processes; they concluded shallow tectonic processes that act more intricately in Dezful embayment, Fars-Arc, and Lorestan side-salients. Similarly, our results show the more intricate temporal distribution of seismicity for the same regions. The fractal study of seismicity temporal distribution is a useful tool for the better understanding of the geodynamic conditions in a region. This approach reveals new seismotectonic aspects of the Zagros region which has not been addressed from this point of view.In this research it is tried to examine the fractal complexity of seismicity temporal dispersion in the Zagros Mountain range. The Wavelet Transform Modoulos Maxima (WTMM) as an innovative strang attractor formalism has been utilized for the multifractal investigation. Earthquakes that occurred from December 2003 to May 2016 have been collected from the master catalog of the International Institute of Earthquake Engineering and Seismology (IIEES). As all events in the master list are reported based on the local magnitude (M<sub>L</sub>), the achieved catalog is already homogeneous. M<sub>L</sub> is saturated for the earthquakes with magnitude greater than 5.5, so they are converted to the moment magnitude (M<sub>W</sub>) using some empirical relations. For a reliable and comprehensive seismicity examination, the Gutenberg–Richter analysis is performed over the cumulative distribution of events, and the minimum magnitude of completeness (M<sub>C</sub>) has been obtained. For M<sub>C</sub> calculation, the maximum curvature method is used and an overall Mc=3.1 is computed for the attained earthquake catalog. To complete the catalog, all events with M<sub>W</sub><M<sub>C</sub> have been removed from the earthquakes list. As the occurrence time is the most reliable seismicity parameter, the time-series are prepared as interevent times between the consecutive earthquakes for the different subzones of the Zagros region. The WTMM technique has been applied to each of the time-series and their fractal characteristics are gaind from the attributes of the related scaling and singularity spectrums. The obtained results revealed that the seismicity is scale invariant; however, its multifractal nature is not constant. There are some differences among the fractal aspects of seismicity temporal changes in the different portions of the belt. Chronological distribution of earthquakes in the simply-folded belt and Dezful embayment are remarkably more complex than the other portions of the Zagros Mountain range. Dezful embayment as an indenter plays an important role on deformation style in the Zagros Mountain. It causes crust materials to escape from the frontal regions toward the Fars-Arc and Lorestan side-salients. Our findings indicate a relatively complex and heterogeneous temporal variation of earthquakes in the salients and Dezful indenter with respect to those in high-Zagros and Izeh frontal subzones. Abadan plain is the quietest subzone seismically and it shows the least amuont of temporal complexity. From the dependency point of view, the seismicity of high-Zagros, Izeh, and Abadan plain has an anticorrelate sharing. On the contrary, Fars-Arc and Lorestan salients have correlated seismic activities and in Dezful embayment the seismicity behaves in a random (stochastic) manner. These findings reveal that the seismicity offers relatively inconsistent configuration in regions with a high-stress concentration and in contrary, earthquakes work dependably in other calm areas. Generally, in the Zagros region independent (scattered) earthquakes are more heterogeneous with respect to the dependent (clustered) seismicity. In other words, the Zagros tectonic setting is such that the independent earthquakes have more intricate temporal spreading with respect to the affiliated temblors. The results of this study are in agreement with Kolahi-Azar and Golriz (2018) examination. In the mentioned work topography complexity has been measured for the different subzones of the Zagros region. Assuming the topography is affected by the superficial tectonic processes; they concluded shallow tectonic processes that act more intricately in Dezful embayment, Fars-Arc, and Lorestan side-salients. Similarly, our results show the more intricate temporal distribution of seismicity for the same regions. The fractal study of seismicity temporal distribution is a useful tool for the better understanding of the geodynamic conditions in a region. This approach reveals new seismotectonic aspects of the Zagros region which has not been addressed from this point of view.Institute of Geophysics, University of TehranJournal of the Earth and Space Physics2538-371X45220190723Application of POCS algorithm for the reconstruction of three-component seismic data in quaternion Fourier domainApplication of POCS algorithm for the reconstruction of three-component seismic data in quaternion Fourier domain2692817096710.22059/jesphys.2019.257996.1007006FAAminEftekhariM.Sc. Graduated, Department of Geophysics, Research and Science Branch, Islamic Azad University, Tehran, IranHamid RezaSiahkoohiProfessor, Department of Earth Physics, Institute of Geophysics, University of Tehran, Tehran, Iran0000-0002-9227-2972Journal Article20180522Three-component (3C) seismic data acquisition method samples seismic wave field at each station along three Cartesian coordinates, simultaneously. Many reservoirs have been discovered and determined by the generation and recording of P waves only, but the P wave alone cannot provide a comprehensive description of the reservoir characteristics. In some studies, S-wave information is required in addition to P-wave information to get a correct estimation from reservoir properties. By the three component seismic acquisition, P and S waves’ information can be recorded simultaneously. More often in seismic surveys, one cannot sample seismic wave field uniformly the along spatial direction due to environment limitations or instrument malfunctions; inevitably we have to use interpolation methods for reconstruction of missing traces. Reconstruction of missing or noisy traces is done using the projection onto convex sets (POCS). The POCS algorithm is a simple algorithm which is suitable for reconstruction of irregularly lost traces in a regular grid using multiple repetitive Fourier transforms. Conventional methods for reconstruction of missing traces in three component acquisition is usually done by implementation of POCS on each component separately, which could damage any subtle features in the record. This research introduces a method to reconstruct all three components at once using the quaternion Fourier transform and Projection onto Convex Sets (QPOCS).Quaternions in mathematics are a commutative numbers system that extend the complex numbers system. As the ordinary complex numbers can be displayed on two dimensions, these numbers can also be displayed on four dimensions. Quaternions were first introduced by William Rowan Hamilton when looking for a way to extend complex numbers to three dimensions. He knew how to sum and multiply three-dimensional numbers, but he was looking for a way to divide these numbers into each other. In 1843, Hamilton discovered that the division of quaternions requires a fourth dimension. Quaternion Algebra is often shown with H (in honor of Hamilton). The two-component data vector representation in the frequency domain can be obtained by putting the real and imaginary parts of each component in the arguments of a quaternion. This method allows operators to apply both components simultaneously. Quaternions are converted to Frequency-wavenumber domain by Quaternion Fourier Transform (QFT) and a single domain spectrum for both components is defined using the polar representation of the Quaternions. Quaternions have other applications in seismic data processing such as computing spectral attributes, multi-component velocity analysis and multi-component deconvolution. The advantage of this method is because of the spectral overlapping of the components in the frequency-wavenumber domain, thus the perpendicularity of input components is preserved (signals are not interconnected) and similarities between components are maintained that helps improve the quality of reconstruction. The coding of this method has been done in MATLAB environment and results of applying the proposed method on 3-component synthetic and real seismic data are compared to that of the POCS algorithm when applied on each component separately. The results of reconstruction using QPOCS algorithm indicate a better quality for reconstructed seismic data and in the output data, the percentage of produced artifacts is lower than that of the POCS algorithm on each component alone.Three-component (3C) seismic data acquisition method samples seismic wave field at each station along three Cartesian coordinates, simultaneously. Many reservoirs have been discovered and determined by the generation and recording of P waves only, but the P wave alone cannot provide a comprehensive description of the reservoir characteristics. In some studies, S-wave information is required in addition to P-wave information to get a correct estimation from reservoir properties. By the three component seismic acquisition, P and S waves’ information can be recorded simultaneously. More often in seismic surveys, one cannot sample seismic wave field uniformly the along spatial direction due to environment limitations or instrument malfunctions; inevitably we have to use interpolation methods for reconstruction of missing traces. Reconstruction of missing or noisy traces is done using the projection onto convex sets (POCS). The POCS algorithm is a simple algorithm which is suitable for reconstruction of irregularly lost traces in a regular grid using multiple repetitive Fourier transforms. Conventional methods for reconstruction of missing traces in three component acquisition is usually done by implementation of POCS on each component separately, which could damage any subtle features in the record. This research introduces a method to reconstruct all three components at once using the quaternion Fourier transform and Projection onto Convex Sets (QPOCS).Quaternions in mathematics are a commutative numbers system that extend the complex numbers system. As the ordinary complex numbers can be displayed on two dimensions, these numbers can also be displayed on four dimensions. Quaternions were first introduced by William Rowan Hamilton when looking for a way to extend complex numbers to three dimensions. He knew how to sum and multiply three-dimensional numbers, but he was looking for a way to divide these numbers into each other. In 1843, Hamilton discovered that the division of quaternions requires a fourth dimension. Quaternion Algebra is often shown with H (in honor of Hamilton). The two-component data vector representation in the frequency domain can be obtained by putting the real and imaginary parts of each component in the arguments of a quaternion. This method allows operators to apply both components simultaneously. Quaternions are converted to Frequency-wavenumber domain by Quaternion Fourier Transform (QFT) and a single domain spectrum for both components is defined using the polar representation of the Quaternions. Quaternions have other applications in seismic data processing such as computing spectral attributes, multi-component velocity analysis and multi-component deconvolution. The advantage of this method is because of the spectral overlapping of the components in the frequency-wavenumber domain, thus the perpendicularity of input components is preserved (signals are not interconnected) and similarities between components are maintained that helps improve the quality of reconstruction. The coding of this method has been done in MATLAB environment and results of applying the proposed method on 3-component synthetic and real seismic data are compared to that of the POCS algorithm when applied on each component separately. The results of reconstruction using QPOCS algorithm indicate a better quality for reconstructed seismic data and in the output data, the percentage of produced artifacts is lower than that of the POCS algorithm on each component alone.Institute of Geophysics, University of TehranJournal of the Earth and Space Physics2538-371X45220190723Detection of the weakness of the vector averaging of magnetization data and a method for treatment of the weaknessDetection of the weakness of the vector averaging of magnetization data and a method for treatment of the weakness2832987098310.22059/jesphys.2019.263079.1007028FAMohammad EhsanHekmatianAssistant Professor, Materials and Nuclear Fuel Research School, Nuclear Science and Technology Research Institute (NSTRI), Atomic Energy Organization of Iran (AEOI), Tehran, IranJournal Article20180820In statistical common population, common or normal distribution is often governed and so that using Gaussian or normal probability density function and arithmetic averaging is appropriate. But if the statistical population has been formed from a number of spatial arbitrary directions, then common or normal distribution is not governed. In this condition Fisher probability density function and vector averaging can be used (Fisher is the name of the scientist who proposed the mentioned density function for the first time). In this function, each direction is shown as a point on a sphere with unit radius. The mentioned function shows the probability of having a particular direction in unit angular area of a particular area that has a definite central direction. This central direction shows the angular difference with the real average direction. In Fisher function, the distribution of the azimuth angles around the real vector average direction is symmetrical. The azimuth and the declination angles are the same and being symmetrical around the their distribution of the real average direction is logical. One of the statistical directional populations is the statistical population of different directions of the magnetization of rocks (Each magnetization direction is specified by two angles. First the angle between the magnetization direction and the surface of the horizon (inclination angle) and second the angle between the magnetization direction projection on the surface of the horizon and the geographic north direction (declination angle)). In this paper after an introduction, both normal and Fisher distributions (the latter is used for directional population) are discussed for better understanding of the difference between normal and directional statistical populations. Then the algorithm for calculating the vector averaging is presented. After that a software having vector averaging ability that is produced in this research is presented and then the vector and arithmetic averages are compared for magnetization data. During this research, it is clear that there is a weakness in the vector averaging and that weakness is that in some conditions the result of the vector averaging is not unique (this non uniqueness is because of the functions used in vector averaging algorithm). For example for calculating the declination angle, the function arc-tangent is used and we know that the result of this function is not unique. For example arctan (0.5637) is equal to both 29.41 and -150.9 degrees). The proposed method for the treatment of this weakness in this research is that, it would be proper to perform an arithmetic averaging beside the vector averaging and by which in the cases of having non unique results for vector averaging, the true result can be detectable (The result of the arithmetic averaging is unique) Between different results of the vector averaging, that result is true which is more similar to the arithmetic averaging. For example if there is a directional population which their declination angles are between -170 to -140 degrees and their arithmetic average is -150.67 degrees and the results of their vector averaging are 29.41 and -150.59 degrees, then the correct vector average is -150.59.In statistical common population, common or normal distribution is often governed and so that using Gaussian or normal probability density function and arithmetic averaging is appropriate. But if the statistical population has been formed from a number of spatial arbitrary directions, then common or normal distribution is not governed. In this condition Fisher probability density function and vector averaging can be used (Fisher is the name of the scientist who proposed the mentioned density function for the first time). In this function, each direction is shown as a point on a sphere with unit radius. The mentioned function shows the probability of having a particular direction in unit angular area of a particular area that has a definite central direction. This central direction shows the angular difference with the real average direction. In Fisher function, the distribution of the azimuth angles around the real vector average direction is symmetrical. The azimuth and the declination angles are the same and being symmetrical around the their distribution of the real average direction is logical. One of the statistical directional populations is the statistical population of different directions of the magnetization of rocks (Each magnetization direction is specified by two angles. First the angle between the magnetization direction and the surface of the horizon (inclination angle) and second the angle between the magnetization direction projection on the surface of the horizon and the geographic north direction (declination angle)). In this paper after an introduction, both normal and Fisher distributions (the latter is used for directional population) are discussed for better understanding of the difference between normal and directional statistical populations. Then the algorithm for calculating the vector averaging is presented. After that a software having vector averaging ability that is produced in this research is presented and then the vector and arithmetic averages are compared for magnetization data. During this research, it is clear that there is a weakness in the vector averaging and that weakness is that in some conditions the result of the vector averaging is not unique (this non uniqueness is because of the functions used in vector averaging algorithm). For example for calculating the declination angle, the function arc-tangent is used and we know that the result of this function is not unique. For example arctan (0.5637) is equal to both 29.41 and -150.9 degrees). The proposed method for the treatment of this weakness in this research is that, it would be proper to perform an arithmetic averaging beside the vector averaging and by which in the cases of having non unique results for vector averaging, the true result can be detectable (The result of the arithmetic averaging is unique) Between different results of the vector averaging, that result is true which is more similar to the arithmetic averaging. For example if there is a directional population which their declination angles are between -170 to -140 degrees and their arithmetic average is -150.67 degrees and the results of their vector averaging are 29.41 and -150.59 degrees, then the correct vector average is -150.59.Institute of Geophysics, University of TehranJournal of the Earth and Space Physics2538-371X45220190723OSL dating of sediments samples from Karun river traces in Khuzestan, SW IranOSL dating of sediments samples from Karun river traces in Khuzestan, SW Iran2993117099010.22059/jesphys.2019.267349.1007049FAMortezaFattahiAssociate Professor, Department of Earth Physics, Institute of Geophysics, University of Tehran, Tehran, Iran0000-0003-2719-5859Kevin P.WoodbridgePh.D. Student, Energy and Environment Institute, Faculty of Science and Engineering, University of Hull, Hull, UKMark D.BatemanProfessor, Faculty of Geography, University of Sheffield, Sheffield, UKJournal Article20181017The River Karun in lowland Khuzestan, SW Iran is influenced by various factors including tectonics, human activities, climate, and relative sea-level changes. Therefore, it is necessary to study these features from different aspects such as geology, geomorphology, paleoclimatology and Archeology. Disentangling these influences can be improved by investigating where river channels incise across active folds to produce river terraces. Determining the age of river terrace deposits has a fundamental role in these studies; especially since average rates of river incision since the time of terrace deposition can be a guide to average rates of tectonic uplift, particularly over longer timescales of thousands or tens of thousands of years where the influences of changes in aggradation and incision due to changes in sediment supply tend to be evened out (Bull, 1991; Burbank and Anderson 2012). <br />River terraces of the Karun river system were found associated with active folds in the Upper Khuzestan Plains. These folds were mostly asymmetric detachment folds and fault bend folds trending approximately NW-SE, with a more steeply dipping fore-limb to the south-west and a more gently dipping back-limb to the north-east (Blanc et al., 2003). <br />Woodbridge (2013) described these river terraces, and assigned each terrace a new name (from a nearby village or fold). As shown in Figure 1, four river terraces were associated with the Naft-e Safid Anticline: the 'Dar Khazineh terrace', the 'Batvand terrace', the 'Naft-e Safid terrace' and the 'Abgah terrace', on the fold fore-limb and back-limb. One river terrace was associated with the Sardarabad Anticline: the 'Kabutarkhan-e Sufla terrace', and one river terrace was associated with the Shushtar Anticline: the 'Kushkak terrace'; both on the fold back-limb Sediment samples were collected from the river terrace deposits and subjected to Optically Stimulated Luminescence (OSL) dating (Woodbridge and Frostick, 2014; Woodbridge et al., 2016). OSL dating was performed in the luminescence laboratory at the University of Sheffield, U.K. Both the palaeodose (De) and the dose rate was determined to derive an OSL age. <br />For De Determination the procedure outlined in Bateman and Catt (1996) was employed. The single aliquot regenerative (SAR) approach (Murray and Wintle 2000), was used for De determination. <br />All the samples showed a weak naturally OSL decay curves. Many aliquots failed to show good growth curves. All aliquots where the recycling ratio exceeded 10% of unity were excluded from further analysis. Thus, only around 10-20 percent of measured aliquots for each sample passed the criteria of the SAR protocol and their De are reported. The most appropriate preheat temperature for each sample was selected using a dose recovery preheat plateau test. This resulted in selection of preheat temperatures of 220 °C for 10 seconds and cutheat of 200 °C for 10 seconds, which were applied to each sample prior to OSL measurement to remove unstable signal generated by laboratory irradiation. <br />Analyst software was used for De determination. All samples demonstrated a high amount of replicate scatter with a large range of De values. Some of the distribution shape may reflect the limited population size of replicates but it also may reflect incomplete bleaching. Typically, poorly bleached sediments retain a significant level of residual signal from previous phases of sedimentary cycling, leading to inherent inaccuracies in the calculation of a palaeodose value. This is difficult to establish with any certainty from OSL data and should be taken in consideration with the site stratigraphy. In principle a well bleached unpost-depositionally disturbed sample should have replicate palaeodose (De) data which is normally distributed (See Bateman et al. 2003, Fig 3). By plotting the replicate data for each sample as a probability density function, some assessment of where older or younger material has been included in the sample measurements can be made. However, by determining the De of aliquots that contains 1000-2000 grains any heterogeneity in De that individual grains have may still be masked. We tried to overcome this problem by using smaller aliquots or at the single grain level. However, for these particular samples the weak OSL signal and low sensitivity to laboratory dose prevented such analysis. <br />In order to calculate an age, different models can be used. Woodbridge and Frostick (2014) and Woodbridge et al. (2016) published the age for each sample based on the mean De value determined by Finite Mixture Modelling or the Central Age Model. This paper provide the technical information behind dating these samples and provide all ages based on different models and without any judgement about partial bleaching, bioturbation or cryoturbation. Ages are quoted in years from the present day (2010) and are presented with one sigma confidence intervals which incorporate systematic uncertainties with the dosimetry data, uncertainties with the palaeomoisture content and errors associated with the De determination.The River Karun in lowland Khuzestan, SW Iran is influenced by various factors including tectonics, human activities, climate, and relative sea-level changes. Therefore, it is necessary to study these features from different aspects such as geology, geomorphology, paleoclimatology and Archeology. Disentangling these influences can be improved by investigating where river channels incise across active folds to produce river terraces. Determining the age of river terrace deposits has a fundamental role in these studies; especially since average rates of river incision since the time of terrace deposition can be a guide to average rates of tectonic uplift, particularly over longer timescales of thousands or tens of thousands of years where the influences of changes in aggradation and incision due to changes in sediment supply tend to be evened out (Bull, 1991; Burbank and Anderson 2012). <br />River terraces of the Karun river system were found associated with active folds in the Upper Khuzestan Plains. These folds were mostly asymmetric detachment folds and fault bend folds trending approximately NW-SE, with a more steeply dipping fore-limb to the south-west and a more gently dipping back-limb to the north-east (Blanc et al., 2003). <br />Woodbridge (2013) described these river terraces, and assigned each terrace a new name (from a nearby village or fold). As shown in Figure 1, four river terraces were associated with the Naft-e Safid Anticline: the 'Dar Khazineh terrace', the 'Batvand terrace', the 'Naft-e Safid terrace' and the 'Abgah terrace', on the fold fore-limb and back-limb. One river terrace was associated with the Sardarabad Anticline: the 'Kabutarkhan-e Sufla terrace', and one river terrace was associated with the Shushtar Anticline: the 'Kushkak terrace'; both on the fold back-limb Sediment samples were collected from the river terrace deposits and subjected to Optically Stimulated Luminescence (OSL) dating (Woodbridge and Frostick, 2014; Woodbridge et al., 2016). OSL dating was performed in the luminescence laboratory at the University of Sheffield, U.K. Both the palaeodose (De) and the dose rate was determined to derive an OSL age. <br />For De Determination the procedure outlined in Bateman and Catt (1996) was employed. The single aliquot regenerative (SAR) approach (Murray and Wintle 2000), was used for De determination. <br />All the samples showed a weak naturally OSL decay curves. Many aliquots failed to show good growth curves. All aliquots where the recycling ratio exceeded 10% of unity were excluded from further analysis. Thus, only around 10-20 percent of measured aliquots for each sample passed the criteria of the SAR protocol and their De are reported. The most appropriate preheat temperature for each sample was selected using a dose recovery preheat plateau test. This resulted in selection of preheat temperatures of 220 °C for 10 seconds and cutheat of 200 °C for 10 seconds, which were applied to each sample prior to OSL measurement to remove unstable signal generated by laboratory irradiation. <br />Analyst software was used for De determination. All samples demonstrated a high amount of replicate scatter with a large range of De values. Some of the distribution shape may reflect the limited population size of replicates but it also may reflect incomplete bleaching. Typically, poorly bleached sediments retain a significant level of residual signal from previous phases of sedimentary cycling, leading to inherent inaccuracies in the calculation of a palaeodose value. This is difficult to establish with any certainty from OSL data and should be taken in consideration with the site stratigraphy. In principle a well bleached unpost-depositionally disturbed sample should have replicate palaeodose (De) data which is normally distributed (See Bateman et al. 2003, Fig 3). By plotting the replicate data for each sample as a probability density function, some assessment of where older or younger material has been included in the sample measurements can be made. However, by determining the De of aliquots that contains 1000-2000 grains any heterogeneity in De that individual grains have may still be masked. We tried to overcome this problem by using smaller aliquots or at the single grain level. However, for these particular samples the weak OSL signal and low sensitivity to laboratory dose prevented such analysis. <br />In order to calculate an age, different models can be used. Woodbridge and Frostick (2014) and Woodbridge et al. (2016) published the age for each sample based on the mean De value determined by Finite Mixture Modelling or the Central Age Model. This paper provide the technical information behind dating these samples and provide all ages based on different models and without any judgement about partial bleaching, bioturbation or cryoturbation. Ages are quoted in years from the present day (2010) and are presented with one sigma confidence intervals which incorporate systematic uncertainties with the dosimetry data, uncertainties with the palaeomoisture content and errors associated with the De determination.Institute of Geophysics, University of TehranJournal of the Earth and Space Physics2538-371X45220190723Static shift estimation in magnetotelluric data using horizontal magnetic tensorStatic shift estimation in magnetotelluric data using horizontal magnetic tensor3133247098610.22059/jesphys.2019.268558.1007055FAEhsanLimooparvar JahromiM.Sc. Graduated, Department of Earth Physics, Institute of Geophysics, University of Tehran, Tehran, IranBanafshehHabibian DehkordiAssistant Professor, Department of Earth Physics, Institute of Geophysics, University of Tehran, Tehran, Iran0000-0002-6847-4291BehroozOskooiAssociate Professor, Department of Earth Physics, Institute of Geophysics, University of Tehran, Tehran, Iran0000-0003-3065-194XJournal Article20181110Interpretation of magnetotelluric data in the presence of galvanic distortions, caused by small-scale near-surface structures, can lead to unreliable results. The simplest manifest of these distortions, while limited only to the changes in the amplitude of electric fields, is vertical displacement of the apparent resistivity curves or static shift phenomenon that constitutes non-determinable part of the distortion matrix. Due to the boundary conditions governing components of the electric and magnetic fields, the occurrence of charge accumulation and therefore the static shift of apparent resistivity curves, affects only TM-mode data in the case of two-dimensional models. Thus, we can use the information available in the TE-mode- impedance phase (tipper and horizontal magnetic data) which are independent of this phenomenon. <br />In this study, geomagnetic transfer functions have been used to estimate this displacement and recover the undistorted TE-mode apparent resistivity based on the Faraday induction law. Ledo et al. (2002) show that tipper data can be used to estimate static shift of magnetotelluric data, if the horizontal variations of the horizontal components of the magnetic field can be ignored. This assumption may be violated in complex situations. We estimate static shift while incorporating such variations and taking into account the horizontal magnetic transfer functions. Estimation of static shift through mathematical methods is only relatively possible and requires the selection of a reference station that has the minimal effect of galvanic distortion. The relations between different components of electric and magnetic fields are integrated and characterized by their mean values. To incorporate the horizontal magnetic tensor, array magnetotelluric data are required, so that components of the magnetic field at the reference and measurement sites are simultaneously provided. <br />Considering two consecutive sites, impedance tensor at one site is written in terms of tipper and horizontal magnetic tensor at that site and the impedance at the adjacent site. By ignoring other types of distortions that can generally exist and using some algebra, the problem of determining the frequency-independent static shift factor becomes a linear fit problem. A set of data points covering different frequency ranges is selected and the quality of their linear fitting is examined through soling procedure. <br />Considering the horizontal variations of the horizontal components of the magnetic field, the method has been applied to two synthetic models. Using two different approaches, the distortions caused by small-scale three-dimensional structures are simulated and added to the model responses. In the first approach, the distortion matrix is considered as the product of four parameters of gain factor, anisotropy, and twist and shear angles in the decomposition model then the distortion simulation is performed by selecting some numerical values of these four parameters and multiplying the resulted distortion matrix by the impedance tensor. In the second approach, some part of the top-layer of the model is replaced by a Gaussian distribution of resistivity with known selected mean and standard deviation. In this way, the effects of various geological processes, such as weathering, erosion- and to some extent- the nature of the deposited sediments are involved. The obtained results confirm that the estimated static shift parameter is more accurate than that of the case in which horizontal magnetic transfer functions are ignored and only the vertical magnetic transfer functions are considered.Interpretation of magnetotelluric data in the presence of galvanic distortions, caused by small-scale near-surface structures, can lead to unreliable results. The simplest manifest of these distortions, while limited only to the changes in the amplitude of electric fields, is vertical displacement of the apparent resistivity curves or static shift phenomenon that constitutes non-determinable part of the distortion matrix. Due to the boundary conditions governing components of the electric and magnetic fields, the occurrence of charge accumulation and therefore the static shift of apparent resistivity curves, affects only TM-mode data in the case of two-dimensional models. Thus, we can use the information available in the TE-mode- impedance phase (tipper and horizontal magnetic data) which are independent of this phenomenon. <br />In this study, geomagnetic transfer functions have been used to estimate this displacement and recover the undistorted TE-mode apparent resistivity based on the Faraday induction law. Ledo et al. (2002) show that tipper data can be used to estimate static shift of magnetotelluric data, if the horizontal variations of the horizontal components of the magnetic field can be ignored. This assumption may be violated in complex situations. We estimate static shift while incorporating such variations and taking into account the horizontal magnetic transfer functions. Estimation of static shift through mathematical methods is only relatively possible and requires the selection of a reference station that has the minimal effect of galvanic distortion. The relations between different components of electric and magnetic fields are integrated and characterized by their mean values. To incorporate the horizontal magnetic tensor, array magnetotelluric data are required, so that components of the magnetic field at the reference and measurement sites are simultaneously provided. <br />Considering two consecutive sites, impedance tensor at one site is written in terms of tipper and horizontal magnetic tensor at that site and the impedance at the adjacent site. By ignoring other types of distortions that can generally exist and using some algebra, the problem of determining the frequency-independent static shift factor becomes a linear fit problem. A set of data points covering different frequency ranges is selected and the quality of their linear fitting is examined through soling procedure. <br />Considering the horizontal variations of the horizontal components of the magnetic field, the method has been applied to two synthetic models. Using two different approaches, the distortions caused by small-scale three-dimensional structures are simulated and added to the model responses. In the first approach, the distortion matrix is considered as the product of four parameters of gain factor, anisotropy, and twist and shear angles in the decomposition model then the distortion simulation is performed by selecting some numerical values of these four parameters and multiplying the resulted distortion matrix by the impedance tensor. In the second approach, some part of the top-layer of the model is replaced by a Gaussian distribution of resistivity with known selected mean and standard deviation. In this way, the effects of various geological processes, such as weathering, erosion- and to some extent- the nature of the deposited sediments are involved. The obtained results confirm that the estimated static shift parameter is more accurate than that of the case in which horizontal magnetic transfer functions are ignored and only the vertical magnetic transfer functions are considered.Institute of Geophysics, University of TehranJournal of the Earth and Space Physics2538-371X45220190723Optimized Estimation of Observation Precisions In Classical Displacement NetworkOptimized Estimation of Observation Precisions In Classical Displacement Network3253427098210.22059/jesphys.2019.268781.1007058FASaeedFarzanehAssistant Professor, Department of Surveying and Geomatics Engineering, Faculty of Engineering, University of Tehran, Tehran, IranKamalParvaziPh.D. Student Department of Surveying and Geomatics Engineering, Faculty of Engineering,
University of Tehran, Tehran, IranJournal Article20181112Any infrastructure such as dams need constant monitoring for the detection of risks of failure and/or to plan civil engineering maintaining work. A recent approach considers precise geodetic instruments and satellite-based geodetic monitoring as a method to estimate potential deformation of such structures. A growing need for a fully automated and continuous monitoring of structural and ground deformations has created new challenges for design and analysis of the monitoring schemes, where multi-sensor geodetic systems can provide essential aid. Combination of different geodetic data helps determining displacements with high precision, hence, the risk of damages is reduced. Corresponding authorities of large man-made structures are faced with the safety problem, as all have aim to reduce risk and cost. Designers try to design large structures to tolerate against different forces like wind, traffic load, temperature, flood, earthquake, land uplift etc. Using geodetic instruments and techniques, we are able to monitor the deformation behavior or deflection in the mentioned structures and eventually provide a structural failure alarm capability (Andersson 2008). <br />It is important to select appropriate sensor and methods to detect the deformation. Slow deforming dams require sub-millimeter to millimeter level accuracy to monitor the displacement and deformation (Lindenbergh et al. 2005). Reaching this level of accuracy is not costly, if geodetic sensors are integrated with other sensors (e.g. geotechnical sensors, and precise total stations, see Hwang et al. 2012). It might be to implement other sensors (e.g. laser scanner and Total Station). Using point clouds data for deformation monitoring is almost new. Gonzalez et al. (2012) studied on point clouds accuracy for applications in civil engineering e.g. deformation monitoring. They showed that the results appear suitable for deformation monitoring, with accuracies less than 1 mm. Bagherbandi et al. (2009) studied on various techniques to find the optimal design of a deformation network using various criteria such as precision, cost and reliability. Better results can be achieved using the control network, provided that an optimal network design is performed for detecting deformations (Kuang 1996). In addition, the methods of geodetic network process can affect the results (Bagherbandi 2016). <br />The aim of this study is primarily to evaluate different deformation monitoring methods and possibilities to physically interpret the deformation and evaluate the risk of failures. In this research, the idea of assigning weights for the observations by least square variance components estimation (LS-VCE) is used (Amiri-Simkooei 2007; Teunissen and Amiri-Simkooei 2008) in order to improve accuracy of adjustment results, which differs from the applied method in Bagherbandi (2016) to determine the variance components. Some issues and parameters should be investigated in LS-VCE such as the effect of variance components estimation on the observations final accuracy, the absolute error ellipsoid estimation, the study of the necessary conditions in a network to achieve higher accuracy and its effect on obtaining real results from the reliability matrix. All results obtained from adjustment by element, LS-VCE, and Tikhonov regularization are compared using a simulated geodetic network and real data. Results from this study provide important information in studying deformation that can be used to interpret the deformation mechanism, which may reduce the risk of potential disasters in large structures. We will evaluate the above-mentioned methods in Jamishan dam in Iran and utilize the geodetic techniques and observations to monitor the deformation of the dam.Any infrastructure such as dams need constant monitoring for the detection of risks of failure and/or to plan civil engineering maintaining work. A recent approach considers precise geodetic instruments and satellite-based geodetic monitoring as a method to estimate potential deformation of such structures. A growing need for a fully automated and continuous monitoring of structural and ground deformations has created new challenges for design and analysis of the monitoring schemes, where multi-sensor geodetic systems can provide essential aid. Combination of different geodetic data helps determining displacements with high precision, hence, the risk of damages is reduced. Corresponding authorities of large man-made structures are faced with the safety problem, as all have aim to reduce risk and cost. Designers try to design large structures to tolerate against different forces like wind, traffic load, temperature, flood, earthquake, land uplift etc. Using geodetic instruments and techniques, we are able to monitor the deformation behavior or deflection in the mentioned structures and eventually provide a structural failure alarm capability (Andersson 2008). <br />It is important to select appropriate sensor and methods to detect the deformation. Slow deforming dams require sub-millimeter to millimeter level accuracy to monitor the displacement and deformation (Lindenbergh et al. 2005). Reaching this level of accuracy is not costly, if geodetic sensors are integrated with other sensors (e.g. geotechnical sensors, and precise total stations, see Hwang et al. 2012). It might be to implement other sensors (e.g. laser scanner and Total Station). Using point clouds data for deformation monitoring is almost new. Gonzalez et al. (2012) studied on point clouds accuracy for applications in civil engineering e.g. deformation monitoring. They showed that the results appear suitable for deformation monitoring, with accuracies less than 1 mm. Bagherbandi et al. (2009) studied on various techniques to find the optimal design of a deformation network using various criteria such as precision, cost and reliability. Better results can be achieved using the control network, provided that an optimal network design is performed for detecting deformations (Kuang 1996). In addition, the methods of geodetic network process can affect the results (Bagherbandi 2016). <br />The aim of this study is primarily to evaluate different deformation monitoring methods and possibilities to physically interpret the deformation and evaluate the risk of failures. In this research, the idea of assigning weights for the observations by least square variance components estimation (LS-VCE) is used (Amiri-Simkooei 2007; Teunissen and Amiri-Simkooei 2008) in order to improve accuracy of adjustment results, which differs from the applied method in Bagherbandi (2016) to determine the variance components. Some issues and parameters should be investigated in LS-VCE such as the effect of variance components estimation on the observations final accuracy, the absolute error ellipsoid estimation, the study of the necessary conditions in a network to achieve higher accuracy and its effect on obtaining real results from the reliability matrix. All results obtained from adjustment by element, LS-VCE, and Tikhonov regularization are compared using a simulated geodetic network and real data. Results from this study provide important information in studying deformation that can be used to interpret the deformation mechanism, which may reduce the risk of potential disasters in large structures. We will evaluate the above-mentioned methods in Jamishan dam in Iran and utilize the geodetic techniques and observations to monitor the deformation of the dam.Institute of Geophysics, University of TehranJournal of the Earth and Space Physics2538-371X45220190723Monitoring and predicting the trend of changes in residential areas using multi-timed images (Case study: Songhor city)Monitoring and predicting the trend of changes in residential areas using multi-timed images (Case study: Songhor city)3433547096910.22059/jesphys.2019.275076.1007084FASaeedNegahbanAssistant Professor, Department of Geography, Faculty of Economics, Management & Social sciences, Shiraz University, Shiraz, Iran0000-0003-0701-5638HamidGanjaeianPh.D. Student, Department of Physical Geography, Faculty of Geography, University of Tehran, Tehran, IranAtrinEbrahimiPh.D. Student, Department of Physical Geography, Faculty of Geography and planning, Tabriz University, Tabriz, IranKamyarEmamiM.Sc. Graduated, Department of Physical Geography, Faculty of Geography, University of Tehran, Tehran, IranJournal Article20190206Considering the ever-increasing changes in land uses and the need for managers and experts to know how changes have taken place in policy and options for solving the existing problems. Detection of changes to determine the trend over time seems necessary. On the other hand, modeling future changes is important for understanding the quality of future changes. Therefore, the full recognition of land use, its past changes and the prediction of future changes plays an important role in the sustainable management of resources. Modeling land use processes is an important tool in optimizing land use and land use planning. One of the models used to predict landslide changes is the model of artificial neural networks and Markov chain analysis. The features of the artificial neural network include the ability to learn and generalize and process information in parallel. <br />Considering the goal of urban development during the years 2000 to 2012, satellite imagery of the years 2000 and 2012 in June has been used. After the preparation of satellite imagery and pre-processing of images, the landuse in the study area for the years 2000 and 2012 has been prepared. Then useing the LCM model landuse change patterns of changes were analyzed. Then, based on the Markov chain model, the potential for changing each use to residential use is measured. This means that each pixel was capable of showing change the image from one land use to another. Then, based on the major changes in the region in the survey, three sub-models of shifting change were identified as transforming pastures into habitat areas, converting agricultural production into settlements, and transforming dryland farming into settlements. After calculating the potential for the transfer of any land use to a settlement using descriptive data, a plan for predicting the use of land for 2025 and 2040 was then provided. <br />Given that the purpose of the present study was to assess the development of residential areas, the extent of changes in these areas were assessed during the years 2000 to 2012. The results indicate that the residential areas increased from 8.3 square kilometers in 2000 to 12.6 square kilometers in 2012, according to the land use map, and mostly changes in the urban area of Songhor area have been made. The results of the assessment of changes indicate that the land use change from irrigated agricultural to residential use during the 12 year period was 1.9 km<sup>2</sup>, which for dryland agriculture it was 0.6 kilometers, Also 1.8 km<sup>2</sup> of rangelands has become residential. The results of this study indicate that the irrigated agricultural lands of the city of Sangar, especially the southern regions and pastures near the urban area, have had most changes. Among the changes in other uses, about 11.5 km<sup>2</sup> of the rangeland has been converted into rainfed farming, and about 12.3 km<sup>2</sup> of land has also become rangelands and also, about 4.7 km<sup>2</sup> of irrigated agricultural has become arable land or Bayer land and about 1.5 km<sup>2</sup> of rangelands has become irrigated agricultural land. <br />The growing population has led to an increase in the number of habitat areas and, as a result, agricultural lands and pastures have undergone changes. The growing trend of settlement development varies from region to region, and in the urban area of Songhor more are moving toward the southern regions of the urban area. Considering the geomorphologic status of the study area, a large part of the range is covered by rangelands. Irrigated agricultural lands which have a significant share, are located on the outskirts of the city of Songhor, which are undergoing further changes. According to the main objective of the research, based on descriptive data such as distance from communication, distance from urban boundaries, elevation and slope, the amount of development of residential areas for 2025 and 2040 is also projected. The results of the forecast indicate that in the case of the growing trend, the development of the settlements will reach about 18.2 km<sup>2</sup> in 2025, and will reach 24.2 km<sup>2</sup> in 2040, due to the high potential of the southern regions of the city of Songhor, the highest rate of development of settlements will be towards these areas. The results indicate that the increasing number of settlements in the city of Songhor will lead to the degradation of high-quality agricultural lands and pastures. If the trend is continued, the irrigated agricultural around of the city of Songhor will reach the lowest level by 2040. Also most of the pastures will also be degraded. Hence, it is necessary to identify areas suitable for the development of a settlement before increasing of rate the destruction occur, so that less prone areas for agriculture and pastures can be degraded.Considering the ever-increasing changes in land uses and the need for managers and experts to know how changes have taken place in policy and options for solving the existing problems. Detection of changes to determine the trend over time seems necessary. On the other hand, modeling future changes is important for understanding the quality of future changes. Therefore, the full recognition of land use, its past changes and the prediction of future changes plays an important role in the sustainable management of resources. Modeling land use processes is an important tool in optimizing land use and land use planning. One of the models used to predict landslide changes is the model of artificial neural networks and Markov chain analysis. The features of the artificial neural network include the ability to learn and generalize and process information in parallel. <br />Considering the goal of urban development during the years 2000 to 2012, satellite imagery of the years 2000 and 2012 in June has been used. After the preparation of satellite imagery and pre-processing of images, the landuse in the study area for the years 2000 and 2012 has been prepared. Then useing the LCM model landuse change patterns of changes were analyzed. Then, based on the Markov chain model, the potential for changing each use to residential use is measured. This means that each pixel was capable of showing change the image from one land use to another. Then, based on the major changes in the region in the survey, three sub-models of shifting change were identified as transforming pastures into habitat areas, converting agricultural production into settlements, and transforming dryland farming into settlements. After calculating the potential for the transfer of any land use to a settlement using descriptive data, a plan for predicting the use of land for 2025 and 2040 was then provided. <br />Given that the purpose of the present study was to assess the development of residential areas, the extent of changes in these areas were assessed during the years 2000 to 2012. The results indicate that the residential areas increased from 8.3 square kilometers in 2000 to 12.6 square kilometers in 2012, according to the land use map, and mostly changes in the urban area of Songhor area have been made. The results of the assessment of changes indicate that the land use change from irrigated agricultural to residential use during the 12 year period was 1.9 km<sup>2</sup>, which for dryland agriculture it was 0.6 kilometers, Also 1.8 km<sup>2</sup> of rangelands has become residential. The results of this study indicate that the irrigated agricultural lands of the city of Sangar, especially the southern regions and pastures near the urban area, have had most changes. Among the changes in other uses, about 11.5 km<sup>2</sup> of the rangeland has been converted into rainfed farming, and about 12.3 km<sup>2</sup> of land has also become rangelands and also, about 4.7 km<sup>2</sup> of irrigated agricultural has become arable land or Bayer land and about 1.5 km<sup>2</sup> of rangelands has become irrigated agricultural land. <br />The growing population has led to an increase in the number of habitat areas and, as a result, agricultural lands and pastures have undergone changes. The growing trend of settlement development varies from region to region, and in the urban area of Songhor more are moving toward the southern regions of the urban area. Considering the geomorphologic status of the study area, a large part of the range is covered by rangelands. Irrigated agricultural lands which have a significant share, are located on the outskirts of the city of Songhor, which are undergoing further changes. According to the main objective of the research, based on descriptive data such as distance from communication, distance from urban boundaries, elevation and slope, the amount of development of residential areas for 2025 and 2040 is also projected. The results of the forecast indicate that in the case of the growing trend, the development of the settlements will reach about 18.2 km<sup>2</sup> in 2025, and will reach 24.2 km<sup>2</sup> in 2040, due to the high potential of the southern regions of the city of Songhor, the highest rate of development of settlements will be towards these areas. The results indicate that the increasing number of settlements in the city of Songhor will lead to the degradation of high-quality agricultural lands and pastures. If the trend is continued, the irrigated agricultural around of the city of Songhor will reach the lowest level by 2040. Also most of the pastures will also be degraded. Hence, it is necessary to identify areas suitable for the development of a settlement before increasing of rate the destruction occur, so that less prone areas for agriculture and pastures can be degraded.Institute of Geophysics, University of TehranJournal of the Earth and Space Physics2538-371X45220190723Snow Cover Changes of Zagros Range in 2001-2016 Using Daily Data of MODISSnow Cover Changes of Zagros Range in 2001-2016 Using Daily Data of MODIS3553717096610.22059/jesphys.2019.256133.1006997FAPaymanMohammadi AhmadmahmoudiM.Sc. Student, Department of Geographical science, University of Hormozgan, Bandaabbas, IranAsadollahKhooraniAssociate Professor, Department of Geographical science, University of Hormozgan, Bandaabbas, Iran0000-0003-0311-9624Journal Article20180508Snow cover (SC) is the main part of global cryosphere that covers about 46 million square kilometer of earth surface. Seasonal SC and snow storage provides the required water of around 1/6 of world population. About 98 percent of global SC is in the northern hemisphere. Southern and southwestern mountain range of Iran, called Zagros is located in the mid-latitudes that is exposed to the northern hemisphere westerlies resulting in high amount of rain and snowfall. This region contains 40 percent of Iran's water storage. Monitoring and evaluation of seasonal SC is essential for managing winter tourism, transportation, environment, land cover, water and irrigation and also essential in runoff modeling. In our study area snow gauges are insufficient and sparse, therefore monitoring spatial changes of Snow Cover Area (SCA) is difficult. MODIS images lunched on TERRA and AQUA platforms are suitable data source for estimating SC, because they provide daily images with 500 meter spatial resolution. The aim of this study is detecte and analyze SC variations and trends in Zagros mountain range using daily images of MODIS in order to determine SC variation. <br />In this research, 348 daily products of MODIS for first days of snowy months (Months with some snow on the ground in Zagros; Dec, Jan, Feb, Mar, Apr, May, Jun) are obtained from National Snow & Ice Data Center (NSIDC) to calculate SCA. Digital Elevation Model (DEM) is used to determining snow line (SRTM 90m from Earthexplorer.usgs.gov) and total monthly precipitation of weather stations above 1500 meter elevation of the same period with SC data is used. Snow line is extracted by overlaying DEM on SCA. Then the snow line pattern is determined as the average of snow line for each day during 2001 to 2016. SCA pattern is also extracted as the average of each day during the same period. In order to calculate trends in SCA and Snow Line Man-Kendal test is used. Man-Kendal test is less sensitive too extreme values of time series. <br />The trends of SC are investigated using Man-Kendal test results show a decreasing trend in January, May and June in 2007-2016. The highest amount of SCA decrease and snow line retrogression for January, May and June are recorded in 1250-4500 m, 2300-4500 m and 3500-4500 m altitude belts, respectively. Based on SC departure from the average, two periods are recognizable: first, from 2001 to 2008, in this period SC is above the average and is more similar to the snow cover pattern of the whole period (2001-2016). Second period from 2009 to 2016, in this period SC is more below the average and is not similar to the snow cover pattern of the whole period (2001-2016). In the second period mean SC reduction compared to the first period in Dec, Jan, Feb, Mar, Apr, May and June are 4, 67.6, 24.5, 21, 10, 47.3 and 48.4 respectively. Mean monthly precipitation changes of weather stations above 1500 meter elevation is similar to the changes of SC changes.Snow cover (SC) is the main part of global cryosphere that covers about 46 million square kilometer of earth surface. Seasonal SC and snow storage provides the required water of around 1/6 of world population. About 98 percent of global SC is in the northern hemisphere. Southern and southwestern mountain range of Iran, called Zagros is located in the mid-latitudes that is exposed to the northern hemisphere westerlies resulting in high amount of rain and snowfall. This region contains 40 percent of Iran's water storage. Monitoring and evaluation of seasonal SC is essential for managing winter tourism, transportation, environment, land cover, water and irrigation and also essential in runoff modeling. In our study area snow gauges are insufficient and sparse, therefore monitoring spatial changes of Snow Cover Area (SCA) is difficult. MODIS images lunched on TERRA and AQUA platforms are suitable data source for estimating SC, because they provide daily images with 500 meter spatial resolution. The aim of this study is detecte and analyze SC variations and trends in Zagros mountain range using daily images of MODIS in order to determine SC variation. <br />In this research, 348 daily products of MODIS for first days of snowy months (Months with some snow on the ground in Zagros; Dec, Jan, Feb, Mar, Apr, May, Jun) are obtained from National Snow & Ice Data Center (NSIDC) to calculate SCA. Digital Elevation Model (DEM) is used to determining snow line (SRTM 90m from Earthexplorer.usgs.gov) and total monthly precipitation of weather stations above 1500 meter elevation of the same period with SC data is used. Snow line is extracted by overlaying DEM on SCA. Then the snow line pattern is determined as the average of snow line for each day during 2001 to 2016. SCA pattern is also extracted as the average of each day during the same period. In order to calculate trends in SCA and Snow Line Man-Kendal test is used. Man-Kendal test is less sensitive too extreme values of time series. <br />The trends of SC are investigated using Man-Kendal test results show a decreasing trend in January, May and June in 2007-2016. The highest amount of SCA decrease and snow line retrogression for January, May and June are recorded in 1250-4500 m, 2300-4500 m and 3500-4500 m altitude belts, respectively. Based on SC departure from the average, two periods are recognizable: first, from 2001 to 2008, in this period SC is above the average and is more similar to the snow cover pattern of the whole period (2001-2016). Second period from 2009 to 2016, in this period SC is more below the average and is not similar to the snow cover pattern of the whole period (2001-2016). In the second period mean SC reduction compared to the first period in Dec, Jan, Feb, Mar, Apr, May and June are 4, 67.6, 24.5, 21, 10, 47.3 and 48.4 respectively. Mean monthly precipitation changes of weather stations above 1500 meter elevation is similar to the changes of SC changes.Institute of Geophysics, University of TehranJournal of the Earth and Space Physics2538-371X45220190723Experimental study of the effect of hydrophilic particles on fog modificationExperimental study of the effect of hydrophilic particles on fog modification3733826915110.22059/jesphys.2019.262718.1007027FAAmirBagheri Mosleh-AbadiM.Sc. Student, Department of Space Physics, Institute of Geophysics, University of Tehran, Tehran, IranAbbas AliAliakbari-BidokhtiProfessor, Department of Space Physics, Institute of Geophysics, University of Tehran, Tehran, Iran0000-0003-4841-2218MaryamGharaylouAssistant Professor, Department of Space Physics, Institute of Geophysics, University of Tehran, Tehran, Iran0000-0002-0014-9409RezaKhalifeM.Sc. Student, Department of Statistics, Faculty of mathematics, statistics & computer Science, University of Tehran, Tehran, IranJournal Article20180818Fog formation reduces the visibility, and low visibility cause problems mainly in airports and roads. Fog modification or even clearance can reduce such problems. Here, in the laboratory the effects of some aerosols on the modification of artificial fog have been investigated. This research is carried out using the method used in the work of Vajaiac et al. (2014), for hydrophilic particles. These particles, as the fog modification factor, can be effective as condensation nuclei by increasing the effective radius of the droplets and reducing the concentrations of the fog droplets in the fog.. We used a glass chamber with a height of 41cm and diameter of 28 cm in which a red laser and a detector were placed on both sides of the chamber. We dissolved water-soluble compounds in water and by pumping the air into water of the chamber, the bubbles are formed and their breaking led to formation of particles. Then the particles are dispersed inside the chamber, while the pressure inside the chamber was increased with an air pump, and so the temperature was rised. Then suddenly or quasi-adiabaticaly, the pressure is released and hence, the air temperature was dropped creating the fog (cloud). The cloud droplets are mainly formed on these particles. With the formation of droplets, the received laser signal was decreased and after falling all of droplets, the received laser signal reached its primary level. Hence, by measuring the level of drop in the signal and its duration some properties to the cloud structure can be found. <br />In this study, we used sodium hydroxide, urea, and salt particles that was compared with background aerosols affecting the cloud. First, the optical depth of fog was calculated and then size distribution of fog droplets was estimated. Also other parameters like the concentration of the fog droplets, liquid water content and fog lifetime were evaluated. The results show that hydrophilic particles as condensation nuclei played a significant role in fog modification. So, in the presence of these particles, the concentration of fog droplets in comparison with that of the background aerosols, was decreased and the fog was diluted significantly. The maximum of fog lifetime was observed in the presence of background aerosols as 63 s, also these particles have had the lowest effective radius of 7.77µm. In the presence of sodium hydroxide particles, the concentration of total droplets was decreased significantly. In addition, the maximum effective radius was 10.48 µm for these particles. Also, fog lifetime was reduced to 31s which suggests that this component for fog modification has a better performance. The area under the curve’s of the size distribution of droplets for salt and urea were nearly identical, that indicated the concentration of the droplets was close to each other. Parasitic structures might appear in curves of the size distribution of droplets that were smoothed. As the time passed all of the effective radius curves showed decreasing trends due to the fact that, the larger droplets can droped out of the laser beam, so gradually the effective radius should be decreased. In comparison with other compounds, when sodium hydroxide particles were used in the cloud chamber, observationally small droplets that moved up an down in the laser beam weren’t present, hence no oscillations were seen in the laser signal. Considering experimental errors, while neglecting the process of coalescence that might occur in the falling stage of the droplets, when the droplets were passing the laser beam, they appeared at larger sizes, which could lead to an overestimation. The signal noises and non-uniform distribution of condensation nuclei could also be considered as experimental errors.Fog formation reduces the visibility, and low visibility cause problems mainly in airports and roads. Fog modification or even clearance can reduce such problems. Here, in the laboratory the effects of some aerosols on the modification of artificial fog have been investigated. This research is carried out using the method used in the work of Vajaiac et al. (2014), for hydrophilic particles. These particles, as the fog modification factor, can be effective as condensation nuclei by increasing the effective radius of the droplets and reducing the concentrations of the fog droplets in the fog.. We used a glass chamber with a height of 41cm and diameter of 28 cm in which a red laser and a detector were placed on both sides of the chamber. We dissolved water-soluble compounds in water and by pumping the air into water of the chamber, the bubbles are formed and their breaking led to formation of particles. Then the particles are dispersed inside the chamber, while the pressure inside the chamber was increased with an air pump, and so the temperature was rised. Then suddenly or quasi-adiabaticaly, the pressure is released and hence, the air temperature was dropped creating the fog (cloud). The cloud droplets are mainly formed on these particles. With the formation of droplets, the received laser signal was decreased and after falling all of droplets, the received laser signal reached its primary level. Hence, by measuring the level of drop in the signal and its duration some properties to the cloud structure can be found. <br />In this study, we used sodium hydroxide, urea, and salt particles that was compared with background aerosols affecting the cloud. First, the optical depth of fog was calculated and then size distribution of fog droplets was estimated. Also other parameters like the concentration of the fog droplets, liquid water content and fog lifetime were evaluated. The results show that hydrophilic particles as condensation nuclei played a significant role in fog modification. So, in the presence of these particles, the concentration of fog droplets in comparison with that of the background aerosols, was decreased and the fog was diluted significantly. The maximum of fog lifetime was observed in the presence of background aerosols as 63 s, also these particles have had the lowest effective radius of 7.77µm. In the presence of sodium hydroxide particles, the concentration of total droplets was decreased significantly. In addition, the maximum effective radius was 10.48 µm for these particles. Also, fog lifetime was reduced to 31s which suggests that this component for fog modification has a better performance. The area under the curve’s of the size distribution of droplets for salt and urea were nearly identical, that indicated the concentration of the droplets was close to each other. Parasitic structures might appear in curves of the size distribution of droplets that were smoothed. As the time passed all of the effective radius curves showed decreasing trends due to the fact that, the larger droplets can droped out of the laser beam, so gradually the effective radius should be decreased. In comparison with other compounds, when sodium hydroxide particles were used in the cloud chamber, observationally small droplets that moved up an down in the laser beam weren’t present, hence no oscillations were seen in the laser signal. Considering experimental errors, while neglecting the process of coalescence that might occur in the falling stage of the droplets, when the droplets were passing the laser beam, they appeared at larger sizes, which could lead to an overestimation. The signal noises and non-uniform distribution of condensation nuclei could also be considered as experimental errors.Institute of Geophysics, University of TehranJournal of the Earth and Space Physics2538-371X45220190723The role of interaction subtropical- polar jet stream on heavy rainfall in the spring and autumn of IranThe role of interaction subtropical- polar jet stream on heavy rainfall in the spring and autumn of Iran3834006915310.22059/jesphys.2019.264764.1007033FAAzarJalilianPh.D. Student, Department of Geography, Razi University, Kermanshah, IranHasanZolfaghariAssociate Professor, Department of Geography, Razi University, Kermanshah, IranJalilSahraeiAssistant Professor, Department of Atmospheric physics, Faculty of Science, Razi University, Kermanshah, IranJournal Article20180918The narrow bands of strong wind (≤30ms<sup>-1</sup>) in the upper levels of the atmosphere called jet stream. The formation of these streams in the maximum gradient zone of the atmosphere (the location of the maximum temperature and energy transition) has caused the establishment and movement of their speed nuclei at the western extremities of the equatorial fringes and on the polar front in the extraterrestrial region and the formation of western margin and polar fronts. This winds blow from west to east, but they often shifts to the north and south of seasonal average position. In this case, jets may be deployed simultaneously in one region and occurrence the interaction between them may affect the heavy rainfall. Based on several definitions that have been presented with heavy rainfall according to different criteria and thresholds, Cumulative precipitation of more than 30 mm in 24 hours, causing severe financial losses to various urban and rural areas, is called heavy rainfall. In This article, shifting and intensity of subtropical - polar front jet stream in the mid-latitude region and that’s effects on heavy precipitation in Iran has been studied. <br />In this study, the effects of interactions patterns subtropical-polar front jet stream on heavy spring and autumn rainfall in Iran (1988-2017) with Atmospheric circulation to environment Approach has been investigated. Therefore after determining the maximum coefficient of variation of this jet streams, in order to determine the effective circulation patterns in these synchronies, by applying factor analysis by PCA on the data of a height of 500 hp, four main components with The explanation of 82% of data diffraction was determined.Then, by using cluster analysis method, the days in each cluster and by intra-group correlation, the representative day of each cluster was determined.After drawing different patterns of simultaneous deployment in the Surfer software, four synchronous patterns were determined in the range of maximum coefficient of variation of jet stream (27-40°N, 36-56 ° E). Finally, the effects of each of these patterns on the occurrence of heavy rainfall were determined by obtaining divergence, vertical velocity, vorticity, relative humidity, specific humidity, temperature and mean sea pressure from the ECMWF center. <br />The results of the study show the effects of three interaction patterns on heavy rainfalls. The effects of direct and indirect jet streams interactions in the first and second patterns, with the change in the thickness of the polar fronts and the interaction of jets in the third model, has been affected on atmospheric conditions and occurrence of heavy rainfall. Despite the direct interactions of the first pattern of jet stream synchronization in the rainfall concentration, the indirect interactions effect of the second pattern on the cyclone circulation boost of the polar jet stream has led to the occurrence of extensive rainfall. Also, the mutual interactions of the jet stream in the third pattern, with the intensification of the frontogenesis process, have been effective in the occurrence of precipitation along these fronts, while the proximity of these opposite turns in the fourth jet streams concurrency pattern is not effective on rainfall, and the occurrence of heavy rainfall has been affected only from the cyclonic circulation of the subtropical jet stream.The narrow bands of strong wind (≤30ms<sup>-1</sup>) in the upper levels of the atmosphere called jet stream. The formation of these streams in the maximum gradient zone of the atmosphere (the location of the maximum temperature and energy transition) has caused the establishment and movement of their speed nuclei at the western extremities of the equatorial fringes and on the polar front in the extraterrestrial region and the formation of western margin and polar fronts. This winds blow from west to east, but they often shifts to the north and south of seasonal average position. In this case, jets may be deployed simultaneously in one region and occurrence the interaction between them may affect the heavy rainfall. Based on several definitions that have been presented with heavy rainfall according to different criteria and thresholds, Cumulative precipitation of more than 30 mm in 24 hours, causing severe financial losses to various urban and rural areas, is called heavy rainfall. In This article, shifting and intensity of subtropical - polar front jet stream in the mid-latitude region and that’s effects on heavy precipitation in Iran has been studied. <br />In this study, the effects of interactions patterns subtropical-polar front jet stream on heavy spring and autumn rainfall in Iran (1988-2017) with Atmospheric circulation to environment Approach has been investigated. Therefore after determining the maximum coefficient of variation of this jet streams, in order to determine the effective circulation patterns in these synchronies, by applying factor analysis by PCA on the data of a height of 500 hp, four main components with The explanation of 82% of data diffraction was determined.Then, by using cluster analysis method, the days in each cluster and by intra-group correlation, the representative day of each cluster was determined.After drawing different patterns of simultaneous deployment in the Surfer software, four synchronous patterns were determined in the range of maximum coefficient of variation of jet stream (27-40°N, 36-56 ° E). Finally, the effects of each of these patterns on the occurrence of heavy rainfall were determined by obtaining divergence, vertical velocity, vorticity, relative humidity, specific humidity, temperature and mean sea pressure from the ECMWF center. <br />The results of the study show the effects of three interaction patterns on heavy rainfalls. The effects of direct and indirect jet streams interactions in the first and second patterns, with the change in the thickness of the polar fronts and the interaction of jets in the third model, has been affected on atmospheric conditions and occurrence of heavy rainfall. Despite the direct interactions of the first pattern of jet stream synchronization in the rainfall concentration, the indirect interactions effect of the second pattern on the cyclone circulation boost of the polar jet stream has led to the occurrence of extensive rainfall. Also, the mutual interactions of the jet stream in the third pattern, with the intensification of the frontogenesis process, have been effective in the occurrence of precipitation along these fronts, while the proximity of these opposite turns in the fourth jet streams concurrency pattern is not effective on rainfall, and the occurrence of heavy rainfall has been affected only from the cyclonic circulation of the subtropical jet stream.Institute of Geophysics, University of TehranJournal of the Earth and Space Physics2538-371X45220190723Dynamical analysis of the effect of tropopause folding on the upper and lower-level frontogenesisDynamical analysis of the effect of tropopause folding on the upper and lower-level frontogenesis4014216915410.22059/jesphys.2019.265143.1007034FAMansoorehAbdollahiM.Sc. Student, Department of Space Physics, Institute of Geophysics, University of Tehran, Tehran, IranFarhangAhmadi-GiviAssociate Professor, Department of Space Physics, Institute of Geophysics, University of Tehran, Tehran, Iran0000-0002-9487-4862MohammadMirzaeiAssistant Professor, Department of Space Physics, Institute of Geophysics, University of Tehran, Tehran, Iran0000-0003-0813-3994Journal Article20180922Tropopause folds are intimately linked to upper level frontogenesis and jet stream dynamics. They play an important role for stratosphere-troposphere exchange, the dynamical coupling of upper and lower tropospheric levels, and for generating severe weather events. In this study, the effect of upper-level positive potential vorticity (PV) anomaly on upper- and lower-level frontogenesis over the Middle East and Iran is investigated. In this regard, first three frontal systems associated with deep tropopause folding and strong fronts were selected by using the ECMWF data with the horizontal resolution of 0.75×0.75 degrees on the latitude and longitude. Then, PV anomaly was removed by replacing the zonal mean of PV in the study area and inverted to obtain the modified fields. To do this, a program package (PV inversion), comprising several different steps, is used which allows to isolate PV elements and then to study their impact on the atmospheric flow field as well as the temperature distribution. In the next step, the weather research and forecasting (WRF) model was applied by using the ECMWF data to perform two simulations with real (unchanged) and modified data, as initial conditions, in two domains with 9km resolution for the inner domain. Finally, the frontogenesis function terms, including deformation, tilting, diabatic heating and vertical frontogenesis were computed, using the WRF outputs. By comparing the results of the two simulations, we can determine the effect of tropopause folding on the frontogenesis function terms in the upper- and lower levels throughout the lifecycle of the fronts. <br />Results show that in the absence of tropopause folding, the horizontal and vertical temperature gradients, horizontal velocity, as well as negative vertical velocity are decreased significantly in the upper and lower levels. Also, positive vertical velocity is increased and its pattern is changed mainly in the lower levels. Generally, large positive values of deformation and vertical frontogenesis terms are collocated well with the gradient of potential temperature (frontal zone), and these terms are declined due to decreasing of horizontal velocity and temperature gradient when the fold of tropopause is removed. The tilting term does not have a fixed pattern in the upper and lower levels, and it follows the vertical velocity pattern. By removing the tropopause folding, vertical velocities are changed, thereby having variable effects on the tilting term and the gradient of potential temperature. The diabatic heating term produced by the release of latent heat intensifies ascending motions, and so affects the gradient of potential temperature (frontogenesis) in the lower levels. This term is also increased in the absence of tropopause folding in the lower levels. The other noticeable point is that the diabatic heating term does not significantly affect the frontogenesis in upper levels, because heating mostly takes place in the lower half of the troposphere far to the east of the upper level front. The total amount of frontogenesis function follows the vertical frontogenesis term in the upper and lower levels, and vertical gradient of potential temperature is generally very large. Although, the amount of the diabatic heating term is larger than the vertical frontogenesis term in the lower levels, but it is limited to a small area. Removal of the tropopause fold causes the total amount of frontogenesis function, similar to the vertical frontogenesis term, to decrease throughout the region of the upper-level front, especially in the downstream of the upper-level trough close to the center of tropopause folding omission. The results of the three cases studied here indicate that frontogenesis function terms are considerably more intense in the presence of tropopause folding. Therefore, it is concluded that tropopause folding has a remarkable positive effect on the formation and intensification of the upper-and lower-level fronts.Tropopause folds are intimately linked to upper level frontogenesis and jet stream dynamics. They play an important role for stratosphere-troposphere exchange, the dynamical coupling of upper and lower tropospheric levels, and for generating severe weather events. In this study, the effect of upper-level positive potential vorticity (PV) anomaly on upper- and lower-level frontogenesis over the Middle East and Iran is investigated. In this regard, first three frontal systems associated with deep tropopause folding and strong fronts were selected by using the ECMWF data with the horizontal resolution of 0.75×0.75 degrees on the latitude and longitude. Then, PV anomaly was removed by replacing the zonal mean of PV in the study area and inverted to obtain the modified fields. To do this, a program package (PV inversion), comprising several different steps, is used which allows to isolate PV elements and then to study their impact on the atmospheric flow field as well as the temperature distribution. In the next step, the weather research and forecasting (WRF) model was applied by using the ECMWF data to perform two simulations with real (unchanged) and modified data, as initial conditions, in two domains with 9km resolution for the inner domain. Finally, the frontogenesis function terms, including deformation, tilting, diabatic heating and vertical frontogenesis were computed, using the WRF outputs. By comparing the results of the two simulations, we can determine the effect of tropopause folding on the frontogenesis function terms in the upper- and lower levels throughout the lifecycle of the fronts. <br />Results show that in the absence of tropopause folding, the horizontal and vertical temperature gradients, horizontal velocity, as well as negative vertical velocity are decreased significantly in the upper and lower levels. Also, positive vertical velocity is increased and its pattern is changed mainly in the lower levels. Generally, large positive values of deformation and vertical frontogenesis terms are collocated well with the gradient of potential temperature (frontal zone), and these terms are declined due to decreasing of horizontal velocity and temperature gradient when the fold of tropopause is removed. The tilting term does not have a fixed pattern in the upper and lower levels, and it follows the vertical velocity pattern. By removing the tropopause folding, vertical velocities are changed, thereby having variable effects on the tilting term and the gradient of potential temperature. The diabatic heating term produced by the release of latent heat intensifies ascending motions, and so affects the gradient of potential temperature (frontogenesis) in the lower levels. This term is also increased in the absence of tropopause folding in the lower levels. The other noticeable point is that the diabatic heating term does not significantly affect the frontogenesis in upper levels, because heating mostly takes place in the lower half of the troposphere far to the east of the upper level front. The total amount of frontogenesis function follows the vertical frontogenesis term in the upper and lower levels, and vertical gradient of potential temperature is generally very large. Although, the amount of the diabatic heating term is larger than the vertical frontogenesis term in the lower levels, but it is limited to a small area. Removal of the tropopause fold causes the total amount of frontogenesis function, similar to the vertical frontogenesis term, to decrease throughout the region of the upper-level front, especially in the downstream of the upper-level trough close to the center of tropopause folding omission. The results of the three cases studied here indicate that frontogenesis function terms are considerably more intense in the presence of tropopause folding. Therefore, it is concluded that tropopause folding has a remarkable positive effect on the formation and intensification of the upper-and lower-level fronts.Institute of Geophysics, University of TehranJournal of the Earth and Space Physics2538-371X45220190723A study of the impacts of the NAO on the relation between the North Atlantic and Mediterranean storm tracks using the NCEP/NCAR and JRA-55 reanalysis dataA study of the impacts of the NAO on the relation between the North Atlantic and Mediterranean storm tracks using the NCEP/NCAR and JRA-55 reanalysis data4234407098710.22059/jesphys.2019.267521.1007050FAAmeneMollasharifiM.Sc. Student, Department of Space Physics, Institute of Geophysics, University of Tehran, Tehran, IranAli RezaMohebalhojehProfessor, Department of Space Physics, Institute of Geophysics, University of Tehran, Tehran, Iran0000-0002-5906-8486FarhangAhmadi-GiviAssociate Professor, Department of Space Physics, Institute of Geophysics, University of Tehran, Tehran, Iran0000-0002-9487-4862Journal Article20181017The weather and climate are affected significantly by the storm tracks, which play an important role in midlatitude atmospheric dynamics. There are indications that the Mediterranean region and its downstream areas like the western Asia are affected by the North Atlantic (NA) storm track. As NAO is the most important factor affecting weather in Europe and the Mediterranean region by altering the location and direction of storms, the NA and Mediterranean storm tracks are influenced by this teleconnection. On the basis of the monthly index of the NAO, the critical positive (negative) months, denoted by NAO<sup>+</sup> (NAO<sup>-</sup>), are defined as the ones with the NAO index greater (less) than the long-term mean of the NAO index, plus (minus) one standard deviation. The aim of this research is to study the impacts of the NAO on the NA and Mediterranean storm tracks from the energetic point of view by the JRA-55 reanalysis dataset and to compare the results with those previously obtained by the NCEP/NCAR reanalysis dataset. To this end, composite maps of the most important terms of the time tendency equation of eddy kinetic energy (EKE), as well as the baroclinic generation term, were computed for the critical positive and negative months of the NAO using both datasets in 1959–2017 period for the winter. <br />Based on the both dataset’s results, although though baroclinicity plays the most important role in increasing EKE in the western NA, it is not very important in Western Europe and the Mediterranean region. Instead, convergence of total energy flux is the dominant factor in producing EKE in the latter areas. Interestingly, the ageostrophic flux has a greater role in strengthening the Mediterranean storm track compared to the NA storm track. <br />Results of the JRA-55 are generally consistent with the results of the NCEP/NCAR, and the two datasets lead to almost similar patterns for energetics but some differences are also observed in terms of magnitude and extension of the main centers. In almost all the fields, the JRA-55 results show higher values with more details. The JRA-55 results display stronger and more extended maxima for the NA and Mediterranean storm tracks in both phases. The difference between the storm tracks in the two datasets, in terms of EKE amplitudes, is more clear in the negative phase of NAO. In NAO<sup>-</sup>, the maxima of the NA storm track is about 10% stronger in the JRA-55 results than that of the NCEP/NCAR. Also, both the northern and southern branches of the NA storm track are more intense with a greater eastward extension in the JRA-55. In a similar way, the Mediterranean storm track obtained from the JRA-55 dataset is also more intense with a greater extension, especially over the east and south of the Mediterranean and its downstream regions like the Middle East, than the corresponding results by the NCEP/NCAR. In the negative phase, the difference between the maximum of divergence over the Mediterranean Sea and the maximum of convergence over the Red Sea is larger in the JRA-55. In NAO<sup>- </sup>this can play an important role in making the strength and spread of the Mediterranean storm track in the JRA-55 greater than those of the NCEP/NCAR. In addition, the impact of NAO phases on the intensity of the divergence (convergence) center over the Mediterranean (Red) Sea is greater when the JRA-55 data are being used. <br />Unlike the previous results on the NA storm track as being stronger in the positive phase of the NAO, results of the JRA-55 dataset do not exhibit a noticeable difference between the two phases. In addition, although being consistent with the NCEP/NCAR results in the way the Mediterranean storm track is related to the NA storm track in each phase, however the JRA-55 results show weaker linkage between the two storm tracks in NAO<sup>-</sup> when compared to previous studies.The weather and climate are affected significantly by the storm tracks, which play an important role in midlatitude atmospheric dynamics. There are indications that the Mediterranean region and its downstream areas like the western Asia are affected by the North Atlantic (NA) storm track. As NAO is the most important factor affecting weather in Europe and the Mediterranean region by altering the location and direction of storms, the NA and Mediterranean storm tracks are influenced by this teleconnection. On the basis of the monthly index of the NAO, the critical positive (negative) months, denoted by NAO<sup>+</sup> (NAO<sup>-</sup>), are defined as the ones with the NAO index greater (less) than the long-term mean of the NAO index, plus (minus) one standard deviation. The aim of this research is to study the impacts of the NAO on the NA and Mediterranean storm tracks from the energetic point of view by the JRA-55 reanalysis dataset and to compare the results with those previously obtained by the NCEP/NCAR reanalysis dataset. To this end, composite maps of the most important terms of the time tendency equation of eddy kinetic energy (EKE), as well as the baroclinic generation term, were computed for the critical positive and negative months of the NAO using both datasets in 1959–2017 period for the winter. <br />Based on the both dataset’s results, although though baroclinicity plays the most important role in increasing EKE in the western NA, it is not very important in Western Europe and the Mediterranean region. Instead, convergence of total energy flux is the dominant factor in producing EKE in the latter areas. Interestingly, the ageostrophic flux has a greater role in strengthening the Mediterranean storm track compared to the NA storm track. <br />Results of the JRA-55 are generally consistent with the results of the NCEP/NCAR, and the two datasets lead to almost similar patterns for energetics but some differences are also observed in terms of magnitude and extension of the main centers. In almost all the fields, the JRA-55 results show higher values with more details. The JRA-55 results display stronger and more extended maxima for the NA and Mediterranean storm tracks in both phases. The difference between the storm tracks in the two datasets, in terms of EKE amplitudes, is more clear in the negative phase of NAO. In NAO<sup>-</sup>, the maxima of the NA storm track is about 10% stronger in the JRA-55 results than that of the NCEP/NCAR. Also, both the northern and southern branches of the NA storm track are more intense with a greater eastward extension in the JRA-55. In a similar way, the Mediterranean storm track obtained from the JRA-55 dataset is also more intense with a greater extension, especially over the east and south of the Mediterranean and its downstream regions like the Middle East, than the corresponding results by the NCEP/NCAR. In the negative phase, the difference between the maximum of divergence over the Mediterranean Sea and the maximum of convergence over the Red Sea is larger in the JRA-55. In NAO<sup>- </sup>this can play an important role in making the strength and spread of the Mediterranean storm track in the JRA-55 greater than those of the NCEP/NCAR. In addition, the impact of NAO phases on the intensity of the divergence (convergence) center over the Mediterranean (Red) Sea is greater when the JRA-55 data are being used. <br />Unlike the previous results on the NA storm track as being stronger in the positive phase of the NAO, results of the JRA-55 dataset do not exhibit a noticeable difference between the two phases. In addition, although being consistent with the NCEP/NCAR results in the way the Mediterranean storm track is related to the NA storm track in each phase, however the JRA-55 results show weaker linkage between the two storm tracks in NAO<sup>-</sup> when compared to previous studies.Institute of Geophysics, University of TehranJournal of the Earth and Space Physics2538-371X45220190723A Case study of WRF model performance to hindcast of 10-m wind and 2-m temperature against the satellite and synoptic stations data over the Gulf of Oman and The Arabian SeaA Case study of WRF model performance to hindcast of 10-m wind and 2-m temperature against the satellite and synoptic stations data over the Gulf of Oman and The Arabian Sea4414587098910.22059/jesphys.2019.267709.1007051FAFatemeBahmanzadePh.D. Student, Department of Space Physics, Institute of Geophysics, University of Tehran, Tehran, IranSarmadGhaderAssociate Professor, Department of Space Physics, Institute of Geophysics, University of Tehran, Tehran, Iran0000-0001-9666-5493Seyed AbbasHaghshenasAssistant Professor, Department of Space Physics, Institute of Geophysics, University of Tehran, Tehran, IranDanielYazgiPost-Doc, Department of Space Physics, Institute of Geophysics, University of Tehran, Tehran, IranJournal Article20181021Reliable and sufficient information of 10-m wind and temperature fields over open seas and near coastlines is a necessary and important data that has impact on many marine activities. Assimilation of NWP models can be used to assess an estimation of these fields. This study reports the performance of the weather research and forecasting (WRF) model to hindcast 10-mwind and temperature fields that were evaluated under two different physical options of planetary boundary layer (PBL) and surface layer (SL) for an area over the Gulf of Oman and the Arabian Sea. The case study includes 16 simulations of 8 different days from WRF model version 3.7.1. The WRF model is configured with two nests. Parent nest has 0.3 degree and the inner nest has 0.1 degree horizontal grid resolution. The grid spacing of the inner domain is almost 11-km. The Lat-Lon (latitude-longitude) method is used as the map projection method. For all domains and all runs 39 terrain following vertical levels are set. The validation of the simulated fields is done considering two observational datasets (the weather stations for 10-m wind and 2-m temperature and satellite instruments just for 10-m wind). Near-surface observations of 2-m temperature and 10-m wind speed and direction are collected from 55 weather stations, located within the chosen area. The measurements from satellite instruments have become an important source of data in the regions that in-situ observations are sparse like seas and oceans, hence observations from two different scatterometers (ASCAT and OSCAT) are also used to evaluate 10-m wind simulations. Moreover, in order to better understand the model performance for different choices of the physical schemes, sensitivity of the model has been investigated. There is plenty of choices for the combination of parameterization schemes available for WRF model; for the current study two configurations are taken from other’s previous published research works. The physical parameterization that used in this study are Revised MM5 and Monin-Obukhov for surface layer and MRF and MYJ for planetary boundary layer. These choices are used to create two different configurations called Phys1 and Phys2. Comparison between winds from satellite scatterometer and simulated winds show very little difference and hence good agreement with observations. Acceptable accuracy has been obtained from statistical analyses. These analyses demonstrate that the maximum average RMSE of wind field is 2.39 m/s, based on results of comparing with ASCAT data and it is 2.37 m/s, based on results of comparing with OSCAT data. The analyses also show that simulation of wind fields have better results over offshore regions than coastlines weather stations. The outcome shows that the simulated 10-m wind present acceptable general skills over the sea. The validation of 2-m temperature presents that the model has a proper estimation about temperature field over the coasts and near coastal station within the simulation domain. The maximum average RMSE of temperature field is 2.6 degrees of centigrade. Finally, without any justification to run WRF for longer periods from a quantitative and qualitative assessment of the results, it can be concluded that for the WRF model has an acceptable performance to simulate 10-m wind and 2-m temperature over the Gulf of Oman and the Arabian Sea. It should be noted that to verify these results for longer periods more similar experiments must be tested.Reliable and sufficient information of 10-m wind and temperature fields over open seas and near coastlines is a necessary and important data that has impact on many marine activities. Assimilation of NWP models can be used to assess an estimation of these fields. This study reports the performance of the weather research and forecasting (WRF) model to hindcast 10-mwind and temperature fields that were evaluated under two different physical options of planetary boundary layer (PBL) and surface layer (SL) for an area over the Gulf of Oman and the Arabian Sea. The case study includes 16 simulations of 8 different days from WRF model version 3.7.1. The WRF model is configured with two nests. Parent nest has 0.3 degree and the inner nest has 0.1 degree horizontal grid resolution. The grid spacing of the inner domain is almost 11-km. The Lat-Lon (latitude-longitude) method is used as the map projection method. For all domains and all runs 39 terrain following vertical levels are set. The validation of the simulated fields is done considering two observational datasets (the weather stations for 10-m wind and 2-m temperature and satellite instruments just for 10-m wind). Near-surface observations of 2-m temperature and 10-m wind speed and direction are collected from 55 weather stations, located within the chosen area. The measurements from satellite instruments have become an important source of data in the regions that in-situ observations are sparse like seas and oceans, hence observations from two different scatterometers (ASCAT and OSCAT) are also used to evaluate 10-m wind simulations. Moreover, in order to better understand the model performance for different choices of the physical schemes, sensitivity of the model has been investigated. There is plenty of choices for the combination of parameterization schemes available for WRF model; for the current study two configurations are taken from other’s previous published research works. The physical parameterization that used in this study are Revised MM5 and Monin-Obukhov for surface layer and MRF and MYJ for planetary boundary layer. These choices are used to create two different configurations called Phys1 and Phys2. Comparison between winds from satellite scatterometer and simulated winds show very little difference and hence good agreement with observations. Acceptable accuracy has been obtained from statistical analyses. These analyses demonstrate that the maximum average RMSE of wind field is 2.39 m/s, based on results of comparing with ASCAT data and it is 2.37 m/s, based on results of comparing with OSCAT data. The analyses also show that simulation of wind fields have better results over offshore regions than coastlines weather stations. The outcome shows that the simulated 10-m wind present acceptable general skills over the sea. The validation of 2-m temperature presents that the model has a proper estimation about temperature field over the coasts and near coastal station within the simulation domain. The maximum average RMSE of temperature field is 2.6 degrees of centigrade. Finally, without any justification to run WRF for longer periods from a quantitative and qualitative assessment of the results, it can be concluded that for the WRF model has an acceptable performance to simulate 10-m wind and 2-m temperature over the Gulf of Oman and the Arabian Sea. It should be noted that to verify these results for longer periods more similar experiments must be tested.Institute of Geophysics, University of TehranJournal of the Earth and Space Physics2538-371X45220190723A synoptic-scale investigation of entropy fluxes during Tropical Cyclone GonuA synoptic-scale investigation of entropy fluxes during Tropical Cyclone Gonu4594727136210.22059/jesphys.2019.272345.1007077FANafisehPegahfarAssistant Professor, Atmospheric Science Center, Iranian National Institute for Oceanography and Atmospheric Science, Tehran, Iran0000-0003-4885-7428Journal Article20181229Tropical cyclones (TC) have been investigated from different points of view. Development of forecast of TC intensity and its track is often the shared purpose of all previous researches. To this aim, various empirical indices and different frameworks, based on various parameters, have been defined to provide deep knowledge of TC dynamics and thermodynamics. In this research, using the thermodynamic parameter of entropy, entropy fluxes (including surface, lateral and vertical fluxes) have been calculated. A theoretical framework based on hypothesized mechanism, introduced by Tang and Emanuel (2010), has been used to calculate the vertical flux of entropy. This ideal framework used a set of rigid assumptions including steadiness, axisymmetry and slantwise neutrality to assess the effects of vertical entropy flux on TC intensity via the possible pathway of downdrafts outside the eyewall. The lateral entropy flux has been computed based on radial component of surface wind. Azimuthal average of lateral entropy flux has been calculated to analyze vertical extension and strength of inflow (in the lower part of boundary layer) and also outflow (in the upper part of troposphere). Here, Tropical Cyclone Gonu (TCG) has been focusedon. TCG, formed at 18:00 UTC 1 June 2007 and decayed on 7 June, passed intensity of Saffir-Simpson Category-5 and affected southern coast (Makran) of Iran. All above parameters have been computed and analyzed during TCG lifetime using (1) Era-Interim reanalysis data (from European Center for Medium Range Weather Forecast) with 0.125 degree horizontal resolution, 12 vertical levels from 1000 to 200 hPa and 6-hour time intervals, and (2) data produced by India Meteorological Department. The variables were used both at the surface and also at pressure levels, the surface values were temperature and humidity (both at 2 m height), wind vector (at 10 m height), mixing ratio and sea level pressure. Synoptic–scale analysis has been done using data in a circular region centered by TCG center with a radius of 500 km. Results of horizontal patterns and time series of radial and tangential components of wind vector indicated that the value of radial component was maximized simultaneously with maximum activity of TCG. At TCG peak activity time, the tangential component had a comparatively minimum value embedded between two relative maximum values. Time series analysis showed that the integrated values of all three parameters of surface, vertical and lateral entropy fluxes experienced their extremum values before TCG reached its maximum intensity. It is worthwhile to be noted that their lead time varied from 6 hours (for surface entropy flux), 18 hours (for lateral entropy flux) to 30 hours (for vertical entropy flux). A comparative analysis between the values of entropy fluxes during TCG and those for Haiyan Tropical Cyclone (TCH, the strongest TC formed over the Pacific Ocean), reported by Pegahfar and Gharaylou (2019), indicated that entropy surface flux and lateral entropy flux during TCG were respectively two and one order of magnitude larger than the related values during TCH. In contrast, TCG experienced entropy vertical flux with two orders of magnitude smaller than that during TCH. Hence it can be concluded that the accumulation of energy helped TCG to travel to the higher latitudes. Moreover, the strongest inflow and outflow occurred before and after TCG maximum intensity, respectively. In a period that TCG reaches category-5 intensity and more, vertical extension of inflow layer was minimized while outflow layer started from the lower levels, comparing with results from the beginning of TCG life cycle. Conclusively, findings of this research showed that surface, vertical and lateral entropy fluxes, even in synoptic scale, have the ability to be served as empirical indices and also need to be focused in theoretical, computational and practical frameworks, for all prognostic purposes of TC intensity.Tropical cyclones (TC) have been investigated from different points of view. Development of forecast of TC intensity and its track is often the shared purpose of all previous researches. To this aim, various empirical indices and different frameworks, based on various parameters, have been defined to provide deep knowledge of TC dynamics and thermodynamics. In this research, using the thermodynamic parameter of entropy, entropy fluxes (including surface, lateral and vertical fluxes) have been calculated. A theoretical framework based on hypothesized mechanism, introduced by Tang and Emanuel (2010), has been used to calculate the vertical flux of entropy. This ideal framework used a set of rigid assumptions including steadiness, axisymmetry and slantwise neutrality to assess the effects of vertical entropy flux on TC intensity via the possible pathway of downdrafts outside the eyewall. The lateral entropy flux has been computed based on radial component of surface wind. Azimuthal average of lateral entropy flux has been calculated to analyze vertical extension and strength of inflow (in the lower part of boundary layer) and also outflow (in the upper part of troposphere). Here, Tropical Cyclone Gonu (TCG) has been focusedon. TCG, formed at 18:00 UTC 1 June 2007 and decayed on 7 June, passed intensity of Saffir-Simpson Category-5 and affected southern coast (Makran) of Iran. All above parameters have been computed and analyzed during TCG lifetime using (1) Era-Interim reanalysis data (from European Center for Medium Range Weather Forecast) with 0.125 degree horizontal resolution, 12 vertical levels from 1000 to 200 hPa and 6-hour time intervals, and (2) data produced by India Meteorological Department. The variables were used both at the surface and also at pressure levels, the surface values were temperature and humidity (both at 2 m height), wind vector (at 10 m height), mixing ratio and sea level pressure. Synoptic–scale analysis has been done using data in a circular region centered by TCG center with a radius of 500 km. Results of horizontal patterns and time series of radial and tangential components of wind vector indicated that the value of radial component was maximized simultaneously with maximum activity of TCG. At TCG peak activity time, the tangential component had a comparatively minimum value embedded between two relative maximum values. Time series analysis showed that the integrated values of all three parameters of surface, vertical and lateral entropy fluxes experienced their extremum values before TCG reached its maximum intensity. It is worthwhile to be noted that their lead time varied from 6 hours (for surface entropy flux), 18 hours (for lateral entropy flux) to 30 hours (for vertical entropy flux). A comparative analysis between the values of entropy fluxes during TCG and those for Haiyan Tropical Cyclone (TCH, the strongest TC formed over the Pacific Ocean), reported by Pegahfar and Gharaylou (2019), indicated that entropy surface flux and lateral entropy flux during TCG were respectively two and one order of magnitude larger than the related values during TCH. In contrast, TCG experienced entropy vertical flux with two orders of magnitude smaller than that during TCH. Hence it can be concluded that the accumulation of energy helped TCG to travel to the higher latitudes. Moreover, the strongest inflow and outflow occurred before and after TCG maximum intensity, respectively. In a period that TCG reaches category-5 intensity and more, vertical extension of inflow layer was minimized while outflow layer started from the lower levels, comparing with results from the beginning of TCG life cycle. Conclusively, findings of this research showed that surface, vertical and lateral entropy fluxes, even in synoptic scale, have the ability to be served as empirical indices and also need to be focused in theoretical, computational and practical frameworks, for all prognostic purposes of TC intensity.Institute of Geophysics, University of TehranJournal of the Earth and Space Physics2538-371X45220190723Hydrological drought response of regulated river flow under the influence of dam reservoir in Ardabil ProvinceHydrological drought response of regulated river flow under the influence of dam reservoir in Ardabil Province4734867096810.22059/jesphys.2019.272671.1007078FAHamedAminiM.Sc. Student, Department of Natural Resources, Faculty of Agriculture and Natural Resources, University of
Mohaghegh Ardabili, Ardabil, IranAbazarEsmali-OuriAssociate Professor, Department of Natural Resources, Faculty of Agriculture and Natural Resources, University of Mohaghegh Ardabili, Ardabil, Iran0000-0001-9227-0419RaoofMostafazadehAssistant Professor, Department of Natural Resources, Faculty of Agriculture and Natural Resources, University of Mohaghegh Ardabili, Ardabil, Iran0000-0002-0401-0260MerajSharariAssistant Professor, Department of Natural Resources, Faculty of Agriculture and Natural Resources, University of
Mohaghegh Ardabili, Ardabil, IranMohsenZabihiPh.D. Student, Department of Watershed management, Faculty of Agriculture and Natural Resources, University of Tarbiat Modares, Tehran, IranJournal Article20190101Due to the importance of rivers in water supply, it is crucial to understand the flow condition and river regime. The river flows are affected by dam construction, climate change, flow diversion, backwater and urban development, which cause various problems, such as drought, flow reduction and environmental issues. Hydraulic structures have great impacts on the hydraulic conditions of rivers downstream. The reservoir dams as the most important hydraulic structures constructed on rivers, have considerable changes on flow regime conditions, therefore, drought and river flow drying area are major impacts of dams on downstream rivers. One of the common types of droughts is the hydrological drought, which is followed by meteorological drought with a delayed time. The most commonly definition of hydrological drought is based on a truncation level. The hydrological drought occurs when the river flow falls below a defined threshold level. The Yamchi and Sabalan dams are two main important dams in Ardabil Province due to the water supply for drinking water and irrigated agricultural lands. <br />Due to the consequences of short and long-term socio-economic and environmental effects of droughts as a natural disaster, it is necessary to study the frequency, severity, and spatial extent of drought spells. Accordingly, the drought characterization plays an important role in macroeconomic planning and decision-making to achieve sustainable development goals. The purpose of this study is to evaluate the effect of the Yamchi and Sabalan dams on the flow regime and hydrological drought occurrence of the downstream rivers, areas. <br />The used data in this study were the daily discharges, which were obtained from four hydrometric stations located upstream and downstream of the dams under the study. In this study, the seasonal threshold of was chosen and drought characteristics of the upstream and downstream stations of each of the dams were evaluated and then the drought severity at different durations was calculated using the Runs theory. Then, the Easy Fit software was employed to calculate the drought intensity at different return periods. The drought severity amounts were plotted over different return periods and the Severity-Duration-Frequency curves that were drawn for the upland and downstream of the constructed dam reservoirs. The percent change in discharge at the upstream and downstream stations of the studied dams were calculated to quantify the degree of influence. <br />According to the results, the highest drought occurrence in both studied dams were observed at 1-month duration and at Pole-almas and Dost-beiglou river gauge stations. It can be concluded that the severity of the observed hydrological droughts increases under the influence of both Yamchi and Sabalan dams as the return period increases. Also, the most severe drought spells occurred is related to the Sabalan dam and its upstream station (Arbabkandi). According to SDF curves, the drought severity increases in higher drought durations in the upstream and downstream stations of the Yamchi dam. According to the results, the percentage of changes in the downstream hydrometery station of the Sabalan dam is positive in almost all water years, which confirms the positive effect of Sabalan Dam on hydrological drought reduction and flow regulation. The highest percentage of changes in both of these dams is for to June, which indicates the difference in discharge values in upstream and downstream of the dams. In other words, this indicates that both of the studied dams, especially the Sabalan Dam, have increased the monthly base flow rate over the downstream reaches. <br />It can be concluded that the construction of the Yamchi Dam caused a decrease in the flow amount in the downstream area and increased the occurrence of hydrological drought events. While, the Sabalan Dam regulated the flow regime of the Garesou river in the dam downstream and reduced the severity of hydrological drought spells and modifies the river flow regime for all months of the year.Due to the importance of rivers in water supply, it is crucial to understand the flow condition and river regime. The river flows are affected by dam construction, climate change, flow diversion, backwater and urban development, which cause various problems, such as drought, flow reduction and environmental issues. Hydraulic structures have great impacts on the hydraulic conditions of rivers downstream. The reservoir dams as the most important hydraulic structures constructed on rivers, have considerable changes on flow regime conditions, therefore, drought and river flow drying area are major impacts of dams on downstream rivers. One of the common types of droughts is the hydrological drought, which is followed by meteorological drought with a delayed time. The most commonly definition of hydrological drought is based on a truncation level. The hydrological drought occurs when the river flow falls below a defined threshold level. The Yamchi and Sabalan dams are two main important dams in Ardabil Province due to the water supply for drinking water and irrigated agricultural lands. <br />Due to the consequences of short and long-term socio-economic and environmental effects of droughts as a natural disaster, it is necessary to study the frequency, severity, and spatial extent of drought spells. Accordingly, the drought characterization plays an important role in macroeconomic planning and decision-making to achieve sustainable development goals. The purpose of this study is to evaluate the effect of the Yamchi and Sabalan dams on the flow regime and hydrological drought occurrence of the downstream rivers, areas. <br />The used data in this study were the daily discharges, which were obtained from four hydrometric stations located upstream and downstream of the dams under the study. In this study, the seasonal threshold of was chosen and drought characteristics of the upstream and downstream stations of each of the dams were evaluated and then the drought severity at different durations was calculated using the Runs theory. Then, the Easy Fit software was employed to calculate the drought intensity at different return periods. The drought severity amounts were plotted over different return periods and the Severity-Duration-Frequency curves that were drawn for the upland and downstream of the constructed dam reservoirs. The percent change in discharge at the upstream and downstream stations of the studied dams were calculated to quantify the degree of influence. <br />According to the results, the highest drought occurrence in both studied dams were observed at 1-month duration and at Pole-almas and Dost-beiglou river gauge stations. It can be concluded that the severity of the observed hydrological droughts increases under the influence of both Yamchi and Sabalan dams as the return period increases. Also, the most severe drought spells occurred is related to the Sabalan dam and its upstream station (Arbabkandi). According to SDF curves, the drought severity increases in higher drought durations in the upstream and downstream stations of the Yamchi dam. According to the results, the percentage of changes in the downstream hydrometery station of the Sabalan dam is positive in almost all water years, which confirms the positive effect of Sabalan Dam on hydrological drought reduction and flow regulation. The highest percentage of changes in both of these dams is for to June, which indicates the difference in discharge values in upstream and downstream of the dams. In other words, this indicates that both of the studied dams, especially the Sabalan Dam, have increased the monthly base flow rate over the downstream reaches. <br />It can be concluded that the construction of the Yamchi Dam caused a decrease in the flow amount in the downstream area and increased the occurrence of hydrological drought events. While, the Sabalan Dam regulated the flow regime of the Garesou river in the dam downstream and reduced the severity of hydrological drought spells and modifies the river flow regime for all months of the year.