Silicones in Chesapeake Bay sediments

Sediment samples from the upper reaches to the mouth of the Chesapeake Bay were analysed for extractable organic silicon (silicone) content. Organic silicon in the sediments ranged from less than 0.2 to some 36 ppm, w/w. In general, silicone tended to accumulate in high depositional areas: the northern bay, and along the shores to the south. While present, anthropogenic inputs in terms of defined source were not clear. It appears that processes in the bay act to rapidly remove silicone from the water column, probably on a function of the flocculation/precipitation of seston, especially in the northern area.     

Wave measurement and modeling in Chesapeake Bay

Three recently measured wind and wave data sets in the northern part of Chesapeake Bay (CB) are presented. Two of the three data sets were collected in late 1995. The third one was collected in July of 1998. The analyzed wind and wave data show that waves were dominated by locally generated, fetch limited young wind seas. Significant wave heights were highly correlated to the local driving wind speeds and the response time of the waves to the winds was about 1 h. We also tested two very different numerical wave models, Simulation of WAves Nearshore (SWAN) and Great Lakes Environmental Research Laboratory (GLERL), to hind-cast the wave conditions against the data sets. Time series model–data comparisons made using SWAN and GLERL showed that both models behaved well in response to a suddenly changing wind. In general, both SWAN and GLERL over-predicted significant wave height; SWAN over-predicted more than GLERL did. SWAN had a larger scatter index and a smaller correlation coefficient for wave height than GLERL had. In addition, both models slightly under-predicted the peak period with a fairly large scatter and low correlation coefficient. SWAN predicted mean wave direction better than GLERL did. Directional wave spectral comparisons between SWAN predictions and the data support these statistical comparisons. The GLERL model was much more computationally efficient for wind wave forecasts in CB. SWAN and GLERL predicted different wave height field distributions for the same winds in deeper water areas of the Bay where data were not available, however. These differences are as yet unresolved.     

Phytoplankton Kinetics in the Chesapeake Bay Eutrophication Model

The CE-Qual-ICM model computes phytoplankton biomass and production as a function of temperature, light, and nutrients. Biomass is computed as carbon while inorganic nitrogen, phosphorus, and silica are considered as nutrients. Model formulations for production, metabolism, predation, nutrient limitation, and light limitation are detailed. Methods of parameter determination and parameter values are presented. Results of model application to a ten-year period in Chesapeake Bay indicate the model provides reasonable representations of observed biomass, nutrient concentrations, and limiting factors. Computed primary production agrees with observed under light-limited conditions. Under strongly nutrient-limited conditions, computed product is less than observed. The production characteristics of the model are similar to behavior reported for several similar models. Process omitted from the model that may account for production shortfalls include variable algal stoichiometry, use of urea as nutrient, and vertical migration by phytoplankton.     

Resonance and sea level variability in Chesapeake Bay

A numerical model is used to determine the resonant period and quality factor Q of Chesapeake Bay and explore physical mechanisms controlling the resonance response in semi-enclosed seas. At the resonant period of 2 days, the mouth-to-head amplitude gain is 1.42 and Q is 0.9, indicating that Chesapeake Bay is a highly dissipative system. The modest amplitude gain results from strong frictional dissipation in shallow water. It is found that the spatial distribution of energy dissipation varies with forcing frequency. While energy at tidal frequencies is dissipated around topographic hotspots distributed throughout the Bay, energy dissipation at subtidal frequencies is mainly concentrated in the shallow-water lower Bay. An analytic calculation shows that the bottom friction parameter is much larger in Chesapeake Bay than in other coastal systems with strong resonance response. The model-predicted amplitude gains and phase changes agree well with the observations at semidiurnal and diurnal tidal frequencies. However, the predicted amplitude gain in the resonant frequency band (34–54 h period) falls below that inferred from band-passed sea level observations. This discrepancy can be attributed to the local wind forcing which amplifies the sea level response in the upper Bay. The model is also used to show that rising sea levels associated with global warming will shift the resonance period of Chesapeake Bay closer to the diurnal tides and thus exacerbate flooding problems by causing an increase in tidal ranges.     

Flow patterns at the Chesapeake Bay entrance

The spatial and temporal variability of water entering and leaving the Chesapeake Bay estuary was determined with a spatial resolution of 75 m. The four cruises during which the observations were made took place under different conditions of freshwater discharge, tidal phase, and wind forcing. The tidal variability of the flows was dominated by the semidiurnal constituents that displayed greatest amplitudes and phase lags near the surface and in the channels that lie at the north and south sides of the entrance. The subtidal variability of the flows was classified into two general scenarios. The first scenario occurred during variable or persistently non-southwesterly winds. Under these conditions there was surface outflow and bottom inflow in the two channels, inflow over the shoal between the two channels, and possible anticyclonic gyre formation over the shoal. The flow pattern in the channels was produced by gravitational circulation and wind forcing. Over the shoal it was caused by tidal rectification and wind forcing. The second scenario occurred during persistently southwesterly winds. The anticyclonic gyre over the shoal vanished suggesting that wind forcing dominated the tidal rectification mechanism over the shoal, while gravitational circulation and wind forcing continued to cause the flows in the channels. In both scenarios, most of the volume exchange took place in the channels.     

Cross-Media Models of the Chesapeake Bay Watershed and Airshed

A continuous, deterministic watershed model of the Chesapeake Bay watershed, linked to an atmospheric deposition model is used to examine nutrient loads to the Chesapeake Bay under different management scenarios. The Hydrologic Simulation Program - Fortran, Version 11 simulation code is used at an hourly time-step for ten years of simulation in the watershed. The Regional Acid Deposition Model simulates management options in reducing atmospheric deposition of nitrogen. Nutrient loads are summed over daily periods and used for loading a simulation of the Chesapeake estuary employing the Chesapeake Bay Estuary Model Package. Averaged over the ten-year simulation, loads are compared for scenarios under 1985 conditions, forecasted conditions in the year 2000, and estimated conditions under a limit of technology scenario. Limit of technology loads are a 50%, 64%, and 42% reduction from the 1985 loads in total nitrogen, total phosphorus, and total suspended solids, respectively. Urban loads, which include point source, on-site wastewater disposal systems, combined sewer overflows, and nonpoint source loads have the highest flux of nutrient loads to the Chesapeake, followed by crop land uses.on assignment from NOAA Air Resources Laboratory     

Coupling Suspension Feeders to the Chesapeake Bay Eutrophication Model

A model predicting suspension-feeding bivalve biomass and its interactions with water quality has been developed and coupled with the Chesapeake Bay Eutrophication Model. This coupling included deposition of filtered particulate matter to the sediments and the recycling of inorganic nutrients back to the water column. Because individual size is a crucial determinant of bivalve filtration and respiration rates, an empirical function, was developed from data, relating computed areal biomass to size, which was then used to adjust these rates during the simulation. Biomass was strongly related to the eutrophication model's predictions of organic and total solids distributions, as well as to bottom water dissolved oxygen. The tight coupling between seasonal organic matter concentration and biomass suggested that food, or the ability of suspension feeders to ingest it given present total solids loadings, is a limiting factor baywide. Hypoxia and anoxia also reduced benthic biomass in affected locations. High site-specific temporal variability observed in the data may contain a large component of spatial patchiness, on scales below which the present estuarine eutrophication model could resolve. Further insights will be needed to incorporate the effects of patchiness, as well as other important spatial and temporal signals, such as predation and recruitment.     

Tidal effects on estuarine circulation and outflow plume in the Chesapeake Bay

The response of the Chesapeake Bay to river discharge under the influence and absence of tide is simulated with a numerical model. Four numerical experiments are examined: (1) response to river discharge only; (2) response to river discharge plus an ambient coastal current along the shelf outside the bay; (3) response to river discharge and tidal forcing; and (4) response to river discharge, tidal forcing, and ambient coastal current. The general salinity distribution in the four cases is similar to observations inside the bay. Observed features, such as low salinity in the western side of the bay, are consistent in model results. Also, a typical estuarine circulation with seaward current in the upper layer and landward current in the lower layer is obtained in the four cases. The two cases without tide produce stronger subtidal currents than the cases with tide owing to greater frictional effects in the cases with tide. Differences in salinity distributions among the four cases appear mostly outside the bay in terms of the outflow plume structure. The two cases without tide produce an upstream (as in a Kelvin wave sense) or northward branch of the outflow plume, while the cases with tide produce an expected downstream or southward plume. Increased friction in the cases with tide changes the vertical structure of outflow at the entrance to the bay and induces large horizontal variations in the exchange flow. Consequently, the outflow from the bay is more influenced by the bottom than in the cases without tide. Therefore, a tendency for a bottom-advected plume appears in the cases with tide, rather than a surface-advected plume, which develops in the cases without tide. Further analysis shows that the tidal current favors a salt balance between the horizontal and vertical advection of salinity around the plume and hinders the upstream expansion of the plume outside the bay.     

Sea level acceleration and variability in the Chesapeake Bay: past trends,future projections,and spatial variations within the Bay

Ocean Dynamics - Fast sea level rise (SLR) is causing a growing risk of flooding to coastal communities around the Chesapeake Bay (hereafter, CB or “the Bay”), but there are also...     

Aerial flux of particulate hydrocarbons to the Chesapeake Bay Estuary

Bulk hydrocarbon deposition rates have been measured over a 15 month period at four stations in south-eastern Virginia surrounding the lower Chesapeake Bay. A nearly linear trend of atmospheric particulate deposition was recorded. Deposition rates at the urban station (195 μg m^?2 day^?1) were aproximately three times greater than those recorded for nonurban and coastal locations (mean value 69 μg m^?2 day^?1). The increased levels at the urban location were attributed to localized source inputs. Anthropogenic hydrocarbons accounted for approximately 50% of the total deposition at all stations. Significant biogenic inputs were indicated by the odd/even n-alkane distribution. A minimum flux to the water surface, based on mean nonurban deposition rates (24 mg yr^?1), indicated an annual particulate hydrocarbon flux of +275 metric tons. Little information is available for the comparison of additional source inputs; however, the data reported here indicate that the aerial deposition of hydrocarbons is of the same order of magnitude as the input from municipal wastewater facilities and accidental discharge and is a potentially significant source of hydrocarbon pollution to the Chesapeake Bay Estuary.     

Effluent trading for water quality management: concept and application to the Chesapeake Bay watershed

This paper examines the present and potential role of effluent trading in water quality management. In particular, it focuses upon the case of the Chesapeake Bay on the east coast of the US, where the implementation of a trading system has been discussed and undertaken. Potential benefits of effluent trading include advantages such as the following: (1) With appropriate monitoring and enforcement, the total pollutant loadings can be kept at or below the prespecified level. (2) New and expanding dischargers can be accommodated, as long as they purchase credits. (3) Abatement costs of pollutants can be reduced. (4) Flexible regulations incorporating trading can reduce the incentive for industries to relocate to areas with less stringent water quality regulation. (5) Broader environmental goals can be addressed, such as wildlife habitat provision and endangered species protection. (6) Preliminary studies with a view to trading-system implementation encourage discussion and dialogue among stakeholders, and positively foster concerted, holistic solutions for maintenance of water bodies.     

Wind effects on the lateral structure of density-driven circulation in Chesapeake Bay

The response of the density-driven circulation in the Chesapeake Bay to wind forcing was studied with numerical experiments. A model of the bay with realistic bathymetry was first applied to produce the density-driven flow under average river discharge and tidal forcing. Subsequently, four spatially uniform wind fields (northeasterly, northwesterly, southwesterly, and southeasterly) were imposed to examine the resulting cross-estuary structure of salinity and flow fields. In general, northeasterly and northwesterly winds intensified the density-driven circulation in the upper and middle reaches of the bay, whereas southeasterly and southwesterly winds weakened it. The response was different in the lower bay, where downwind flow from the upper and middle reaches of the bay competed with onshore/offshore coastal flows. Wind remote effects were dominant, over local effects, on volume transports through the bay entrance. However, local effects were more influential in establishing the sea-level slopes that drove subtidal flows and salinity fields in most of the bay. The effect of vertical stratification on wind-induced flows was also investigated by switching it off. The absence of stratification allowed development of Ekman layers that reached depths of the same order as the water depth. Consequently, bathymetric effects became influential on the homogeneous flow structure causing the wind-induced flow inside the bay to show a marked transverse structure: downwind over the shallow areas and upwind in the channels. In the presence of stratification, Ekman layers became shallower and the wind-induced currents showed weaker transverse structure than those that developed in the absence of stratification. In essence, the wind-driven flows were horizontally sheared under weak stratification and vertically sheared under stratified conditions.     

Evaluation of butyltin compounds in Maryland waters of Chesapeake Bay

A monthly sampling programme for dibutyltin (DBT), tributyltin (TBT) and tetrabutyltin (TTBT) was initiated for a period of one year (July 1985–June 1986) in the Maryland waters of Chesapeake Bay. Concentrations of the above butyltin species were evaluated in the microlayer and water column of eight sampling stations representing two small and two large marinas, a large harbour, two major river systems and a heavily used shipping channel. DBT concentrations in the microlayer were generally higher in the four marinas when compared with the other stations. The highest DBT concentration reported in the microlayer was 1156 ng l⁻¹. Mean microlayer TBT concentrations ranged from 54–310 ng l⁻¹ in the four marinas. Three TBT concentrations ranging from 1049–1171 ng l⁻¹ were reported in the microlayer of the marinas. TBT concentrations of 41 and 29 mg l⁻¹ were detected in the microlayer of a heavily used shipping channel (C & D Canal) during May and June. TTBT concentrations were not detected in the microlayer at most stations during the 12 month sampling period.Mean DBT concentrations in the water column ranged from 23–145 ng l⁻¹ in the four marinas. DBT concentrations in the water column of the other stations were < 35 ng l⁻¹. Mean water column concentrations of TBT ranged from 51–408 ng l⁻¹ in all four marinas. Peak concentrations of TBT were reported in May and June for the various marinas. The highest TBT concentration reported in the water column was 998 ng l⁻¹. TBT concentrations of 20–24 ng l⁻¹ were reported in one of the river systems (Potomac River). TTBT concentrations were not detected in the water column at most of the stations.     

Evaluation of weather forecast systems for storm surge modeling in the Chesapeake Bay

Accurate forecast of sea-level heights in coastal areas depends, among other factors, upon a reliable coupling of a meteorological forecast system to a hydrodynamic and wave system. This study evaluates the predictive skills of the coupled circulation and wind-wave model system (ADCIRC+SWAN) for simulating storm tides in the Chesapeake Bay, forced by six different products: (1) Global Forecast System (GFS), (2) Climate Forecast System (CFS) version 2, (3) North American Mesoscale Forecast System (NAM), (4) Rapid Refresh (RAP), (5) European Center for Medium-Range Weather Forecasts (ECMWF), and (6) the Atlantic hurricane database (HURDAT2). This evaluation is based on the hindcasting of four events: Irene (2011), Sandy (2012), Joaquin (2015), and Jonas (2016). By comparing the simulated water levels to observations at 13 monitoring stations, we have found that the ADCIR+SWAN System forced by the following: (1) the HURDAT2-based system exhibited the weakest statistical skills owing to a noteworthy overprediction of the simulated wind speed; (2) the ECMWF, RAP, and NAM products captured the moment of the peak and moderately its magnitude during all storms, with a correlation coefficient ranging between 0.98 and 0.77; (3) the CFS system exhibited the worst averaged root-mean-square difference (excepting HURDAT2); (4) the GFS system (the lowest horizontal resolution product tested) resulted in a clear underprediction of the maximum water elevation. Overall, the simulations forced by NAM and ECMWF systems induced the most accurate results best accuracy to support water level forecasting in the Chesapeake Bay during both tropical and extra-tropical storms.     

基于FVCOM的太湖梅梁湾夏季水温、溶解氧模拟及其影响机制初探

水体中的溶解氧是表征水生生态系统健康与否的重要参数之一.本研究基于太湖2008年8月16 20日的风速、风向、短波辐射等气象场资料以及实测的相关水质参量,利用FVCOM(即非结构化网格有限体积近海海洋模型)模式对太湖梅梁湾三维水温以及水体中溶解氧进行模拟,模拟结果与实测值基本吻合,水温的验证回归方程为y=1.02x,R²为0.690;溶解氧的R²为0.760.同时对溶解氧浓度的时空分布,梅梁湾溶解氧的"源"和"汇"及其贡献进行了分析.结果表明:太阳辐射、风速是影响水温日成层现象的重要因子;受水温和光照的影响,夏季梅梁湾的溶解氧存在垂直差异,呈现出"双峰双谷"的日变化特征;浮游植物光合作用制氧是水中溶解氧的最重要来源,水下光衰减直接控制着初级生产力的垂直分布;浮游植物呼吸及死亡是溶解氧的最大消耗者,余下依次为底泥耗氧、碳化需氧、细菌呼吸耗氧和硝化作用耗氧.     

Validation and Application of the Second Generation Three Dimensional Hydrodynamic Model of Chesapeake Bay

The validation and subsequent application of the current three-dimensional numerical hydrodynamic model of Chesapeake Bay is presented. The numerical model solves conservation equations for water mass, momentum, salinity, and heat on a boundary-fitted grid in the horizontal plane and a Cartesian z-grid in the vertical. A generalized ADI finite difference scheme is employed in conjunction with mode splitting technique, solving external and the internal modes. The 10-year boundary conditions including tide, slinity, temperature, wind, heat exchange coefficient, river and non-point source flows were constructed. Model validation was accomplished by demonstrating the model's ability to reproduce observed data over time scales ranging from tidal to seasonal periods. The major parameters compared include tidal elevation, intra-tidal and residual velocities, salinity, temperature, stratification, and flux calculated through the Bay mouth.After validation, the model was applied to simulate bay hydrodynamics for the 10 years of 1985–94. These results were used to drive the three-dimensional water quality model of Chesapeake Bay, which is discussed in a companion paper.     

湖泊中溶解氧极大值之成因

本文对位于台湾南部两个次高海拔湖泊－大鬼湖及万山神池进行研究,并试图探讨此二无明显进,出水口之封闭型湖泊中,水体溶解氧垂直分布出现极大值之成因,此二次高海拔湖泊地处偏远,人烟罕至,因此较不受人为干扰,为研究湖泊水体中种种自然作用的良好对象,大鬼湖平均水深约14.8m,最深处约40m,降冬季外,水体均有分层现象,1988年夏季资料显示,水深16m以下水体趋于无氧状态,且于此深度以上的溶氧饱和值,均接近当地的大气饱和值（约78％）,经各种资料推断,此极大值的成因,除水团乃于春季时留下主溶氧值之外,应综合了季节增温效应下使表水向下温合的物理作用（尤其是山区明显的口,夜温差所引起的表水冷却向下混合作用,且其混合深度随季节增温而逐步变浅）,消耗溶解氧的生物作用及跃层存在等影响因素,而非单纯的物理或生物作用所造成,1991年4月万山神池观测资料显示,湖水平均深度约8m,最深可达14m ,其氧饱和程序分布在80％－104％之间,表水接近近当地饱和值（78％）,而于1.5m处往下增加,于2m处有溶解氧及和程度极大值,比当时大气饱合值高约20％,此极大值深度与叶绿素a极大值深度相吻合,主要由生物之光合作用造成,此外,于此深度pH值亦有明显增加现象,更证实了生物作用的存在,此二湖环境背景相似,且水泥中均出现溶氧极大值,但二者间极大值却有不同的成因,经曲两不同成因的比较,将可提供许多相关研究资料。     

Occurrence of Irgarol 1051 and its major metabolite in Maryland waters of Chesapeake Bay

Irgarol and its major metabolite (GS26575) were measured in Maryland waters of Chesapeake Bay: (1) in and near 10 marinas, a mainstem Bay site and two Severn River locations during a general survey in July and December of 2002; (2) at various sites in the Port Annapolis Marina and the Severn River area during March of 2002 before the boating season began; and (3) during July (peak boating season) in the same Port Annapolis Marina and Severn River sites area during both an ebb and flood tide. Irgarol concentrations ranged from 1.82 ng/l at the mid-Bay site to 585 ng/l in Port Annapolis marina during the July and December general survey. An Irgarol 90th centile of 239 ng/l was reported for the 10 marina sites, two Severn River sites and one mainstem site sampled during the general survey conducted in July and December. Temporal analysis of all pooled data showed that 90th centiles were over seven times higher in July when compared to December. A comparison of Irgarol concentrations at 12 sites in the Port Annapolis marina and Severn River area during both an ebb and flood tide in July showed no consistent trend with tidal cycle by site although significant reductions in concentrations were reported with distance from the three Port Annapolis marina sites. Ecological risk from Irgarol exposure was judged to be low for most Chesapeake Bay sites sampled. Possible exceptions were Port Annapolis marina, Severn River sites in close proximity to this marina and Chesapeake Harbor marina where Irgarol concentrations exceeded a conservative effects threshold during the peak boating season in July. Ecological risk from GS26575 exposure was low for all sites.