Submarine Groundwater Discharge Research Articles

The interplay of freshwater inputs and catchment geology in regulating seawater chemistry in Irish coastal areas

Open ocean carbonate chemistry is altered by the dissolution of atmospheric CO2 in seawater. Up to 40% of anthropogenic CO2 emissions have dissolved in the surface ocean since the pre-industrial era, driving changes in marine carbonate pools and promoting ocean acidification. Under open ocean conditions, Total Alkalinity (TA) generally relates with salinity and temperature due to the conservative nature of its constituents. In coastal areas, however, land-ocean interactions may greatly contribute to making TA far less predictable, since freshwater inputs can affect the chemistry of coastal water masses by increasing TA and Dissolved Inorganic Carbon (DIC) or, alternatively, lowering them through dilution. Here we analysed and compared coastal carbonate chemistry dynamics in selected coastal areas of Ireland, in order to assess whether rivers and their catchment geology can influence coastal water chemistry and to verify the occurrence of local ocean acidification processes. Data on TA, DIC, salinity, temperature, and nutrients (total nitrogen, phosphate, and silicate) collected during several surveys along the Irish coastline by The Marine Institute Foras Na Mara (MIFNM) from 2009 to 2018 were analysed against those available on the main watercourses feeding the selected coastal areas. The relationships among the different variables were scrutinized, also considering the geological characteristics of river catchments. Results showed a clear linkage between coastal water chemistry and local freshwater inputs, whereas changes in the geology of catchment areas exerted only a secondary influence. An in-depth scrutiny of the Galway Bay sector performed by comparing MIFNM data with those collected during an oceanographic cruise carried out in October 2017 further indicated remarkable effects of the interaction between river inputs and submarine groundwater discharges on seasonal variations in the water chemistry of the area. Future studies, integrating chemical data across the land-ocean boundary, may provide additional, useful information on the factors actually affecting the variation in water carbonates in coastal areas. Elucidating how land-ocean interactions influence the coastal carbonate system may ultimately provide a key for a more advanced understanding of the resilience of coastal areas to ocean acidification.

Characterizing 226Ra and its daughters in coastal zone groundwater of a typical human-activity affected bay: occurrence, safety, and source evaluation.

Due to their extremely toxic properties, 226Ra and it daughters (222Rn, 210Pb, and 210Po) in drinking groundwater require monitoring. Recent studies have reported exceptionally high levels of natural 210Po (up to 10,000Bq/m3), 226Ra, and 222Rn isotopes in groundwater. This study aims to provide background data on 226Ra and its daughter radionuclides in the typical agricultural-industrial Dongshan Bay (DSB) before the construction of Zhangzhou Nuclear Power Plant (Zhangzhou NPP). The measurement results indicate that no abnormally high activities of 210Po and 210Pb were detected in the investigated wells. Strong positive correlations between 210Pb and 210Po, as well as between 222Rn and 210Pb activities, suggest that the origins of 210Pb and 210Po in groundwater are strongly influenced by the decay of the parent radionuclides 222Rn and 210Pb, respectively. In the DSB coastal zone groundwater, significant deficiencies of 210Po relative to 210Pb and 210Pb relative to 222Rn were observed, providing further evidence that 210Po and 210Pb are also effectively scavenged due to their geochemical properties (specifically particle affinity) within the groundwater-aquifer system. A systematic comparison among all relevant water bodies in the DSB revealed that the activity concentrations of 210Pb and 210Po in groundwater were the highest, except for rainwater. Based on the evaluation of 210Pb sources, the results imply that submarine groundwater discharge (SGD) is an important pathway for transferring radionuclides (such as 210Pb) from land to the nearshore marine environment, even though the study area has a lower 210Pb background groundwater. By considering all the 210Pb's sources in the DSB, we found low 210Pb background groundwater discharge still needs to be taken into account for small-scale bays. This is because SGD was calculated to be one of the most important 210Pb sources in the bay during observation season. Regardless of whether the system is in a normal state or a nuclear accident emergency state, greater attention should be paid to the groundwater discharge of radionuclides into the ocean.

Dissolved carbon dynamics and exchange in a high permeability beach aquifer

The biogeochemical transformation of dissolved carbon in tidal beaches affects the coastal carbon cycle and budget. Yet, the influence of groundwater flow on the behavior of dissolved carbon and CO2 flux across the different interfaces in these dynamic beach systems remains poorly understood, especially in deep beach aquifers. In this study, we investigated the biogeochemical reactions and origins of dissolved carbon and quantitatively evaluated CO2 outgassing driven by submarine groundwater discharge (SGD) in a high permeability beach aquifer (9–13 m depth) located in the South China Sea shelf. Results revealed that both dissolved inorganic carbon (DIC) and total alkalinity (TA) exhibited non-conservative additions across the salinity gradient. The excess DIC and TA most likely resulted from organic carbon decomposition in the brackish and saline zones. Aerobic respiration (upper saline plume zone), anaerobic reactions, and tidal mixing were identified as the dominant processes controlling DIC and TA dynamics. The mean CO2 outgassing flux across the beach sediment-air interface was estimated to be 1.05 g m−2 d−1, significantly higher than the CO2 flux across the sea-air interface. The outgassing pattern of groundwater CO2 from beach aquifer appeared to be driven by tidal pumping. Considering the combined effects of SGD-derived dissolved carbon and nutrients, we suggest that SGD yielded a net CO2 flux, contributing to 45 % of the CO2 flux across the sea-air interface. SGD directly and indirectly delivered substantial CO2 from the beach aquifer to the atmosphere, approximately (4.89–9.21) × 103 tons of CO2 per year along the 76.2 km long sandy coast. These findings demonstrate that intertidal beach aquifers, as active biogeochemical reactors have the strong potential to intensify CO2 emissions to the atmosphere through both beach sediment-air and sea-air interfaces. This study provides a better understanding of the links between carbon biogeochemical cycles and hydrologic processes in the coastal zones.

Submarine groundwater discharge into a semi-enclosed coastal bay of the southern Baltic Sea: A multi-method approach

The present study aims to understand the impact of submarine groundwater discharge (SGD) on a coastal area with different lithology and degrees of SGD. Sampling campaigns took place in Puck Bay and the Gulf of Gdańsk, southern Baltic Sea encompassing years between 2009 and 2021. The methodological approach combined geophysical characterization of the surface sediments with detailed spatial and temporal (isotope) biogeochemical investigations of pore and surface waters, and was supported by nearshore groundwater and river surveys. Acoustic investigations identified areas of disturbance that may indicate zones of preferential SGD release. The composition of porewater and the differences in the bay's surface waters disclosed SGD as common phenomenon in the study area. Regional SGD was estimated through a radium mass balance. Local estimation of SGD, based on porewater profiles, revealed highest SGD fluxes at the sandy shoreline, but relatively low elemental fluxes. Though SGD was low at the muddy sites corresponding elemental fluxes of nutrients and dissolved carbon exceeded those determined at the sandy sites due to intense diagenesis in the top sediments. SGD appears to be sourced from different freshwater endmembers; however, diagenesis in surface sediments substantially modified the composition of the mixed solutions that are finally discharged to coastal waters. Overall, this study provides a better understanding of the SGD dynamics in the region by a multi-approach and emphasizes the need to understand the processes occurring at the sediment-water interface when estimating SGD.

Nitrogen and carbon cycling and relationships to radium behavior in porewater and surface water: Insight from a dry year sampling in a hypersaline estuary

Biogeochemical transformations within highly saline subterranean estuaries (STE) dramatically affect solute cycling, resulting in submarine groundwater discharge (SGD) with distinct chemical signatures. The study hypothesizes that biogeochemical processes within hypersaline bay porewaters (PW) simultaneously affect nitrogen, carbon, and radium cycling. We measured radium isotopes (226Ra, 224Ra, and 223Ra), nutrients (dissolved inorganic nitrogen [DIN: NH4+ + NO2− + NO3−], HPO42− [DIP], HSiO3− [DSi], dissolved organic carbon [DOC]), total alkalinity (TA), dissolved inorganic carbon (DIC), stable isotopes, and major cations in PW and surface water (SW) of Baffin Bay, a well-mixed, semi-enclosed estuary along the semiarid northwestern Gulf of Mexico coast, over three seasons in a characteristically dry year. This study's findings show a concurrent increase in NH4+, DIP, DSi, and TA/DIC with reduced metal species (e.g., Mn and Fe) and Ra during the hot and dry seasons, particularly in PW, under increasingly reducing conditions. Principal component analyses (PCA) suggest these increases are primarily driven by dissimilatory nitrate/nitrite reduction to ammonium (DNRA) and dissolution of lithogenic particles and biogenic CaCO3, modulated by organic matter degradation or remineralization. While more significant terrestrial groundwater inputs may contribute to solutes and Ra supply in the STE, the biogeochemically induced variability in solute concentrations in PW primarily drives larger SGD-derived fluxes, particularly notable in hot months. During a typically dry year, these fluxes, estimated as the average of 226Ra and 223Ra mass balance models (e.g., July/November fluxes in Mmol∙d−1: 0.093/0.092 of NO3−; 0.2/0.02 of NO2−; 72/16 of NH4+; 72.2/18 of DIN; 1.5/0.2 of HPO42−; 20/9 of HSiO3−; 42/37 of DOC; 503/399 of TA; 582/431 of DIC) are orders of magnitude (∼4 for DIN and DIC, ∼3 for DIP, DSi, and DOC, and ∼2 for TA) greater than surface runoff inputs. These substantial SGD inputs likely sustain phytoplankton growth and potentially fuel harmful algal blooms while countering estuarine acidification.

Time-Lapse Geophysical Measurements for Monitoring Coastal Groundwater Dynamics in an Unconfined Aquifer.

The coastal zone, which is the interface between land and sea, is hydrodynamically very active due to the complex interactions of various hydrological controls and variable-density fluids. These forces vary over time, resulting in a state of dynamic equilibrium in the system. The major hydrological processes in coastal aquifer systems are salt water intrusion and submarine groundwater discharge, which are interdependent. Monitoring these complex processes is crucial for sustainable coastal zone management but poses a significant research challenge. In this study, we demonstrate the effectiveness of non-invasive geophysical techniques, specifically the time-lapse electrical resistivity imaging method, in conjunction with groundwater monitoring, for monitoring coastal groundwater dynamics in an unconfined aquifer at varying time scales and hydrogeological settings present at formerly glaciated sites worldwide. We generated two-dimensional baseline salt water intrusion maps for the test site, located on the coast of Rhode Island, USA. The time-lapse electrical resistivity survey method enables the rapid estimation of fresh groundwater discharge. Our approach offers insight into the mechanisms and seasonably variable salt water-freshwater interactions in unconfined heterogeneous aquifers. Although the results are site-specific, their implications are broad and may stimulate other studies related to sea to land pollution (sea water intrusion) and land to sea pollution (groundwater discharge) in heterogeneous coastal aquifer settings.

Coastline changes under natural and anthropogenic drivers in a macro-tidal estuary between 2000-2020

Coastline changes in estuarine areas can result from a combination of natural processes such as erosion, sedimentation, and sea-level rise, as well as human activities, including urbanization and infrastructure development. These changes have the potential to affect the local environment, including submarine groundwater discharge, wetlands, and navigation routes. The Yalu River Estuary (YRE), situated on the border between China and North Korea, has been experiencing significant changes in its coastline over recent years. This study aims to investigate the coastline dynamics in the YRE from 2000 to 2020. The study employs Landsat 5/7/8 satellite data and proposes a modified Normalized Difference Water Index (NDWI) to accurately delineate the coast boundary, particularly in areas with extensive tidal flats like the YRE. The research findings indicate that from 2010 to 2020, significant changes occurred in the YRE shoreline, with erosion being the dominant trend. Human activities and alternations in hydrological conditions are important factors affecting the YRE coastline changes, contributing to the formation of distinctive spatiotemporal patterns. An extreme flooding event in the year 2010 also altered the inner estuarine coastline in the YRE, indicating the impact of strong natural drivers. Findings from this study provide a comprehensive understanding of the evolving coastal environment, considering natural and anthropogenic drivers, and highlight the importance of continuous monitoring in a region of ecological and geopolitical significance.

Groundwater dynamics in intertidal zones of beach aquifers during salt-fingering flow

At aquifer-sea interfaces, submarine groundwater discharge (SGD) and related chemical loads are closely related to groundwater environment in the intertidal zone of coastal aquifers. Although unstable salinity distribution in the intertidal zone associated with salt-fingering flow has been known as a natural phenomenon for some coastal aquifers, a quantitative understanding of groundwater dynamics in the intertidal zone during salt-fingering flow is still lacking. This study investigated the movement features of salt fingers and the response laws of freshwater discharge zone (FDZ) in the intertidal zone, considering the fixed and seasonal inland conditions. The results revealed that during salt-fingering flow, the finger was repeatedly generated and detached from the intertidal zone. Horizontally, the finger moved toward the sea, and the finger speed first increased at the development stage, then remained unchanged at the separation stage, and finally decreased at the mergence stage. The movement speed of salt fingers was closely related to ambient freshwater flow in the FDZ. In the vertical direction, the finger first permeated downward with a decreasing speed, and then the upward discharge of fresh groundwater in the FDZ region with an increased flow velocity pushed the finger upwards with an increasing speed. In contrast to the results in a stable intertidal regime, the width of FDZ at the aquifer-sea interface fluctuated over time, and the variation trends possessed a well agreement with the flux of fresh groundwater discharge. This suggested that the width of FDZ at the aquifer-sea interface can be used as an indicator of fresh groundwater discharge, even in an unstable regime, to provide a basis for field monitoring. The intensity of salt-fingering flow varied over time under the condition of seasonal freshwater inflow. The mean duration of salt fingers decreased slightly with the increase of seasonal freshwater periods due to a decrease in the interface perturbation for the intertidal environment. These findings lay the foundation for gaining a better understanding of the potential implications of salt-fingering flow on SGD and related solute transports.

Transport of nitrate-containing groundwater to coastal areas through buried tunnel valleys, Denmark

Nitrogen impact on the aquatic environment, including coastal areas, is too high in many countries worldwide, particularly in areas with intensive agriculture. Efficient mitigation initiatives demand that important pathways and the fate of nitrate in the hydrological cycle are known. In this study, we focus on groundwater nitrate contamination in two near-shore catchment areas in north-west Denmark. Groundwater in the area is mainly located in buried tunnel valleys, which are subsurface structures eroded by meltwater during Pleistocene glaciations in former glaciated areas. Groundwater samples from the aquifers inside the buried valleys reveal the presence of up to 120 mg/l nitrate down to 10 m below sea level and about 1 km down from the stream outlet towards the coast. We interpret the complex tunnel-valley infill to be responsible for the spatial heterogeneity of the groundwater geochemistry, where sandy geological windows create localised hydraulic pathways and complex redox structures. Groundwater and stream water chemistry in the study area clearly demonstrate the role of groundwater in nitrate transport within the catchment as well as the direct pathway to the coast bypassing the stream and riverine systems. Our results show that the buried tunnel valleys potentially contribute to submarine groundwater discharge and therefore could be responsible for a hitherto unaccounted input of nitrogen to the marine environment.

Permafrost and groundwater interaction: current state and future perspective

This study reviews the available and published knowledge of the interactions between permafrost and groundwater. In its content, the paper focuses mainly on groundwater recharge and discharge in the Arctic and the Qinghai-Tibet Plateau. The study revealed that the geochemical composition of groundwater is site-specific and varies significantly within the depth of the aquifers reflecting the water-rock interactions and related geological history. All reviewed studies clearly indicated that the permafrost thaw causes an increase in groundwater discharge on land. Furthermore, progressing climate warming is likely to accelerate permafrost degradation and thus enhance hydrological connectivity due to increased subpermafrost groundwater flow through talik channels and higher suprapermafrost groundwater flow. In the case of submarine groundwater discharge (SGD), permafrost thaw can either reinforce or reduce SGD, depending on how much pressure changes affecting the aquifers will be caused by the loss of permafrost. Finally, this comprehensive assessment allowed also for identifying the lack of long-term and interdisciplinary in situ measurements that could be used in sophisticated computational simulations characterizing the current status and predicting groundwater flow and permafrost dynamics in the future warmer climate.

Coral Reveals a Transition in Sources of Seawater Vanadium From Submarine Groundwater Discharge to Upwelling in the Northern South China Sea

AbstractIn recent decades, vanadium (V) levels in the environment have dramatically increased and will continue to rise alarmingly in the future. Conventional methods for assessing oceanic V typically only involved the collection of seawater and sediments, which significantly limited our understanding of the dynamics of oceanic V and its relationship to environment and climate. Here, the interannual temporal variations in surface seawater V levels in the northern South China Sea, were reconstructed utilizing a 137‐year record of V/Ca ratios from Porites lutea coral from Xiaodonghai Bay off Sanya. We observed that anthropogenic activities had a limited impact on the V levels of surface seawater in the study area, whereas V variability exhibited significant disparities across different time periods. Precipitation‐related submarine groundwater discharge (SGD) recharge was the primary driver of V variation prior to 1990, with tropical cyclones (TCs) playing a detectable role in total precipitation during this period. Conversely, changes in surface seawater V levels after 1990 were predominantly linked to increasingly intense upwelling driven by summer winds. The recharge processes in unconfined groundwater systems within the study area might generate favorable oxic and alkaline conditions, which facilitated mineral dissolution and/or desorption reactions that elevated V concentrations along the groundwater flow path, which were recorded in coral skeletal V/Ca ratios. This study confirmed the seawater V was dominated by SGD and upwelling near Sanya, providing essential evidence for comprehending the geochemical behavior of V and the impact of ocean‐atmospheric climate interactions on long‐term variations in seawater V levels.

Comparing nearshore and embayment scale assessments of submarine groundwater discharge: Significance of offshore groundwater discharge as a nutrient pathway

Submarine groundwater discharge (SGD) can influence biogeochemical cycles in coastal seas by delivering nutrients from the seafloor. Comparison between the nearshore and embayment scale assessments of SGD against river water discharge would be crucial for understanding biogeochemical impacts on the coastal seas because the discharge pattern (non-point or point pathway) is different. Here, we quantified SGD contribution to rivers in nutrient budgets at two scales within a coastal embayment (Obama Bay, Japan) by mass balance models of radon and radium isotopes. We then compared the SGD contribution between the two scales by the meta-analysis for regional data sets conducted in nearshore and embayment scales. The estimated SGD rates in the nearshore and embayment scales in the bay were 7.8 cm d−1 and 20.0 cm d−1, indicating that offshore SGD was more significant than nearshore. The ratios of nutrient fluxes derived from SGD to rivers (SGD:River) in the nearshore scale were 1.7 for dissolved inorganic nitrogen (DIN), 3.0 for phosphorus (DIP), and 0.5 for silica (DSi), while those in the embayment scale increased to 10.4 for DIN, 18.5 for DIP, and 3.9 for DSi. This result indicates that SGD-derived nutrients become more important at larger spatial scales. Meta-analysis revealed that the difference in the contribution of SGD to rivers was affected by the seafloor size and there was no significant difference in SGD rates between nearshore and embayment scale studies. However, our regional study shows the site-specific pattern that SGD rates in the embayment scale were higher than those in the nearshore scale. Overall, we clarified that SGD can be a crucial nutrient pathway for coastal embayments regardless of the spatial scales and contribute to coastal nutrient biogeochemistry in more offshore areas.

Intercomparison of marine tracer and catchment-based submarine groundwater discharge estimates for a karst aquifer

Submarine groundwater discharge (SGD) is a major source of solutes to the ocean. Consequently, the accuracy of estimates of the magnitude and seasonal variability of SGD is critical for the management of coastal marine ecosystems. However, very few intercomparison studies between different conceptual approaches have been made to assess the reliability and consistency of SGD flow estimates while simultaneously gauging the applicability of underlying models and budget assumptions in real-world scenarios. To date, intercomparison efforts have focused on homogeneous aquifers, and less so on highly heterogeneous aquifers such as karst systems. Here we compare the results of tracer-based (salt mass balance, tidal prism and radon mass balances) and catchment-based methods (water balance, semi-distributed karst network model) to assess SGD at a data-rich coastal location. All approaches tested here reveal a positive relationship between SGD fluxes and groundwater level measured close to the preferential groundwater flow pathways in the phreatic karst aquifer (e.g. water-saturated conduits or fracture network). We use this finding to derive site-specific relationships between groundwater level and SGD flow rates. We derive these relationships from least square fitting of Darcian and Non-Darcian models and regression analysis. The relationships validate the findings derived independently from hydraulic/hydrogeological modelling of the complex karst networks. By combining various estimates made with distinct assumptions, this method produces robust estimates of the SGD seasonal variability, essential to develop optimal management strategies of SGD resources in coastal karstic regions. Where direct measurements of SGD are not possible, monitoring of coastal groundwater level close to preferential phreatic groundwater pathways can be thus used as a useful indicator to compare past groundwater flow measurements and derive site-specific relationships, particularly in settings where large flow rate changes are present such as karst aquifers.

Anthropogenic As pollution mediated by submarine groundwater discharge in a marine ranch

Arsenic (As) pollution, is a global problem, threatening human health and ecological security, especially in the bay environment with dense population and human activities. Among potential pathways of As into the bay, submarine groundwater discharge (SGD) has not received adequate attention due to its invisibility. We determined As and 222Rn activity concentrations in different water mass. Spatial variation of dissolved As concentration in the groundwater was large and attributed to the adjacent local industries. By combining 222Rn mass balance modeling with As concentrations measured, the SGD-derived As fluxes was conservatively estimated to be 1310 kg As d−1 and 5880 kg As d−1 in the dry and wet seasons, respectively. The migration of arsenic may be enhanced by rainfall and dissolved carbon. The amount of SGD derived As input to the bay was greater than the total combined As input from river discharge, atmosphere, sewage drainage, and diffusion from sediment.

Effect of submarine groundwater discharge on nutrient distribution and eutrophication in Liaodong Bay, China

Driven by the anthropogenic activities associated with coastal settlements, eutrophication has become a global issue. Submarine groundwater discharge (SGD) is a significant continuous pathway for transporting nutrients from land to coastal waters, but its influence on eutrophication in Liaodong Bay (LDB) has received limited attention. In this study, radium isotopes and nutrient data from coastal waters were analyzed to evaluate the SGD flux and its implications for potential eutrophication in LDB. We found that the mean concentrations of dissolved inorganic nitrogen (DIN), phosphorous (DIP), and silicate (DSi) in groundwater were higher than those of seawater and river water. Based on 223Ra and 228Ra mass balance models, the SGD fluxes were estimated to be (0.53–2.03) × 109 m3/d, of which the fresh SGD accounted for 4 %–15 %. SGD is a vital invisible source of nutrients, contributing more than 79 % of the total inputs of DIN, DIP, and DSi into LDB. With high DIN/DIP ratios (average=85.8) and large nutrient inputs, SGD may significantly drive the phosphorus limitation and eutrophication in LDB. This study shows that SGD-derived nutrient fluxes should be considered in the assessment of water eutrophication for the formulation of future environmental management protocols in coastal systems.

Submarine Groundwater Discharge Exceeds River Inputs as a Source of Nutrients to the Great Barrier Reef.

Rivers are often assumed to be the main source of nutrients triggering eutrophication in the Great Barrier Reef (GBR). However, existing nutrient budgets suggest a major missing source of nitrogen and phosphorus sustaining primary production. Here, we used radium isotopes to resolve submarine groundwater discharge (SGD)-derived, shelf-scale nutrient inputs to the GBR. The total SGD was ∼10-15 times greater than average river inputs, with nearshore groundwater discharge accounting for ∼30% of this. Total SGD accounted for >30% of all known dissolved inorganic N and >60% of inorganic P inputs and exceeded regional river inputs. However, SGD was only a small proportion of the nutrients necessary to sustain primary productivity, suggesting that internal recycling processes still dominate the nutrient budget. With millions of dollars spent managing surface water nutrient inputs to reef systems globally, we argue for a shift in the focus of management to safeguard reefs from the impacts of excess nutrients.

Sustained upward groundwater discharge through salt marsh tidal creeks

AbstractSalt marshes can export considerable nutrients and carbon to the ocean through submarine groundwater discharge (SGD). However, the complicated SGD processes in salt marshes remain poorly understood. Here, we first report the phenomenon of numerous highly saline artesian springs found in a salt marsh system of East China. Multiple methods including time‐series thermal monitoring, isotope signatures, and high‐resolution electrical resistivity tomography were combined to determine their origin and trajectory. Strong evidence suggests that these springs keep discharging even during high tide and represent a long‐term re‐distribution process of the ancient marine water trapped in the unconfined aquifer. This new pattern of spring‐derived groundwater flow indicates a hidden SGD pathway and has significant implications for studies concerning SGD‐derived fluxes in similar multi‐aquifer‐aquitard coastal systems.

Complex Effects of Tides on Coastal Groundwater Revealed by High‐Resolution Integrated Flow Modeling

AbstractRiver deltas typically have high population density and support a wide range of intensive and prosperous socioeconomic activities. The hydrological processes in these regions are complex, primarily due to the interactions among the river, aquifer, and sea. However, a systematic and quantitative elaboration of the river‐aquifer‐sea interactions is still lacking. Here, we developed an integrated hydrological flow model for the Pearl River Delta (PRD), which contains the world's largest urban area in both size and population, to gain a deeper understanding of the complexities in the river‐aquifer‐sea interactions. The model performance was validated and cross‐checked via observations at gauging stations and independent remote‐sensing products (e.g., soil moisture, evapotranspiration, and total water storage anomalies). Based on the 10‐year simulation results (2004–2013), the major findings of this study are as follows: (a) accurate representation of the tidal effect is important not only for simulating short‐term flow dynamics but also for capturing the characteristics of long‐term hydrological fluxes and states; (b) the flow‐model‐computed average groundwater discharge rate per unit length of the coastline for the PRD is 3.01 m3/d/m, which is comparable with those derived from water budget approaches but 1–2 orders of magnitude lower than the total submarine groundwater discharge (SGD) estimated by using isotope tracer‐based methods; (c) the temporal variation of SGD is controlled by tidal forcing on an hourly time scale, but by terrestrial hydrological processes on monthly and annual time scales; and (d) an integrated hydrological flow model can be used to identify distinct and large subsurface zones sensitive to tidal fluctuations, quantifying the pivotal role of ocean tides in shaping the coastal groundwater system. This study represents a first step in using an integrated hydrological model to explore river‐aquifer‐sea interactions and their effects on the regional groundwater system simultaneously driven by meteorological and tidal forcings.

Microplastic pollution in the groundwater under a bedrock island in the South China sea

Groundwater is the only freshwater resource on islands. Research on microplastic pollution in groundwater on islands is scarce. This study is the first to explore microplastic pollution in the groundwater under a bedrock island (Dawanshan Island) located in the South China Sea. The influence of hydrogeological factors on the distribution, source, and ageing features of microplastics in the groundwater were investigated. Despite the small scale of industrial and agricultural activities on the island, the amount of microplastics in the groundwater ranged from 34 to 64 particles/L, with over 80% of the microplastics being polyester fibres with diameters smaller than 2 mm, which is comparable to those in coastal cities. These microplastics were originated from inland plastic usage, rather than from the surrounding sea, which was confirmed by the lack of seawater intrusion on the island. Owing to the low permeability of granite, microplastics were mainly distributed in the water of the loose layer of porous sediment, and their quantity decreased with depth. In addition, the abundance of microplastics in pore groundwater increased with an increase in the velocity of groundwater flow. The severity of microplastic pollution in the groundwater increased with an increase and decrease in the content of total dissolved solids and dissolved oxygen, respectively. The microplastics originated from plastic waste disposed of on the island, rather than from seawater intrusion. Also, through groundwater infiltration into exposed soil at recharge areas, artificial wells at residential areas, and water exchange with surface water at valley areas. Microplastics buried in the groundwater aged faster along the migration path of the groundwater. These microplastics threaten the safety of people and plants on the island through exposure resulting from the extraction of groundwater for irrigation, while they endanger marine life through submarine groundwater discharge.

Development of a high-resolution marine ecosystem model for predicting the combined impacts of ocean acidification and deoxygenation

An approach was developed to help evaluate and predict the combined effects of ocean acidification and deoxygenation on calcifying organisms along the coast of Japan. The Coastal and Regional Ocean COmmunity (CROCO) modeling system was set up to couple the Regional Ocean Modeling System (ROMS) to the Pelagic Interaction Scheme for Carbon and Ecosystem Studies (PISCES) biogeochemical model and used to reproduce physical and biochemical processes in the area around Miyako Bay, Iwate Prefecture, Japan. Future scenario cases were also set up, which used initial and boundary conditions based on Future Ocean Regional Projection (FORP) simulations. Present day simulations were able to reproduce the general features of observed physical and biochemical parameters, except for some rapid decreases in salinity, pH and aragonite saturation state (Ωarag). This suggests that more local factors which have not been introduced into the model, such as submarine groundwater discharge, may be involved, or that river inputs may be underestimated. Results of the future projections suggest a significant impact of global warming and ocean acidification on calcifying organisms for the worst case of climate change under the Representative Concentration Pathway (RCP) 8.5 scenario. In particular, it is feared that values of Ωaragwould approach the critical level for calcifying organisms (Ωarag< 1.1) throughout the year, under which decreased larval shell lengths and malformation have been observed experimentally for the locally grownHaliotis discus hannai(Ezo Abalone) species. However, these findings may not be true for a different coastal locality, and this study highlights and continues to stress the importance of developing model setups capable of incorporating both regional and local factors affecting ocean acidification and deoxygenation.