Numerical modeling of the effects of water flow,sediment transport and vegetation growth on the spatiotemporal patterning of the ridge and slough landscape of the Everglades wetland

A numerical model has been developed to simulate the spatiotemporal patterning of the ridge and slough landscape in wetlands, characterized by crests (ridges) and valleys (sloughs) that are elongated parallel to the direction of water flow. The model formulation consists of governing equations for integrated surface water and groundwater flow, sediment transport, and soil accretion, as well as litter production by vegetation growth. The model simulations show how the spatial pattern self-organizes over time with the generation of ridges and sloughs through sediment deposition and erosion driven by the water flow field. The spatial and temporal distributions of the water depth, flow rates and sediment transport processes are caused by differential flow due to vegetation and topography heterogeneities. The model was parameterized with values that are representative of the Everglades wetland in the southern portion of the Florida peninsula in the USA. Model simulation sensitivity was tested with respect to numerical grid size, lateral vegetation growth and the rate of litter production. The characteristic wavelengths of the pattern in the directions along and perpendicular to flow that are simulated with this model develop over time into ridge and slough shapes that resemble field observations. Also, the simulated elevation differences between the ridges and sloughs are of the same order of those typically found in the field. The width of ridges and sloughs was found to be controlled by a lateral vegetation growth distance parameter in a simplified formulation of vegetation growth, which complements earlier modeling results in which a differential peat accretion mechanism alone did not reproduce observations of ridge and slough lateral wavelengths. The results of this work suggest that ridge and slough patterning occurs as a result of vegetation's ability to grow laterally, enhancing sediment deposition in ridge areas, balanced by increased sediment erosion in slough areas to satisfy flow continuity. The interplay between sediment transport, water flow and vegetation and soil dynamic processes needs to be explored further through detailed field experiments, using a model formulation such as the one developed in this work to guide data collection and interpretation. This should be one of the focus areas of future investigations of pattern formation and stability in ridge and slough areas.     

Optimization of photovoltaic solar power plant locations in northern Chile

The optimization of photovoltaic solar power plants location in Atacama Desert, Chile, is presented in this study. The study considers three objectives: (1) Find sites with the highest solar energy potential, (2) determine sites with the least impact on the environment, and (3) locate the areas which produce small social impact. To solve this task, multi-criteria decision analyses (MCDAs) such as analytical hierarchy process and ordered weighted averaging were applied in a GIS environment. In addition, survey results of social impacts were analyzed and included into the decision-making process, including landscape values. The most suitable sites for solar energy projects were found near roads and power lines throughout the study area. Large suitable areas were found also from central valley from Arica and Parinacota to the north edge of Atacama region. In Atacama region, most suitable sites were found in the Andes. On the contrary, Andes were also found to have high environmental values and scenically valuable landscapes. Moderate and low suitability were found on the coast, especially in Atacama region. Factors such as slope and distance to power lines and roads influenced largely the sensitivity analysis. Area of high suitability increased by 15% when distance to roads was excluded and 18% when distance to power lines or slope was removed. MCDA-GIS method was found to be useful and applicable to the optimization of solar power plant locations in northern Chile.     

Final stages of N-body star cluster encounters

We performed numerical simulations of star cluster encounters with Hernquist's treecode on a CRAY YMP-2E computer. We used different initial conditions (relative positions and velocities, cluster sizes, masses and concentration degrees) with the total number of particles per simulation ranging from 4608 to 20 480. Long-term interaction stages (up to 1 Gyr) when the pair coalesces into a single cluster are compared with isolated LMC clusters. Evidence is found that, when seen in a favourable plane, these resulting clusters show elliptical shapes as a result of the disruption of one of the companions. These elliptical shapes are essentially time-independent, but they do depend on the initial structural parameters of the pair components. We also analysed the fraction of stars that are ejected to the field by the interaction. We found that this fraction can be almost 50 per cent for the disrupted cluster. These simulations can represent a possible mechanism with which to explain the ellipticity observed in several star clusters in the Magellanic Clouds.     

Biological,physical and chemical properties at the Subtropical Shelf Front Zone in the SW Atlantic Continental Shelf

The physical aspects of the Subtropical Shelf Front (STSF) for the Southwest Atlantic Continental Shelf were previously described. However, only scarce data on the biology of the front is available in the literature. The main goal of this paper is to describe the physical, chemical and biological properties of the STSF found in winter 2003 and summer 2004. A cross-section was established at the historically determined location of the STSF. Nine stations were sampled in winter and seven in summer. Each section included a series of conductivity-temperature-depth (CTD) stations where water samples from selected depths were filtered for nutrient determination. Surface samples were taken for chlorophyll a (Chl-a) determination and plankton net tows carried out above and below the pycnocline. Results revealed that winter was marked by an inner-shelf salinity front and that the STSF was located on the mid-shelf. The low salinity waters in the inner-shelf indicated a strong influence of freshwater, with high silicate (72 μM), suspended matter (45 mg l⁻¹), phosphate (2.70 μM) and low nitrate (1.0 μM) levels. Total dissolved nitrogen was relatively high (22.98 μM), probably due to the elevated levels of organic compound contribution close to the continental margin. Surface Chl-a concentration decreased from coastal well-mixed waters, where values up to 8.0 mg m⁻³ were registered, to offshore waters. Towards the open ocean, high subsurface nutrients values were observed, probably associated to South Atlantic Central Waters (SACW). Zooplankton and ichthyoplankton abundance followed the same trend; three different groups associated to the inner-, mid- and outer-shelf region were identified. During summer, diluted waters extended over the shelf to join the STSF in the upper layer; the concentration of inorganic nutrients decreased in shallow waters; however, high values were observed between 40 and 60 m and in deep offshore waters. Surface Chl-a ranged 0.07–1.5 mg m⁻³; winter levels were higher. Three groups of zoo and ichthyoplankton, separated by the STSF, were also identified. Results of the study performed suggest that the influence of freshwater was stronger during winter and that abundance distribution of Chl-a, copepods and ichthyoplankton was related to the Plata Plume Waters (PPW), rather than to the presence of the STSF. During summer, when the presence of freshwater decreases, plankton interactions seem to take place in the STSF.     

Allogenic and autogenic effects on mangrove dynamics from the Ceará Mirim River,north‐eastern Brazil,during the middle and late Holocene

It is possible that climate changes and sea level fluctuations (allogenic processes) are and will cause major changes in mangrove dynamics. However, other driving forces may be significantly affecting this system. Distinguishing allogenic and autogenic influence on mangroves is a challenging question, because mechanisms related to the natural dynamics of depositional environments (autogenic processes) have strong influences on the establishment and degradation of mangroves. Thus, impacts on mangroves caused by autogenic processes may be erroneously attributed to allogenic mechanisms. Therefore, it is imperative to identify the ‘fingerprint’ of global changes in modern mangrove dynamics. In order to characterize the influence of these forces on mangroves, this work has used geomorphology and vegetation maps integrated with sedimentological and palynological data, radiocarbon dating, as well as δ¹³C, δ¹⁵N and C/N from sedimentary organic matter. The inter‐proxy analyses reveal an estuarine influence with mangrove development along the Ceará Mirim River, north‐eastern Brazil, since ~6920 cal yr bp , after the post‐glacial sea level rise. Relative sea level (RSL) has been stable during the middle and late Holocene. Mangrove establishment along this fluvial valley begins at about 6920 cal yr bp , caused by the sea‐level stabilization, an allogenic influence. However, after its establishment, wetland dynamics were mainly controlled by autogenic factors, related to channel migrations, instead of allogenic process. Some influence of sea‐level and climate changes on mangrove dynamics in this estuarine channel have been weakened by more intense tidal channels activities. Therefore, the expansion and contraction of mangrove areas along the estuary of the Ceará Mirim River since 6920 cal yr bp has been mainly influenced by channel dynamics that regulate the accretion and erosion of mangrove substrates. Copyright © 2018 John Wiley & Sons, Ltd.     

Sediment geochemistry of the urban Lake Paulo Gorski

The objective of this study was to evaluate the concentration and distribution of heavy metals in the sediments of Paulo Gorski Lake, as well as the metals’ bioavailability and potential ecological risk, and to define the anthropogenic and natural heavy metal contributions to the lake. The chemical elements calcium (Ca), cadmium (Cd), cobalt (Co), chromium (Cr), copper (Cu), iron (Fe), potassium (K), magnesium (Mg), manganese (Mn), sodium (Na), nickel (Ni), lead (Pb), and zinc (Zn) were quantified by flame atomic absorption spectrophotometry with two extraction methods to quantify the bioavailable and non-bioavailable fractions. The data were evaluated using multivariate statistics and sediment quality indices. All sediment collection points (S1, S2, S3, S4 and S5) are different in terms of the concentration of heavy metals, except for S4 and S5, which were statistically equal. The bioavailable fraction of the elements in the sediment follows the sequence Pb>Cu>Mn>Zn>Ni>Cr>phosphorus (P) for all points. The elements Co, Cr, Pb, and Zn showed moderate to considerable contamination at all points. Only points S3 and S5 had moderate ecological risk. Urbanization has been affecting Paulo Gorski Lake via the input of chemical elements, especially Co and Pb. The points most affected by heavy metal contamination are S3 and S5 when the sedimentological sensitivity factor is considered. The lake has high hydrodynamics, causing some of the contaminants that enter the system to leave it, leading to potential negative impacts downstream.     

Three‐dimensional numerical modeling of the Bulle effect: the nonlinear distribution of near‐bed sediment at fluvial diversions

The Bulle effect is a phenomenon in which a disproportionately higher amount of near‐bed sediment load at a fluvial diversion moves into the diverted channel, even for cases in which the proportion of water (with respect to the main flow) entering the diversion channel is relatively small. This phenomenon has wide‐ranging implications for both engineered and natural systems: from efficient design of channels to redirect water and sediment for reclaiming sinking deltas, designing navigational channels that do not need frequent dredging, to morphological evolution of river bifurcations. The first ever, and one of the most extensive set of experiments conducted to explore this phenomenon, were conducted by Bulle in 1926 . In the current study the experiments conducted by Bulle have been simulated using an open‐source, free‐surface finite‐element‐based hydrodynamic solver. The main objectives were to explore to what extent the complex phenomenon of the Bulle effect at the scale of a laboratory experiment can be simulated accurately using Reynolds‐averaged Navier–Stokes (RANS)‐based hydrodynamic solver, and to understand the details of the hydrodynamics that Bulle could not analyze through his experiments. The hydrodynamics captured by the simulations were found to match the observations made by Bulle through his experiments, and the distributions of sediment at the diversion predicted by the numerical simulations were found to match the general trend observed in the laboratory experiments. The results from the numerical simulations were also compared with existing one‐dimensional models for sediment distribution at bifurcations, and the three‐dimensional numerical model was found to perform appreciably better. This is expected due to the complex flow features at the diversion, which can only be captured satisfactorily using a three‐dimensional hydrodynamic model. Copyright © 2017 John Wiley & Sons, Ltd.     

Effects of levee cover strength on flood mapping in the case of levee breach due to overtopping

The reliability of a levee system is a crucial factor in flood risk management. In this study we present a probabilistic methodology to assess the effects of levee cover strength on levee failure probability, triggering time, flood propagation and consequent impacts on population and assets. A method for determining fragility curves is used in combination with the results of a one-dimensional hydrodynamic model to estimate the conditional probability of levee failure in each river section. Then, a levee breach model is applied to calculate the possible flood hydrographs, and for each breach scenario a two-dimensional hydrodynamic model is used to estimate flood hazard (flood extent and timing, maximum water depths) and flood impacts (economic damage and affected population) in the areas at risk along the river reach. We show an application for levee overtopping and different flood scenarios for a 98 km reach of the lower Po River in Italy. The results show how different design solutions for the levee cover can influence the probability of levee failure and the consequent flood scenarios. In particular, good grass cover strength can significantly delay levee failure and reduce maximum flood depths in the flood-prone areas, thus helping the implementation of flood risk management actions.

EDITOR D. Koutsoyiannis

ASSOCIATE EDITOR A. Viglione     

Does Poisson’s downward continuation give physically meaningful results?

The downward continuation (DWC) of the gravity anomalies from the Earth’s surface to the geoid is still probably the most problematic step in the precise geoid determination. It is this step that motivates the quasi-geoid users to opt for Molodenskij’s rather than Stokes’s theory. The reason for this is that the DWC is perceived as suffering from two major flaws: first, a physically meaningful DWC technique requires the knowledge of the irregular topographical density; second, the Poisson DWC, which is the only physically meaningful technique we know, presents itself mathematically in the form of Fredholm integral equation of the 1st kind. As Fredholm integral equations are often numerically ill-conditioned, this makes some people believe that the DWC problem is physically ill-posed. According to a revered French mathematician Hadamard, the DWC problem is physically well-posed and as such gives always a finite and unique solution. The necessity of knowing the topographical density is, of course, a real problem but one that is being solved with an ever increasing accuracy; so sooner or later it will allow us to determine the geoid with the centimetre accuracy.     

Length scales and statistical characteristics of outer bank roughness for large elongate meander bends: The influence of bank material properties,floodplain vegetation and flow inundation

This article explores the length scales and statistical characteristics of form roughness along the outer banks of two elongate bends on a large meandering river through investigation of topographic variability of the bank face. The analysis also examines how roughness varies over the vertical height of the banks and when the banks are exposed subaerially and inundated during flood stage. Detailed data on the topography of the outer banks were obtained subaerially using terrestrial LiDAR during low flow conditions and subaqueously using multibeam echo sounding (MBES) during near‐bankfull conditions. The contributions of various length scales of topographic irregularity to roughness for subaerial conditions were evaluated for different elevation contours on the bank faces using Hilbert–Huang Transform (HHT) spectral analysis. Statistical characteristics for discrete areas on the bank faces were determined by calculating the root‐mean‐square of normal distances from a triangulated irregular network (TIN) surface. Results of the HHT analysis show that the characteristics of roughness along bank faces composed primarily of non‐cohesive sediment, and eroding into cropland, vary with bank elevation and exhibit a dominant range of roughness length scales (~15–50 m). However, bank faces composed predominantly of cohesive material and carved into a forested floodplain have relatively uniform topographic roughness characteristics over the vertical extent of the bank face and do not exhibit a dominant roughness length scale or range of length scales. Additionally, comparison between local surface roughness for subaerial versus subaqueous conditions shows that roughness decreases considerably when the banks are submerged, most likely because of the removal of vegetation and eradication of small‐scale erosional features in non‐cohesive bank materials by flow along the bank face. Thus, roughness appears to be linked to the hydraulic conditions affecting the bank, at least relative to conditions that develop when banks are exposed subaerially. Copyright © 2017 John Wiley & Sons, Ltd.