Pub Date : 2023-10-23DOI: 10.3390/geosciences13100318
Evangelos Machairas, Emmanouil A. Varouchakis
Mining activities are prolific worldwide in light of the perpetual production of metal. The high need for metal materials in human life necessitates the development of mining operations, especially in places characterized as being highly enriched in metal ions. After the separation of beneficial and non-beneficial materials, industrial enrichment mechanisms take place to increase metal output. These mechanisms, known as metallurgical procedures, produce a vast volume of mining/metallurgical waste (MMW) at final disposal sites. MMW’s composition usually includes metal filings in low-pH site conditions. Thus, the environmental pollution hazard is high unless sustainable methods are implemented to reduce both heavy and toxic metals’ concentration in MMW at every disposal site. The scope of this review is to determine how cost–benefit analysis (CBA) and risk assessment (RA) could contribute positively to (a) the environmental effect of MMW reduction, (b) decreasing the environmental rehabilitation cost, and (c) research into economically sustainable methods of recovering metal from MMW.
{"title":"Cost–Benefit Analysis and Risk Assessment for Mining Activities in Terms of Circular Economy and Their Environmental Impact","authors":"Evangelos Machairas, Emmanouil A. Varouchakis","doi":"10.3390/geosciences13100318","DOIUrl":"https://doi.org/10.3390/geosciences13100318","url":null,"abstract":"Mining activities are prolific worldwide in light of the perpetual production of metal. The high need for metal materials in human life necessitates the development of mining operations, especially in places characterized as being highly enriched in metal ions. After the separation of beneficial and non-beneficial materials, industrial enrichment mechanisms take place to increase metal output. These mechanisms, known as metallurgical procedures, produce a vast volume of mining/metallurgical waste (MMW) at final disposal sites. MMW’s composition usually includes metal filings in low-pH site conditions. Thus, the environmental pollution hazard is high unless sustainable methods are implemented to reduce both heavy and toxic metals’ concentration in MMW at every disposal site. The scope of this review is to determine how cost–benefit analysis (CBA) and risk assessment (RA) could contribute positively to (a) the environmental effect of MMW reduction, (b) decreasing the environmental rehabilitation cost, and (c) research into economically sustainable methods of recovering metal from MMW.","PeriodicalId":38189,"journal":{"name":"Geosciences (Switzerland)","volume":"94 1","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"2023-10-23","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"135366867","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2023-10-20DOI: 10.3390/geosciences13100317
Ernesto Ausilio, Maria Giovanna Durante, Paolo Zimmaro
Geosynthetic-reinforced soil structures are often used to support shallow foundations of various infrastructure systems including bridges, railways, and highways. When such infrastructures are located in seismic areas, their performance is linked to the seismic bearing capacity of the foundation. Various approaches can be used to calculate this quantity such as analytical solutions and advanced numerical models. Building upon a robust upper bound limit analysis, we created a database comprising 732 samples. The database was then used to train and test a model based on a random forest machine learning algorithm. The trained random forest model was used to develop a publicly available web application that can be readily used by researchers and practitioners. The model considers the following input factors: (1) the ratio of the distance of the foundation from the edge and the width of the foundation (D/B), (2) the slope angle (β), (3) the horizontal seismic intensity coefficient (kh), and (4) the dimensionless geosynthetic factor, which accounts for the tensile strength of the geosynthetic. Leveraging the model developed in this study, we show that the most important features to predict the seismic bearing capacity of strip footings positioned on the crest of geosynthetic-reinforced soil structures are D/B and kh.
{"title":"On the Potential of Using Random Forest Models to Estimate the Seismic Bearing Capacity of Strip Footings Positioned on the Crest of Geosynthetic-Reinforced Soil Structures","authors":"Ernesto Ausilio, Maria Giovanna Durante, Paolo Zimmaro","doi":"10.3390/geosciences13100317","DOIUrl":"https://doi.org/10.3390/geosciences13100317","url":null,"abstract":"Geosynthetic-reinforced soil structures are often used to support shallow foundations of various infrastructure systems including bridges, railways, and highways. When such infrastructures are located in seismic areas, their performance is linked to the seismic bearing capacity of the foundation. Various approaches can be used to calculate this quantity such as analytical solutions and advanced numerical models. Building upon a robust upper bound limit analysis, we created a database comprising 732 samples. The database was then used to train and test a model based on a random forest machine learning algorithm. The trained random forest model was used to develop a publicly available web application that can be readily used by researchers and practitioners. The model considers the following input factors: (1) the ratio of the distance of the foundation from the edge and the width of the foundation (D/B), (2) the slope angle (β), (3) the horizontal seismic intensity coefficient (kh), and (4) the dimensionless geosynthetic factor, which accounts for the tensile strength of the geosynthetic. Leveraging the model developed in this study, we show that the most important features to predict the seismic bearing capacity of strip footings positioned on the crest of geosynthetic-reinforced soil structures are D/B and kh.","PeriodicalId":38189,"journal":{"name":"Geosciences (Switzerland)","volume":"2 1","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"2023-10-20","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"135569444","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2023-10-19DOI: 10.3390/geosciences13100316
Evgeny Chuvilin, Dinara Davletshina, Boris Bukhanov, Sergey Grebenkin, Elena Pankratova
High-latitude permafrost, including hydrate-bearing frozen ground, changes its properties in response to natural climate change and to impacts from petroleum production. Of special interest is the behavior of thermal conductivity, one of the key parameters that control the thermal processes in permafrost containing gas hydrate accumulations. Thermal conductivity variations under pressure and temperature changes were studied in the laboratory through physical modeling using sand sampled from gas-bearing permafrost of the Yamal Peninsula (northern West Siberia, Russia). When gas pressure drops to below equilibrium at a constant negative temperature (about −6 °C), the thermal conductivity of the samples first becomes a few percent to 10% lower as a result of cracking and then increases as pore gas hydrate dissociates and converts to water and then to ice. The range of thermal conductivity variations has several controls: pore gas pressure, hydrate saturation, rate of hydrate dissociation, and amount of additionally formed pore ice. In general, hydrate dissociation can cause up to 20% thermal conductivity decrease in frozen hydrate-bearing sand. As the samples are heated to positive temperatures, their thermal conductivity decreases by a magnitude depending on residual contents of pore gas hydrate and ice: the decrease reaches ~30% at 20–40% hydrate saturation. The thermal conductivity decrease in hydrate-free saline frozen sand is proportional to the salinity and can become ~40% lower at a salinity of 0.14%. The behavior of thermal conductivity in frozen hydrate-bearing sediments under a pressure drop below the equilibrium and a temperature increase to above 0 °C is explained in a model of pore space changes based on the experimental results.
{"title":"Thermal Conductivity Variations in Frozen Hydrate-Bearing Sand upon Heating and Dissociation of Pore Gas Hydrate","authors":"Evgeny Chuvilin, Dinara Davletshina, Boris Bukhanov, Sergey Grebenkin, Elena Pankratova","doi":"10.3390/geosciences13100316","DOIUrl":"https://doi.org/10.3390/geosciences13100316","url":null,"abstract":"High-latitude permafrost, including hydrate-bearing frozen ground, changes its properties in response to natural climate change and to impacts from petroleum production. Of special interest is the behavior of thermal conductivity, one of the key parameters that control the thermal processes in permafrost containing gas hydrate accumulations. Thermal conductivity variations under pressure and temperature changes were studied in the laboratory through physical modeling using sand sampled from gas-bearing permafrost of the Yamal Peninsula (northern West Siberia, Russia). When gas pressure drops to below equilibrium at a constant negative temperature (about −6 °C), the thermal conductivity of the samples first becomes a few percent to 10% lower as a result of cracking and then increases as pore gas hydrate dissociates and converts to water and then to ice. The range of thermal conductivity variations has several controls: pore gas pressure, hydrate saturation, rate of hydrate dissociation, and amount of additionally formed pore ice. In general, hydrate dissociation can cause up to 20% thermal conductivity decrease in frozen hydrate-bearing sand. As the samples are heated to positive temperatures, their thermal conductivity decreases by a magnitude depending on residual contents of pore gas hydrate and ice: the decrease reaches ~30% at 20–40% hydrate saturation. The thermal conductivity decrease in hydrate-free saline frozen sand is proportional to the salinity and can become ~40% lower at a salinity of 0.14%. The behavior of thermal conductivity in frozen hydrate-bearing sediments under a pressure drop below the equilibrium and a temperature increase to above 0 °C is explained in a model of pore space changes based on the experimental results.","PeriodicalId":38189,"journal":{"name":"Geosciences (Switzerland)","volume":"7 1","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"2023-10-19","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"135779337","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2023-10-19DOI: 10.3390/geosciences13100315
Miguel Ángel Álvarez-Vázquez, Elena De Uña-Álvarez, Alexandra María Ramírez-Pérez, Esther de Blas, Ricardo Prego
Sediment compositions and enrichment patterns are investigated in an urban reach of a bedrock river, the Miño River passing through Ourense City, Spain. This study focuses on the trace element distribution in different fractions to gain insights into trace element enrichment. To assess enrichment, a context-specific approach was employed, based on the mean, the standard deviation of the estimated background, and the empirical rule, avoiding the pitfalls of general and arbitrary thresholds. Notably, the <0.063 mm and <2 mm fractions showed differential accumulation patterns. Both fractions serve to detect enrichments that can be indicative of contamination, but they measure different things, the maturity of sediments and postdepositional processes being key factors in understanding the sediment composition and enrichments. These findings also highlight the role of rock cavities, particularly those hosting permanent deposits, as traps for trace elements and their potential significance in assessing environmental enrichment. This work contributes to understanding sediment compositions and enrichment dynamics in bedrock rivers. It also underscores the significance of considering site-specific approaches for enrichment assessment and the necessity for further research to unravel the mechanisms driving differential accumulation within distinct depositional environments.
{"title":"Distinctive Accumulation Patterns of Trace Elements in Sediments of Bedrock Rivers (Miño River, NW Iberian Peninsula)","authors":"Miguel Ángel Álvarez-Vázquez, Elena De Uña-Álvarez, Alexandra María Ramírez-Pérez, Esther de Blas, Ricardo Prego","doi":"10.3390/geosciences13100315","DOIUrl":"https://doi.org/10.3390/geosciences13100315","url":null,"abstract":"Sediment compositions and enrichment patterns are investigated in an urban reach of a bedrock river, the Miño River passing through Ourense City, Spain. This study focuses on the trace element distribution in different fractions to gain insights into trace element enrichment. To assess enrichment, a context-specific approach was employed, based on the mean, the standard deviation of the estimated background, and the empirical rule, avoiding the pitfalls of general and arbitrary thresholds. Notably, the <0.063 mm and <2 mm fractions showed differential accumulation patterns. Both fractions serve to detect enrichments that can be indicative of contamination, but they measure different things, the maturity of sediments and postdepositional processes being key factors in understanding the sediment composition and enrichments. These findings also highlight the role of rock cavities, particularly those hosting permanent deposits, as traps for trace elements and their potential significance in assessing environmental enrichment. This work contributes to understanding sediment compositions and enrichment dynamics in bedrock rivers. It also underscores the significance of considering site-specific approaches for enrichment assessment and the necessity for further research to unravel the mechanisms driving differential accumulation within distinct depositional environments.","PeriodicalId":38189,"journal":{"name":"Geosciences (Switzerland)","volume":"78 1","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"2023-10-19","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"135730719","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2023-10-18DOI: 10.3390/geosciences13100314
Ken’ichi Koshimizu, Taro Uchida
After multiple simultaneous landslides caused by heavy rainfall, expanding landslides continue to occur for a certain duration. Evaluation of the influencing period of sediment yield due to expanding landslides is vital for comprehensive sediment management of the basin. In this study, we investigated a region with a low frequency of heavy rainfall that has not received its due level of attention until now. Consequently, the transition of expanding landslides depends on the transition of the number of remaining landslides, based on the difference in the frequency of heavy rainfall. Furthermore, the transition of expanding landslides depends on the maximum daily rainfall after the landslides. These findings indicate that “the number of remaining landslides” and “maximum daily rainfall after a landslide” are related factors that determine the period during which expanding landslides frequently occur. An estimation formula based on elapsed time was developed to calculate the number of remaining landslides. An empirical formula for the number of expanding landslides was obtained by multiplying the function of the daily maximum rainfall after the landslide by the estimation formula for the number of remaining landslides. The developed empirical formula can be used effectively for evaluation during periods when rainfall-induced landslides are subject to subsequent expansion.
{"title":"Time-Series Variation of Landslide Expansion in Areas with a Low Frequency of Heavy Rainfall","authors":"Ken’ichi Koshimizu, Taro Uchida","doi":"10.3390/geosciences13100314","DOIUrl":"https://doi.org/10.3390/geosciences13100314","url":null,"abstract":"After multiple simultaneous landslides caused by heavy rainfall, expanding landslides continue to occur for a certain duration. Evaluation of the influencing period of sediment yield due to expanding landslides is vital for comprehensive sediment management of the basin. In this study, we investigated a region with a low frequency of heavy rainfall that has not received its due level of attention until now. Consequently, the transition of expanding landslides depends on the transition of the number of remaining landslides, based on the difference in the frequency of heavy rainfall. Furthermore, the transition of expanding landslides depends on the maximum daily rainfall after the landslides. These findings indicate that “the number of remaining landslides” and “maximum daily rainfall after a landslide” are related factors that determine the period during which expanding landslides frequently occur. An estimation formula based on elapsed time was developed to calculate the number of remaining landslides. An empirical formula for the number of expanding landslides was obtained by multiplying the function of the daily maximum rainfall after the landslide by the estimation formula for the number of remaining landslides. The developed empirical formula can be used effectively for evaluation during periods when rainfall-induced landslides are subject to subsequent expansion.","PeriodicalId":38189,"journal":{"name":"Geosciences (Switzerland)","volume":"554 1","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"2023-10-18","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"135883400","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2023-10-17DOI: 10.3390/geosciences13100311
Zuhal Akyurek, Semih Kuter, Çağrı H. Karaman, Berkay Akpınar
Understanding the distribution, patterns, and characteristics of snowfall and snow cover within a given region over extended periods is important. Snow climatology provides valuable insights into the seasonal and long-term variations in snowfall, helping researchers and meteorologists understand the impacts of climate change on snow accumulation, melt rates, and snowmelt runoff. In this study, in order to understand the spatial and temporal variation in snow cover in Turkey, the temporal and spatial dynamics of snow cover in the country were analyzed during the latest and longest period from 1970 to 2022 using ERA5-Land reanalysis product. It is aimed (1) to show snow-covered area (SCA), snow duration, and snow depth trends over the country; (2) to examine the altitudinal difference of snow phenology response to climate change; and (3) to evaluate the Snow Cover Frequency Maps from MODIS Snow Cover Products with the reanalysis snow depth data. It is found that the “false snow” mapping problem still exists in the MOD10C1_CGF Snow Cover Frequency maps over Turkey, especially in the melting period. We found that an increasing trend of 0.4 °C/decade and snow duration have a decreasing trend due to the early melting between 1970 and 2022. This trend is even more noticeable at elevations below 2000 m. Another important finding is the decreasing trend in snow duration at altitudes below 500 m, indicating a shift from snow to rain for precipitation types.
{"title":"Understanding the Snow Cover Climatology over Turkey from ERA5-Land Reanalysis Data and MODIS Snow Cover Frequency Product","authors":"Zuhal Akyurek, Semih Kuter, Çağrı H. Karaman, Berkay Akpınar","doi":"10.3390/geosciences13100311","DOIUrl":"https://doi.org/10.3390/geosciences13100311","url":null,"abstract":"Understanding the distribution, patterns, and characteristics of snowfall and snow cover within a given region over extended periods is important. Snow climatology provides valuable insights into the seasonal and long-term variations in snowfall, helping researchers and meteorologists understand the impacts of climate change on snow accumulation, melt rates, and snowmelt runoff. In this study, in order to understand the spatial and temporal variation in snow cover in Turkey, the temporal and spatial dynamics of snow cover in the country were analyzed during the latest and longest period from 1970 to 2022 using ERA5-Land reanalysis product. It is aimed (1) to show snow-covered area (SCA), snow duration, and snow depth trends over the country; (2) to examine the altitudinal difference of snow phenology response to climate change; and (3) to evaluate the Snow Cover Frequency Maps from MODIS Snow Cover Products with the reanalysis snow depth data. It is found that the “false snow” mapping problem still exists in the MOD10C1_CGF Snow Cover Frequency maps over Turkey, especially in the melting period. We found that an increasing trend of 0.4 °C/decade and snow duration have a decreasing trend due to the early melting between 1970 and 2022. This trend is even more noticeable at elevations below 2000 m. Another important finding is the decreasing trend in snow duration at altitudes below 500 m, indicating a shift from snow to rain for precipitation types.","PeriodicalId":38189,"journal":{"name":"Geosciences (Switzerland)","volume":"15 1","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"2023-10-17","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"136033142","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2023-10-17DOI: 10.3390/geosciences13100313
Luigi Massaro, Giovanni Forte, Melania De Falco, Antonio Santo
The morphoevolution of coastal areas is due to the interactions of multiple continental and marine processes that define a highly dynamic environment. These processes can occur as rapid catastrophic events (e.g., landslides, storms, and coastal land use) or as slower continuous processes (i.e., wave, tidal, and current actions), creating a multi-hazard scenario. Maronti Bay (Ischia Island, Southern Italy) can be classified as a pocket beach that represents an important tourist and environmental area for the island, although it has been historically affected by slope instability, sea cliff recession, and coastal erosion. In this study, the historical morphoevolution of the shoreline was analysed by means of a dataset of aerial photographs and cartographic information available in the literature over a 25-year period. Furthermore, the role of cliff recession and its impact on the beach was also explored, as in recent years, the stability condition of the area was worsened by the occurrence of a remarkable landslide in 2019. The latter was reactivated following a cloudburst on the 26th of November 2022 that affected the whole Island and was analysed with the Dem of Difference technique. It provided an estimate of the mobilised volumes and showed how the erosion and deposition areas were distributed and modified by wave action. The insights from this research can be valuable in developing mitigation strategies and protective measures to safeguard the surrounding environment and ensure the safety of residents and tourists in this multi-hazard environment.
{"title":"Geomorphological Evolution of Volcanic Cliffs in Coastal Areas: The Case of Maronti Bay (Ischia Island)","authors":"Luigi Massaro, Giovanni Forte, Melania De Falco, Antonio Santo","doi":"10.3390/geosciences13100313","DOIUrl":"https://doi.org/10.3390/geosciences13100313","url":null,"abstract":"The morphoevolution of coastal areas is due to the interactions of multiple continental and marine processes that define a highly dynamic environment. These processes can occur as rapid catastrophic events (e.g., landslides, storms, and coastal land use) or as slower continuous processes (i.e., wave, tidal, and current actions), creating a multi-hazard scenario. Maronti Bay (Ischia Island, Southern Italy) can be classified as a pocket beach that represents an important tourist and environmental area for the island, although it has been historically affected by slope instability, sea cliff recession, and coastal erosion. In this study, the historical morphoevolution of the shoreline was analysed by means of a dataset of aerial photographs and cartographic information available in the literature over a 25-year period. Furthermore, the role of cliff recession and its impact on the beach was also explored, as in recent years, the stability condition of the area was worsened by the occurrence of a remarkable landslide in 2019. The latter was reactivated following a cloudburst on the 26th of November 2022 that affected the whole Island and was analysed with the Dem of Difference technique. It provided an estimate of the mobilised volumes and showed how the erosion and deposition areas were distributed and modified by wave action. The insights from this research can be valuable in developing mitigation strategies and protective measures to safeguard the surrounding environment and ensure the safety of residents and tourists in this multi-hazard environment.","PeriodicalId":38189,"journal":{"name":"Geosciences (Switzerland)","volume":"32 1","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"2023-10-17","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"135994695","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2023-10-17DOI: 10.3390/geosciences13100312
De Nyago Tafen, Barbara Kutchko, Mehrdad Massoudi
Cement fractures represent preferential leakage pathways in abandoned wells upon exposure to a CO2-rich fluid. Understanding fracture alteration resulting from geochemical reactions is critical for assessing well integrity in CO2 storage. This paper describes a mathematical model used to investigate the physical and the chemical changes in cement properties when CO2-saturated water is injected into a wellbore. This study examines the flow of a solution of CO2-saturated water in a two-dimensional fractured cement. In this approach, a micro-continuum equation based on the Darcy–Brinkman–Stokes (DBS) equation is used as the momentum balance equation; in addition, reactive transport equations are used to study the coupled processes of reactant transport and geochemical reactions, and the model for cement porosity alteration and fracture enhancement. This paper focuses on the effects of cement porosity, fracture aperture size, and surface roughness. Mineral dissolution and precipitation mechanisms are also considered. Our simulations show that smaller initial fracture apertures tend to a high mineral precipitation self-sealing. However, a complete sealing of the fracture is not observed due to the continuous flow of CO2-saturated water. The calcite precipitation mechanism of a rough fracture (random zigzag shape) differs from that of a smooth/flat fracture surface.
{"title":"On the Flow of CO2-Saturated Water in a Cement Fracture","authors":"De Nyago Tafen, Barbara Kutchko, Mehrdad Massoudi","doi":"10.3390/geosciences13100312","DOIUrl":"https://doi.org/10.3390/geosciences13100312","url":null,"abstract":"Cement fractures represent preferential leakage pathways in abandoned wells upon exposure to a CO2-rich fluid. Understanding fracture alteration resulting from geochemical reactions is critical for assessing well integrity in CO2 storage. This paper describes a mathematical model used to investigate the physical and the chemical changes in cement properties when CO2-saturated water is injected into a wellbore. This study examines the flow of a solution of CO2-saturated water in a two-dimensional fractured cement. In this approach, a micro-continuum equation based on the Darcy–Brinkman–Stokes (DBS) equation is used as the momentum balance equation; in addition, reactive transport equations are used to study the coupled processes of reactant transport and geochemical reactions, and the model for cement porosity alteration and fracture enhancement. This paper focuses on the effects of cement porosity, fracture aperture size, and surface roughness. Mineral dissolution and precipitation mechanisms are also considered. Our simulations show that smaller initial fracture apertures tend to a high mineral precipitation self-sealing. However, a complete sealing of the fracture is not observed due to the continuous flow of CO2-saturated water. The calcite precipitation mechanism of a rough fracture (random zigzag shape) differs from that of a smooth/flat fracture surface.","PeriodicalId":38189,"journal":{"name":"Geosciences (Switzerland)","volume":"29 1","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"2023-10-17","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"136032512","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2023-10-16DOI: 10.3390/geosciences13100310
Hanna N. Vaidya, Robert D. Breininger, Marisela Madrid, Steven Lazarus, Nezamoddin N. Kachouie
Within the last century, the global sea level has risen between 16 and 21 cm and will likely accelerate into the future. Projections from the Intergovernmental Panel on Climate Change (IPCC) show the global mean sea level (GMSL) rise may increase to up to 1 m (1000 mm) by 2100. The primary cause of the sea level rise can be attributed to climate change through the thermal expansion of seawater and the recession of glaciers from melting. Because of the complexity of the climate and environmental systems, it is very difficult to accurately predict the increase in sea level. The latest estimate of GMSL rise is about 3 mm/year, but as GMSL is a global measure, it may not represent local sea level changes. It is essential to obtain tailored estimates of sea level rise in coastline Florida, as the state is strongly impacted by the global sea level rise. The goal of this study is to model the sea level in coastal Florida using climate factors. Hence, water temperature, water salinity, sea surface height anomalies (SSHA), and El Niño southern oscillation (ENSO) 3.4 index were considered to predict coastal Florida sea level. The sea level changes across coastal Florida were modeled using both multiple regression as a broadly used parametric model and the generalized additive model (GAM), which is a nonparametric method. The local rates and variances of sea surface height anomalies (SSHA) were analyzed and compared to regional and global measurements. The identified optimal model to explain and predict sea level was a GAM with the year, global and regional (adjacent basins) SSHA, local water temperature and salinity, and ENSO as predictors. All predictors including global SSHA, regional SSHA, water temperature, water salinity, ENSO, and the year were identified to have a positive impact on the sea level and can help to explain the variations in the sea level in coastal Florida. Particularly, the global and regional SSHA and the year are important factors to predict sea level changes.
{"title":"Generalized Additive Models for Predicting Sea Level Rise in Coastal Florida","authors":"Hanna N. Vaidya, Robert D. Breininger, Marisela Madrid, Steven Lazarus, Nezamoddin N. Kachouie","doi":"10.3390/geosciences13100310","DOIUrl":"https://doi.org/10.3390/geosciences13100310","url":null,"abstract":"Within the last century, the global sea level has risen between 16 and 21 cm and will likely accelerate into the future. Projections from the Intergovernmental Panel on Climate Change (IPCC) show the global mean sea level (GMSL) rise may increase to up to 1 m (1000 mm) by 2100. The primary cause of the sea level rise can be attributed to climate change through the thermal expansion of seawater and the recession of glaciers from melting. Because of the complexity of the climate and environmental systems, it is very difficult to accurately predict the increase in sea level. The latest estimate of GMSL rise is about 3 mm/year, but as GMSL is a global measure, it may not represent local sea level changes. It is essential to obtain tailored estimates of sea level rise in coastline Florida, as the state is strongly impacted by the global sea level rise. The goal of this study is to model the sea level in coastal Florida using climate factors. Hence, water temperature, water salinity, sea surface height anomalies (SSHA), and El Niño southern oscillation (ENSO) 3.4 index were considered to predict coastal Florida sea level. The sea level changes across coastal Florida were modeled using both multiple regression as a broadly used parametric model and the generalized additive model (GAM), which is a nonparametric method. The local rates and variances of sea surface height anomalies (SSHA) were analyzed and compared to regional and global measurements. The identified optimal model to explain and predict sea level was a GAM with the year, global and regional (adjacent basins) SSHA, local water temperature and salinity, and ENSO as predictors. All predictors including global SSHA, regional SSHA, water temperature, water salinity, ENSO, and the year were identified to have a positive impact on the sea level and can help to explain the variations in the sea level in coastal Florida. Particularly, the global and regional SSHA and the year are important factors to predict sea level changes.","PeriodicalId":38189,"journal":{"name":"Geosciences (Switzerland)","volume":"57 1","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"2023-10-16","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"136114223","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2023-10-16DOI: 10.3390/geosciences13100309
Lyudmila N. Morozova, Sergey G. Skublov, Dmitry R. Zozulya, Pavel A. Serov, Elena S. Borisenko, Anna N. Solovjova, Alexandra K. Gavrilchik
Beryl is both an accessory and a rock-forming mineral in pegmatites that contain beryl, making it a major source of Be. Beryl-bearing pegmatites of the Shongui deposit, located in the Kola province of the Northeastern Fennoscandian Shield, hold beryl with a yellowish-greenish color. An investigation into the chemical composition of this beryl from pegmatite dike No. 7 has been performed for the first time via the secondary ion mass spectrometry (SIMS) technique, and the chemical composition of the beryl-bearing pegmatites has been analyzed for the first time by the inductively coupled plasma mass spectrometry (ICP-MS) method. These pegmatites have high concentrations (ppm) of Be (11.8), Li (30.9), Rb (482), Nb (50.3), Ta (14.6), Cs (66.8), and Mn (283) and low concentrations of Sr, Y, Ba, rare earth elements (REE), Zr, and Th. In the Shongui pegmatite field, concentrations of Be, Li, Rb, Cs, Nb, Ta, and Mn increase from barren to beryl-bearing pegmatites, whereas concentrations of Ba, Sr, Y, and REE decline. Rb/Ba, Rb/Sr, and Zr/Hf ratios, showing the fractionation degree, change from the barren to beryl-bearing pegmatites: Rb/Ba and Rb/Sr increase from 111 and 0.46 to 1365 and 8.06, respectively, and Zr/Hf decreases from 18.9 to 14.5. The chemical composition of beryl from the Shongui deposit is unique. This mineral has a concentration of 25,300 ppm of alkalis (Li, Cs, K, Rb, Na) and the average Li, Ce, and Na content is 4430, 5000, and 15,400 ppm, respectively. According to its chemical composition, the Shongui beryl belongs to the Li-Cs-Na type, a type that is not recognized in the available classifications. It is supposed that this beryl was mainly crystallized in the magmatic stage rather than in any hydrothermal and metasomatic stages. Two beryl groups have been distinguished in beryl-bearing pegmatite dike No. 7: beryl from the intermediate zone (Brl-I) and beryl from the core zone (Brl-II). These beryls are concluded to have crystallized in the following order: Brl-I and then Brl-II. Compared with Brl-I, Brl-II is depleted in Cs, Na, Cl, and H2O and is enriched in Fe and Mn. The Fe/Mn ratio varies from 9.18 to 16.50 in these beryls and their yellowish-greenish shades are thought to be driven by a large amount of Fe compared to Mn.
{"title":"Li-Cs-Na-Rich Beryl from Beryl-Bearing Pegmatite Dike No. 7 of the Shongui Deposit, Kola Province, Russia","authors":"Lyudmila N. Morozova, Sergey G. Skublov, Dmitry R. Zozulya, Pavel A. Serov, Elena S. Borisenko, Anna N. Solovjova, Alexandra K. Gavrilchik","doi":"10.3390/geosciences13100309","DOIUrl":"https://doi.org/10.3390/geosciences13100309","url":null,"abstract":"Beryl is both an accessory and a rock-forming mineral in pegmatites that contain beryl, making it a major source of Be. Beryl-bearing pegmatites of the Shongui deposit, located in the Kola province of the Northeastern Fennoscandian Shield, hold beryl with a yellowish-greenish color. An investigation into the chemical composition of this beryl from pegmatite dike No. 7 has been performed for the first time via the secondary ion mass spectrometry (SIMS) technique, and the chemical composition of the beryl-bearing pegmatites has been analyzed for the first time by the inductively coupled plasma mass spectrometry (ICP-MS) method. These pegmatites have high concentrations (ppm) of Be (11.8), Li (30.9), Rb (482), Nb (50.3), Ta (14.6), Cs (66.8), and Mn (283) and low concentrations of Sr, Y, Ba, rare earth elements (REE), Zr, and Th. In the Shongui pegmatite field, concentrations of Be, Li, Rb, Cs, Nb, Ta, and Mn increase from barren to beryl-bearing pegmatites, whereas concentrations of Ba, Sr, Y, and REE decline. Rb/Ba, Rb/Sr, and Zr/Hf ratios, showing the fractionation degree, change from the barren to beryl-bearing pegmatites: Rb/Ba and Rb/Sr increase from 111 and 0.46 to 1365 and 8.06, respectively, and Zr/Hf decreases from 18.9 to 14.5. The chemical composition of beryl from the Shongui deposit is unique. This mineral has a concentration of 25,300 ppm of alkalis (Li, Cs, K, Rb, Na) and the average Li, Ce, and Na content is 4430, 5000, and 15,400 ppm, respectively. According to its chemical composition, the Shongui beryl belongs to the Li-Cs-Na type, a type that is not recognized in the available classifications. It is supposed that this beryl was mainly crystallized in the magmatic stage rather than in any hydrothermal and metasomatic stages. Two beryl groups have been distinguished in beryl-bearing pegmatite dike No. 7: beryl from the intermediate zone (Brl-I) and beryl from the core zone (Brl-II). These beryls are concluded to have crystallized in the following order: Brl-I and then Brl-II. Compared with Brl-I, Brl-II is depleted in Cs, Na, Cl, and H2O and is enriched in Fe and Mn. The Fe/Mn ratio varies from 9.18 to 16.50 in these beryls and their yellowish-greenish shades are thought to be driven by a large amount of Fe compared to Mn.","PeriodicalId":38189,"journal":{"name":"Geosciences (Switzerland)","volume":"28 1","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"2023-10-16","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"136114345","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}