Convective Front Research Articles

A Novel Method for Diagnosing Land–Atmosphere Coupling Sensitivity in a Single-Column Model

Abstract The response of boundary layer properties and cloudiness to changes in surface evaporative fraction (EF) is investigated in a single-column model to quantify the locally coupled impact of subgrid surface variations on the atmosphere during summer. Sensitive coupling days are defined when the model atmosphere exhibits large variations across a range of EFs centered on the analyzed value. Coupling sensitivity exists as both positive feedback (cloudiness increases with EF) and negative feedback (clouds increase with decreasing EF) regimes. The positive regime manifests in shallow convection situations, which are capped by a strengthened inversion and subsidence, restricting the vertical extent of convection to just above the boundary layer. Surfaces with larger EF (greater surface latent heat flux) can inject more moisture into the vertically confined system, lowering the cloud base and an increasing cloud liquid water path (LWP). Negative feedback regimes tend to manifest when large-scale deep convection, such as from mesoscale convective systems and fronts, is advected through the domain, where convection strengthens over surfaces with a lower EF (greater surface sensible heat flux). The invigoration of these systems by the land surface leads to an increase in LWP through strengthened updrafts and stronger coupling between the boundary layer and the free atmosphere. These results apply in the absence of heterogeneity-induced mesoscale circulations, providing a one-dimensional dynamical perspective on the effect of surface heterogeneity. This study provides a framework of intermediate complexity, lying between parcel theory and high-resolution coupled land–atmosphere modeling, and therefore isolates the relevant first-order processes in land–atmosphere interactions. Significance Statement Cloud formation, distribution, and other properties may be sensitive to heterogeneous surfaces depending on the strength and location of such heterogeneities and the background atmospheric state. This may drive differences in the cloud population depending on which part of the domain one is located. This may also lead to mesoscale circulations, which may strengthen or weaken this effect. Currently, climate models act on scales (∼100 km) that are too large to explicitly represent these processes, which are strongest at smaller scales (around 5–40 km). Therefore, subgrid-scale heterogeneity is neglected, and any predictability and model fidelity it may provide is lost. We use a simple model to diagnose sensitivity of the local atmosphere to surface variations meant to represent possible subgrid heterogeneity, providing a first-order estimate of its effect. We conclude that preferentially sensitive atmospheric states exist that lead to positive and/or negative feedback between land and atmosphere. This information is valuable to future climate model parameterizations aimed at improving the representation of these feedbacks.

Interacting density fronts in saturated brines cooled from above

When a saturated brine layer is cooled from above, both a convective temperature front as well as a front of sedimenting salt crystals can form. We employ direct numerical simulations to investigate the evolution and interaction of these two density fronts. Depending on the ratio of the temperature front velocity and the crystal settling velocity, which is governed by a dimensionless parameter in the form of a Rayleigh number, we find that either two separate fronts exist for all times, two initially separate fronts combine into a single front after some time or a single front exists at all times. We furthermore propose approximate scaling laws for the propagation of the thermal and crystal fronts in each regime and compare them with the simulation data, with generally good agreement.

Convective Density Current Circulations That Modulated Meso-γ Surface Winds near the Yarnell Hill Fire

On 30 June 2013, 19 Granite Mountain Hotshots firefighters were killed fighting a wildfire near Yarnell in the mountains of Central Arizona. They succumbed when the wildfire, driven by erratic winds, blocked their escape route and overran their location. A previous study is extended to simulate and analyze the downscale organization of convective circulations that redirected the wildfire, which started from the scale of the Rossby Wave Breaking over North America to a convective gust front that redirected the wildfire, trapping the firefighters. Five stages are found: Stage I, the initial deep prolonged gust front; Stage II, a front-to-rear jet and its ascending motions that organized high-based convection; Stage III, high-based dry microburst-induced downdrafts organized initially by ascending flow in Stage II that transported mass and entropy to the surface; Stage IV; multiple meso-γ-scale high centers and confluence zones formed that encompassed the firefighters’ location, which established a favorable environment leading to Stage V, canyon-scale circulations formed surrounding the fire. The atmosphere thus transitioned from supporting a deep and long-lived convective density current to elevated dry microbursts with mass and wind outflow into a canyon, redirecting the ongoing wildfire.

Transitions between convective reaction fronts in a Poiseuille flow

Transitions between convective reaction fronts in a Poiseuille flow

Explaining prolonged fluctuations in light curves of classical novae via modelling

ABSTRACT Fluctuations during a prolonged maximum have been observed in several nova eruptions, although it is not clear, and cannot be deduced directly from observations, whether the phenomenon is an actual physical reaction to some mechanism originating in the erupting white dwarf, whether it is occurring in the expanding ejected shell, or whether it is a form of interaction with the red dwarf companion. A handful of erupting nova models are investigated in this work, in order to assess the possibility of this sort of feature being an actual part of the eruption itself. The results explain that the mechanism that may produce these fluctuations is the repeated approach and recession of the convective front from the surface. The efficiency of this mechanism, being dependent on the mass of the WD (white dwarf) envelope and the time-scale of the nova cycle, favours low-mass WDs and long accretion phases.

The Thermal Imbalances Recorded at the NE Rift during the 2012 Explosive Activity at the South East Cone (Mt. Etna, Italy)

Mild thermal anomalies are sensitive to change in the advection processes in a volcanic system. A mild thermal anomaly, near the top of the North-East Rift of Mt. Etna (Italy), has been monitored from January 2010 to September 2012 by means of four temperature sensors buried in the shallow ground. The pulses of the convective circulation have been tracked and the diffuse heat flux has been evaluated. The positive pulses of the convective front reflected the local increases of volcanic degassing; conversely, the negative pulses showed the contraction of the convective front emerging through the North-East Rift. The steam condensation depth fluctuated below the monitoring site, from depths of a couple of meters to more than 30 meters, while the New South-East crater was erupting. The data hourly recorded, relative to the 2012 eruptive period, were compared to the radiant energy released by the paroxysms. We registered a dramatic decrease in the diffuse heat flux several hours before the onset of the two most energetic paroxysms (12 and 23 April). Thereafter, the convective front (the steam condensation depth) showed many negative pulses, reaching the deepest recorded levels. Thermal transients could be one of the early signals, possibly heralding transitions in the dynamic equilibrium conditions.

Absolute vs Convective Instabilities and Front Propagation in Lipid Membrane Tubes.

We analyze the stability of biological membrane tubes, with and without a base flow of lipids. Membrane dynamics are completely specified by two dimensionless numbers: the well-known Föppl-von Kármán number Γ and the recently introduced Scriven-Love number SL, respectively quantifying the base tension and base flow speed. For unstable tubes, the growth rate of a local perturbation depends only on Γ, whereas SL governs the absolute versus convective nature of the instability. Furthermore, nonlinear simulations of unstable tubes reveal an initially localized disturbance result in propagating fronts, which leave a thin atrophied tube in their wake. Depending on the value of Γ, the thin tube is connected to the unperturbed regions via oscillatory or monotonic shape transitions-reminiscent of recent experimental observations on the retraction and atrophy of axons. We elucidate our findings through a weakly nonlinear analysis, which shows membrane dynamics may be approximated by a model of the class of extended Fisher-Kolmogorov equations. Our study sheds light on the pattern selection mechanism in axonal shapes by recognizing the existence of two Lifshitz points, at which the front dynamics undergo steady-to-oscillatory bifurcations.

The Multi-Scale Dynamics Organizing a Favorable Environment for Convective Density Currents That Redirected the Yarnell Hill Fire

The deadly shift of the Yarnell Hill, Arizona wildfire was associated with an environment exhibiting gusty wind patterns in response to organized convectively driven circulations. The observed synoptic (>2500 km) through meso-β (approximately 100 km) scale precursor environment that organized a mid-upper tropospheric cross-mountain mesoscale jet streak circulation and upslope thermally direct flow was examined. Numerical simulations and observations indicated that both circulations played a key role in focusing the upper-level divergence, ascent, downdraft potential, vertical wind shear favoring mobile convective gust fronts, and a microburst. This sequence was initiated at the synoptic scale by a cyclonic Rossby Wave Break (RWB) 72 h prior, followed by an anticyclonic RWB. These RWBs combined to produce a mid-continent baroclinic trough with two short waves ushering in cooler air with the amplifying polar jet. Cool air advection with the second trough and surface heating across the Intermountain West (IW) combined to increase the mesoscale pressure gradient, forcing a mid-upper tropospheric subsynoptic jet around the periphery of the upstream ridge over Southern Utah and Northern New Mexico. Convection was triggered by an unbalanced secondary jetlet circulation within the subsynoptic jet in association with a low-level upslope flow accompanying a mountain plains solenoidal circulation above the Mogollon Rim (MR) and downstream mountains.

Capillarity and phase-mobility of a hydrocarbon gas–liquid system

When oil fields fall during their lifetime below the bubble point gas comes out of solution. The key questions are at which saturation the gas becomes mobile (“critical gas saturation”) and what the gas mobility is, because mobile gas reduces the production of oil significantly. The traditional view is that the gas phase becomes mobile once gas bubbles grow or expand to a size where they connect and form a percolating path. For typical 3D porous media the saturation corresponding to this percolation limit is on the order of 20%. However, significant literature report on gas mobility below lower limits of percolation thresholdsi.e.below 0.1%. A direct experimental insight for that is lacking because laboratory measurements are notoriously difficult since the formation of gas bubbles below the bubble point includes thermodynamic and kinetic aspects, and the pressure decline rates achievable in laboratory experiments are orders of magnitude higher than the decline rates in the field. Here we study the nucleation and transport of gas coming out of solutionin-situin 3D rock using X-ray computed micro tomography which allows direct visualization of the nucleation kinetics and connectivity of gas. We use either propane or a propane–decane mixture as model system and conduct pressure depletion in absence of flow finding that – consistent with the literature – observation of the bubble point in the porous medium is decreased and becomes pressure decline rate dependent because of the bubble nucleation kinetics. That occurs in single-component systems and in hydrocarbon mixtures. Pressure depletion in absence of flow results in critical gas saturations between 20 and 30% which is consistent with typical percolation thresholds in 3D porous structures. That does not explain experimentally observed critical gas saturations significantly below 20%. Also, the respective pore level fluid occupancy where pores are filled with either gas or liquid phase but not partially with both as in normal 2-phase immiscible systems rather diminishes connectivity of gas and liquid phases. This observation indicates that likely other mechanisms play a role in establishing gas mobility at saturations significantly below 20%. Experiments under flow conditions, where gas is injected near the bubble point suggest that diffusion may significantly contribute to the transport of gas and may even be the dominant transport mechanism at field relevant flow rates. The consequence of diffusive transport are compositional gradients where locally the composition is such gas nucleation may occur. That would lead to a disconnected but mobile gas distribution ahead of the convective front. Furthermore, diffusive exchange leads to ripening and anti-ripening effects which influences the distribution for which we see evidence in pressure depletion experiments but not so much at low rate gas injection. Respective relative permeability computed from the imaged fluid distributions using a lattice Boltzmann approach show distinctly different behavior between pressure depletion and flowing conditions. These findings suggest that capillarity in a gas–liquid hydrocarbon mixture is far more complex than in a 2-phase immiscible system. Capillarity is coupled to phase behavior thermodynamics and kinetics on a fast time scale and diffusion-dominated mechanisms such as ripening and anti-ripening effects at a slow time scale. While the consequences for the current experimental and field modelling approaches are not yet fully clear, this shows that more research is needed to fully understand these effects and their implications.

Improving well productivity – Ways to reduce the lag between the diffusive and convective time of flight in shale wells

This short communication explains how the key to improve shale well productivity lies in reducing the gap between the diffusive and convective time of flight. Diffusive time of flight in a hydrocarbon reservoir controls which reservoir regions will be affected by the pressure transient at a certain time after the onset of production. The convective time of flight is a measure of tracer front advance as well as for fluid withdrawal rates due to a well system draining the reservoir. The difference between the advance rate of a diffusive pressure front and a convective tracer front has been long recognized. However, only recently has it become clear that a thorough grasp of the two concepts is crucial for better decisions about fracture and well spacing in low permeability reservoirs. Although the diffusive time of flight is highlighted in the traditional pressure depletion models, the convective time of flight has a more direct impact on well productivity because it controls the advance rate, to the well system, of hydrocarbon fluids to be produced. The two concepts, diffusive and convective time of flight, are briefly highlighted below, with emphasis on their relevance for decision-making in shale field development solutions.

Experimental and numerical study on melting process of paraffin in a vertical annular cylinder

The melting process of paraffin wax in a central heating vertical annular cylinder was simulated by a numerical model based on the enthalpy-porosity method. The numerical results were validated against the experimental results. The melting experiment was performed in a thermostatic water tank which can provide good initial conditions and constant temperature boundary conditions. The electrical heating rod located in the central of the cylinder was supposed to provide heating at a constant power. Good agreements between experimental and numerical heat transfer data were achieved. Natural convection and melt front interface were well predicted by simulation. However, considerable differences appeared due to mushy zone constant in the model. A recommended value of mushy zone constant was obtained by comparison to the experiment. The influence of thermal conductivity coefficient was also analyzed numerically.

Internal sub-sonic burning during an explosion viewed via dynamic X-ray radiography

We observe internal convective and conductive burn front propagation and solid consumption subsequent to thermal ignition for plastic bonded formulations of the solid organic secondary explosives octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine and 1,3,5-triamino-2,4,6-trinitrobenzene. This work describes x-ray radiographic diagnostics enabling the study of solid density in a fully encased explosive during internal burning subsequent to ignition. The result of this study is the ability to directly observe and measure rates of energy release during a thermal explosion.

Optical Hilbert diagnostics of convective structures and phase transition in a horizontal layer of supercooled water

Methods of Hilbert optics are used to visualize evolution of convective structures and wave front of crystallization in a horizontal layer of supercooled water bounded by thermally stabilized plane surfaces. Measured phase velocities and shapes of the crystallization front approximated with the aid of the Bezier curves are presented. The spatiotemporal quasi-periodicity of the phase-velocity field can be related to oscillatory modes of the crystallization front as a dynamic surface of the water–ice interphase transition.

Fabrication of PMDA-ODA hollow fibers with regular cross-section morphologies and study on the formation mechanism

Since PMDA-ODA polyimide is insoluble, its prepolymer (PMDA-ODA polyamic acid) has to be used to prepare asymmetric membranes which are then converted to polyimide using the imidization method. However, the phase inversion process of polyamic acid (PAA) in water is slow, which causes difficulties in the process of spinning PMDA-ODA hollow fiber (HF). In this paper, a variety of spinning parameters, including the PAA concentration, the content of ethanol additive, the bore fluid composition, the dope and bore flow rates, the HF take-up speed, the air gap distance, the temperature and size of coagulation bath have been adjusted to investigate their impacts on the HF cross-section morphologies. Based on our observations, the guidelines of spinning PMDA-ODA HF with regularly round inner and outer surfaces are given. That is, first, increase the phase inversion rate of the PAA dope; second, increase the residence period of the HF in the water coagulant; and third, maintain a relatively low coagulant temperature. Moreover, the viscosities, the grown rates of the precipitation and convection fronts of different PAA dopes, and the ternary phase diagram have been measured. The formation mechanism of the cross-section morphologies of the PAA HFs have been well explained using the abovementioned results.

Tropospheric Turbulence over the Tropical Open Ocean: Role of Gravity Waves

Abstract A large set of soundings obtained in the Indian Ocean during three field campaigns is used to provide statistical characteristics of tropospheric turbulence and its link with gravity wave (GW) activity. The Thorpe method is used to diagnose turbulent regions of a few hundred meters depth. Above the mixed layer, turbulence frequency varies from ~10% in the lower troposphere up to ~30% around 12-km height. GWs are captured by their signature in horizontal wind, normalized temperature, and balloon vertical ascent rate. These parameters emphasize different parts of the wave spectrum from longer to shorter vertical wavelengths. Composites are constructed in order to reveal the vertical structure of the waves and their link with turbulence. The relatively longer-wavelength GWs described by their signature in temperature (GWTs) are more active in the lower troposphere, where they are associated with clear variations in moisture. Turbulence is then associated with minimum static stability and vertical shear, stressing the importance of the former and the possibility of convective instability. Conversely, the short waves described by their signature in balloon ascent rate (GWws) are detected primarily in the upper troposphere, and their turbulence is associated with a vertical shear maximum, suggesting the importance of dynamic instability. Furthermore, GWws appear to be linked with local convection, whereas GWTs are more active in suppressed and dry phases in particular of the Madden–Julian oscillation. These waves may be associated with remote sources, such as organized convection or local fronts, such as those associated with dry-air intrusions.

On the Dynamical Control of the Mesosphere–Lower Thermosphere by the Lower and Middle Atmosphere

Abstract The NCAR Whole Atmosphere Community Climate Model (WACCM) is used to investigate the dynamical influence of the lower and middle atmosphere on the upper mesosphere and lower thermosphere. In simulations using a methodology adapted from the “specified dynamics” (nudged) version of the model, horizontal winds and temperature over part of the vertical range of the atmosphere are relaxed toward results from a previous simulation that serves as the true simulation, equivalent to meteorological analysis. In the upper mesosphere, the magnitude of the divergence of the constrained simulations from the true simulation depends on the vertical extent and frequency of the data used for nudging the model and grows with altitude. The simulations quantify the error growth of the model dynamical fields when data and forcing terms are known exactly and there are no model biases. The error growth rate and the ultimate discrepancy between the nudged and true fields depend strongly on the method used for representing gravity wave drag. The largest error growth occurs when the gravity wave parameterization uses interactive wave sources that depend on convective activity or fronts. Errors are reduced when the same parameterization is used with smoothly varying specified wave sources. The smallest errors are seen when the parameterized gravity wave drag is replaced by linear Rayleigh friction damping on the wind speed. These comparisons demonstrate the role of gravity waves in transporting the variability of the troposphere into the mesosphere and lower thermosphere.

Downdraft outflows: climatological potential to influence fire behaviour

Sudden wind shifts caused by atmospheric gust fronts can lead to firefighter entrapments and fatalities. In this study, we describe the physical processes involved in the related phenomena of convective downdrafts, gust fronts and downbursts. We focus on the dominant process, evaporative cooling in a dry surface layer, as characterised by the measure known as downdraft maximum available potential energy (DMAPE). We present a climatological analysis of DMAPE for the coterminous United States, developed from the Climate Forecast System Reanalysis data for the period 1979–2008. Diurnal and seasonal patterns are described. We conclude with a discussion of the implications and limitations of the analysis and DMAPE as an indicator of gust front strength or downburst occurrence.

Influence of Marangoni flows on the dynamics of isothermal A + B → C reaction fronts.

The nonlinear dynamics of A + B → C fronts is analyzed both numerically and theoretically in the presence of Marangoni flows, i.e., convective motions driven by surface tension gradients. We consider horizontal aqueous solutions where the three species A, B, and C can affect the surface tension of the solution, thereby driving Marangoni flows. The resulting dynamics is studied by numerically integrating the incompressible Navier-Stokes equations coupled to reaction-diffusion-convection (RDC) equations for the three chemical species. We show that the dynamics of the front cannot be predicted solely on the basis of the one-dimensional reaction-diffusion profiles as is the case for buoyancy-driven convection around such fronts. We relate this observation to the structure of Marangoni flows which lead to more complex and exotic dynamics. We find in particular the surprising possibility of a reversal of the front propagation direction in time for some sets of Marangoni numbers, quantifying the influence of each chemical species concentration on the solution surface tension. We explain this reversal analytically and propose a new classification of the convective effects on A + B → C reaction fronts as a function of the Marangoni numbers. The influence of the layer thickness on the RDC dynamics is also presented. Those results emphasize the importance of flow symmetry properties when studying convective front dynamics in a given geometry.

An explicit approach to capture diffusive effects in finite water-content method for solving vadose zone flow

Summary Vadose zone flow problems are usually solved from the Richards equation. Solution to the Richards equation is generally challenging because the hydraulic conductivity and diffusivity in the equation are strongly non-linear functions of water content. The finite water-content method was proposed as an alternative general solution method of the vadose zone flow problem for infiltration, falling slugs, and vadose zone response to water table dynamics based on discretizing the water content domain into numerous bins instead of the traditional spatial discretization. In this study, we develop an improved approach to the original finite water-content method (referred to as TO method hereinafter) that better simulates diffusive effects but retains the robustness of the TO method. The approach treats advection and diffusion separately and considers diffusion on a bin by bin basis. After discretizing into water content bins, we treat the conductivity and diffusivity in individual bins as water content dependent constant evaluated at given water content corresponding to each bin. For each bin, we can solve the flow equations analytically since the hydraulic conductivity and diffusivity can be treated as a constant. We then develop solutions for each bin to determine the diffusive water amounts at each time step. The water amount ahead of the convective front for each bin is redistributed among water content bins to account for diffusive effects. The application of developed solution is straightforward only involving algebraic manipulations at each time step. The method can mainly improve water content profiles, but has no significant difference for the total infiltration rate and cumulative infiltration compared to the TO method. Although the method separately deals with advection and diffusion, it can account for the coupling effects of advection and diffusion reasonably well.

Properties of small cirrus ice crystals from commercial aircraft measurements and implications for flight operations

Measurements of cloud ice crystal size distributions have been made by a backscatter cloud probe (BCP) mounted on five commercial airliners flying international routes that cross five continents. Bulk cloud parameters were also derived from the size distributions. As of 31 December 2014, a total of 4399 flights had accumulated data from 665 hours in more than 19 000 cirrus clouds larger than 5 km in length. The BCP measures the equivalent optical diameter (EOD) of individual crystals in the 5–90 µm range from which size distributions are derived and recorded every 4 seconds. The cirrus cloud property database, an ongoing development stemming from these measurements, registers the total crystal number and mass concentration, effective and median volume diameters and extinction coefficients derived from the size distribution. This information is accompanied by the environmental temperature, pressure, aircraft position, date and time of each sample. The seasonal variations of the cirrus cloud properties measured from 2012 to 2014 are determined for six geographic regions in the tropics and extratropics. Number concentrations range from a few per litre for thin cirrus to several hundreds of thousands for heavy cirrus. Temperatures range from 205 to 250 K and effective radii from 12 to 20 µm. A comparison of the regional and seasonal number and mass size distributions, and the bulk microphysical properties derived from them, demonstrates that cirrus properties cannot be easily parameterised by temperature or by latitude. The seasonal changes in the size distributions from the extratropical Atlantic and Eurasian air routes are distinctly different, showing shifts from mono-modal to bi-modal spectra out of phase with one another. This phase difference may be linked to the timing of deep convection and cold fronts that lead to the cirrus formation. Likewise, the size spectra of cirrus over the tropical Atlantic and Eastern Brazil differ from each other although they were measured in adjoining regions. The cirrus crystals in the maritime continental tropical region over Malaysia form tri-modal spectra that are not found in any of the other regions measured by the IAGOS aircraft so far, a feature that is possibly linked to biomass burning or dust. Frequent measurements of ice crystal concentrations greater than 1×105 L−1, often accompanied by anomalously warm temperature and erratic airspeed readings, suggest that aircraft often experience conditions that affect their sensors. This new instrument, if used operationally, has the potential of providing real-time and valuable information to assist in flight operations as well as providing real-time information for along-track nowcasting.