Turbulent and numerical mixing in a salt wedge estuary : dependence on grid resolution, bottom roughness, and turbulence closure

Author Posting. © American Geophysical Union, 2017. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Journal of Geophysical Research: Oceans 122 (2017): 692–712, doi:10.1002/2016JC011738.The Connecticut River is a tidal salt wedge estuary, where advection of sharp salinity gradients through channel constrictions and over steeply sloping bathymetry leads to spatially heterogeneous stratification and mixing. A 3-D unstructured grid finite-volume hydrodynamic model (FVCOM) was evaluated against shipboard and moored observations, and mixing by both the turbulent closure and numerical diffusion were calculated. Excessive numerical mixing in regions with strong velocities, sharp salinity gradients, and steep bathymetry reduced model skill for salinity. Model calibration was improved by optimizing both the bottom roughness (z0), based on comparison with the barotropic tidal propagation, and the mixing threshold in the turbulence closure (steady state Richardson number, Rist), based on comparison with salinity. Whereas a large body of evidence supports a value of Rist ∼ 0.25, model skill for salinity improved with Rist ∼ 0.1. With Rist = 0.25, numerical mixing contributed about 1/2 the total mixing, while with Rist = 0.10 it accounted for ∼2/3, but salinity structure was more accurately reproduced. The combined contributions of numerical and turbulent mixing were quantitatively consistent with high-resolution measurements of turbulent mixing. A coarser grid had increased numerical mixing, requiring further reductions in turbulent mixing and greater bed friction to optimize skill. The optimal Rist for the fine grid case was closer to 0.25 than for the coarse grid, suggesting that additional grid refinement might correspond with Rist approaching the theoretical limit. Numerical mixing is rarely assessed in realistic models, but comparisons with high-resolution observations in this study suggest it is an important factor.NSF Grant Number: OCE 0926427; ONR Grant Number: N00014-08-1-11152017-07-2

A nonhydrostatic version of FVCOM : 2. Mechanistic study of tidally generated nonlinear internal waves in Massachusetts Bay

Author Posting. © American Geophysical Union, 2010. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Journal of Geophysical Research 115 (2010): C12049, doi:10.1029/2010JC006331.The generation, propagation, and dissipation processes of large-amplitude nonlinear internal waves in Massachusetts Bay during the stratified season were examined using the nonhydrostatic Finite-Volume Coastal Ocean Model (FVCOM-NH). The model reproduced well the characteristics of the high-frequency internal waves observed in Massachusetts Bay in August 1998. The model experiments suggested that internal waves over Stellwagen Bank are generated by the interaction of tidal currents with steep bottom topography through a process of forming a large-density front on the western slope of the bank by the release of an initial density perturbation near ebb-flood transition, nonlinear steepening of the density front into a deep density depression, and disintegrating of the density depression into a wave train. Earth's rotation tends to transfer the cross-bank tidal kinetic energy into the along-bank direction and thus reduces the intensity of the density perturbation at ebb-flood transition and density depression in the flood period. The internal wave packet propagates as a leading edge feature of the internal tidal wave, and the faster propagation speed of the high-frequency internal waves in Massachusetts Bay is caused by Earth's rotation. The model experiments suggested that bottom friction can significantly influence the cross-bank scale of the density perturbation and thus the density depression during wave generation and the dissipation during the wave's shoaling. Inclusion of vertical mixing using the Mellor-Yamada level 2.5 turbulence closure model had only a marginal effect on wave evolution. The model results support the internal wave theory proposed by Lee and Beardsley (1974) but are in disagreement with the lee-wave mechanism proposed by Maxworthy (1979).This research was supported by NOAA g r a n t s DOC/NOAA/NA04NMF4720332 and DOC/NOAA/ NA05NMF4721131, U.S. GLOBEC Northwest Atlantic/Georges Bank Program NSF grants (OCE‐0606928, OCE‐0712903, OCE‐0732084, OCE‐0726851, OCE0814505), and MIT Sea Grant funds (2006‐RC‐103 and 2010‐R/RC‐116), NOAA NERACOOS Program for the UMASSD team and the Smith Chair in Coastal Oceanography, and NOAA grant (NA‐17RJ1223) for R.C. Beardsley. C. Chen’s contribution is also supported by Shanghai Ocean University under grants A‐2302‐10‐0003 and 09320503700 and the State Key Laboratory for Estuarine and Coastal Research, East China Normal University

Connectivity of the bay scallop (Argopecten irradians) in Buzzards Bay, Massachusetts, U.S.A.

Author Posting. © The Author(s), 2015. This is the author's version of the work. It is posted here by permission of John Wiley & Sons for personal use, not for redistribution. The definitive version was published in Fisheries Oceanography 24 (2015): 364-382, doi:10.1111/fog.12114.The harvest of bay scallops (Argopecten irradians) from Buzzards Bay, Massachusetts, USA undergoes large interannual fluctuations, varying by more than an order of magnitude in successive years. To investigate the extent to which these fluctuations may be due to yearly variations in the transport of scallop larvae from spawning areas to suitable juvenile habitat (settlement zones), a high-resolution hydrodynamic model was used to drive an individual-based model of scallop larval transport. Model results revealed that scallop spawning in Buzzards Bay occurs during a time when nearshore bay currents were principally directed up-bay in response to a persistent southwesterly sea breeze. This nearshore flow results in substantial transport of larvae from lower-bay spawning areas to settlement zones further up-bay. Averaged over the entire bay, the spawning-to-settlement zone connectivity exhibits little interannual variation. However, connectivities between individual spawning and settlement zones vary by up to an order of magnitude. The model results identified spawning areas that have the greatest probability of transporting larvae to juvenile habitat. Because managers may aim to increase scallop populations either locally or broadly, the high-connectivity spawning areas were divided into: 1) high larval retention and relatively little larval transport to adjoining settlement areas, 2) both significant larval retention and transport to more distant settlement areas, and 3) little larval retention but significant transport to distant settlement areas.This project was supported by the Woods Hole Sea Grant through award NA10OAR4170083. All modeling computations were made on the University of Massachusetts at Dartmouth’s (UMD’s) GPGPU cluster, which was acquired with support from NSF award CNS-0959382 and AFOSR DURIP award FA9550-10-1-0354.2016-07-1

A wet/dry point treatment method of FVCOM, part I: stability experiments

© The Author(s), 2022. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Chen, C., Qi, J., Liu, H., Beardsley, R., Lin, H., & Cowles, G. A wet/dry point treatment method of FVCOM, part I: stability experiments. Journal of Marine Science and Engineering, 10(7), (2022): 896, https://doi.org/10.3390/jmse10070896.A 3-dimensional wet/dry point treatment method was developed for the unstructured-grid Finite-Volume Community Ocean Model (FVCOM). Analytical equations were derived to examine discretized errors that occurred during the flooding/drying process by the wet/dry point treatment for the flooding/drying process. Numerical experiments were carried out for an idealized estuary, including the inter-tidal zone. The model results show that if the ratio of internal to external mode time steps (Isplit) is appropriately selected, FVCOM was capable of simulating the flooding/drying process with sufficient accuracy to ensure the mass conservation. The up-bound limit of Isplit was restricted by the bathymetric slope of the inter-tidal zone, external mode time step, horizontal/vertical resolution, and amplitude of tidal forcing at the open boundary, as well as the thickness of the viscous layer specified in the model. Criteria for time steps via these parameters were derived from these experiments, which provide a helpful guide in selectingIsplit for applying FVCOM to realistic geometric estuaries.This research was funded by the Georgia Sea Grant College Program under grant numbers NA26RG0373 and NA66RG0282, the Georgia DNR grants 024409-01 and 026450-01, the NSF Georges Bank/Northwest Atlantic GLOBEC program under grant number NSF-OCE 02-27679, and the SMAST fishery program under the NASA grant number NAG 13-02042

Positive interspecific associations consistent with social information use shape juvenile fish assemblages

Social information obtained from heterospecifics can enhance individual fitness by reducing environmental uncertainty, making it an important driver of mixed-species grouping behavior. Heterospecific groups are well documented among fishes, yet are notably more prevalent among juveniles than more advanced life stages, implying that the adaptive value of joining other species is greater during this developmental period. We propose this phenomenon can be explained by the heightened ecological relevance of heterospecifically produced cues pertaining to predation risk and or resources, as body-size uniformity inherent in early ontogeny yields greater overlap in predator and prey guild membership across juveniles of disparate taxa. To evaluate the putative role of information in shaping juvenile fish assemblages, we employed a joint species distribution model (JSDM), identifying nonrandom relationships among fishes collected in 785 seine hauls within the shallow littoral zones of a subtropical island. After accounting for species-environment relationships, which explained 39% of observed covariation in the abundance of 11 taxa, we detected high rates of positive association (84% of significant correlations) predominantly between mutual foraging guild members, consistent with assemblage patterns predicted to evolve under widespread interspecific information use. Affiliations occurred primarily between species characterized by neutral (i.e., noninteracting) or negative (i.e., predator-prey) relationships in later life stages, supporting the notion that heightened niche overlap due to body size homogeneity acted to increase the pertinence of information among juveniles. Taxa exerted varying degrees of influence on assemblage structure; however Eucinostomus spp., a gregarious generalist with exceptional information-production potential, had an effect several times that of all other species combined, further evidencing the likely role of information in motivating observed relationships. Co-occurrence and qualitative behavioral data inferred from remote underwater video surveys reinforced these conclusions. Collectively, these results suggest that positive interactions linked to information exchange can be among the principal factors organizing juvenile fish assemblages at local scales, highlighting the role of ontogeny in mediating the relevance and exploitation of information across speciesFunding from the Bonefish & Tarpon Trust, UMass Intercampus Marine Science Graduate Program, and the Department of Environmental Conservation at the University of Massachusetts Amherst. A. Danylchuk was supported the National Institute of Food & Agriculture, U.S. Department of Agriculture, the Massachusetts Agricultural Experiment Station, and Department of Environmental Conservation and is also a Bonefish & Tarpon Trust Research Fellow

A new high-resolution unstructured grid finite volume Arctic Ocean model (AO-FVCOM) : an application for tidal studies

Author Posting. © American Geophysical Union, 2009. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Journal of Geophysical Research 114 (2009): C08017, doi:10.1029/2008JC004941.A spherical coordinate version of the unstructured grid 3-D FVCOM (finite volume coastal ocean model) has been applied to the Arctic Ocean to simulate tides with a horizontal resolution ranging from 1 km in the near-coastal areas to 15 km in the deep ocean. By accurately resolving the irregular coastlines and bathymetry in the Arctic Ocean coastal regions, this model reproduces the diurnal (K1 and O1) and semidiurnal (M2 and S2) tidal wave dynamics and captures the complex tidal structure along the coast, particularly in the narrow straits of the Canadian Archipelago. The simulated tidal parameters (harmonic constituents of sea surface elevation and currents) agree well with the available observational data. High-resolution meshes over the continental shelf and slope capture the detailed spatial structure of topographic trapped shelf waves, which are quite energetic along the Greenland, Siberia, and Spitsbergen continental slope and shelf break areas. Water stratification influences the vertical distribution of tidal currents but not the water transport and thus tidal elevation. The comparison with previous finite difference models suggests that horizontal resolution and geometric fitting are two prerequisites to simulate realistically the tidal energy flux in the Arctic Ocean, particularly in the Canadian Archipelago.This research was supported by the NSF Office of Polar Programs through grants OPP ARC-0712903, ARC- 0732084, and ARC-0804029 for C. Chen, G. Gao, and G. Cowles; OPP ARC-0804010 and ARC-0712848 for A. Proshutinsky; OPP ANT-0523223, ARC0712848, NOAA Cooperative Agreement NA17RJ1223 (409) and the WHOI Smith Chair for R. C. Beardsley. J. Qi was supported by the SMAST fishery program under NOAA grants NA04NMF4720332 and NA05NMF4721131. The spherical coordinate version of FVCOM was developed with initial funds from NSF grants OCE-0606928 and OCE- 0726851. Gao was also supported by the Chinese NSF Arctic Ocean grant under contract 40476007

FVCOM validation experiments : comparisons with ROMS for three idealized barotropic test problems

Author Posting. © American Geophysical Union, 2008. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Journal of Geophysical Research 113 (2008): C07042, doi:10.1029/2007JC004557.The unstructured-grid Finite-Volume Coastal Ocean Model (FVCOM) is evaluated using three idealized benchmark test problems: the Rossby equatorial soliton, the hydraulic jump, and the three-dimensional barotropic wind-driven basin. These test cases examine the properties of numerical dispersion and damping, the performance of the nonlinear advection scheme for supercritical flow conditions, and the accuracy of the implicit vertical viscosity scheme in barotropic settings, respectively. It is demonstrated that FVCOM provides overall a second-order spatial accuracy for the vertically averaged equations (i.e., external mode), and with increasing grid resolution the model-computed solutions show a fast convergence toward the analytic solutions regardless of the particular triangulation method. Examples are provided to illustrate the ability of FVCOM to facilitate local grid refinement and speed up computation. Comparisons are also made between FVCOM and the structured-grid Regional Ocean Modeling System (ROMS) for these test cases. For the linear problem in a simple rectangular domain, i.e., the wind-driven basin case, the performance of the two models is quite similar. For the nonlinear case, such as the Rossby equatorial soliton, the second-order advection scheme used in FVCOM is almost as accurate as the fourth-order advection scheme implemented in ROMS if the horizontal resolution is relatively high. FVCOM has taken advantage of the new development in computational fluid dynamics in resolving flow problems containing discontinuities. One salient feature illustrated by the three-dimensional barotropic wind-driven basin case is that FVCOM and ROMS simulations show different responses to the refinement of grid size in the horizontal and in the vertical.For this work, H. Huang and G. Cowles were supported by the Massachusetts Marine Fisheries Institute (MFI) through NOAA grants DOC/NOAA/NA04NMF4720332 and DOC/ NOAA/NA05NMF472113; C. Chen was supported by NSF grants (OCE0234545, OCE0606928, OCE0712903, OCE0732084, and OCE0726851), NOAA grants (NA160P2323, NA06RG0029, and NA960P0113), MIT Sea grant (2006-RC-103), and Georgia Sea grant (NA26RG0373 and NA66RG0282); C. Winant was supported through NSF grant OCE-0726673; R. Beardsley was supported through NSF OCE—0227679 and the WHOI Smith Chair; K. Hedstrom was supported through NASA grant NAG13– 03021 and the Arctic Region Supercomputing Center; and D. Haidvogel was supported through grants ONR N00014- 03-1-0683 and NSF OCE 043557

Physical mechanisms for the offshore detachment of the Changjiang Diluted Water in the East China Sea

Author Posting. © American Geophysical Union, 2008. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Journal of Geophysical Research 113 (2008): C02002, doi:10.1029/2006JC003994.Physical mechanisms for the summertime offshore detachment of the Changjiang Diluted Water (CDW) into the East China Sea are examined using the high-resolution, unstructured-grid, Finite-Volume Coastal Ocean Model (FVCOM). The model results suggest that isolated low salinity water lens detected west of Cheju Island can be formed by (1) a large-scale adjustment of the flow field to the Changjiang discharge and (2) the detachment of anticyclonic eddies as a result of baroclinic instability of the CDW front. Adding the Changjiang discharge intensifies the clockwise vorticity of the subsurface current (originating from the Taiwan Warm Current) flowing along the 50-m isobath and thus drives the low-salinity water in the northern coastal area of the Changjiang mouth offshore over a submerged plateau that extends toward Cheju Island. Given a model horizontal resolution of less than 1.0 km, the CDW front becomes baroclinically unstable and forms a chain of anticyclonic and cyclonic eddies. The offshore detachment of anticyclonic eddies can carry the CDW offshore. This process is enhanced under northward winds as a result of the spatially nonuniform interaction of wind-induced Ekman flow and eddy-generated frontal density currents. Characteristics of the model-predicted eddy field are consistent with previous theoretical studies of baroclinic instability of buoyancy-driven coastal density currents and existing satellite imagery. The plume stability is controlled by the horizontal Ekman number. In the Changjiang, this number is much smaller than the criterion suggested by a theoretical analysis.The development of FVCOM is supported by the Massachusetts Fisheries Institute through NOAA grants DOC/ NOAA/NA04NMF4720332 and DOC/NOAA/NA05NMF4721131 and also the U.S. GLOBEC Northwest Atlantic/Georges Bank program through NSF grants OCE-0234545 and OCE-0227679, NOAA grant NA160P2323 and ONR subcontract grant from Woods Hole Oceanographic Institution. P. Ding is supported by the Chinese National Key Basic Research Project grant 2002CB412403. X. Mao is supported by the National Natural Science Foundation of China (NSFC) grant 40576079

Application and comparison of Kalman filters for coastal ocean problems : an experiment with FVCOM

Author Posting. © American Geophysical Union, 2009. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Journal of Geophysical Research 114 (2009): C05011, doi:10.1029/2007JC004548.Twin experiments were made to compare the reduced rank Kalman filter (RRKF), ensemble Kalman filter (EnKF), and ensemble square-root Kalman filter (EnSKF) for coastal ocean problems in three idealized regimes: a flat bottom circular shelf driven by tidal forcing at the open boundary; an linear slope continental shelf with river discharge; and a rectangular estuary with tidal flushing intertidal zones and freshwater discharge. The hydrodynamics model used in this study is the unstructured grid Finite-Volume Coastal Ocean Model (FVCOM). Comparison results show that the success of the data assimilation method depends on sampling location, assimilation methods (univariate or multivariate covariance approaches), and the nature of the dynamical system. In general, for these applications, EnKF and EnSKF work better than RRKF, especially for time-dependent cases with large perturbations. In EnKF and EnSKF, multivariate covariance approaches should be used in assimilation to avoid the appearance of unrealistic numerical oscillations. Because the coastal ocean features multiscale dynamics in time and space, a case-by-case approach should be used to determine the most effective and most reliable data assimilation method for different dynamical systems.P. Malanotte-Rizzoli and J. Wei were supported by the Office of Naval Research (ONR grant N00014-06-1- 0290); C. Chen and Q. Xu were supported by the U.S. GLOBEC/Georges Bank program (through NSF grants OCE-0234545, OCE-0227679, OCE- 0606928, OCE-0712903, OCE-0726851, and OCE-0814505 and NOAA grant NA-16OP2323), the NSF Arctic research grants ARC0712903, ARC0732084, and ARC0804029, and URI Sea Grant R/P-061; P. Xue was supported through the MIT Sea Grant 2006-RC-103; Z. Lai, J. Qi, and G. Cowles were supported through the Massachusetts Marine Fisheries Institute (NOAA grants NA04NMF4720332 and NA05NMF4721131); and R. Beardsley was supported through U.S. GLOBEC/Georges Bank NSF grant OCE-02227679, MIT Sea Grant NA06OAR1700019, and the WHOI Smith Chair in Coastal Oceanography

Ernst Freund as Precursor of the Rational Study of Corporate Law

Gindis, David, Ernst Freund as Precursor of the Rational Study of Corporate Law (October 27, 2017). Journal of Institutional Economics, Forthcoming. Available at SSRN: https://ssrn.com/abstract=2905547, doi: https://dx.doi.org/10.2139/ssrn.2905547The rise of large business corporations in the late 19th century compelled many American observers to admit that the nature of the corporation had yet to be understood. Published in this context, Ernst Freund's little-known The Legal Nature of Corporations (1897) was an original attempt to come to terms with a new legal and economic reality. But it can also be described, to paraphrase Oliver Wendell Holmes, as the earliest example of the rational study of corporate law. The paper shows that Freund had the intuitions of an institutional economist, and engaged in what today would be called comparative institutional analysis. Remarkably, his argument that the corporate form secures property against insider defection and against outsiders anticipated recent work on entity shielding and capital lock-in, and can be read as an early contribution to what today would be called the theory of the firm.Peer reviewe