Radar interferometry time series analysis of Mashhad subsidence

A large agricultural area located in 20 km north of the city of Mashhad in the north-east of Iran is subject to land subsidence. The subsidence rate was achieved in a couple of sparse points by precise leveling between 1995 and 2005, and continuous GPS measurements obtained from 2005 to 2006. In order to study the temporal behavior of the deformation in high spatial resolution, the small baseline subset (SBAS) algorithm was used to generate the interferometric SAR time series analysis. Time series analysis was performed using 19 interferograms calculated from 12 ENVISAT ASAR data spanning between 2003 and 2006. The time series results exhibited that the area is subsiding continuously without a significant seasonal effect. Mean LOS deformation velocity map obtained from time series analysis demonstrated a considerable subsidence rate up to 24 (cm/yr). In order to evaluate the time series analysis results, continuous GPS measurements as a geodetic approach were applied. The comparisons showed a great agreement between interferometry results and geodetic technique. Moreover, the information of various piezometric wells distributed in the area corresponding to 1995 to 2005 showed a significant decline in water table up to 20 meters. The correlation between the piezometric information and the surface deformation at well’s locations showed that the subsidence occurrence in Mashhad is due to the excess groundwater withdrawal.     

Assessment of changes in urban green spaces of Mashad city using satellite data

Green spaces play important functions in urban environments. Reducing air pollution, providing shade and habitat for arboreal birds, producing oxygen, providing shelter against winds, recreational and aesthetic qualities and architectural applications are the main functions of urban green spaces. With the rapid change of urban area in Mashad city during the past decades, green spaces have been fragmented and dispersed causing impairment and dysfunction of these important urban elements. The objective of this study was to detect changes in extent and pattern of green areas of Mashad city and to analyze the results in terms of landscape ecology principles and functioning of the green spaces. In this research, we classified a Landsat TM and an IRS LISS-III image belonging to the years 1987 and 2006, respectively. We then used a post-classification comparison to determine the changes in green space areas of Mashad city during the 19 years covered by the images. Then, we applied landscape ecology calculations to derive metrics that quantified pattern of the changes in the green areas. The results showed that during 19 years from 1987, a significant decrease had occurred in the extent of urban green spaces with a concomitant fragmentation resulting in downgrading and destruction of the functions and services these areas provide. We conclude that the general quality of life in the central parts of the city has been diminished. We also state that a combination of remote sensing image classification, landscape metrics assessment and vegetation indices can provide a tool for assessing life quality and its trend for urban areas.     

Isoparametric Boundary Element Method Applied to the Rectangular and Delta Hydrofoils near the Free Surface

In this paper, the hydrodynamic characteristics and flow field around rectangular and delta hydrofoils, moving with a constant speed beneath the free surface are numerically studied by means of isoparametric boundary element method (IBEM). The quantities (source and dipole strengths) and the geometry of the elements are represented by a linear distribution. Two types of three-dimensional hydrofoils (rectangular and delta) are selected with NACA4412 and symmetric Joukowski sections. Some numerical results of pressure distribution, lift, wave-making drag coefficients and velocity field around the hydrofoils are presented. Also, the wave pattern due to moving hydrofoil is predicted at different operational conditions. Comparisons are made between computational results obtained through this method and those from the experimental measurements and other numerical results which reveal good agreement.     

Conjunction of 2D and 3D modified flow solvers for simulating spatio-temporal wind induced hydrodynamics in the Caspian Sea

The main objective of this study is the simulation of flow dynamics in the deep parts of the Caspian Sea, in which the southern and middle deep regions are surrounded by considerable areas of shallow zones. To simulate spatio-temporal wind induced hydrodynamics in deep waters, a conjunctive numerical model consisting of a 2D depth average model and a 3D pseudo compressible model is proposed. The 2D model is applied to determine time dependent free surface oscillations as well as the surface velocity patterns and is conjunct to the 3D flow solver for computing three-dimensional velocity and pressure fields which coverage to steady state for the top boundary condition. The modified 2D and 3D sets of equations are conjunct considering interface shear stresses. Both sets of 2D and 3D equations are solved on unstructured triangular and tetrahedral meshes using the Galerkin Finite Volume Method. The conjunctive model is utilized to investigate the deep currents affected by wind, Coriolis forces and the river inflow conditions of the Caspian Sea. In this study, the simulation of flow field due to major winds as well as transient winds in the Caspian Sea during a period of 6 hours in the winter season has been conducted and the numerical results for water surface level are then compared to the 2D numerical results.     

Hybrid frequency–time domain models for dynamic response analysis of marine structures

Time-domain models of marine structures based on frequency domain data are usually built upon the Cummins equation. This type of model is a vector integro-differential equation which involves convolution terms. These convolution terms are not convenient for analysis and design of motion control systems. In addition, these models are not efficient with respect to simulation time, and ease of implementation in standard simulation packages. For these reasons, different methods have been proposed in the literature as approximate alternative representations of the convolutions. Because the convolution is a linear operation, different approaches can be followed to obtain an approximately equivalent linear system in the form of either transfer function or state-space models. This process involves the use of system identification, and several options are available depending on how the identification problem is posed. This raises the question whether one method is better than the others. This paper therefore has three objectives. The first objective is to revisit some of the methods for replacing the convolutions, which have been reported in different areas of analysis of marine systems: hydrodynamics, wave energy conversion, and motion control systems. The second objective is to compare the different methods in terms of complexity and performance. For this purpose, a model for the response in the vertical plane of a modern containership is considered. The third objective is to describe the implementation of the resulting model in the standard simulation environment Matlab/Simulink.     

Reaching law based sliding mode control for a frame structure under seismic load

Implementation of efficient vibration control schemes for seismically excited structures is becoming more and more important in recent years. In this study, the influence of different control schemes on the dynamic performance of a frame structure excited by El Centro wave, with an emphasis on reaching law based control strategies, is examined. Reaching law refers to the reachable problem and criteria for the sliding state of a control system. Three reaching laws are designed to present different sliding mode control strategies by incorporating a state space model that describes structural dynamic characteristics of a frame structure. Both intact and damaged structures are studied by using the aforementioned control strategies. The influence of different structural damage extents, control locations and reaching law based control methods are further investigated. The results show that the structure can be well controlled using the sliding mode strategy when the induced structural damage extent does not exceed the standard percentage for considering the structure was damaged, which is 20% reduction in structure stiffness, as reported in the literature. The control effectiveness is more satisfactory if the control location is the same as the direction of external excitation. Furthermore, to study the chattering phenomenon of the sliding mode control method, approximation and detail components extracted from the phase plots of the sliding mode control system are compared via wavelet transform at different scales. The results show that for the same type of control law, the system behaves with similar chattering phenomenon.     

Shoreline spatial and temporal response to natural and human effects in Boujagh National Park,Iran

Shoreline variation and river deltas are among the most dynamic systems in marine environments. The related different variations in spatial and temporal scales play significant roles in land planning and different management applications. Modeling the dynamics of seashore of Boujagh National Park (BNP) which is located on the southern coast of the Caspian Sea in the Sefidrud Delta (SD), considering natural and anthropogenic factors, was the main objective of the current study. To achieve this goal, a combination of remote sensing data, historical data, and numerical simulations was utilized. The BNP covers an area of 3,270 ha and includes two international wetlands, Boujagh and Kiashahr. In earlier periods, this area faced severe morphological changes whereas recently its shoreline has experienced gradual variations. Accordingly, at the first stage, the shoreline variation from 2006 to 2017 was extracted by processing and classifying Operational Land Imager (OLI) and Thematic Mapper (TM) images from Landsat satellites using the Maximum Likelihood approach. In the second stage, the two dimensional MIKE21 model was utilized to identify wave and coastal current patterns and parameters for the year 2015. Morphologically, the results showed that, the shoreline of the BNP is affected by several natural and anthropogenic factors. Seaward advancement of the shoreline occurred in zones A (east zone) and C (west zone) due to Caspian Sea Level drop and sedimentation while retreating occurred at Zone B (north zone) influenced by wave and current patterns and reduction of the Sefidrud River flows. Also, the results imply that maintaining the existing conditions results in the disappearance of a considerable part of the ecological area in the BNP. Hence, to manage and preserve the coastline of the BNP complying with the current anthropogenic and natural factors, it is vital to take necessary management measures.     

Bridging the gap between numerical solutions of travel time distributions and analytical storage selection functions

Recent advancements in analytical solutions to quantify water and solute travel time distributions (TTDs) and the related StorAge Selection (SAS) functions synthesize catchment complexity into a simplified, lumped representation. Although these analytical approaches are efficient in application, they require rarely available long‐term and high‐frequency hydrochemical data for parameter estimation. Alternatively, integrated hydrologic models coupled to Lagrangian particle‐tracking approaches can directly simulate age under different catchment geometries and complexity, but at a greater computational expense. Here, we bridge the two approaches, using a physically based model to explore the uncertainty in the estimation of the SAS function shape. In particular, we study the influence of subsurface heterogeneity, interactions between distinct flow domains (i.e., the vadose zone and saturated groundwater), diversity of flow pathways, and recharge rate on the shape of TTDs and the SAS functions. We use an integrated hydrology model, ParFlow, linked with a particle‐tracking model, SLIM, to compute transient residence times (or ages) at every cell in the domain, facilitating a direct characterization of the SAS function. Steady‐state results reveal that the SAS function shape shows a wide range of variation with respect to the variability in the structure of subsurface heterogeneity. Ensembles of spatially correlated realizations of hydraulic conductivity indicate that the SAS functions in the saturated groundwater have an overall weak tendency toward sampling younger ages, whereas the vadose zone gives a strong preference for older ages. We further show that the influence of recharge rate on the TTD is tightly dependent on the variability of subsurface hydraulic conductivity.     

A bridge between dual porosity and multiscale models of heterogeneous deformable porous media

In this paper, a multiscale homogenization approach is developed for fully coupled saturated porous media to represent the idealized sugar cube model, which is generally employed in fractured porous media on the basis of dual porosity models. In this manner, an extended version of the Hill-Mandel theory that incorporates the microdynamic effects into the multiscale analysis is presented, and the concept of the deformable dual porosity model is demonstrated. Numerical simulations are performed employing the multiscale analysis and dual porosity model, and the results are compared with the direct numerical simulation through 2 numerical examples. Finally, a combined multiscale-dual porosity technique is introduced by employing a bridge between these 2 techniques as an alternative approach that reduces the computational cost of numerical simulation in modeling of heterogeneous deformable porous media.