Mesh adaptation on the sphere using optimal transport and the numerical solution of a Monge-Ampère type equation

An equation of Monge-Ampère type has, for the first time, been solved numerically on the surface of the sphere in order to generate optimally transported (OT) meshes, equidistributed with respect to a monitor function. Optimal transport generates meshes that keep the same connectivity as the original mesh, making them suitable for r-adaptive simulations, in which the equations of motion can be solved in a moving frame of reference in order to avoid mapping the solution between old and new meshes and to avoid load balancing problems on parallel computers. The semi-implicit solution of the Monge-Ampère type equation involves a new linearisation of the Hessian term, and exponential maps are used to map from old to new meshes on the sphere. The determinant of the Hessian is evaluated as the change in volume between old and new mesh cells, rather than using numerical approximations to the gradients. OT meshes are generated to compare with centroidal Voronoi tesselations on the sphere and are found to have advantages and disadvantages; OT equidistribution is more accurate, the number of iterations to convergence is independent of the mesh size, face skewness is reduced and the connectivity does not change. However anisotropy is higher and the OT meshes are non-orthogonal. It is shown that optimal transport on the sphere leads to meshes that do not tangle. However, tangling can be introduced by numerical errors in calculating the gradient of the mesh potential. Methods for alleviating this problem are explored. Finally, OT meshes are generated using observed precipitation as a monitor function, in order to demonstrate the potential power of the technique.

Conditioning and preconditioning of the optimal state estimation problem

Optimal state estimation is a method that requires minimising a weighted, nonlinear, least-squares objective function in order to obtain the best estimate of the current state of a dynamical system. Often the minimisation is non-trivial due to the large scale of the problem, the relative sparsity of the observations and the nonlinearity of the objective function. To simplify the problem the solution is often found via a sequence of linearised objective functions. The condition number of the Hessian of the linearised problem is an important indicator of the convergence rate of the minimisation and the expected accuracy of the solution. In the standard formulation the convergence is slow, indicating an ill-conditioned objective function. A transformation to different variables is often used to ameliorate the conditioning of the Hessian by changing, or preconditioning, the Hessian. There is only sparse information in the literature for describing the causes of ill-conditioning of the optimal state estimation problem and explaining the effect of preconditioning on the condition number. This paper derives descriptive theoretical bounds on the condition number of both the unpreconditioned and preconditioned system in order to better understand the conditioning of the problem. We use these bounds to explain why the standard objective function is often ill-conditioned and why a standard preconditioning reduces the condition number. We also use the bounds on the preconditioned Hessian to understand the main factors that affect the conditioning of the system. We illustrate the results with simple numerical experiments.

Transatlantic flight times and climate change

Aircraft do not fly through a vacuum, but through an atmosphere whose meteorological characteristics are changing because of global warming. The impacts of aviation on climate change have long been recognised, but the impacts of climate change on aviation have only recently begun to emerge. These impacts include intensified turbulence and increased take-off weight restrictions. Here we investigate the influence of climate change on flight routes and journey times. We feed synthetic atmospheric wind fields generated from climate model simulations into a routing algorithm of the type used operationally by flight planners. We focus on transatlantic flights between London and New York, and how they change when the atmospheric concentration of carbon dioxide is doubled. We find that a strengthening of the prevailing jet-stream winds causes eastbound flights to significantly shorten and westbound flights to significantly lengthen in all seasons. Eastbound and westbound crossings in winter become approximately twice as likely to take under 5 h 20 min and over 7 h 00 min, respectively. For reasons that are explained using a conceptual model, the eastbound shortening and westbound lengthening do not cancel out, causing round-trip journey times to increase. Even assuming no future growth in aviation, the extrapolation of our results to all transatlantic traffic suggests that aircraft will collectively be airborne for an extra 2000 h each year, burning an extra 7.2 million gallons of jet fuel at a cost of US$ 22 million, and emitting an extra 70 million kg of carbon dioxide, which is equivalent to the annual emissions of 7100 average British homes. Our results provide further evidence of the two-way interaction between aviation and climate change.

Ecuador Paraiso Escondido virus, a new flavivirus isolated from New World sand flies in Ecuador, is the first representative of a novel clade in the genus flavivirus

A new flavivirus, Ecuador Paraiso Escondido virus (EPEV), named after the village where it was discovered, was isolated from sand flies (Psathyromyia abonnenci, formerly Lutzomyia abonnenci) that are unique to the New World. This represents the first sand fly-borne flavivirus identified in the New World. EPEV exhibited a typical flavivirus genome organization. Nevertheless, the maximum pairwise amino acid sequence identity with currently recognized flaviviruses was 52.8%. Phylogenetic analysis of the complete coding sequence showed that EPEV represents a distinct clade which diverged from a lineage that was ancestral to the nonvectored flaviviruses Entebbe bat virus, Yokose virus, and Sokoluk virus and also the Aedes-associated mosquito-borne flaviviruses, which include yellow fever virus, Sepik virus, Saboya virus, and others. EPEV replicated in C6/36 mosquito cells, yielding high infectious titers, but failed to reproduce either in vertebrate cell lines (Vero, BHK, SW13, and XTC cells) or in suckling mouse brains. This surprising result, which appears to eliminate an association with vertebrate hosts in the life cycle of EPEV, is discussed in the context of the evolutionary origins of EPEV in the New World.The flaviviruses are rarely (if ever) vectored by sand fly species, at least in the Old World. We have identified the first representative of a sand fly-associated flavivirus, Ecuador Paraiso Escondido virus (EPEV), in the New World. EPEV constitutes a novel clade according to current knowledge of the flaviviruses. Phylogenetic analysis of the virus genome showed that EPEV roots the Aedes-associated mosquito-borne flaviviruses, including yellow fever virus. In light of this new discovery, the New World origin of EPEV is discussed together with that of the other flaviviruses.

Impact of the North Atlantic Oscillation on transatlantic flight routes and clear-air turbulence

The variation of wind-optimal transatlantic flight routes and their turbulence potential is investigated to understand how upper-level winds and large-scale flow patterns can affect the efficiency and safety of long-haul flights. In this study, the wind-optimal routes (WORs) that minimize the total flight time by considering wind variations are modeled for flights between John F. Kennedy International Airport (JFK) in New York, New York, and Heathrow Airport (LHR) in London, United Kingdom, during two distinct winter periods of abnormally high and low phases of North Atlantic Oscillation (NAO) teleconnection patterns. Eastbound WORs approximate the JFK–LHR great circle (GC) route following northerly shifted jets in the +NAO period. Those WORs deviate southward following southerly shifted jets during the −NAO period, because eastbound WORs fly closely to the prevailing westerly jets to maximize tailwinds. Westbound WORs, however, spread meridionally to avoid the jets near the GC in the +NAO period to minimize headwinds. In the −NAO period, westbound WORs are north of the GC because of the southerly shifted jets. Consequently, eastbound WORs are faster but have higher probabilities of encountering clear-air turbulence than westbound ones, because eastbound WORs are close to the jet streams, especially near the cyclonic shear side of the jets in the northern (southern) part of the GC in the +NAO (−NAO) period. This study suggests how predicted teleconnection weather patterns can be used for long-haul strategic flight planning, ultimately contributing to minimizing aviation’s impact on the environment