Geological Engineering Commons^™

Constructing A Complete Landslide Inventory Dataset For The 2018 Monsoon Disaster In Kerala, India, For Land Use Change Analysis, Lina Hao, Rajaneesh A., Cees Van Westen, Sajinkumar K. S., Tapas Ranjan Martha, Pankaj Jaiswal, Brian G. Mcadoo

Michigan Tech Publications, Part 1

Event-based landslide inventories are important for analyzing the relationship between the intensity of the trigger (e.g., rainfall, earthquake) and the density of the landslides in a particular area as a basis for the estimation of the landslide probability and the conversion of susceptibility maps into hazard maps required for risk assessment. They are also crucial for the establishment of local rainfall thresholds that are the basis of early warning systems and for evaluating which land use and land cover changes are related to landslide occurrence. The completeness and accuracy of event-based landslide inventories are crucial aspects to derive reliable results …

Using Tandem-X To Measure Pyroclastic Flow Source Location, Thickness And Volume: Application To The 3rd June 2018 Eruption Of Fuego Volcano, Guatemala, F. Albino, J. Biggs, R. Escobar-Wolf, A. Naismith, M. Watson, J. C. Phillips, G. A. Chigna Marroquin

Michigan Tech Publications, Part 1

The estimation of the volume of volcanic flows during an ongoing eruption is challenging but this information is crucial for improving risk assessment and for forecasting future events. Although previous studies have shown the ability of TanDEM-X satellite data to derive the thickness and the volume of lava flow fields during effusive eruptions, the method has not been explored yet for pyroclastic flows. Using bi-static interferometry, we produce TanDEM-X DEM on Fuego volcano (Guatemala) to measure the significant topographic changes caused by the 3rd June 2018 eruption, which destroyed the town of San Miguel Los Lotes. We estimate the volume …

Effect Of Geometry And Fluid Viscosity On Dynamics Of Fluid-Filled Cracks: Insights From Analog Experimental Observations, Haitao Cao, Ezequiel F. Medici, Gregory P. Waite, Roohollah Askari

Michigan Tech Publications, Part 1

Fluid-filled volumes in geological systems can change the local stress field in the host rock and may induce brittle deformation as well as crack propagation. Although the mechanisms relating fluid pressure perturbations and seismicity have been widely studied, the fluid-solid interaction inside the crack of a host rock is still not well understood. An analog experimental model of fluid intrusion in cracks between planar layers has been developed to study stress conditions at the margins and tips. A combined high-speed shadowgraph and a photoelasticity imaging system is used to visualize the fluid dynamics and induced stresses on the solid matrix. …

Design And Development Of A Machine Vision System Using Artificial Neural Network-Based Algorithm For Automated Coal Characterization, Amit Kumar Gorai, Simit Raval, Ashok Kumar Patel, Snehamoy Chatterjee, Tarini Gautam

Michigan Tech Publications, Part 1

Coal is heterogeneous in nature, and thus the characterization of coal is essential before its use for a specific purpose. Thus, the current study aims to develop a machine vision system for automated coal characterizations. The model was calibrated using 80 image samples that are captured for different coal samples in different angles. All the images were captured in RGB color space and converted into five other color spaces (HSI, CMYK, Lab, xyz, Gray) for feature extraction. The intensity component image of HSI color space was further transformed into four frequency components (discrete cosine transform, discrete wavelet transform, discrete Fourier …

Global Aeolian Dust Variations And Trends: A Revisit Of Dust Event And Visibility Observations From Surface Weather Stations, Xin Xi

Michigan Tech Publications, Part 1

This study revisits the use of horizontal visibility and manually reported present weather (ww) records from the NOAA Integrated Surface Database (ISD) for characterizing the aeolian dust variability and recent trends over the globe and three largest source regions (North Africa, Middle East, and East Asia). Due to its qualitative nature, ww is combined with visibility to derive a new variable, VI, which has higher correlations with the dust emission and burden from satellite observations and global aerosol reanalyses than does the dust event frequency (FR) derived from ww only. Both FR and VI capture …

Utilizing Hyperspectral Remote Sensing For Soil Gradation, Jordan Ewing, T. Oommen, Paramsothy Jayakumar, Russell Alger

Michigan Tech Publications, Part 1

Soil gradation is an important characteristic for soil mechanics. Traditionally soil gradation is performed by sieve analysis using a sample from the field. In this research, we are interested in the application of hyperspectral remote sensing to characterize soil gradation. The specific objective of this work is to explore the application of hyperspectral remote sensing to be used as an alternative to traditional soil gradation estimation. The advantage of such an approach is that it would provide the soil gradation without having to obtain a field sample. This work will examine five different soil types from the Keweenaw Research Center …

Flow Simulation Considering Adsorption Boundary Layer Based On Digital Rock And Finite Element Method, Yong Fei Yang, Ke Wang, Qian Fei Lv, Roohollah Askari, Qing Yan Mei, Jun Yao, Jie Xin Hou, Kai Zhang, Ai Fen Li, Chen Chen Wang

Michigan Tech Publications, Part 1

Due to the low permeability of tight reservoirs, throats play a significant role in controlling fluid flow. Although many studies have been conducted to investigate fluid flow in throats in the microscale domain, comparatively fewer works have been devoted to study the effect of adsorption boundary layer (ABL) in throats based on the digital rock method. By considering an ABL, we investigate its effects on fluid flow. We build digital rock model based on computed tomography technology. Then, microscopic pore structures are extracted with watershed segmentation and pore geometries are meshed through Delaunay triangulation approach. Finally, using the meshed digital …

Paleomagnetic Evidence For Modern-Like Plate Motion Velocities At 3.2 Ga, Alec R. Brenner, Roger R. Fu, David A.D. Evans, Aleksey V. Smirnov, Raisa Trubko, Ian R. Rose

Michigan Tech Publications, Part 1

The mode and rates of tectonic processes and lithospheric growth during the Archean [4.0 to 2.5 billion years (Ga) ago] are subjects of considerable debate. Paleomagnetism may contribute to the discussion by quantifying past plate velocities. We report a paleomagnetic pole for the ~3180 million year (Ma) old Honeyeater Basalt of the East Pilbara Craton, Western Australia, supported by a positive fold test and micromagnetic imaging. Comparison of the 44°±15° Honeyeater Basalt paleolatitude with previously reported paleolatitudes requires that the average latitudinal drift rate of the East Pilbara was ≥2.5 cm/year during the ~170 Ma preceding 3180 Ma ago, a …

Ambient Noise Seismic Interferometry Through Surface Seismic Data Acquisition Vs Zero-Offset Vertical Seismic Profiling (Vsp), Erdi Apatay

Dissertations, Master's Theses and Master's Reports

Ambient noise seismic interferometry (ANSI) is a method that uses passive noise created by unknown sources such as local or teleseismic earthquakes, waves, wind, fracturing and fault movement in a reservoir to obtain an image of subsurface. ANSI has been used in some various disciplines, including earthquake engineering, time-lapse monitoring, borehole seismology and seismic exploration for years. Compared to conventional active seismic surveys, it is advantageous in terms of cost-efficiency as no active source such dynamite or vibroseis is used. Although many promising results have been obtained from surface wave ANSI, body wave ANSI is still challenging. Such a challenge …

Physical Modeling Of The Krauklis Waves: Insights From Two Experimental Apparatuses, Haitao Cao

Dissertations, Master's Theses and Master's Reports

Fractures play an important role in the geological related processes such as hydraulic fracturing, water-water disposal, and volcanic earthquake. Seismic waves can provide useful information from fractures at a relatively low cost. In particular, the acoustic property of fractures containing magmatic or hydrothermal fluids can provide useful information about the fracture size and the fluid composition within the fracture. For instance, in volcanology, the resonant frequency of long-period events that are linked to crack interface waves is used to obtain fluid properties of cracks in magmatic systems. However, in order to rely on seismic data, they should be precisely characterized …

Developing Innovative Spectral And Machine Learning Methods For Mineral And Lithological Classification Using Multi-Sensor Datasets, Chandan Kumar

Dissertations, Master's Theses and Master's Reports

The sustainable exploration of mineral resources plays a significant role in the economic development of any nation. The lithological maps and surface mineral distribution can be vital baseline data to narrow down the geochemical and geophysical analysis potential areas. This study developed innovative spectral and Machine Learning (ML) methods for mineral and lithological classification. Multi-sensor datasets such as Airborne Visible/Infrared Imaging Spectrometer-Next Generation (AVIRIS-NG), Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER), Advanced Land Observing (ALOS) Phased Array type L-band Synthetic Aperture Radar (PALSAR), Sentinel-1, and Digital Elevation Model (DEM) were utilized. The study mapped the hydrothermal alteration minerals derived …

The Overexploitation Of The Valle De Querétaro Aquifer And Its Impact In Small Peri-Urban Communities, Querétaro, México, Kelsey Kirkland

Dissertations, Master's Theses and Master's Reports

The Valle de Querétaro Aquifer is the only viable local water source satisfying the domestic, agricultural, and industrial water needs of the Querétaro Valley. Severe depletion of the groundwater source has had significant consequences throughout the region, especially in the peri-urban communities of Santo Niño de Praga, Tlacote el Bajo, and La Palma. Historically, residents depended on aquifer-fed freshwater springs (known as Los Tajos) to meet their basic and productive needs. Spring production ceased between 10 and 15 ago and investigating the causes formed the basis of this research.

Analysis of available qualitative and quantitative data was used to identify …

Rain Generated Lahars Prior To The 2018 Catastrophic Eruption Of Fuego Volcano, Guatemala, Claudia Buondonno

Dissertations, Master's Theses and Master's Reports

Fuego volcano is one of the most active and hazardous volcanoes in the world. It is located in the northern part of the Central American Volcanic Arc in Guatemala and its activity can be characterized by long term, low-level background activity, and sporadic larger explosive eruptions. Its historical observations of eruptions date back to 1531, but it has been erupting vigorously since 2002 with major activity throughout 2018, producing three main eruptions in February, June and November. Its almost persistent activity generates major ashfalls, pyroclastic flows, lava flows; when heavy rains mobilize its deposits, they can form damaging lahars. Phenomena, …

Application Of Radar And Optical Satelite Imagery Data In Landslide Potential Mapping Of Shesheer Sub-Catchment In Iran, A. Imanian, M. H. Tangestani, A. Asadi

Michigan Tech Publications, Part 1

One of the serious dangers which threatens human communities especially those who are living in mountainous areas is the occurrence of landslides. Therefore, determination of the areas with potential for landslide events is very important for avoiding establishment of residential areas or industrial facilities. The aim of this study is to provide a landslide potential map in Sheshpeer sub-catchment, Iran, using an integration of remote sensing (RS) techniques and geographical information systems (GIS). Compared with the traditional approaches, these techniques are very fast, inexpensive and trustworthy in landslide mapping. For this purpose, we collected and produced seven data layers using …

Investigation Of Spectral Characteristics Of Carbonate Rocks – A Case Study On Posht Moleh Mount In Iran, A. Imanian, M. H. Tangestani, A. Asadi

Michigan Tech Publications, Part 1

Recent developments in the image processing approaches and the availability of multi and/or hyper spectral remote sensing data with high spectral, spatial and temporal resolutions have made remote sensing technique of great interest in investigations of geological sciences. One of the biggest advantage of the application of remote sensing in geology is recognizing the type of unknown rocks and minerals. In this study, an investigation on spectral features of carbonate rocks (i.e. calcite, dolomite, and dolomitized calcite) were done in terms of main absorptions, the reasons of those absorptions and comparison of these absorption with Johns Hopkins University (JHU) spectral …

A New Discrete Wavelength Backscattered Ultraviolet Algorithm For Consistent Volcanic So2 Retrievals From Multiple Satellite Missions, Bradford L. Fisher, Nickolay A. Krotkov, Pawan K. Bhartia, Can Li, Simon Carn, Eric Hughes, Peter J. T. Leonard

Michigan Tech Publications, Part 1

This paper describes a new discrete wavelength algorithm developed for retrieving volcanic sulfur dioxide (SO₂) vertical column density (VCD) from UV observing satellites. The Multi-Satellite SO₂ algorithm (MS_SO2) simultaneously retrieves column densities of sulfur dioxide, ozone, and Lambertian effective reflectivity (LER) and its spectral dependence. It is used operationally to process measurements from the heritage Total Ozone Mapping Spectrometer (TOMS) onboard NASA's Nimbus-7 satellite (N7/TOMS: 1978–1993) and from the current Earth Polychromatic Imaging Camera (EPIC) onboard Deep Space Climate Observatory (DSCOVR: 2015–ongoing) from the Earth–Sun Lagrange (L1) orbit. Results from MS_SO2 algorithm for several volcanic cases were …

Volcano Monitoring From Space Using High-Cadence Planet Cubesat Images Applied To Fuego Volcano, Guatemala, Anna Aldeghi, Simon Carn, Rudiger Escobar-Wolf, Gianluca Groppelli

Michigan Tech Publications, Part 1

Fuego volcano (Guatemala) is one of the most active and hazardous volcanoes in the world. Its persistent activity generates lava flows, pyroclastic density currents (PDCs), and lahars that threaten the surrounding areas and produce frequent morphological change. Fuego’s eruption deposits are often rapidly eroded or remobilized by heavy rains and its constant activity and inaccessible terrain makes ground-based assessment of recent eruptive deposits very challenging. Earth-orbiting satellites can provide unique observations of volcanoes during eruptive activity, when ground-based techniques may be too hazardous, and also during inter-eruptive phases, but have typically been hindered by relatively low spatial and temporal resolution. …

Analysis Of An Updated Paleointensity Database (Qpi‐Pint) For 65–200 Ma: Implications For The Long‐Term History Of Dipole Moment Through The Mesozoic, E. V. Kulakov, C. J. Sprain, P. V. Doubrovine, A. V. Smirnov, G. A. Paterson, L. Hawkins, L. Fairchild, E. J. Piispa, A. J. Biggin

Michigan Tech Publications, Part 1

The global paleointensity database for 65–200 Ma was analyzed using a modified suite of paleointensity quality criteria (Q_PI) such that the likely reliability of measurements is assessed objectively and as consistently as possible across the diverse data set. This interval was chosen because of dramatic extremes of geomagnetic polarity reversal frequency ranging from greater than 10 reversals per million years in the Jurassic hyperactivity period (155–171 Ma) to effectively zero during the Cretaceous Normal Superchron (CNS; 84–126 Ma). Various attempts to establish a relationship between the strength of Earth's magnetic field and the reversal frequency have been made …

Thermal Conduction Simulation Based On Reconstruction Digital Rocks With Respect To Fractures, Haiyuan Yang, Li Zhang, Ronghe Liu, Xianli Wen, Yongfei Yang, Roohollah Askari, Et. Al.

Michigan Tech Publications, Part 1

Effective thermal conductivity (ETC), as a necessary parameter in the thermal properties of rock, is affected by the pore structure and the thermal conduction conditions. To evaluate the effect of fractures and saturated fluids on sandstone’s thermal conductivity, we simulated thermal conduction along three orthogonal (X, Y, and Z) directions under air- and water-saturated conditions on reconstructed digital rocks with different fractures. The results show that the temperature distribution is separated by the fracture. The significant difference between the thermal conductivities of solid and fluid is the primary factor influencing the temperature distribution, and the thermal conduction mainly depends on …

Emittance Spectroscopy And Broadband Thermal Remote Sensing Applied To Phosphorite And Its Utility In Geoexploration: A Study In The Parts Of Rajasthan, India, Arindam Guha, Yasushi Yamaguchi, Snehamoy Chatterjee, Komai Rani, Kumranchat Vinod Kumar

Michigan Tech Publications, Part 1

The contrast in the emissivity spectra of phosphorite and associated carbonate rock can be used as a guide to delineate phosphorite within dolomite. The thermal emissivity spectrum of phosphorite is characterized by a strong doublet emissivity feature with their absorption minima at 9 µm and 9.5 µm; whereas, host rock dolomite has relatively subdued emissivity minima at ~9 µm. Using the contrast in the emissivity spectra of phosphorite and dolomite, data obtained by the thermal bands of Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) sensor were processed to delineate phosphorite within dolomite. A decorrelation stretched ASTER radiance composite could …

Evidence Of Instability In Previously-Mapped Landslides As Measured Using Gps, Optical, And Sar Data Between 2007 And 2017: A Case Study In The Portuguese Bend Landslide Complex, California, El Hachemi Y. Bouali, Thomas Oommen, Rüdiger Escobar-Wolf

Michigan Tech Publications, Part 1

Velocity dictates the destructive potential of a landslide. A combination of synthetic aperture radar (SAR), optical, and GPS data were used to maximize spatial and temporal coverage to monitor continuously-moving portions of the Portuguese Bend landslide complex on the Palos Verdes Peninsula in Southern California. Forty SAR images from the COSMO-SkyMed satellite, acquired between 19 July 2012 and 27 September 2014, were processed using Persistent Scatterer Interferometry (PSI). Eight optical images from the WorldView-2 satellite, acquired between 20 February 2011 and 16 February 2016, were processed using the Co-registration of Optically Sensed Images and Correlation (COSI-Corr) technique. Displacement measurements were …

Satellite-Based Assessment Of The August 2018 Flood In Parts Of Kerala, India, C. L. Vishnu, K. S. Sajinkumar, Thomas Oommen, R. A. Coffman, K. P. Thrivikramji, V. R. Rani, S. Keerthy

Michigan Tech Publications, Part 1

From 1 June to 29 August 2018, Kerala, a state in southwestern India, recorded 36% excess rainfall than normal levels, leading to widespread floods and landslides events and resulting in 445 deaths. In this study, satellite-based data were used to map the flood inundation in the districts of Thrissur, Ernakulam, Alappuzha, Idukki and Kottayam. Specifically, flood delineation was enabled with Sentinel-1A radar data of 21 August 2018 and was compared with an average pre-flood, water-cover map based on Modified Normalized Difference Water Index (MNDWI) that was developed using a January and February 2018 Sentinel-2A dataset. A 90% increase in water …

Assessment Of Post-Wildfire Debris Flow Occurrence Using Classifier Tree, Priscilla Addison, Thomas Oommen, Qiuying Sha

Michigan Tech Publications, Part 1

Besides the dangers of an actively burning wildfire, a plethora of other hazardous consequences can occur afterwards. Debris flows are among the most hazardous of these, being known to cause fatalities and extensive damage to infrastructure. Although debris flows are not exclusive to fire affected areas, a wildfire can increase a location’s susceptibility by stripping its protective covers like vegetation and introducing destabilizing factors such as ash filling soil pores to increase runoff potential. Due to the associated dangers, researchers are developing statistical models to isolate susceptible locations. Existing models predominantly employ the logistic regression algorithm; however, previous studies have …

Investigation Of Long Term Pillar Stability In An Abandoned Gypsum Mine Based On Dissolution Models, Mohammadhossein Sadeghiamirshahidi

Dissertations, Master's Theses and Master's Reports

The abandoned Domtar Mine in Grand Rapids, MI underlies Interstate I-196. The mine operated for over 140 years, mining a 30-meter-deep, four-meter thick gypsum seam. In 2000, the mine was abandoned with the removal of the mine’s dewatering pumps allowing groundwater to flood the mine. The initial flooding resulted in saturating the mine pillars along with some amount of pillar dissolution. Overtime, gypsum dissolution would cease as the water becomes saturated with the gypsum dissolution products. The mine, however, is located adjacent to the Grand River and has groundwater moving through the mine resulting in the potential for continuing dissolution. …

Revisiting The Paleomagnetism Of The Neoarchean Uauá Mafic Dyke Swarm, Brazil: Implications For Archean Supercratons, J. Salminen, E. P. Oliveira, Elisa J. Piispa, Aleksey Smirnov, R. I. F. Trindade

Michigan Tech Publications, Part 1

The original connections of Archean cratons are becoming traceable due to an increasing amount of paleomagnetic data and refined magmatic barcodes. The Uauá block of the northern São Francisco craton may represent a fragment of a major Archean craton. Here, we report new paleomagnetic data from the 2.62 Ga Uauá tholeiitic mafic dyke swarm of the Uauá block in the northern São Francisco craton, Eastern Brazil. Our paleomagnetic results confirm the earlier results for these units, but our interpretation differs. We suggest that the obtained characteristic remanent magnetization for the 2.62 Ga swarm is of primary origin, supported by a …

Climbing The Crustal Ladder: Magma Storage-Depth Evolution During A Volcanic Flare-Up, Guilherme A. R. Gualda, Darren M. Gravley, Michelle Connor, Brooke Hollmann, Ayla S. Pamukcu, Florence Begue, Mark S. Ghiorso, Chad Deering

Michigan Tech Publications, Part 1

Very large eruptions (>50 km³) and supereruptions (>450 km³) reveal Earth’s capacity to produce and store enormous quantities (>1000 km³) of crystal-poor, eruptible magma in the shallow crust. We explore the interplay between crustal evolution and volcanism during a volcanic flare-up in the Taupo Volcanic Zone (TVZ, New Zealand) using a combination of quartz-feldspar-melt equilibration pressures and time scales of quartz crystallization. Over the course of the flare-up, crystallization depths became progressively shallower, showing the gradual conditioning of the crust. Yet, quartz crystallization times were invariably very short (<100 years), demonstrating that very large reservoirs of eruptible magma were transient crustal features. We conclude that the dynamic nature of the TVZ crust favored magma eruption over storage. Episodic tapping of eruptible magmas likely prevented a supereruption. Instead, multiple very large bodies of eruptible magma were assembled and erupted in decadal time scales.

Landsliding Near Enguri Dam (Caucasus, Georgia) And Possible Seismoectonic Effects, Alessandro Tibaldi, Paolo Oppizzi, John S. Gierke, Thomas Oommen, Nino Tsereteli, Davit Odilavadze

Michigan Tech Publications, Part 1

The Enguri dam and water reservoir, nested in southwestern Caucasus (Republic of Georgia), are surrounded by steep mountain slopes. At a distance of 2.5 km from the dam, a mountain ridge along the reservoir is affected by active deformations with a double vergence. The western slope, directly facing the reservoir, has deformations that involve a subaerial area of 1.2 km2. The head scarp interests the main Jvari-Khaishi-Mestia road with offset of man-made features that indicate slip rates of 2-9 cm/y. Static, pseudostatic and Newmark numerical analyses, based on field and seismological data, suggest different unstable rock volumes basing on the …

Slalom: An All-Surface Snow Water Path Retrieval Algorithm For The Gpm Microwave Imager, Jean-François Rysman, Giulia Panegrossi, Paolo Sanò, Anna Cinzia Marra, Stefano Dietrich, Lisa Milani, Mark S. Kulie

Michigan Tech Publications, Part 1

This paper describes a new algorithm that is able to detect snowfall and retrieve the associated snow water path (SWP), for any surface type, using the Global Precipitation Measurement (GPM) Microwave Imager (GMI). The algorithm is tuned and evaluated against coincident observations of the Cloud Profiling Radar (CPR) onboard CloudSat. It is composed of three modules for (i) snowfall detection, (ii) supercooled droplet detection and (iii) SWP retrieval. This algorithm takes into account environmental conditions to retrieve SWP and does not rely on any surface classification scheme. The snowfall detection module is able to detect 83% of snowfall events including …

Sensitivity Of Atmospheric Aerosol Scavenging To Precipitation Intensity And Frequency In The Context Of Global Climate Change, Pei Hou, Shiliang Wu, Jessica L. Mccarty, Yang Gao

Michigan Tech Publications, Part 1

Wet deposition driven by precipitation is an important sink for atmospheric aerosols and soluble gases. We investigate the sensitivity of atmospheric aerosol lifetimes to precipitation intensity and frequency in the context of global climate change. Our sensitivity model simulations, through some simplified perturbations to precipitation in the GEOS-Chem model, show that the removal efficiency and hence the atmospheric lifetime of aerosols have significantly higher sensitivities to precipitation frequencies than to precipitation intensities, indicating that the same amount of precipitation may lead to different removal efficiencies of atmospheric aerosols. Combining the long-term trends of precipitation patterns for various regions with the …

The Ozone Monitoring Instrument: Overview Of 14 Years In Space, Pieternel F. Levelt, Joanna Joiner, Johanna Tamminen, J. Pepijn Veefkind, Pawan K. Bhartia, Simon Carn, Et. Al.

Michigan Tech Publications, Part 1

This overview paper highlights the successes of the Ozone Monitoring Instrument (OMI) on board the Aura satellite spanning a period of nearly 14 years. Data from OMI has been used in a wide range of applications and research resulting in many new findings. Due to its unprecedented spatial resolution, in combination with daily global coverage, OMI plays a unique role in measuring trace gases important for the ozone layer, air quality, and climate change. With the operational very fast delivery (VFD; direct readout) and near real-time (NRT) availability of the data, OMI also plays an important role in the development …