Clogging in horizontal subsurface flow constructed wetlands

Horizontal Subsurface Flow Treatment Wetlands (HSSF TWs) are used by Severn Trent Water as a low-cost tertiary wastewater treatment for rural locations. Experience has shown that clogging is a major operational problem that reduces HSSF TW lifetime. Clogging is caused by an accumulation of secondary wastewater solids from upstream processes and decomposing leaf litter. Clogging occurs as a sludge layer where wastewater is loaded on the surface of the bed at the inlet. Severn Trent systems receive relatively high hydraulic loading rates, which causes overland flow and reduces the ability to mineralise surface sludge accumulations. A novel apparatus and method, the Aston Permeameter, was created to measure hydraulic conductivity in situ. Accuracy is ±30 %, which was considered adequate given that conductivity in clogged systems varies by several orders of magnitude. The Aston Permeameter was used to perform 20 separate tests on 13 different HSSF TWs in the UK and the US. The minimum conductivity measured was 0.03 m/d at Fenny Compton (compared with 5,000 m/d clean conductivity), which was caused by an accumulation of construction fines in one part of the bed. Most systems displayed a 2 to 3 order of magnitude variation in conductivity in each dimension. Statistically significant transverse variations in conductivity were found in 70% of the systems. Clogging at the inlet and outlet was generally highest where flow enters the influent distribution and exits the effluent collection system, respectively. Surface conductivity was lower in systems with dense vegetation because plant canopies reduce surface evapotranspiration and decelerate sludge mineralisation. An equation was derived to describe how the water table profile is influenced by overland flow, spatial variations in conductivity and clogging. The equation is calibrated using a single parameter, the Clog Factor (CF), which represents the equivalent loss of porosity that would reproduce measured conductivity according to the Kozeny-Carman Equation. The CF varies from 0 for ideal conditions to 1 for completely clogged conditions. Minimum CF was 0.54 for a system that had recently been refurbished, which represents the deviation from ideal conditions due to characteristics of non-ideal media such as particle size distribution and morphology. Maximum CF was 0.90 for a 15 year old system that exhibited sludge accumulation and overland flow across the majority of the bed. A Finite Element Model of a 15 m long HSSF TW was used to indicate how hydraulics and hydrodynamics vary as CF increases. It was found that as CF increases from 0.55 to 0.65 the subsurface wetted area increases, which causes mean hydraulic residence time to increase from 0.16 days to 0.18 days. As CF increases from 0.65 to 0.90, the extent of overland flow increases from 1.8 m to 13.1 m, which reduces hydraulic efficiency from 37 % to 12 % and reduces mean residence time to 0.08 days.

Thomas, Howard, O'Neal, Don and Ghertman, Michael (Eds). Strategy, structure and style. Chichester: John Wiley, 1997

Book review

Surface and subsurface Vickers indentation cracks in SiC, Si3N4, and sialon ceramics

The nature of subsurface cracks formed under and around Vickers hardness indentations is often assumed rather than identified. Subsurface cracks in four engineering ceramics are revealed using a penetrant technique, and flaw dimensions are recorded. The resulting data are used to investigate several aspects of indentation cracking, such as crack shape, functional relationships between indentation load and flaw dimensions, and the performance of indentation fracture toughness equations. An R curve is constructed for each of the materials. © 1995 The Institute of Materials.

A comparison of in situ constant and falling head permeameter tests to assess the distribution of clogging within horizontal subsurface flow constructed wetlands

Clogging is the main operational problem associated with horizontal subsurface flow constructed wetlands (HSSF CWs). The measurement of saturated hydraulic conductivity has proven to be a suitable technique to assess clogging within HSSF CWs. The vertical and horizontal distribution of hydraulic conductivity was assessed in two full-scale HSSF CWs by using two different in situ permeameter methods (falling head (FH) and constant head (CH) methods). Horizontal hydraulic conductivity profiles showed that both methods are correlated by a power function (FH= CH 0.7821, r 2=0.76) within the recorded range of hydraulic conductivities (0-70 m/day). However, the FH method provided lower values of hydraulic conductivity than the CH method (one to three times lower). Despite discrepancies between the magnitudes of reported readings, the relative distribution of clogging obtained via both methods was similar. Therefore, both methods are useful when exploring the general distribution of clogging and, specially, the assessment of clogged areas originated from preferential flow paths within full-scale HSSF CWs. Discrepancy between methods (either in magnitude and pattern) aroused from the vertical hydraulic conductivity profiles under highly clogged conditions. It is believed this can be attributed to procedural differences between the methods, such as the method of permeameter insertion (twisting versus hammering). Results from both methods suggest that clogging develops along the shortest distance between water input and output. Results also evidence that the design and maintenance of inlet distributors and outlet collectors appear to have a great influence on the pattern of clogging, and hence the asset lifetime of HSSF CWs. © Springer Science+Business Media B.V. 2011.

Clogging in subsurface-flow treatment wetlands:measurement, modeling and management

This paper reviews the state of the art in measuring, modeling, and managing clogging in subsurface-flow treatment wetlands. Methods for measuring in situ hydraulic conductivity in treatment wetlands are now available, which provide valuable insight into assessing and evaluating the extent of clogging. These results, paired with the information from more traditional approaches (e.g., tracer testing and composition of the clog matter) are being incorporated into the latest treatment wetland models. Recent finite element analysis models can now simulate clogging development in subsurface-flow treatment wetlands with reasonable accuracy. Various management strategies have been developed to extend the life of clogged treatment wetlands, including gravel excavation and/or washing, chemical treatment, and application of earthworms. These strategies are compared and available cost information is reported. © 2012 Elsevier Ltd.

Clogging in subsurface-flow treatment wetlands:occurrence and contributing factors

Clogging is a major operational and maintenance issue associated with the use of subsurface flow wetlands for wastewater treatment, and can ultimately limit the lifetime of the system. This review considers over two decades of accumulated knowledge regarding clogging in both vertical and horizontal subsurface flow treatment wetlands. The various physical, chemical and biological factors responsible for clogging are identified and discussed. The occurrence of clogging is placed into the context of various design and operational parameters such as wastewater characteristics, upstream treatment processes, intermittent or continuous operation, influent distribution, and media type. This information is then used to describe how clogging develops within, and subsequently impacts, common variants of subsurface flow treatment wetland typically used in the U.S., U.K., France and Germany. Comparison of these systems emphasized that both hydraulic loading rate and solids loading rate need to be considered when designing systems to operate robustly, i.e. hydraulic overloading makes horizontal-flow tertiary treatment systems in the U.K. more susceptible to clogging problems than vertical-flow primary treatment systems in France. Future research should focus on elucidating the underlying mechanisms of clogging as they relate to the design, operation, and maintenance of subsurface flow treatment wetlands. © 2010 Elsevier B.V.

A finite element approach to modelling the hydrological regime in horizontal subsurface flow constructed wetlands for wastewater treatment

A Finite Element Analysis (FEA) model is used to explore the relationship between clogging and hydraulics that occurs in Horizontal Subsurface Flow Treatment Wetlands (HSSF TWs) in the United Kingdom (UK). Clogging is assumed to be caused by particle transport and an existing single collector efficiency model is implemented to describe this behaviour. The flow model was validated against HSSF TW survey results obtained from the literature. The model successfully simulated the influence of overland flow on hydrodynamics, and the interaction between vertical flow through the low permeability surface layer and the horizontal flow of the saturated water table. The clogging model described the development of clogging within the system but under-predicted the extent of clogging which occurred over 15 years. This is because important clogging mechanisms were not considered by the model, such as biomass growth and vegetation establishment. The model showed the usefulness of FEA for linking hydraulic and clogging phenomenon in HSSF TWs and could be extended to include treatment processes. © 2011 Springer Science+Business Media B.V.

The L-31E Surface Water Rediversion Project Final Report: Implementation, Results, and Recommendations

Throughout the Biscayne Bay watershed, existing coastal wetland communities have been cut off from sheet flow for decades. With the expectation that reconnection of these wetlands to upstream water sources would alter existing hydrologic conditions and recreate a more natural sheet flow to Biscayne National Park, a demonstration project on freshwater rediversion was undertaken. The objectives of the project were to document the effects of freshwater diversion on: (a) swamp and nearshore water chemistry and hydrology; (b) soil development processes; (c) macrophyte and benthic algal community composition, structure and production; (d) abundance of epiphytic and epibenthic invertebrates; (e) zonation, production, and phenology of primary producers in the nearshore environment, and (f) exchanges of nutrients and particulates between nearshore and mangrove ecosystems.

Integrated surface-ground water modeling in wetlands with improved methods to simulate vegetative resistance to flow

In topographically flat wetlands, where shallow water table and conductive soil may develop as a result of wet and dry seasons, the connection between surface water and groundwater is not only present, but perhaps the key factor dominating the magnitude and direction of water flux. Due to their complex characteristics, modeling waterflow through wetlands using more realistic process formulations (integrated surface-ground water and vegetative resistance) is an actual necessity. This dissertation focused on developing an integrated surface – subsurface hydrologic simulation numerical model by programming and testing the coupling of the USGS MODFLOW-2005 Groundwater Flow Process (GWF) package (USGS, 2005) with the 2D surface water routing model: FLO-2D (O’Brien et al., 1993). The coupling included the necessary procedures to numerically integrate and verify both models as a single computational software system that will heretofore be referred to as WHIMFLO-2D (Wetlands Hydrology Integrated Model). An improved physical formulation of flow resistance through vegetation in shallow waters based on the concept of drag force was also implemented for the simulations of floodplains, while the use of the classical methods (e.g., Manning, Chezy, Darcy-Weisbach) to calculate flow resistance has been maintained for the canals and deeper waters. A preliminary demonstration exercise WHIMFLO-2D in an existing field site was developed for the Loxahatchee Impoundment Landscape Assessment (LILA), an 80 acre area, located at the Arthur R. Marshall Loxahatchee National Wild Life Refuge in Boynton Beach, Florida. After applying a number of simplifying assumptions, results have illustrated the ability of the model to simulate the hydrology of a wetland. In this illustrative case, a comparison between measured and simulated stages level showed an average error of 0.31% with a maximum error of 2.8%. Comparison of measured and simulated groundwater head levels showed an average error of 0.18% with a maximum of 2.9%. The coupling of FLO-2D model with MODFLOW-2005 model and the incorporation of the dynamic effect of flow resistance due to vegetation performed in the new modeling tool WHIMFLO-2D is an important contribution to the field of numerical modeling of hydrologic flow in wetlands.

Tracking rates of ecotone migration due to salt-water encroachment using fossil mollusks in coastal South Florida

We determined the rate of migration of coastal vegetation zones in response to salt-water encroachment through paleoecological analysis of mollusks in 36 sediment cores taken along transects perpendicular to the coast in a 5.5 km2 band of coastal wetlands in southeast Florida. Five vegetation zones, separated by distinct ecotones, included freshwater swamp forest, freshwater marsh, and dwarf, transitional and fringing mangrove forest. Vegetation composition, soil depth and organic matter content, porewater salinity and the contemporary mollusk community were determined at 226 sites to establish the salinity preferences of the mollusk fauna. Calibration models allowed accurate inference of salinity and vegetation type from fossil mollusk assemblages in chronologically calibrated sediments. Most sediments were shallow (20–130 cm) permitting coarse-scale temporal inferences for three zones: an upper peat layer (zone 1) representing the last 30–70 years, a mixed peat-marl layer (zone 2) representing the previous ca. 150–250 years and a basal section (zone 3) of ranging from 310 to 2990 YBP. Modern peat accretion rates averaged 3.1 mm yr)1 while subsurface marl accreted more slowly at 0.8 mm yr)1. Salinity and vegetation type for zone 1 show a steep gradient with freshwater communities being confined west of a north–south drainage canal constructed in 1960. Inferences for zone 2 (pre-drainage) suggest that freshwater marshes and associated forest units covered 90% of the area, with mangrove forests only present along the peripheral coastline. During the entire pre-drainage history, salinity in the entire area was maintained below a mean of 2 ppt and only small pockets of mangroves were present; currently, salinity averages 13.2 ppt and mangroves occupy 95% of the wetland. Over 3 km2 of freshwater wetland vegetation type have been lost from this basin due to salt-water encroachment, estimated from the mollusk-inferred migration rate of freshwater vegetation of 3.1 m yr)1 for the last 70 years (compared to 0.14 m yr)1 for the pre-drainage period). This rapid rate of encroachment is driven by sea-level rise and freshwater diversion. Plans for rehydrating these basins with freshwater will require high-magnitude re-diversion to counteract locally high rates of sea-level rise.

Rapid responses of vegetation to hydrological changes in Taylor Slough, Everglades National Park, Florida, USA

We analyzed the dynamics of freshwater marsh vegetation of Taylor Slough in eastern Everglades National Park for the 1979 to 2003 period, focusing on cover of individual plant species and on cover and composition of marsh communities in areas potentially influenced by a canal pump station (‘‘S332’’) and its successor station (‘‘S332D’’). Vegetation change analysis incorporated the hydrologic record at these sites for three intervals: pre-S332 (1961–1980), S332 (1980–1999), post-S332 (1999–2002). During S332 and post-S332 intervals, water level in Taylor Slough was affected by operations of S332 and S332D. To relate vegetation change to plot-level hydrological conditions in Taylor Slough, we developed a weighted averaging regression and calibration model (WA) using data from the marl prairies of Everglades National Park and Big Cypress National Preserve. We examined vegetation pattern along five transects. Transects 1–3 were established in 1979 south of the water delivery structures, and were influenced by their operations. Transects 4 and 5 were established in 1997, the latter west of these structures and possibly under their influence. Transect 4 was established in the northern drainage basin of Taylor Slough, beyond the likely zones of influence of S332 and S332D. The composition of all three southern transects changed similarly after 1979. Where muhly grass (Muhlenbergia capillaris var. filipes) was once dominant, sawgrass (Cladium jamaicense), replaced it, while where sawgrass initially predominated, hydric species such as spikerush (Eleocharis cellulosa Torr.) overtook it. Most of the changes in species dominance in Transects 1–3 occurred after 1992, were mostly in place by 1995–1996, and continued through 1999, indicating how rapidly vegetation in seasonal Everglades marshes can respond to hydrological modifications. During the post-S332 period, these long-term trends began reversing. In the two northern transects, total cover and dominance of both muhly grass and sawgrass increased from 1997 to 2003. Thus, during the 1990’s, vegetation composition south of S332 became more like that of long hydroperiod marshes, but afterward it partially returned to its 1979 condition, i.e., a community characteristic of less prolonged flooding. In contrast, the vegetation change along the two northern transects since 1997 showed little relationship to hydrologic status.

Hydrologic measurements and implications for tree island formation within Everglades National Park

Tree islands in the Shark River Slough of the Everglades National Park (ENP), in the southern state of Florida in the United States, are part of a wetland system of densely vegetated ridges interspersed within relatively open sloughs. Human alteration of this system has had dramatic negative effects on the landscape of the region and restoration efforts will require adjusting the hydrology of the region to assure the preservation of these important ecologic features. The primary objectives of this study were to document the hydrology in the vicinity of tree islands in ENP by measuring velocities in time and space and by characterizing suspended sediments. The results of such measurements were interpreted with respect to factors that may limit tree island growth. The measurements were conducted in the vicinity of three tree islands known as Black Hammock (BH), Gumbo Limbo (GL), and an unnamed island that was named for this study as Satin Leaf (SL). Acoustical Doppler Velocity (ADV) meters were used for measuring the low velocities of the Everglades water flow. Properties of suspended sediments were characterized through measurements of particle size distribution, turbidity, concentration and particle density. Mean velocities observed at each of the tree islands varied from 0.9 to 1.4 cm/s. Slightly higher mean velocities were observed during the wet season (1.2–1.6 cm/s) versus the dry season (0.8–1.3 cm/s). Maximum velocities of more than 4 cm/s were measured in areas of Cladium jamaicense die-off and at the hardwood hammock (head) of the islands. At the island’s head, water is channelized around obstructions such as tree trunks in relatively rapid flow, which may limit the lateral extent of tree island growth. Channelization is facilitated by shade from the tree canopy, which limits the growth of underwater vegetation thereby minimizing the resistance to flow and limiting sediment deposition. Suspended sediment concentrations were low (0.5–1.5 mg/L) at all study sites and were primarily of organic origin. The mean particle size of the suspended sediments was 3 μm with a distribution that was exponential. Critical velocities needed to cause re-suspension of these particles were estimated to be above the actual velocities observed. Sediment transport within the water column appears to be at a near steady state during the conditions evaluated with low rates of sediment loss balanced by presumably the release of equivalent quantities of particles of organic origin. Existing hydrologic conditions do not appear to transport sufficient suspended sediments to result in the formation of tree islands. Of interest would be to collect hydrologic and sediment transport data during extreme hydrologic events to determine if enough sediment is transported under these conditions to promote sufficient sediment accumulations.

Developing performance measures of mangrove wetlands using simulation models of hydrology, nutrient biogeochemistry, and community dynamics

The goal of mangrove restoration projects should be to improve community structure and ecosystem function of degraded coastal landscapes. This requires the ability to forecast how mangrove structure and function will respond to prescribed changes in site conditions including hydrology, topography, and geophysical energies. There are global, regional, and local factors that can explain gradients of regulators (e.g., salinity, sulfides), resources (nutrients, light, water), and hydroperiod (frequency, duration of flooding) that collectively account for stressors that result in diverse patterns of mangrove properties across a variety of environmental settings. Simulation models of hydrology, nutrient biogeochemistry, and vegetation dynamics have been developed to forecast patterns in mangroves in the Florida Coastal Everglades. These models provide insight to mangrove response to specific restoration alternatives, testing causal mechanisms of system degradation. We propose that these models can also assist in selecting performance measures for monitoring programs that evaluate project effectiveness. This selection process in turn improves model development and calibration for forecasting mangrove response to restoration alternatives. Hydrologic performance measures include soil regulators, particularly soil salinity, surface topography of mangrove landscape, and hydroperiod, including both the frequency and duration of flooding. Estuarine performance measures should include salinity of the bay, tidal amplitude, and conditions of fresh water discharge (included in the salinity value). The most important performance measures from the mangrove biogeochemistry model should include soil resources (bulk density, total nitrogen, and phosphorus) and soil accretion. Mangrove ecology performance measures should include forest dimension analysis (transects and/or plots), sapling recruitment, leaf area index, and faunal relationships. Estuarine ecology performance measures should include the habitat function of mangroves, which can be evaluated with growth rate of key species, habitat suitability analysis, isotope abundance of indicator species, and bird census. The list of performance measures can be modified according to the model output that is used to define the scientific goals during the restoration planning process that reflect specific goals of the project.

Development of an integrated surface and subsurface model of Everglades National Park

An integrated surface-subsurface hydrological model of Everglades National Park (ENP) was developed using MIKE SHE and MIKE 11 modeling software. The model has a resolution of 400 meters, covers approximately 1050 square miles of ENP, includes 110 miles of drainage canals with a variety of hydraulic structures, and processes hydrological information, such as evapotranspiration, precipitation, groundwater levels, canal discharges and levels, and operational schedules. Calibration was based on time series and probability of exceedance for water levels and discharges in the years 1987 through 1997. Model verification was then completed for the period of 1998 through 2005. Parameter sensitivity in uncertainty analysis showed that the model was most sensitive to the hydraulic conductivity of the regional Surficial Aquifer System, the Manning's roughness coefficient, and the leakage coefficient, which defines the canal-subsurface interaction. The model offers an enhanced predictive capability, compared to other models currently available, to simulate the flow regime in ENP and to forecast the impact of topography, water flows, and modifying operation schedules.