Storm track response to ocean fronts in a global high-resolution climate model

Synoptic atmospheric eddies are affected by lower tropospheric air-temperature gradients and by turbulent heat fluxes from the surface. In this study we examine how ocean fronts affect these quantities and hence the storm tracks. We focus on two midlatitude regions where ocean fronts lie close to the storm tracks: the north-west Atlantic and the Southern Ocean. An atmospheric climate model of reasonably high resolution (~50 km) is applied in a climate-length (60 year) simulation in order to obtain stable statistics. Simulations with frontal structure in the sea surface temperature (SST) in one of the regions are compared against simulations with globally smoothed SST. We show that in both regions the ocean fronts have a strong influence on the transient eddy heat and moisture fluxes, not just in the boundary layer, but also in the free troposphere. Local differences in these quantities between the simulations reach 20–40 % of the maximum values in the simulation with smoothed SST. Averaged over the entire region of the storm track over the ocean the corresponding differences are 10–20 %. The effect on the transient eddy meridional wind variance is strong in the boundary layer but relatively weak above that. The potential mechanisms by which the ocean fronts influence the storm tracks are discussed, and our results are compared against previous studies with regional models, Aquaplanet models, and coarse resolution coupled models.     

A Parameterization of Third-order Moments for the Dry Convective Boundary Layer

We describe one-dimensional (1D) simulations of the countergradient zone of mean potential temperature observed in the convective boundary layer (CBL). The method takes into account the third-order moments (TOMs) in a turbulent scheme of relatively low order, using the turbulent kinetic energy equation but without prognostic equations for other second-order moments. The countergradient term is formally linked to the third-order moments and , and a simple parameterization of these TOMs is proposed. It is validated for several cases of a dry CBL, using large-eddy simulations that have been realized from the MESO-NH model. The analysis of the simulations shows that TOMs are responsible for the inversion of the sign of in the higher part of the CBL, and budget analysis shows that the main terms responsible for turbulent fluxes and variances are now well reproduced.     

Effects of a Fence on Pollutant Dispersion in a Boundary Layer Exposed to a Rural-to-Urban Transition

Simultaneous particle-image velocimetry and laser-induced fluorescence combined with large-eddy simulations are used to investigate the flow and pollutant dispersion behaviour in a rural-to-urban roughness transition. The urban roughness is characterized by an array of cubical obstacles in an aligned arrangement. A plane fence is added one obstacle height h upstream of the urban roughness elements, with three different fence heights considered. A smooth-wall turbulent boundary layer with a depth of 10h is used as the approaching flow, and a passive tracer is released from a uniform line source 1h upstream of the fence. A shear layer is formed at the top of the fence, which increases in strength for the higher fence cases, resulting in a deeper internal boundary layer (IBL). It is found that the mean flow for the rural-to-urban transition can be described by means of a mixing-length model provided that the transitional effects are accounted for. The mixing-length formulation for sparse urban canopies, as found in the literature, is extended to take into account the blockage effect in dense canopies. Additionally, the average mean concentration field is found to scale with the IBL depth and the bulk velocity in the IBL.     

Pollutant Dispersion in Boundary Layers Exposed to Rural-to-Urban Transitions: Varying the Spanwise Length Scale of the Roughness

Both large-eddy simulations (LES) and water-tunnel experiments, using simultaneous stereoscopic particle image velocimetry and laser-induced fluorescence, have been used to investigate pollutant dispersion mechanisms in regions where the surface changes from rural to urban roughness. The urban roughness was characterized by an array of rectangular obstacles in an in-line arrangement. The streamwise length scale of the roughness was kept constant, while the spanwise length scale was varied by varying the obstacle aspect ratio l / h between 1 and 8, where l is the spanwise dimension of the obstacles and h is the height of the obstacles. Additionally, the case of two-dimensional roughness (riblets) was considered in LES. A smooth-wall turbulent boundary layer of depth 10h was used as the approaching flow, and a line source of passive tracer was placed 2h upstream of the urban canopy. The experimental and numerical results show good agreement, while minor discrepancies are readily explained. It is found that for \(l/h=2\) the drag induced by the urban canopy is largest of all considered cases, and is caused by a large-scale secondary flow. In addition, due to the roughness transition the vertical advective pollutant flux is the main ventilation mechanism in the first three streets. Furthermore, by means of linear stochastic estimation the mean flow structure is identified that is responsible for street-canyon ventilation for the sixth street and onwards. Moreover, it is shown that the vertical length scale of this structure increases with increasing aspect ratio of the obstacles in the canopy, while the streamwise length scale does not show a similar trend.     

Benthic respiration and stoichiometry of regenerated nutrients in lake sediments with Dreissena polymorpha

The effect of patchy colonies of the invasive zebra mussel (Dreissena polymorpha) on sedimentary processes was investigated in a mesotrophic lake (Plateliai Lake, Lithuania). Benthic fluxes of O₂, TCO₂, CH₄, Mn²⁺, Fe²⁺, N₂, the inorganic forms of N, Si and P and dissolved organic C and N were quantified by dark incubations of sediments cores, with and without D. polymorpha. Individual mussels also were incubated for metabolic measurements. Sediments with D. polymorpha had significantly higher O₂ and TCO₂ fluxes and displayed higher rates of denitrification. The presence of mussels also resulted in higher regeneration of P and N (mostly as ammonium) while the effect on Si was not significant. However, likely due to the low zebra mussel biomass (57.2 ± 25.3 g_SFDW m^?2), biodeposition has not changed the ratio between anaerobic and total respiration. Methane and reduced metals fluxes were in fact similar in the presence and absence of D. polymorpha. Incubations of mussels without sediments confirmed that bivalve metabolism was the main driver of benthic respiration and nutrient recycling. Nitrate production suggested the presence of nitrifiers associated with the molluscs. The main outcome of this study was that zebra mussels alter the quantity and the stoichiometry of nutrients regenerated by the benthic compartment. The enhancement of nitrogen loss via denitrification, by a factor of 1.5, was much less pronounced than the increase in ammonium recycling rate, stimulated by a factor of 33. Negligible PO₄ ^3? fluxes in bare sediments (?3.4 ± 6.8 μmol m^?2 h^?1) increased in the presence of mussels and considerable amounts of this nutrient (69.6 ± 29.4 μmol m^?2 h^?1) were mobilized to the water column. Further research should address other nutrient sources to the lake to verify whether altered rates and stoichiometry of benthic regeneration can affect primary producer community composition and activity.     

Differential herbivory of invasive algae by native fish in the Mediterranean Sea

The potential role of generalist herbivores to serve as a source of biotic resistance against algal invasion in marine ecosystems has been poorly examined. The present study investigates the capacity of Mediterranean herbivorous fishes to consume three of the most invasive seaweeds of the Western Mediterranean (Caulerpa racemosa, Lophocladia lallemandii and Womersleyella setacea) and examines vertical and temporal variations of such consumption. Our results show that although fish feed throughout the depth gradient examined (5–35 m), they concentrate in shallow waters, and can consume high amounts of C. racemosa. Such high ingestion of C. racemosa does not appear to be random, since this alga is consistently chosen when offered in pairs with several native species. Conversely, L. lallemandii and W. setacea are barely eaten by fish even though they can be very abundant in the field throughout the year. Our results suggest that fish could be an important controlling agent that has been overlooked in temperate marine invasions, and they may be able to provide certain resistance to C. racemosa invasion. In contrast, they are unlikely to exert any important control effects on L. lallemandii or W. setacea.     

An Efficient System for Creating Synthetic InSAR Images from Simulations

In this work we visualize tsunami and earthquake simulation results with graphics hardware acceleration. The rapid improvement in the computational power of graphics hardware and its programmability has made general computation on Graphics Processing Units (GPUs) very compelling. We generate Synthetic InSAR images using GPUs. Interference phenomena have formed the underlying theory for Interferometric Synthetic Aperture Radar (InSAR) in unveiling dynamical Earth movements. In our approach light path differences are defined by the surface values to be visualized. These path differences then modulate the lighting intensity to generate the interference patterns. We can interactively visualize surface deformation patterns by leveraging the computational power of GPUs. Our visualization method is applied to simulations of rupture fault displacements during the tsunamogenic earthquake events, which are vital to understanding the subsequent wave propagation. We also integrate the visualization results into Google Earth virtual globe to provide the geological context of the visualized regions.     

Orthogonal Variance Structures in Lake Water Quality Data and Their Use for Geo-chemical Classification of Dimictic, Glacial/Boreal Lakes

The accumulating volumes of data collected within environmental monitoring programs facilitate the use of exploratory statistical methods of data analysis as a supplement to traditional methods of characterizing lake water quality. When principal component analysis and multidimensional scaling are applied to a matrix containing approximately 24000 samples of lake water quality variables pH, alkalinity, conductivity, hardness, color, Secchi depth and total phosphorus concentration, it is found that the total matrix variance can be approximately reproduced in an orthogonal two-dimensional base with transformations of hardness and color as best principal component representatives. This base is suggested as an empirical lake classification standard where the variance structure of subset lake populations (such as single lakes) can be referenced to the water quality standard of the generic population. Since the principal axes of the base exclusively contain inorganic and organic related variables respectively, the combined inorganic/organic characteristics of the lake can be expressed with the hardness and color variables alone. With the data matrix being large enough to produce high significance levels, and with variable ranges wide enough to represent a majority of dimictic, glacial/boreal lakes, the analysis results should be valid in many lakes throughout the world.     

Experimental Investigation of Turbulent Momentum Transfer in a Neutral Boundary Layer over a Rough Surface

The turbulent characteristics of the neutral boundary layer developing over rough surfaces are not well predicted with operational weather-forecasting models. The problem is attributed to inadequate mixing-length models, to the anisotropy of the flow and to a lack of controlled experimental data against which to validate numerical studies. Therefore, in order to address directly the modelling difficulties for the development of a neutral boundary layer over rough surfaces, and to investigate the turbulent momentum transfer of such a layer, a set of hydraulic flume experiments were carried out. In the experiments, the mean and turbulent quantities were measured by a particle image velocimetry (PIV) technique. The measured velocity variances and fluxes \({(\overline{{u_{i}^{\prime}}{u_{j}^{\prime}}})}\) in longitudinal vertical planes allowed the vertical and longitudinal gradients (?/?z and ?/?x) of the mean and turbulent quantities (fluxes, variances and third-order moments) to be evaluated and the terms of the evolution equations for ?e/?t, \({\partial \overline{u^{\prime 2}}/\partial t}\), \({\partial \overline{w^{\prime 2}}/\partial t}\) and \({\partial \overline{{u^{\prime}}{w^{\prime}}}/\partial t}\) to be quantified, where e is the turbulent kinetic energy. The results show that the pressure-correlation terms allow the turbulent energy to be transferred equitably from \({\overline{{u^{\prime}}^{2}}}\) to \({\overline{{w^{\prime}}^{2}}}\). It appears that the repartition between the constitutive terms of the budget of e, \({\overline{{u^{\prime}}^{2}}}\), \({\overline{{w^{\prime}}^{2}}}\) and \({\overline{{u^{\prime}}{w^{\prime}}}}\) is not significantly affected by the development of the rough neutral boundary layer. For the whole evolution, the transfers of energy are governed by the same terms that are also very similar to the smooth-wall case. The PIV measurements also allowed the spatial integral scales to be computed directly and to be compared with the dissipative and mixing length scales, which were also computed from the data.     

Improved Parameter Resolution with Markov Chain Monte Carlo Simulation of Different Aquifer Tests

Hydraulic testing for aquifer characterization at contaminated sites often includes tests of short duration and of different types, such as slug tests and pumping tests, conducted at different phases of investigation. Tests conducted on a well cluster installed in a single aquifer can be combined in aggregate inverse analysis using an analytical model for groundwater flow near a test well. A genetic algorithm performs parallel search of the parameter space and provides starting parameter values for a Markov chain Monte Carlo simulation to estimate the parameter distribution. This sequence of inverse methods avoids guessing of the initial parameter vector and the often encountered difficult convergence of gradient-based methods and estimates the parameter covariance matrix from a distribution rather than from a single point in the parameter space. Combination of different tests improves the resolution of the estimated aquifer properties and allows an assessment of the uniformity of the aquifer. Estimated parameter correlations and standard deviations are used as relative metrics to distinguish well resolved and poorly resolved parameters. The methodology is demonstrated on example field tests in unconfined and leaky aquifers.