The advection of high-resolution tracers by low-resolution winds

The usefulness of any simulation of atmospheric tracers using low-resolution winds relies on both the dominance of large spatial scales in the strain and time dependence that results in a cascade in tracer scales. Here, a quantitative study on the accuracy of such tracer studies is made using the contour advection technique. It is shown that, although contour stretching rates are very insensitive to the spatial truncation of the wind field, the displacement errors in filament position are sensitive. A knowledge of displacement characteristics is essential if Lagrangian simulations are to be used for the inference of airmass origin. A quantitative lower estimate is obtained for the tracer scale factor (TSF): the ratio of the smallest resolved scale in the advecting wind field to the smallest “trustworthy” scale in the tracer field. For a baroclinic wave life cycle the TSF = 6.1 ± 0.3 while for the Northern Hemisphere wintertime lower stratosphere the TSF = 5.5 ± 0.5, when using the most stringent definition of the trustworthy scale. The similarity in the TSF for the two flows is striking and an explanation is discussed in terms of the activity of potential vorticity (PV) filaments. Uncertainty in contour initialization is investigated for the stratospheric case. The effect of smoothing initial contours is to introduce a spinup time, after which wind field truncation errors take over from initialization errors (2–3 days). It is also shown that false detail from the proliferation of finescale filaments limits the useful lifetime of such contour advection simulations to 3σ−1 days, where σ is the filament thinning rate, unless filaments narrower than the trustworthy scale are removed by contour surgery. In addition, PV analysis error and diabatic effects are so strong that only PV filaments wider than 50 km are at all believable, even for very high-resolution winds. The minimum wind field resolution required to accurately simulate filaments down to the erosion scale in the stratosphere (given an initial contour) is estimated and the implications for the modeling of atmospheric chemistry are briefly discussed.

A closer look at chaotic advection in the stratosphere: part I: geometric structure

The relevance of chaotic advection to stratospheric mixing and transport is addressed in the context of (i) a numerical model of forced shallow-water flow on the sphere, and (ii) a middle-atmosphere general circulation model. It is argued that chaotic advection applies to both these models if there is suitable large-scale spatial structure in the velocity field and if the velocity field is temporally quasi-regular. This spatial structure is manifested in the form of “cat’s eyes” in the surf zone, such as are commonly seen in numerical simulations of Rossby wave critical layers; by analogy with the heteroclinic structure of a temporally aperiodic chaotic system the cat’s eyes may be thought of as an “organizing structure” for mixing and transport in the surf zone. When this organizing structure exists, Eulerian and Lagrangian autocorrelations of the velocity derivatives indicate that velocity derivatives decorrelate more rapidly along particle trajectories than at fixed spatial locations (i.e., the velocity field is temporally quasi-regular). This phenomenon is referred to as Lagrangian random strain.

A closer look at chaotic advection in the stratosphere: part II: statistical diagnostics

Statistical diagnostics of mixing and transport are computed for a numerical model of forced shallow-water flow on the sphere and a middle-atmosphere general circulation model. In particular, particle dispersion statistics, transport fluxes, Liapunov exponents (probability density functions and ensemble averages), and tracer concentration statistics are considered. It is shown that the behavior of the diagnostics is in accord with that of kinematic chaotic advection models so long as stochasticity is sufficiently weak. Comparisons with random-strain theory are made.

A traceable physical calibration of the vertical advection-diffusion equation for modelling ocean heat uptake

The classic vertical advection-diffusion (VAD) balance is a central concept in studying the ocean heat budget, in particular in simple climate models (SCMs). Here we present a new framework to calibrate the parameters of the VAD equation to the vertical ocean heat balance of two fully-coupled climate models that is traceable to the models’ circulation as well as to vertical mixing and diffusion processes. Based on temperature diagnostics, we derive an effective vertical velocity w∗ and turbulent diffusivity k∗ for each individual physical process. In steady-state, we find that the residual vertical velocity and diffusivity change sign in mid-depth, highlighting the different regional contributions of isopycnal and diapycnal diffusion in balancing the models’ residual advection and vertical mixing. We quantify the impacts of the time-evolution of the effective quantities under a transient 1%CO2 simulation and make the link to the parameters of currently employed SCMs.

Robust future changes in temperature variability under greenhouse gas forcing and the relationship with thermal advection

Recent temperature extremes have highlighted the importance of assessing projected changes in the variability of temperature as well as the mean. A large fraction of present day temperature variance is associated with thermal advection, as anomalous winds blow across the land-sea temperature contrast for instance. Models project robust heterogeneity in the 21st century warming pattern under greenhouse gas forcing, resulting in land-sea temperature contrasts increasing in summer and decreasing in winter, and the pole-to-equator temperature gradient weakening in winter. In this study, future monthly variability changes in the 17 member ensemble ESSENCE are assessed. In winter, variability in midlatitudes decreases while in very high latitudes and the tropics it increases. In summer, variability increases over most land areas and in the tropics, with decreasing variability in high latitude oceans. Multiple regression analysis is used to determine the contributions to variability changes from changing temperature gradients and circulation patterns. Thermal advection is found to be of particular importance in the northern hemisphere winter midlatitudes, where the change in mean state temperature gradients alone could account for over half the projected changes. Changes in thermal advection are also found to be important in summer in Europe and coastal areas, although less so than in winter. Comparison with CMIP5 data shows that the midlatitude changes in variability are robust across large regions, particularly high northern latitudes in winter and mid northern latitudes in summer.

Simulation results and applications of an advection bounded scheme to practical flows

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Simulation results and applications of an advection bounded scheme to practical flows

This paper reports experiments on the use of a recently introduced advection bounded upwinding scheme, namely TOPUS (Computers & Fluids 57 (2012) 208-224), for flows of practical interest. The numerical results are compared against analytical, numerical and experimental data and show good agreement with them. It is concluded that the TOPUS scheme is a competent, powerful and generic scheme for complex flow phenomena.

Flux and resident injection in gaseous advection experiments

The occurrence of gaseous pollutants in soils has stimulated many experimental activities, including forced ventilation in the field as well as laboratory transport experiments with gases. The dispersion coefficient in advective-dispersive gas phase transport is often dominated by molecular diffusion, which leads to a large overall dispersivity gamma. Under such conditions it is important to distinguish between flux and resident modes of solute injection and detection. The influence of the inlet type oil the macroscopic injection mode was tested in two series of column experiments with gases at different mean flow velocities nu. First we compared infinite resident and flux injections, and second, semi-infinite resident and flux injections. It is shown that the macroscopically apparent injection condition depends on the geometry of the inlet section. A reduction of the cross-sectional area of the inlet relative to that of the column is very effective in excluding the diffusive solute input, thus allowing us to use the solutions for a flux Injection also at rather low mean flow velocities nu. If the whole cross section of a column is exposed to a large reservoir like that of ambient air, a semi-infinite resident injection is established, which can be distinguished from a flux injection even at relatively high velocities nu, depending on the mechanical dispersivity of the porous medium.

Permeable coral reef sediment dissolution driven by elevated pCO2 and pore water advection

Ocean acidification (OA) is expected to drive the transition of coral reef ecosystems from net calcium carbonate (CaCO3) precipitating to net dissolving within the next century. Although permeable sediments represent the largest reservoir of CaCO3 in coral reefs, the dissolution of shallow CaCO3 sands under future pCO2 levels has not been measured under natural conditions. In situ, advective chamber incubations under elevated pCO2 (~800 µatm) shifted the sediments from net precipitating to net dissolving. Pore water advection more than doubled dissolution rates (1.10 g CaCO3/m**2/day) when compared to diffusive conditions (0.42 g CaCO3/m**2 /day). Sediment dissolution could reduce net ecosystem calcification rates of the Heron Island lagoon by 8% within the next century, which is equivalent to a 25% reduction in the global average calcification rate of coral lagoons. The dissolution of CaCO3 sediments needs to be taken into account in order to address how OA will impact the net accretion of coral reefs under future predicted increases in CO2.

Robust spatially resolved pressure measurements using MRI with novel buoyant advection-free preparations of stable microbubbles in polysaccharide gels

MRI of fluids containing lipid coated microbubbles has been shown to be an effective toot for measuring the local fluid pressure. However, the intrinsically buoyant nature of these microbubbles precludes lengthy measurements due to their vertical migration under gravity and pressure-induced coalescence. A novel preparation is presented which is shown to minimize both these effects for at least 25 min. By using a 2% polysaccharide gel base with a small concentration of glycerol and 1,2-distearoyl-sn-glycero-3-phosphocholine coated gas microbubbles, MR measurements are made for pressures between 0.95 and 1.44 bar. The signal drifts due to migration and amalgamation are shown to be minimized for such an experiment whilst yielding very high NMR sensitivities up to 38% signal change per bar.

Summertime Influences of Tidal Energy Advection on the Surface Energy Balance in a Mangrove Forest

Mangrove forests are ecosystems susceptible to changing water levels and temperatures due to climate change as well as perturbations resulting from tropical storms. Numerical models can be used to project mangrove forest responses to regional and global environmental changes, and the reliability of these models depends on surface energy balance closure. However, for tidal ecosystems, the surface energy balance is complex because the energy transport associated with tidal activity remains poorly understood. This study aimed to quantify impacts of tidal flows on energy dynamics within a mangrove ecosystem. To address the research objective, an intensive 10-day study was conducted in a mangrove forest located along the Shark River in the Everglades National Park, FL, USA. Forest–atmosphere turbulent exchanges of energy were quantified with an eddy covariance system installed on a 30-m-tall flux tower. Energy transport associated with tidal activity was calculated based on a coupled mass and energy balance approach. The mass balance included tidal flows and accumulation of water on the forest floor. The energy balance included temporal changes in enthalpy, resulting from tidal flows and temperature changes in the water column. By serving as a net sink or a source of available energy, flood waters reduced the impact of high radiational loads on the mangrove forest. Also, the regression slope of available energy versus sink terms increased from 0.730 to 0.754 and from 0.798 to 0.857, including total enthalpy change in the water column in the surface energy balance for 30-min periods and daily daytime sums, respectively. Results indicated that tidal inundation provides an important mechanism for heat removal and that tidal exchange should be considered in surface energy budgets of coastal ecosystems. Results also demonstrated the importance of including tidal energy advection in mangrove biophysical models that are used for predicting ecosystem response to changing climate and regional freshwater management practices.

Neodymium isotopic characterization of Ross Sea Bottom Water and its advection through the southern South Pacific, POLARSTERN cruise ANT-XXVI/2

Since the inception of the international GEOTRACES program, studies investigating the distribution of trace elements and their isotopes in the global ocean have significantly increased. In spite of this large-scale effort, the distribution of neodymium isotopes (143Nd/144Nd) and concentrations ([Nd]) in the high latitude south Pacific is still understudied. Here we report dissolved Nd isotopes and concentrations from 11 vertical water column profiles from the south Pacific between South America and New Zealand. Results suggest that Ross Sea Bottom Water (RSBW) is represented by an epsilon-Nd value of ~ -7, and is thus more radiogenic than Circumpolar Deep Water (epsilon-Nd ~ -8). RSBW and its characteristic epsilon-Nd signature can be traced far into the SE Pacific until progressive mixing with ambient Lower Circumpolar Deep water (LCDW) dilutes this signal north of the Antarctic Polar Front (APF). The SW-NE trending Pacific-Antarctic Ridge restricts the advection of RSBW into the SW Pacific, where bottom water density, salinity, and epsilon-Nd values of -9 indicate the presence of bottom waters of an origin different from the Ross Sea. Neodymium concentrations show low surface concentrations and a linear increase with depth north of the Polar Front. South of the APF, surface [Nd] is high and increases with depth but remains almost constant below ~1000 m. This vertical and spatial [Nd] pattern follows the southward shoaling density surfaces of the Southern Ocean frontal system and hence suggests supply of Nd to the upper ocean through upwelling of Nd-rich deep water. Low particle abundance dominated by reduced opal production and seasonal sea ice cover likely contributes to the maintenance of the high upper ocean [Nd] south of the APF. The reported data highlights the use of Nd isotopes as a water mass tracer in the Southern Ocean, with the potential for paleocenaographic reconstructions, and contributes to an improved understanding of Nd biogeochemistry.

Discontinuous Galerkin Methods for Advection-Diffusion-Reaction Problems on Anisotropically Refined Meshes

In this paper we consider the a posteriori and a priori error analysis of discontinuous Galerkin interior penalty methods for second-order partial differential equations with nonnegative characteristic form on anisotropically refined computational meshes. In particular, we discuss the question of error estimation for linear target functionals, such as the outflow flux and the local average of the solution. Based on our a posteriori error bound we design and implement the corresponding adaptive algorithm to ensure reliable and efficient control of the error in the prescribed functional to within a given tolerance. This involves exploiting both local isotropic and anisotropic mesh refinement. The theoretical results are illustrated by a series of numerical experiments.