Modeling the Effects of Tidal Energy Extraction on Estuarine Hydrodynamics in a Stratified Estuary

A 3-D coastal ocean model with a tidal turbine module was used in this paper to study the effects of tidal energy extraction on temperature and salinity stratification and density-driven two-layer estuarine circulation. Numerical experiments with various turbine array configurations were carried out to investigate the changes in tidally averaged temperature, salinity, and velocity profiles in an idealized stratified estuary that connects to coastal water through a narrow tidal channel. The model was driven by tides, river inflow, and sea surface heat flux. To represent the realistic size of commercial tidal farms, model simulations were conducted based on a small percentage (less than 10 %) of the total number of turbines that would generate the maximum extractable energy in the system. Model results show that extraction of tidal in-stream energy will increase the vertical mixing and decrease the stratification in the estuary. Installation of in-stream tidal farm will cause a phase lag in tidal wave, which leads to large differences in tidal currents between baseline and tidal farm conditions. Extraction of tidal energy in an estuarine system has stronger impact on the tidally averaged salinity, temperature, and velocity in the surface layer than the bottom layer even though the turbine hub height is close to the bottom. Finally, model results also indicate that extraction of tidal energy weakens the two-layer estuarine circulation, especially during neap tides when tidal mixing is weakest and energy extraction is smallest.     

钱塘江河口段盐水入侵的时空变化及预测模型

杭州城市供水85%取自钱塘江河口段,取水水质在枯水大潮期都不同程度地受到盐水入侵的威胁,分析钱塘江河口盐水入侵时空变化及研制二维数值预测模型对保障城市供水安全十分必要。根据钱塘江河口段实测水文氯度资料,分析了强潮作用下盐水入侵的时空变化特征;据此构建考虑斜压作用的二维水流、盐度输移的耦合数学模型,计算格式采用守恒性较好的有限体积法;在模型验证的基础上,数值分析了径流和潮汐对钱塘江河口段盐水入侵的影响,结果表明河口段的盐水入侵明显地受径流和潮汐的影响,据此可通过增大上游新安江水库的下泄流量抑制盐水入侵上溯以减小取水口氯度及超标时间,确保用水安全。     

Numerical modeling of tide-surge interaction along Orissa coast of India

The Orissa coast of India is one of the most vulnerable regions of extreme sea levels associated with severe tropical cyclones. There was extensive loss of life and property due to the October 1999 super cyclone, which devastated large part of the Orissa coast. The shallow nature of the head bay, presence of a large number of deltas formed by major rivers of Orissa such as Mahanadi and Dhamra, and high tidal range are responsible for storm surge flooding in the region. Specifically, rising and falling tidal phases influence the height, duration, and arrival time of peak surge along the coast. The objective of the present study is to evaluate the tide-surge interaction during the 1999 Orissa cyclone by using nonlinear vertically integrated numerical models. The pure tidal solution for the head bay region of the Bay of Bengal provides the initial condition for the fine resolution nested grid Orissa model. However, the feedback from the Orissa model does not affect the head bay model as the study provides a one-way interaction. Numerical experiments are performed to study the tide-surge interaction by considering various relative phases of the tidal waves with the surge-wave produced by 1999 Orissa cyclone. The comparison, although utilizing only the limited estimates of tidal data, appears adequate to assert that the principal features are reproduced correctly.     

Patterns of estuarine use by juvenile English sole (Parophrys vetulus) and Dungeness crab (Cancer magister)

Extensive trawl surveys were conducted in two large estuaries (Grays Harbor and Willapa Bay) on the Washington coast during 1983–1987, and in adjacent areas of the open coast. These surveys have shown that both English sole and Dungeness crab rely heavily on these estuaries as nursery areas, although the pattern of utilization differs substantially. Juvenile migration patterns can show substantial interannual variability and can only be delineated by concurrent surveys in both coastal and estuarine areas, conducted over a period of several years. English sole eggs and Dungeness crab larvae are released in coastal waters. Larvae of both species transform to the benthic stage in both coastal and estuarine areas, but most English sole eventually migrate into the estuaries during the first year of life, even if initial settlement is along the open coast. By the time English sole have attained a length of 55 mm (TL), most of them are found in estuaries. English sole begin emigrating from the estuaries at about 75 mm, and few remain there during the second year of life. In contrast, Dungeness crab appear to remain in the area of initial settlement throughout the first year of life. Growth is substantially faster in estuaries where 0+ crab reach a mean size of about 40 mm carapace width (CW) by September, with those off the coast are only about 14 mm CW. Juveniles remain in the area of settlement over their first winter but, in contrast to English sole, most coastal 1+ crab immigrate to estuaries to join siblings that settled there the previous year. By September of the second year, crab at about 100 mm CW emigrate to the open coast where they reach maturity. Advantages to juvenile stages that reside in estuaries are discussed in terms of accelerated growth at higher temperatures and potentially greater food supplies than found nearshore along the coast.     

Modeling diagnosis of suspended sediment transport in tidal estuarine system

A three-dimensional, time-dependent hydrodynamic and suspended sediment transport model was performed and applied to the Danshuei River estuarine system and adjacent coastal sea in northern Taiwan. The model was validated with observed time-series salinity in 2001, and with salinity and suspended sediment distributions in 2002. The predicted results quantitatively agreed with the measured data. A local turbidity maximum was found in the bottom water of the Kuan-Du station. The validated model then was conducted with no salinity gradient, no sediment supply from the sediment bed, wind stress, and different freshwater discharges from upstream boundaries to comprehend the influences on suspended sediment dynamics in the Danshuei River estuarine system. The results reveal that concentrations of the turbidity maximum simulated without salinity gradient are higher than those of the turbidity maximum simulated with salinity gradient at the Kuan-Du station. Without bottom resuspension process, the estuarine turbidity maximum zone at the Kuan-Du station vanishes. This suggests that bottom sediment resuspension is a very important sediment source to the formation of estuarine turbidity maximum. The wind stress with northeast and southwest directions may contribute to decrease the suspended sediment concentration. When the freshwater discharges increase at the upstream boundaries, the limits of salt intrusion pushes downriver toward river mouth. Suspended sediment concentrations increase at the upriver reaches in the Danshuei River to Tahan Stream, while decrease at Kuan-Du station.     

上海复合极端风暴洪水淹没模拟

河口三角洲和沿海城市面临着台风、暴雨、高潮位和上游下泄洪水叠加的“四碰头”复合极端风暴洪水的严重威胁。构建了大气-海洋-陆地相耦合的一体化数值模拟系统,实现了上海市“风”“暴”“潮”“洪”多灾种复合情景的极端洪涝淹没模拟,并验证了耦合方法的有效性,为复合风暴洪水的一体化模拟提供了一套可行的数值模拟方法。在9711台风影响下,模拟了1998年堤防升级改造后淹没面积(水深>0.2m)比改造前减少了62%,表明沿海沿江堤防设施建设在上海市防台防汛中起着关键性的作用。复合极端风暴洪水的有效模拟可为财产保险和未来市政规划提供参考。     

The influence of episodic events on transport of striped bass eggs to the estuarine turbidity maximum nursery area

The estuarine turbidity maximum (ETM) is an important nursery area for anadromous fish where early-life stages can be retained in high prey concentrations and favorable salinities. Episodic freshwater flow and wind events could influence the transport of striped bass (Morone saxatilis) eggs to the ETM. This hypothesis was evaluated with regression analysis of observational data and with a coupled biological-physical model of a semi-idealized upper Chesapeake Bay driven by observed wind and freshwater flow. A particle-tracking model was constructed within a numerical circulation model (Princeton Ocean Model) to simulate the transport of fish eggs in a 3-dimensional flow field. Particles with the sinking speed of striped bass eggs were released up-estuary of the salt front in both 2-d event-scale and 60-d seasonal-scale scenarios. In event scenarios, egg-like particles with observed specific gravities (densities) of striped bass eggs were transported to the optimum ETM nursery area after 2 d, the striped bass egg-stage duration. Wind events and pulses in river discharge decreased the number of egg-like particles transported to the ETM area by 20.9% and 13.2%, respectively, compared to nonevent conditions. In seasonal scenarios, particle delivery to the ETM depended upon the timing of the release of egg-like particles. The number of particles transported to the ETM area decreased when particles were released before and during wind and river pulse events. Particle delivery to the ETM area was enhanced when the salt front was moving up-estuary after river pulse events and as base river flow receded over the spawning season. Model results suggest that the timing of striped bass spawning in relation to pulsed events may have a negative (before or during events) or positive (after river flow events) effect on egg transport. Spawning after river flow events may promote early-stage survival by taking advantage of improved transport, enhanced turbidity refuge, and elevated prey production that may occur after river pulse events. In multiple regression analysis of observed data, mean spring freshwater flow rates and the number of pulsed freshwater flow events during the striped bass spawning season explained 71% of the variability in striped bass juvenile abundance in upper Chesapeake Bay from 1986 to 2002. Positive parameter estimates for these effects support the hypothesis that pulsed freshwater flow events, coupled with spawning after the events, may enhance striped bass early-stage survival. Results suggest that episodic events may have an important role in controlling fish recruitment.     

Toward a unified theory of tidally-averaged estuarine salinity structure

Equations are developed for the tidally-averaged, width-averaged estuarine salinity and circulation in a rectangular estuary. Width and depth may vary along the length of the channel, as may coefficients of vertical turbulent mixing and along channel diffusion. The system is reduced to a single first-order, nonlinear, ordinary differential equation governing the section-averaged salinity. A technique for specifying the ocean boundary condition is given, and solutions are found by numerical integration. Under different assumptions for the diffusion it is possible to reproduce the few existing analytical solutions, in particular the Hansen and Rattray (1965) Central Regime solution, and Chatwin's (1976) solution. The mathematical framework allows easy comparison of the results of different channel geometries and mixing coefficients. Of particular interest is the along-channel distribution of the diffusive fraction of up-estuary salt flux. It is shown that the Hansen and Rattray solution is always diffusion-dominated near the mouth. A theory is presented for estimating the diffusion coefficient within a tidal excursion of the mouth. It is shown that the resulting rapid along-channel increase of diffusion may explain some observed patterns of salinity structure: a decrease in both stratification and along-channel salinity gradient near the mouth. The theory is applied to the Delaware Estuary and Northern San Francisco Bay, and shows reasonable agreement with observed sensitivities of salt intrusion distance to river flow.     

珠江口磨刀门枯季水文特征及河口动力过程

根据磨刀门2003年12月9~15日的大、中潮同步水文观测资料,分析了磨刀门枯季的潮汐、潮流、余流、悬移质含沙量、盐度等水文特征,并对枯季河口动力过程,如咸淡水混合、河口射流等进行了初步研究。在枯季由于径流较弱,潮流成为主要动力。表层由于受径流和风的影响余流基本上沿河道走向向下游,中层以下有稳定的向上的余流存在。枯季磨刀门含沙量较小(<1 kg/m3),盐度在平面上和垂向上均有一定变化。磨刀门枯季咸淡水混合类型为缓混合型,各站盐度分层参数均在0.01~1.0。从实测流速的分布情况来看,河口下层有反向的水流,存在明显的因密度差而形成的密度环流。根据枯季实测资料计算所得的密度弗劳德数,磨刀门枯季以浮力射流为主。     

A Bay–River Coupled Model for Storm Surge Prediction Along the Andhra Coast of India

The paper describes a two-dimensional bay–river coupled numerical model for storm surges along the Andhra coast of India. The effect of the Krishna and Godavari rivers on the surge development is analysed. A comparative study of the surge generated by a tropical cyclone with and without the inclusion of rivers is done in detail. Three cyclones that struck the Andhra coast in November 1977, May 1990 and November 1996 were used for the simulation studies. It is found that the idealized model without a river overestimates the sea-level elevation as compared to a more realistic bay–river coupled model. The temporal variation of surge values at the mouth of the rivers is also studied for all three cyclone cases. It is found that the effect of the presence of rivers depends on the strength of the cyclone, its point of landfall and the location of the rivers with respect to the landfall point.     

Temporal and spatial variations in estuarine sedimentation

Sedimentation in estuaries is normally viewed as a continuous process controlled by inputs of fluvial and marine materials, which migrate in response to tidal circulation. Short-term examination yields reproducible and measurable rates of sedimentation; measurement of annual cycles reveals that sediment may be eroded from a site normally subject to deposition. Likewise annual migration patterns of turbidity maximum zones permits deposition of fine sediments at different parts of the estuary according to season. Study of long-term variations in the positions of sand banks shows that advances of the tidal flat margins may occur suddenly in response to catastrophic events; events that result in mid-channel banks becoming attached to the outer margins of the tidal flat. An example is also given of a channel-margin sand bank known to have repeatedly migrated round an elliptical path, taking approximately 25 yr to regain its original position. Long-term migration of bay head beach bars driven into the Eden Estuary and moving slowly landward are shown to create areas of temporary shelter for salt marshes, which are then over-run as the bar continues toward the head of the estuary. Several examples of anthropogenic influences on estuarine sediments are noted, with major deposits of coarse gravels attributed to known shipping trade within the Tay Estuary. Finally attention is drawn to the effects of deposition of coal mine wastes on the Fife coast: the erosion and redeposition of the wastes to form coastal platforms, their removal once mining ceased (all within a period of 90 yr), and the associated mining-induced subsidence of the coastal cliffs.     

Pathways of manganese in an open estuarine system

The distribution of Mn was examined in the bottom sediments and water column (suspended paniculate matter) of the Laurentian Trough. Gulf of St. Lawrence. A characteristic profile of Mn with depth in the sediment consisted of a Mn-enriched surface oxidized zone, less than 20 mm thick, and a Mn-depleted subsurface reducing zone. A subsurface Mn maximum occurred within the oxidized zone. Below this maximum the concentration dropped sharply to nearly constant residual levels in the reducing zone. The accumulating estuarine sediments are deficient in Mn compared to the river input of suspended matter and are definitely not the ultimate sink for manganese. Manganese escapes from the sediment by diffusion and resuspension, forming Mn-enriched, fine-grained particles which are flushed out in the estuarine circulation. 5.0 × 10⁹gyr^?1 of Mn, or 50% more than the river input of dissolved Mn. are exported to the open ocean. In spite of the efficient mobilization and export of Mn, the quantity exported is a small fraction (0.2%) of the total flux to the deep-sea sediments. This is related to the low levels of paniculate matter transported by the St. Lawrence River. The export phénomenon, however, is probably true of many coastal regions of muddy sediments and thus has interesting implications for the oceanic budget of Mn.     

Sedimentary structures of tidal flats: A journey from coast to inner estuarine region of eastern India

Sedimentary structures of some coastal tropical tidal flats of the east coast of India, and inner estuarine tidal point bars located at 30 to 50 kilometers inland from the coast, have been extensively studied under varying seasonal conditions. The results reveal that physical features such as flaser bedding, herringbone cross-bedding, lenticular bedding, and mud/silt couplets are common to both the environments. In fact, flaser bedding and lenticular bedding are more common in the point bar facies during the monsoon months than in the coastal tidal flat environments. Interference ripples, though common in both the environments, show different architectural patterns for different environmental domains. Interference ripples with thread-like secondary set overriding the earlier ripple-form, resembling wrinkle marks, are the typical features in estuarine point bars near the high water region. Because structures which are so far considered as key structures for near-coastal tidal flats are common to both the environments, caution should be exercised for deciphering palaeo-environments, particularly for Proterozoic rocks, where one has to depend only on physical sedimentary structures.     

Nitrogen loads to estuaries: Using loading models to assess the effectiveness of management options to restore estuarine water quality

Nitrogen (N) loading to estuaries has become a major concern for coastal planners. As urban development on coastal watershed continues, estuaries and bays are becoming more eutrophic, and cascading effects are being felt at every trophic level. Managers and stakeholders need to have a suite of effective management tools that can be applied to coastal watersheds to minimize the effects of eutrophication. We applied an N loading model and an estuarine loading model to examine the effectiveness of a suite of potential management options that could be implemented in Waquoit Bay, Cape Cod, Massachusetts. This estuarine system is a case study in which we can explore the relative potential effectiveness of decreasing inputs from wastewater and fertilizer-derived N, diverting nitrogenous runoff from impervious surfaces, altering zoning ordinances, preserving forested tracts of land as well as freshwater and saltwater wetlands, harvesting macroalgae, dredging estuary channels, and exterminating waterfowl. From a combination of simulation results, assessment of the magnitude of loads from different sources, and through different land covers, and the additional consideration of feasibility we identified management options with high, intermediate, and low potential effectiveness. Improvement of septic system performance, use of zoning regulations, preservation of forested tracts and freshwater bodies, and conservation of salt marshes emerged as the most promising avenues to manage N loads in our system. Installation of wastewater treatment plants, controlling fertilizer use, and harvesting macroalgae would potentially have intermediate success. Diversion of runoff from impervious surfaces, dredging, and extermination of waterfowl show little promise at reducing N loads. These conclusions potentially set priorities for decision-makers charged with the management of Waquoit Bay. The same procedures applied to another watershed-estuary system with different land covers and different estuarine features may differ. Evaluation studies like this need to be done for any particular site, since the watershed-estuary coupling and the loads delivered to the receiving estuary could differ. The Waquoit Bay case study provides an example of a protocol that leads to identification of the most promising management options.     

Evidence for millennial-scale climatic events in the sedimentary infilling of a macrotidal estuarine system,the Seine estuary (NW France)

High-resolution sedimentological and rock magnetic analyses from sediment cores collected in the Seine estuary record changes in coastal sedimentary dynamics linked to climatic variations during the late Holocene. Using AMS ¹⁴C and paleomagnetic data we present a first attempt in developing a reliable age model on macrotidal estuarine archives, with a decadal resolution. Correlations between sedimentary successions from the outer Seine estuary document the main sedimentary infilling phases of the system during the last 3000 years. Between 3000 and 1150 cal. BP sedimentary patterns reveal that sequence deposition and preservation are predominantly controlled by marine and tidal hydrodynamics while severe storm events are recorded at ca. 2700 and 1250 cal. BP in the outermost estuary. Conversely, the Medieval Warm Period (MWP; 900–1200 AD) is characterized by a drastic waning of the influence of marine hydrodynamics on sediment preservation. Pronounced episodes of Seine river floods indicate a much stronger impact of continental inputs on sedimentary patterns during this period. The onset of the Little Ice Age marks a diminishing influence of continental inputs while tidal and open marine hydrodynamics again exerted a primary control on the sedimentary evolution of the system during 1200–2003 AD. Coastal sedimentary dynamics as preserved within sedimentary successions appear to have been largely influenced by changes in storminess during the last 3000 years. We have matched the preservation of MWP Seine river floods, as revealed by sedimentological and rock magnetic proxy data, to a prolonged interval of weakened storminess in Normandy permitting the preservation of estuarine flood deposits within a context of reduced coastal erosion in northern Europe. The preservation of sedimentary successions in the Seine estuary is therefore maximal when climate conditions resembled those of the preferred low phase of the NAO on multidecadal timescales such as during 800–1200 AD (MWP). In contrast, increased removal and transport of estuarine sediments occur when winter storm activity greatly intensified over northwestern France. We report four prominent centennial-scale periods of stronger storminess, occurring with a pacing of ～1500 years, which are likely to be related to the last four Bond's Holocene cold events. Our results documenting a close link between coastal sedimentary dynamics, millennial-scale variations in Holocene climate and North Atlantic atmospheric circulation are fairly consistent with other records from Scandinavia, central Greenland and southern Europe.     

A review of recent developments in estuarine scalar flux estimation

The purpose of this contribution is to review recent developments in calculation of estuarine scalar fluxes, to suggest avenues for future improvement, and to place the idea of flux calculation in a broader physical and biogeochemical context. A scalar flux through an estuarine cross section is the product of normal velocity and scalar concentration, sectionally integrated and tidally averaged. These may vary on interannual, reasonal, tidal monthly, and event time scales. Formulation of scalar fluxes in terms of an integral scalar conservation expression shows that they may be determined either through “direct” means (measurement of velocity and concentration) or by “indirect” inference (from changes in scalar, inventory and source/sink terms). Direct determination of net flux at a cross section has a long and generally discouraging history in estuarine oceanography. It has proven difficult to extract statistically significant net (tidally averaged) fluxes from much larger flood and ebb transports, and the best mathematical representation of flux mechanisms is unclear. Observations further suggest that both lateral and vertical variations in scalar transport through estuarine cross sections are large, while estuarine circulation theory has focused on two-dimensional analyses that treatment either vertical or lateral variations but not both. Indirect estimates of net fluxes by determination of the other relevant terms in an integral scalar conservation balance may be the best means of determining scalar import-export in systems with residence times long relative to periods of tidal monthly fluctuations. But this method offers, little insight into the interaction of circulation modes and scalar fluxes, little help in verifying predictive models, and may also be difficult to apply in some circumstances. Thus, the need to understand, measure, and predict anthropogenic influences on transport or carbon, nutrient, suspended matter, trace metals, and other substances across the land-margin brings a renewed urgency to the issue of how to best carry out estuarine scalar flux determination. An interdisciplinary experiment is suggested to test present understanding, available instrument, and numerical models.     

Oceanography of the U.S. Pacific Northwest Coastal Ocean and estuaries with application to coastal ecology

Ocean processes are generally large scale on the U.S. Pacific Northwest coast; this is true of both seasonal variations and event-scale upwelling-downwelling fluctuations., which are highly energetic. Coastal upwelling supplies most of the macronutrients available for production, although the intensity of upwelling-favorable wind forcing increases southward while primary production and chlorophyll are higher in the north, off the Washington coast. This discrepancy could be related to several mesoscale features: the wider, more gently sloping shelf to the north, the existence of numerous submarine canyons to the north, the availability of Columbia River plume water and sediment north of the river mouth, and the existence of a semi-permanent eddy offshore of the Strait of Juan de Fuca. We suggest that these features have important effects on the magnitude and timing of macronutrient or micronutrient delivery to the plankton. These features are potentially important as well to transport pathways and residence times of planktonic larvae and to the development of harmful algal blooms. The coastal plain estuaries, with the exception of the Columbia River, are relatively small, with large tidal forcing and highly seasonal direct river inputs that are low to negligible during the growing season. Primary production in these estuaries is likely controlled not by river-driven stratification but by coastal upwelling and exchange with the ocean. Both baroclinic mechanisms (the gravitational circulation) and barotropic ones (lateral stirring by tide and, possibly, wind) contribute to this exchange. Because estuarine hydrography and ecology are so dominated by ocean signals, the coastal estuaries, like the coastal ocean, are largely synchronous on seasonal and event time scales, though, intrusions of the Columbia River plume can cause strong asymmetries between Washington and Oregon estuaries especially during spring downwelling conditions. Water property correlation increases between spring and summer as wind forcing becomes more spatially coherent along the coast. Estuarine habitat is structure not only, by large scale forcing but also by fine scale processes in the extensive intertidal zone, such as by solar heating or differential advection by tidal, curents.     

Simulation of water levels and extent of coastal inundation due to a cyclonic storm along the east coast of India

The devastation due to storm surge flooding caused by extreme wind waves generated by the cyclones is a severe apprehension along the coastal regions of India. In order to coexist with nature’s destructive forces in any vulnerable coastal areas, numerical ocean models are considered today as an essential tool to predict the sea level rise and associated inland extent of flooding that could be generated by a cyclonic storm crossing any coastal stretch. For this purpose, the advanced 2D depth-integrated (ADCIRC-2DDI) circulation model based on finite-element formulation is configured for the simulation of surges and water levels along the east coast of India. The model is integrated using wind stress forcing, representative of 1989, 1996, and 2000 cyclones, which crossed different parts of the east coast of India. Using the long-term inventory of cyclone database, synthesized tracks are deduced for vulnerable coastal districts of Tamil Nadu. Return periods are also computed for the intensity and frequency of cyclones for each coastal district. Considering the importance of Kalpakkam region, extreme water levels are computed based on a 50-year return period data, for the generation of storm surges, induced water levels, and extent of inland inundation. Based on experimental evidence, it is advocated that this region could be inundated/affected by a storm with a threshold pressure drop of 66 hpa. Also it is noticed that the horizontal extent of inland inundation ranges between 1 and 1.5 km associated with the peak surge. Another severe cyclonic storm in Tamil Nadu (November 2000 cyclone), which made landfall approximately 20 km south of Cuddalore, has been chosen to simulate surges and water levels. Two severe cyclonic storms that hit Andhra coast during 1989 and 1996, which made landfall near Kavali and Kakinada, respectively, are also considered and computed run-up heights and associated water levels. The simulations exhibit a good agreement with available observations from the different sources on storm surges and associated inundation caused by these respective storms. It is believed that this study would help the coastal authorities to develop a short- and long-term disaster management, mitigation plan, and emergency response in the event of storm surge flooding.     

A conceptual and numerical model for groundwater management: a case study on a coastal aquifer in southern Tuscany, Italy

Ongoing hydrogeological research aims to develop a correct management model for the Plio-Pleistocene multi-aquifer system of the Albegna River coastal plain (southern Tuscany, Italy); overexploitation of this aquifer for irrigation and tourism has caused seawater intrusion. The conceptual model is based on field and laboratory data collected during the 1995–2003 period. Meteoric infiltration and flows from the adjoining carbonate aquifer recharge the aquifer. Natural outflow occurs through a diffuse flow into the sea and river; artificial outflow occurs through intensive extraction of groundwater from wells. Water exchanges in the aquifer occur naturally (leakage, closing of aquitard) and artificially (multiscreened wells). The aquifer was represented by a three-dimensional finite element model using the FEFLOW numerical code. The model was calibrated for steady-state and transient conditions by matching computed and measured piezometric levels (February 1995–February 1996). The model helped establish that seawater intrusion is essentially due to withdrawals near the coast during the irrigation season and that it occurs above all in the Osa-Albegna sector, as well as along the river that at times feeds the aquifer. The effects of hypothetical aquifer exploitation were assessed in terms of water budget and hydraulic head evolution.