CMS-Flow:Salinity Calculation
Salinity Calculations in the Coastal Modeling System
Introduction
Salinity refers to the salt content of water, with values ranging from 0 for fresh water to 31-35 ppt (parts per thousand) for ocean water. In water bodies with poor mixing and limited water exchange, or those experiencing high evaporation, salinity can be greater and lead to formation of brine. Table 1 presents typical values and nomenclature for describing the degree of water salinity:
In coastal zones and estuaries, both temporal and spatial variations in salinity are controlled by changes in circulation, waves, tides, precipitation, evaporation, and freshwater inflows. These changes in salinity can have major effects on water density and water stratification, changing circulation patterns. Dynamic behavior of suspended sediment can be controlled by the salinity-driven flow and mixing. Any sustained changes to salinity can directly change the aggregation and consolidation of cohesive sediment as well (Nicholson and O’Connor 1986). Salinity can also alter the water chemistry that is closely related to marine organisms. Distribution and abundance of marine life will change water turbidity and define water quality in coastal and estuarine systems. Modifications of coastal inlets, such as channel deepening and widening and rehabilitation or extension of jetties may alter the salinity distribution within the estuary.
Salinity Calculations in CMS
The CMS calculates the salinity field based on the following 2D salinity conservation equation:
(1) |
where S is depth-averaged salinity; d is total water depth, qx and qy are flow per unit width in the x- and y-axis direction, respectively; Kx and Ky are the salt diffusion coefficients in the corresponding x- and y-axis directions, and P and E are precipitation and evaporation in m/year, respectively. Equation (1) represents the horizontal fluxes of salt in water bodies and is balanced by exchanges of salt via diffusive fluxes. Major processes contributing to the salinity are freshwater inflows from rivers, vertical fluxes of freshwater by precipitation and evaporation at the water surface, and groundwater fluxes, which can be specified as the surface and bottom boundary conditions in the equation.
Model Assumptions
CMS-Flow is presently capable of 2D salinity computations in both the explicit and implicit solvers. 3D representation of salinity in the CMS, discussed herein, is being tested. The simulation of salinity can often require a 3D solution due to the presence of vertical salinity gradients that can significantly influence flow. It is therefore important to understand the limitations of 2D salinity simulations, and apply them only when the assumptions inherent in 2D simulations are valid. Typically, 2D salinity calculations are valid when the salinity is well mixed over the water column. These conditions are usually met for shallow bays with open exchanges to the ocean or gulf, dominant tidal signals, and sufficient wind energy to provide the vertical mixing. Also, the assumption of sufficient energy to mix over the water column is valid under storm conditions, even for deeper water bodies. Finally, when the exchange with the open sea is restricted by an inlet, the tidal range is an important indicator of vertical mixing conditions. For low tide ranges, significant vertical stratification can occur, even in shallow bays and estuaries, especially when the winds are calm. Pritchard (1955) and Cameron and Pritchard (1963) have classified estuaries using stratification and salinity distribution as the governing criteria, and these classifications can be used for guidance in applying the 2D simulations.
The lateral mixing for salinity in the CMS Flow model is the same as the lateral mixing in the momentum equations.
Initial and Boundary Conditions
The initial salinity is specified as a constant in the whole domain. The value of the constant is specified in the SMS 10.1 interface. Inflow salinity concentrations are applied at specified salinity boundary cell strings. Salinity cell strings are specified in the same manner as the hydrodynamic boundary cells strings.
Numerical Methods
The salinity transport equation is solved with an explicit, finite volume method. The advection term is discretized with upwind scheme, and the diffusion term is discretized with the standard central difference scheme.
Additional Information
This wiki technical note was prepared and funded under the Coastal Inlet Research Program (CIRP) and was written by Dr. Honghai Li (Honghai.Li@usace.army.mil, voice: 601-634-2840, fax: 601-634-3080) of the U.S. Army Engineer Research and Development Center (ERDC), Coastal and Hydraulics Laboratory (CHL), Dr. Christopher W. Reed (Chris_Reed@URSCorp.com) of URS Corporation, and Mitchell E. Brown (Mitchell.E.Brown@usace.army.mil) of ERDC, CHL. Alejandro Sanchez of CIRP provided hydrodynamic model information for the Humboldt Bay example. The CIRP Program Manager, Dr. Julie D. Rosati (Julie.D.Rosati@usace.army.mil), the assistant Program Manager, Dr. Nicholas C. Kraus, the Chief of the Coastal Engineering Branch at CHL, Dr. Jeffrey P. Waters, and Dr. Lihwa Lin, the Coastal Engineering Branch, reviewed this CHETN. Files for the example may be obtained by contacting the author.
References
- Aquaveo. 2010. SMS: XY Series Files (*.xys). http://www.xmswiki.com/xms/SMS:XY_Series_Files_(*.xys).
- Buttolph, A. M., C. W. Reed, N. C. Kraus, N. Ono, M. Larson, B. Camenen, H. Hanson, T. Wamsley, and A.
K. Zundel, A. K. 2006. Two-dimensional depth-averaged circulation model CMS-M2D: Version 3.0, Report 2, sediment transport and morphology change. Coastal and Hydraulics Laboratory Technical Report ERDC/CHL-TR-06-7. Vicksburg, MS: U.S. Army Engineer Research and Development Center.
- Cameron, W. M. and D. W. Pritchard. 1963. Estuaries. In M. N. Hill (editor): The Sea vol. 2, John Wiley and Sons, New York, 306 - 324.
- CeNCOOS. 2010. Central & Northern California Ocean Observing System. http://www.cencoos.org/.
- Cooke, W. H. III, K. Grala, and C. L. Wax. 2008. A Method for estimating pan evaporation for inland and coastal regions of the southeastern U. S. Southeastern Geographer, 48(2): 149–171.
- Demirbilek, Z., K. J. Connell, N. J. MacDonald, and A. K. Zundel. 2008. Particle Tracking Model in the SMS10: IV. Link to Coastal Modeling System. Coastal and Hydraulics Engineering Technical Note ERDC/CHL CHETN-IV-71. Vicksburg, MS: U.S. Army Engineer Research and Development Center. http://chl.erdc.usace.army.mil/chetn.
- Lin, L., Z. Demirbilek, H. Mase, J. Zheng, and F. Yamada. 2008. CMS-Wave: A nearshore spectral wave processes model for coastal inlets and navigation projects. Coastal and Hydraulics Laboratory Technical Report ERDC/CHL-TR-08-13. Vicksburg, MS: U.S. Army Engineer Research and Development Center.
- Moon, I. 2005. Impact of a coupled ocean wave-tide-circulation system on coastal modeling. Ocean Modeling, 8: 203-236.
- Pritchard, D. W. 1955. Estuarine circulation patterns. Proceedings of the American Society of Civil Engineers 81, no 717, 1 – 11.
- URS. 2010. USACE White Ditch evaluation and design – hydrodynamic and salinity transport modeling. Tallahassee, FL.