The streaming instability for solid particles in protoplanetary disks is re-examined assuming the familiar alpha ($\alpha$) model for isotropic turbulence. Turbulence always reduces the growth rates of the streaming instability relative to values calculated for globally laminar disks. While for small values of the turbulence parameter, $\alpha < 10^{-5}$, the wavelengths of the fastest-growing disturbances are small fractions of the local gas vertical scale height $H$, we find that for moderate values of the turbulence parameter, i.e., $\alpha \sim 10^{-5}-10^{-3}$, the lengthscales of maximally growing disturbances shift toward larger scales, approaching $H$. At these moderate turbulent intensities and for local particle to gas mass density ratios $\epsilon < 0.5$, the vertical scales of the most unstable modes begin to exceed the corresponding radial scales so that the instability appears in the form of vertically oriented sheets extending well beyond the particle scale height. We find that for hydrodynamical turbulent disk models reported in the literature, with $\alpha = 4\times 10^{-5} - 5\times 10^{-4}$, together with state of the art global evolution models of particle growth, the streaming instability is predicted to be viable within a narrow triangular patch of $\alpha$--$\tau_s$ parameter space centered on Stokes numbers, $\tau_s \sim 0.01$ and $\alpha \sim 4\times 10^{-5}$ and, further, exhibits growth rates on the order of several hundred to thousands of orbit times for disks with 1 percent ($Z= 0.01$) cosmic solids abundance or metallicity. Our results are consistent with, and place in context, published numerical studies of streaming instabilities.

The Kuiper Belt is a distant region of the Solar System. On 1 January 2019, the New Horizons spacecraft flew close to (486958) 2014 MU69, a Cold Classical Kuiper Belt Object, a class of objects that have never been heated by the Sun and are therefore well preserved since their formation. Here we describe initial results from these encounter observations. MU69 is a bi-lobed contact binary with a flattened shape, discrete geological units, and noticeable albedo heterogeneity. However, there is little surface color and compositional heterogeneity. No evidence for satellites, ring or dust structures, gas coma, or solar wind interactions was detected. By origin MU69 appears consistent with pebble cloud collapse followed by a low velocity merger of its two lobes.

The Cold Classical Kuiper Belt, a class of small bodies in undisturbed orbits beyond Neptune, are primitive objects preserving information about Solar System formation. The New Horizons spacecraft flew past one of these objects, the 36 km long contact binary (486958) Arrokoth (2014 MU69), in January 2019. Images from the flyby show that Arrokoth has no detectable rings, and no satellites (larger than 180 meters diameter) within a radius of 8000 km, and has a lightly-cratered smooth surface with complex geological features, unlike those on previously visited Solar System bodies. The density of impact craters indicates the surface dates from the formation of the Solar System. The two lobes of the contact binary have closely aligned poles and equators, constraining their accretion mechanism.

The New Horizons spacecraft's encounter with the cold classical Kuiper belt object (486958) Arrokoth (formerly 2014 MU69) revealed a contact-binary planetesimal. We investigate how it formed, finding it is the product of a gentle, low-speed merger in the early Solar System. Its two lenticular lobes suggest low-velocity accumulation of numerous smaller planetesimals within a gravitationally collapsing, solid particle cloud. The geometric alignment of the lobes indicates the lobes were a co-orbiting binary that experienced angular momentum loss and subsequent merger, possibly due to dynamical friction and collisions within the cloud or later gas drag. Arrokoth's contact-binary shape was preserved by the benign dynamical and collisional environment of the cold classical Kuiper belt, and so informs the accretion processes that operated in the early Solar System.

The outer Solar System object (486958) Arrokoth (provisional designation 2014 MU$_{69}$) has been largely undisturbed since its formation. We study its surface composition using data collected by the New Horizons spacecraft. Methanol ice is present along with organic material, which may have formed through radiation of simple molecules. Water ice was not detected. This composition indicates hydrogenation of carbon monoxide-rich ice and/ or energetic processing of methane condensed on water ice grains in the cold, outer edge of the early Solar System. There are only small regional variations in color and spectra across the surface, suggesting Arrokoth formed from a homogeneous or well-mixed reservoir of solids. Microwave thermal emission from the winter night side is consistent with a mean brightness temperature of 29$\pm$5 K.

Rosetta observations of 67P/Churyumov-Gerasimenko (67P) reveal that most changes occur in the fallback-generated smooth terrains, vast deposits of granular material blanketing the comet's northern hemisphere. These changes express themselves both morphologically and spectrally across the nucleus, yet we lack a model that describes their formation and evolution. Here we present a self-consistent model that thoroughly explains the activity and mass loss from Hapi's smooth terrains. Our model predicts the removal of dust via re-radiated solar insolation localized within depression scarps that are substantially more ice-rich than previously expected. We couple our model with numerous Rosetta observations to thoroughly capture the seasonal erosion of Hapi's smooth terrains, where local scarp retreat gradually removes the uppermost dusty mantle. As sublimation-regolith interactions occur on rocky planets, comets, icy moons and KBOs, our coupled model and observations provide a foundation for future understanding of the myriad of sublimation-carved worlds.

The deep nitrogen-covered Sputnik Planitia (SP; informal name) basin on Pluto is located very close to the longitude of Pluto's tidal axis[1] and may be an impact feature [2], by analogy with other large basins in the solar system[3,4]. Reorientation[5-7] due to tidal and rotational torques can explain SP's location, but requires it to be a positive gravity anomaly[7], despite its negative topography. Here we argue that if SP formed via impact and if Pluto possesses a subsurface ocean, a positive gravity anomaly would naturally result because of shell thinning and ocean uplift, followed by later modest N2 deposition. Without a subsurface ocean a positive gravity anomaly requires an implausibly thick N2 layer (greater than 40 km). A rigid, conductive ice shell is required to prolong such an ocean's lifetime to the present day[8] and maintain ocean uplift. Because N2 deposition is latitude-dependent[9], nitrogen loading and reorientation may have exhibited complex feedbacks[7].

The vast, deep, volatile-ice-filled basin informally named Sputnik Planum is central to Pluto's geological activity[1,2]. Composed of molecular nitrogen, methane, and carbon monoxide ices[3], but dominated by N2-ice, this ice layer is organized into cells or polygons, typically ~10-40 km across, that resemble the surface manifestation of solid state convection[1,2]. Here we report, based on available rheological measurements[4], that solid layers of N2 ice approximately greater than 1 km thick should convect for estimated present-day heat flow conditions on Pluto. More importantly, we show numerically that convective overturn in a several-km-thick layer of solid nitrogen can explain the great lateral width of the cells. The temperature dependence of N2-ice viscosity implies that the SP ice layer convects in the so-called sluggish lid regime[5], a unique convective mode heretofore not definitively observed in the Solar System. Average surface horizontal velocities of a few cm/yr imply surface transport or renewal times of ~500,000 years, well under the 10 Myr upper limit crater retention age for Sputnik Planum[2]. Similar convective surface renewal may also occur on other dwarf planets in the Kuiper belt, which may help explain the high albedos of some of them.

Modern astronomy has finally been able to observe protoplanetary disks in reasonable resolution and detail, unveiling the processes happening during planet formation. These observed processes are understood under the framework of disk-planet interaction, a process studied analytically and modeled numerically for over 40 years. Long a theoreticians' game, the wealth of observational data has been allowing for increasingly stringent tests of the theoretical models. Modeling efforts are crucial to support the interpretation of direct imaging analyses, not just for potential detections but also to put meaningful upper limits on mass accretion rates and other physical quantities in current and future large-scale surveys. This white paper addresses the questions of what efforts on the computational side are required in the next decade to advance our theoretical understanding, explain the observational data, and guide new observations. We identified the nature of accretion, ab initio planet formation, early evolution, and circumplanetary disks as major fields of interest in computational planet formation. We recommend that modelers relax the approximations of alpha-viscosity and isothermal equations of state, on the grounds that these models use flawed assumptions, even if they give good visual qualitative agreement with observations. We similarly recommend that population synthesis move away from 1D hydrodynamics. The computational resources to reach these goals should be developed during the next decade, through improvements in algorithms and the hardware for hybrid CPU/GPU clusters. Coupled with high angular resolution and great line sensitivity in ground based interferometers, ELTs and JWST, these advances in computational efforts should allow for large strides in the field in the next decade.

A unique feature of Pluto's large satellite Charon is its dark red northern polar cap. Similar colours on Pluto's surface have been attributed to organic macromolecules produced by energetic radiation processing of hydrocarbons. The polar location of this material on Charon implicates the temperature extremes that result from Charon's high obliquity and long seasons. The escape of Pluto's atmosphere provides a potential feed stock for production of complex chemistry. Gas from Pluto that is transiently cold-trapped and processed at Charon's winter pole was proposed as an explanation on the basis of an image of Charon's northern hemisphere, but not modelled quantitatively. Here we report images of the southern hemisphere illuminated by Pluto-shine and also images taken during the approach phase showing the northern polar cap over a range of longitudes. We model the surface thermal environment on Charon, the supply and temporary cold-trapping of material escaping from Pluto, and, while cold-trapped, its photolytic processing into more complex and less volatile molecules. The model results are consistent with the proposed mechanism producing the observed colour pattern on Charon.

New Horizons observations suggest that CH4 on Pluto has a complex history, involving reservoirs of different composition, thickness and stability controlled by volatile processes occurring on different timescales. In order to interpret these observations, we use a Pluto volatile transport model able to simulate the cycles of N2 and CH4 ices over millions of years. By assuming fixed solid mixing ratios, we explore how changes in surface albedos, emissivities and thermal inertias impact volatile transport. This work is therefore a direct and natural continuation of the work by Bertrand et al. (2018), which only explored the N2 cycles. Results show that bright CH4 deposits can create cold traps for N2 ice outside Sputnik Planitia, leading to a strong coupling between the N2 and CH4 cycles. Depending on the assumed albedo for CH4 ice, the model predicts CH4 ice accumulation (1) at the same equatorial latitudes where the Bladed Terrain Deposits are observed, supporting the idea that these CH4-rich deposits are massive and perennial, or (2) at mid-latitudes (25{\deg}N-70{\deg}N), forming a thick mantle which is consistent with New Horizons observations. In our simulations, both CH4 ice reservoirs are not in an equilibrium state and either one can dominate the other over long timescales, depending on the assumptions made for the CH4 albedo. This suggests that long-term volatile transport exists between the observed reservoirs. The model also reproduces the formation of N2 deposits at mid-latitudes and in the equatorial depressions surrounding the Bladed Terrain, as observed by New Horizons. At the poles, only seasonal CH4 and N2 deposits are obtained in Pluto's current orbital configuration. Finally, we show that Pluto's atmosphere always contained, over the last astronomical cycles, enough gaseous CH4 to absorb most of the incoming Lyman-flux.

This review examines recent theoretical developments in our understanding of turbulence in cold, non-magnetically active, planetesimal forming regions of protoplanetary disks which we refer to throughout as "Ohmic zones". We give a brief background introduction to the subject of disk turbulence followed by a terse pedagogical review of the phenomenology of hydrodynamic turbulence. The equations governing the dynamics of cold astrophysical disks are given and basic flow states are described. We discuss the Solberg-Hoiland conditions required for stability, and the three recently identified turbulence generating mechanisms possibly active in protoplanetary disk Ohmic zones, namely, (i) the Vertical Shear Instability, (ii) The Convective Overstability and (iii) the Zombie Vortex Instability. We summarize the properties of these processes, identify their limitations and discuss where and under what conditions these processes are active in protoplanetary disk Ohmic zones.

The New Horizons encounter with the cold classical Kuiper Belt object (KBO) 2014 MU69 (informally named 'Ultima Thule,' hereafter Ultima) on 1 January 2019 will be the first time a spacecraft has ever closely observed one of the free-orbiting small denizens of the Kuiper Belt. Related to but not thought to have formed in the same region of the Solar System as the comets that been explored so far, it will also be the largest, most distant, and most primitive body yet visited by spacecraft. In this letter we begin with a brief overview of cold classical KBOs, of which Ultima is a prime example. We give a short preview of our encounter plans. We note what is currently known about Ultima from earth-based observations. We then review our expectations and capabilities to evaluate Ultima's composition, surface geology, structure, near space environment, small moons, rings, and the search for activity.

Observations of Pluto's surface made by the New Horizons spacecraft indicates present-day nitrogen ice glaciation in and around the basin known as Sputnik Planum. Motivated by these observations, we have developed an evolutionary glacial flow model of solid nitrogen ice taking into account its published thermophysical and rheologies properties. This model assumes that glacial ice layers flow laminarly and have low aspect ratios which permits a vertically integrated mathematical formulation. We assess the conditions for the validity of laminar nitrogen ice motion by revisiting the problem of the onset of solid-state buoyant convection of nitrogen ice for a variety of bottom thermal boundary conditions. Subject to uncertainties in nitrogen ice rheology, nitrogen ice layers are estimated to flow laminarly for thicknesses less than 400-1000 meters. The resulting mass-flux formulation for when the nitrogen ice flows as a laminar dry glacier is characterized by an Arrhenius-Glen functional form. The flow model developed is used here to qualitatively answer some questions motivated by observed glacial flow features found on Sputnik Planum. We find that the wavy transverse dark features found along the northern shoreline of Sputnik Planum may be a transitory imprint of shallow topography just beneath the ice surface suggesting the possibility that a major shoreward flow event happened relatively recently within the last few hundred years. Model results also support the interpretation that the prominent darkened features resembling flow lobes observed along the eastern shoreline of the Sputnik Planum basin may be a result of wet nitrogen glacial ice flowing into the basin from the pitted highlands of eastern Tombaugh Regio.

Pluto's landscape is shaped by the cycles of the volatile ices covering its surface. In particular, the Sputnik Planitia (SP) ice sheet displays a large diversity of terrains, with bright and dark plains, pits, topographic depressions and evidences of recent and past glacial flows. Outside SP, New Horizons also revealed numerous N2 ice deposits, in Tombaugh Regio and at mid-northern latitudes. These observations suggest a complex history involving volatile and glacial processes on different timescales. We present numerical simulations of volatile transport on Pluto performed with a model able to simulate the N2 cycle over millions of years (Myrs), taking into account the changes of obliquity and orbital parameters as experienced by Pluto. Results show that over one obliquity cycle, the latitudes of SP between 25{\deg}S-30{\deg}N are dominated by N2 condensation, while the latitudes between 30-50{\deg}N are dominated by N2 sublimation. We find that a net amount of 1 km of ice has sublimed at the northern edge of SP during the last 2 Myrs. By comparing these results with the observed geology of SP, we can relate the formation of the pits and the brightness of the ice to the ice flux occurring at the annual timescale, while the glacial flows at its eastern edge and the erosion of the water ice mountains all around the ice sheet are related to the astronomical timescale. We also perform simulations with a glacial flow scheme which shows that SP is currently at its minimum extent. We also explore the stability of N2 ice outside SP. Results show that it is not stable at the poles but rather in the equatorial regions, in particular in depressions, where thick deposits may persist over tens of Myrs, before being trapped in SP. Finally, another key result is that the minimum and maximum surface pressures obtained over the simulated Myrs remain in the range of mm-Pa and Pa, respectively.

A linear analysis of the zombie vortex instability is performed in a stratified shearing sheet setting for three model barotropic shear flows. The linear analysis is done by utilizing a Green's function formulation to resolve the critical layers of the associated normal-mode problem. The instability is the result of a resonant interaction between a Rossby wave and a gravity wave which we refer to as Z-modes. The associated critical layer is the location where the Doppler shifted frequency of a distant Rossby wave equals the local Brunt-Vaisala frequency. The minimum required Rossby number for instability, Ro= 0.2, is confirmed for parameter values reported in the literature. It is also found that the shear layer supports the instability in the limit where stratification vanishes. The zombie vortex instability is examined in a jet model, finding that the instability can occur for Ro= 0.05. Nonlinear vorticity forcing due to unstable Z-modes is shown to result in the creation of a jet flow at the critical layer emerging as the result of the competition between the vertical lifting of perturbation radial vorticity and the radial transport of perturbation vertical vorticity. We find that the picture of this instability leading to a form of nonlinearly driven self-replicating pattern of creation and destruction is warranted: a parent jet spawns a growing child jet at associated critical layers. A mature child jet creates a next generation of child jets at associated critical layers while simultaneously destroying its parent jet via the Rossby wave instability.

Approach images taken by the LORRI imaging system during the New Horizons spacecraft encounter have been used to determine the mean radii and shapes of Pluto and Charon. The primary observations are limb locations derived using three independent approaches. The resulting mean radii of Pluto and Charon are 1188.3 +/- 1.6 km and 606.0 +/- 1.0 km, respectively (2-sigma). The corresponding densities are 1854 +/- 11 kg m-3 and 1701 +/- 33 kg m-3 (2-sigma). The Charon radius value is consistent with previous Earth-based occultation estimates. The Pluto radius estimate is consistent with solar occultation measurements performed by the ALICE and Fine Sun Sensor instruments on New Horizons. Neither Pluto nor Charon show any evidence for tidal/rotational distortions; upper bounds on the oblateness are <0.6% and <0.5%, respectively.

NASA's New Horizons spacecraft has revealed the complex geology of Pluto and Charon. Pluto's encounter hemisphere shows ongoing surface geological activity centered on a vast basin containing a thick layer of volatile ices that appears to be involved in convection and advection, with a crater retention age no greater than $\approx$10 Ma. Surrounding terrains show active glacial flow, apparent transport and rotation of large buoyant water-ice crustal blocks, and pitting, likely by sublimation erosion and/or collapse. More enigmatic features include tall mounds with central depressions that are conceivably cryovolcanic, and ridges with complex bladed textures. Pluto also has ancient cratered terrains up to ~4 Ga old that are extensionally fractured and extensively mantled and perhaps eroded by glacial or other processes. Charon does not appear to be currently active, but experienced major extensional tectonism and resurfacing (probably cryovolcanic) nearly 4 billion years ago. Impact crater populations on Pluto and Charon are not consistent with the steepest proposed impactor size-frequency distributions proposed for the Kuiper belt.

The dynamical response of edge waves under the influence of self-gravity is examined in an idealized two-dimensional model of a proto-stellar disc, characterized in steady state as a rotating vertically infinite cylinder of fluid with constant density except for a single density interface at some radius r0. The fluid in basic state is prescribed to rotate with a Keplerian profile $\Omega_k(r)\sim r^{-3/2}$ modified by some additional azimuthal sheared flow. A linear analysis shows that there are two azimuthally propagating edge waves, kin to the familiar Rossby waves and surface gravity waves in terrestrial studies, which move opposite to one another with respect to the local basic state rotation rate at the interface. Instability only occurs if the radial pressure gradient is opposite to that of the density jump (unstably stratified) where self-gravity acts as a wave stabilizer irrespective of the stratification of the system. The propagation properties of the waves are discussed in detail in the language of vorticity edge waves. The roles of both Boussinesq and non-Boussinesq effects upon the stability and propagation of these waves with and without the inclusion of self-gravity are then quantified. The dynamics involved with self-gravity non- Boussinesq effect is shown to be a source of vorticity production where there is a jump in the basic state density, in addition, self-gravity also alters the dynamics via the radial main pressure gradient, which is a Boussinesq effect . Further applications of these mechanical insights are presented in the conclusion including the ways in which multiple density jumps or gaps may or may not be stable.

We present the results of 2D and 3D hydrodynamic simulations of idealized protoplanetary discs that examine the formation and evolution of vortices by the vertical shear instability (VSI). In agreement with recent work, we find that discs with radially decreasing temperature profiles and short thermal relaxation time-scales, are subject to the axisymmetric VSI. In three dimensions, the resulting velocity perturbations give rise to quasi-axisymmetric potential vorticity perturbations that break-up into discrete vortices, in a manner that is reminiscent of the Rossby wave instability. Discs with very short thermal evolution time-scales (i.e. {\tau}<0.1 local orbit periods) develop strong vorticity perturbations that roll up into vortices that have small aspect ratios ({\chi}<2) and short lifetimes (~ a few orbits). Longer thermal time-scales give rise to vortices with larger aspect ratios (6<{\chi}<10), and lifetimes that depend on the entropy gradient. A steeply decreasing entropy profile leads to vortex lifetimes that exceed the simulation run times of hundreds of orbital periods. Vortex lifetimes in discs with positive or weakly decreasing entropy profiles are much shorter, being 10s of orbits at most, suggesting that the subcritical baroclinic instability plays an important role in sustaining vortices against destruction through the elliptical instability. Applied to the outer regions of protoplanetary discs, where the VSI is most likely to occur, our results suggest that vortices formed by the VSI are likely to be short lived structures.

The Pluto system was recently explored by NASA's New Horizons spacecraft, making closest approach on 14 July 2015. Pluto's surface displays diverse landforms, terrain ages, albedos, colors, and composition gradients. Evidence is found for a water-ice crust, geologically young surface units, surface ice convection, wind streaks, volatile transport, and glacial flow. Pluto's atmosphere is highly extended, with trace hydrocarbons, a global haze layer, and a surface pressure near 10 microbars. Pluto's diverse surface geology and long-term activity raise fundamental questions about how small planets remain active many billions of years after formation. Pluto's large moon Charon displays tectonics and evidence for a heterogeneous crustal composition, its north pole displays puzzling dark terrain. Small satellites Hydra and Nix have higher albedos than expected.

The Vertical Shear Instability is one of two known mechanisms potentially active in the so-called dead zones of protoplanetary accretion disks. A recent analysis indicates that a subset of unstable modes shows unbounded growth - both as resolution is increased and when the nominal lid of the atmosphere is extended, possibly indicating ill-posedness in previous attempts of linear analysis. The reduced equations governing the instability are revisited and the generated solutions are examined using both the previously assumed separable forms and an improved non-separable solution form that is herewith introduced. Analyzing the reduced equations using the separable form shows that, while the low-order body modes have converged eigenvalues and eigenfunctions (as both the vertical boundaries of the atmosphere are extended and with increased radial resolution), it is also confirmed that the corresponding high-order body modes and the surface modes do indeed show unbounded growth rates. However, the energy contained in both the higher-order body modes and surface modes diminishes precipitously due to the disk's Gaussian density profile. Most of the energy of the instability is contained in the low-order modes. An inseparable solution form is introduced which filters out the inconsequential surface modes leaving only body modes (both low and high-order ones). The analysis predicts a fastest growing mode with a specific radial length scale. The growth rates associated with the fundamental corrugation and breathing modes matches the growth and length scales observed in previous nonlinear studies of the instability.

We compare five general circulation models (GCMs) which have been recently used to study hot extrasolar planet atmospheres (BOB, CAM, IGCM, MITgcm, and PEQMOD), under three test cases useful for assessing model convergence and accuracy. Such a broad, detailed intercomparison has not been performed thus far for extrasolar planets study. The models considered all solve the traditional primitive equations, but employ different numerical algorithms or grids (e.g., pseudospectral and finite volume, with the latter separately in longitude-latitude and `cubed-sphere' grids). The test cases are chosen to cleanly address specific aspects of the behaviors typically reported in hot extrasolar planet simulations: 1) steady-state, 2) nonlinearly evolving baroclinic wave, and 3) response to fast timescale thermal relaxation. When initialized with a steady jet, all models maintain the steadiness, as they should -- except MITgcm in cubed-sphere grid. A very good agreement is obtained for a baroclinic wave evolving from an initial instability in pseudospectral models (only). However, exact numerical convergence is still not achieved across the pseudospectral models: amplitudes and phases are observably different. When subject to a typical `hot-Jupiter'-like forcing, all five models show quantitatively different behavior -- although qualitatively similar, time-variable, quadrupole-dominated flows are produced. Hence, as have been advocated in several past studies, specific quantitative predictions (such as the location of large vortices and hot regions) by GCMs should be viewed with caution. Overall, in the tests considered here, pseudospectral models in pressure coordinate (PEBOB and PEQMOD) perform the best and MITgcm in cubed-sphere grid performs the worst.

The dynamics contained in magnetized layers of exoplanet atmospheres are important to understand in order to characterize what observational signatures they may provide for future observations. It is important to develop a framework to begin studying and learning the physical processes possible under those conditions and what, if any, features contained in them may be observed in future observation missions. The aims of this study is to formally derive, from scaling arguments, a manageable reduced set of equations for analysis, i.e. a magnetic formulation of the equations of quasigeostrophy appropriate for a multi-layer atmosphere. The main goal is to provide a simpler theoretical platform to explore the dynamics possible within confined magnetized layers of exoplanet atmospheres. We primarily use scaling arguments to derive the reduced equations of "magnetoquasigeostrophy" which assumes dynamics to take place in an atmospheric layer which is vertically thin compared to its horizontal scales. The derived equation set retains features existing in standard shallow-water magnetohydrodynamic equations but are absent in more classical derivations of the quasi-geostrophic limit, namely, the non-divergence of the in-plane components of the magnetic field. We liken this non-divergence of the in-plane magnetic fields as indicative of a quantity whose behaviour mimics a two-dimensional "pseudo"-magnetic monopole source. We also find, using the same scaling argument procedures, appropriate limits of the fundamental parameters of the system which yield reduced equations describing the flow dynamics primarily characterized by magnetostrophic balance. The standard scaling arguments employed here show how traditional magnetized quasigeostrophic equations connect to their magnetized shallow water forms. The equations derived are amenable to analysis using well-known techniques.

A new non-dissipative mechanism is proposed for the saturation of the axisymmetric magnetorotational (MRI) instability in thin Keplerian disks that are subject to an axial magnetic field. That mechanism relies on the energy transfer from the MRI to stable magnetosonic (MS) waves. Such mode interaction is enabled due to the vertical stratification of the disk that results in the discretization of its MRI spectrum, as well as by applying the appropriate boundary conditions. A second order Duffing-like amplitude equation for the initially unstable MRI modes is derived. The solutions of that equation exhibit bursty nonlinear oscillations with a constant amplitude that signifies the saturation level of the MRI. Those results are verified by a direct numerical solution of the full nonlinear reduced set of thin disk magnetohydrodynamics equations.

(Abridged) We analyse the stability and evolution of power-law accretion disc models. These have midplane densities that follow radial power-laws, and have either temperature or entropy distributions that are power-law functions of cylindrical radius. We employ two different hydrodynamic codes to perform 2D-axisymmetric and 3D simulations that examine the long-term evolution of the disc models as a function of the power-law indices of the temperature or entropy, the thermal relaxation time of the fluid, and the viscosity. We present a stability analysis of the problem that we use to interpret the simulation results. We find that disc models whose temperature or entropy profiles cause the equilibrium angular velocity to vary with height are unstable to the growth of modes with wavenumber ratios |k_R/k_Z| >> 1 when the thermodynamic response of the fluid is isothermal, or the thermal evolution time is comparable to or shorter than the local dynamical time scale. These discs are subject to the Goldreich-Schubert-Fricke (GSF) or `vertical shear' linear instability. Development of the instability involves excitation of vertical breathing and corrugation modes in the disc, with the corrugation modes in particular being a feature of the nonlinear saturated state. Instability operates when the dimensionless disc kinematic viscosity nu < 10^{-6} (Reynolds numbers Re>H c_s/nu > 2500). In 3D the instability generates a quasi-turbulent flow, and the Reynolds stress produces a fluctuating effective viscosity coefficient whose mean value reaches alpha ~ 6 x 10^{-4} by the end of the simulation. The vertical shear instability in disc models which include realistic thermal physics has yet to be examined. Should it occur, however, our results suggest that it will have significant consequences for their internal dynamics, transport properties, and observational appearance.

We extend exploration of potential vorticity instabilities in cold astrophysical disks whose mean states are baroclinic. In particular, we seek to demonstrate the potential existence of traditional baroclinic instabilities of meteorological studies in a simplified two-layer Philips Disk Model. Each disk layer is of constant but differing densities. The resulting mean azimuthal velocity profile shows a variation in the vertical direction implying that the system is baroclinic in the mean state. The stability of the system is treated in the context of disk shallow water theory wherein azimuthal disturbances are much longer than the corresponding radial or vertical scales. The normal-mode problem is solved numerically using two different methods. The results of a symmetric single layer barotropic model is considered and it is found that instability persists for models in which the potential vorticity profiles are not symmetric, consistent with previous results. The instaiblity is interpreted in terms of interacting Rossby waves. For a two layer system in which the flow is fundamentally baroclinic we report here that instability takes on the form of mixed barotropic-baroclinic type: instability occurs but it qualitatively follows the pattern of instability found in the barotropic models. Instability arises because of the phase locking and interaction of the Rossby waves between the two layers. The strength of the instability weakens as the density contrast between layers increases. (For full abstract see article.)

The Rossby wave instability in astrophysical disks is as a potentially important mechanism for driving angular momentum transport in disks. We aim to understand this instability in an approximate three-dimensional disk model environment which we assume to be a single homentropic annular layer we analyze using disk shallow-water theory. We consider the normal mode stability analysis of two kinds of radial profiles of the mean potential vorticity: The first type is a single step and the second kind is a symmetrical step of finite width describing either a localized depression or peak of the mean potential vorticity. For single potential vorticity steps we find there is no instability. There is no instability when the symmetric step is a localized peak. However, the Rossby wave instability occurs when the symmetrical step profile is a depression, which, in turn, corresponds to localized peaks in the mean enthalpy profile. This is in qualitative agreement with previous two-dimensional investigations of the instability. For all potential vorticity depressions, instability occurs for regions narrower than some maximum radial length scale. We interpret the instability as resulting from the interaction of at least two Rossby edgewaves. We identify the Rossby wave instability in the restricted three-dimensional framework of disk shallow water theory. Additional examinations of generalized barotropic flows are needed. Viewing disk vortical instabilities from the conceptual perspective of interacting edgewaves can be useful.

The short time evolution of three dimensional small perturbations is studied. Exhibiting spectral asymptotic stability, thin discs are nonetheless shown to host intensive hydrodynamical activity in the shape of non modal growth of initial small perturbations. Two mechanisms that lead to such behavior are identified and studied, namely, non-resonant excitation of vertically confined sound waves by stable planar inertia-coriolis modes that results in linear growth with time, as well as resonant coupling of those two modes that leads to a quadratic growth of the initial perturbations. It is further speculated that the non modal growth can give rise to secondary strato-rotational instabilities and thus lead to a new route to turbulence generation in thin discs.

Simulations of astrophysical disks in the shearing box that are subject to the magnetorotational instability (MRI) show that activity appears to be reduced as the magnetic Prandtl number (Pm) is lowered. On the other hand, calculations for laboratory experiments show that saturation is achieved through modification of the background shear for Pm << 1. Guided by the results of calculations appropriate for laboratory experiments when Pm is very low, the axisymmetric stability of inviscid disturbances in a shearing box model immersed in a constant vertical background magnetic field is considered under a variety of shear profiles and boundary conditions in order to evaluate the hypothesis that modifications of the shear bring about saturation of the instability. It is found that the emergence and stability of the MRI is sensitive to the boundary conditions adopted. Channel modes do not appear to be stabilized through any modification of the background shear whose average remains Keplerian. However, systems that have non-penetrative boundaries can saturate the MRI through modification of the background shear. Conceptually equating the qualitative results from laboratory experiments to the conditions in a disk may therefore be misleading.

Thin viscous Keplerian accretion disks are considered asymptotically stable, even though they can show significant dynamic activity on short timescales. In this paper the dynamics of non-axisymmetric hydrodynamical disturbances of disks are investigated analytically building upon the steady state three-dimensional structure and evolution of axisymmetric perturbations explored in previous work. Assuming a polytropic equation of state solutions are found by means of an asymptotic expansion in the small parameter measuring the ratio of the disk thickness to characteristic radius. In-depth analysis shows that every perturbation that disturbs the radial velocity induces significant transient growth in the (acoustic) energy of the evolving disturbance. This effect is most evident in the density and vertical velocity. The transient growth observed is tied to the non-separable nature of the solutions where, in particular, pattern evolution is controlled by a similarity variable composed of the radial coordinate and time. This leads to growing winding perturbations that display successive radial peaks and troughs. We argue that these transient non-axisymmetric structures may precipitate secondary instabilities which, consequently, may be a critical element for a new alternative picture of turbulence arousal in non-magnetized astrophysical disks.

This paper discusses new perspectives and approaches to the problem of disk dynamics where, in this study, we focus on the effects of viscous instabilities influenced by boundary effects. The Boussinesq approximation of the viscous large shearing box equations is analyzed in which the azimuthal length scale of the disturbance is much larger than the radial and vertical scales. We examine the stability of a non-axisymmetric potential vorticity mode, i.e. a PV-anomaly. in a configuration in which buoyant convection and the strato-rotational instability do not to operate. We consider a series of boundary conditions which show the PV-anomaly to be unstable both on a finite and semi-infinite radial domains. We find these conditions leading to an instability which is the disk analog of Tollmien-Schlichting waves. When the viscosity is weak, evidence of the instability is most pronounced by the emergence of a vortex sheet at the critical layer located away from the boundary where the instability is generated. For some boundary conditions a necessary criterion for the onset of instability for vertical wavelengths that are a sizable fraction of the layer's thickness and when the viscosity is small is that the appropriate Froude number of the flow be greater than one. This instability persists if more realistic boundary conditions are applied, although the criterion on the Froude number is more complicated. The unstable waves studied here share qualitative features to the instability seen in rotating Blasius boundary layers. The implications of these results are discussed. An overall new strategy for exploring and interpreting disk instability mechanisms is also suggested.

A scaling argument is presented that leads to a shallow water theory of non-axisymmetric disturbances in annular sections of thin Keplerian disks. To develop a theoretical construction that will aid in physically understanding the relationship of known two-dimensional vortex dynamics to their three-dimensional counterparts in Keplerian disks. Using asymptotic scaling arguments varicose disturbances of a Keplerian disk are considered on radial and vertical scales consistent with the height of the disk while the azimuthal scales are the full $2\pi$ angular extent of the disk. The scalings lead to dynamics which are radially geostrophic and vertically hydrostatic. It follows that a potential vorticity quantity emerges and is shown to be conserved in a Lagrangian sense. Uniform potential vorticity linear solutions are explored and the theory is shown to contain an incarnation of the strato-rotational instability under channel flow conditions. Linearized solutions of a single defect on an infinite domain is developed and is shown to support a propagating Rossby edgewave. Linear non-uniform potential vorticity solutions are also developed and are shown to be similar in some respects to the dynamics of strictly two-dimensional inviscid flows. Based on the framework of this theory, arguments based on geophysical notions are presented to support the assertion that the strato-rotational instability is in a generic class of barotropic/baroclinic potential vorticity instabilities. Extensions of this formalism are also proposed. The shallow water formulation achieved by the asymptotic theory developed here opens a new approach to studying disk dynamics.

The hydrodynamic response of the inviscid small shearing box model of a midplane section of a rotationally supported astrophysical disk is examined. An energy functional ${\cal E}$ is formulated for the general nonlinear problem. It is found that the fate of disturbances is related to the conservation of this quantity which, in turn, depends on the boundary conditions utilized: ${\cal E}$ is conserved for channel boundary conditions while it is not conserved in general for shearing box conditions. Linearized disturbances subject to channel boundary conditions have normal-modes described by Bessel Functions and are qualitatively governed by a quantity $\Sigma$ which is a measure of the ratio between the azimuthal and vertical wavelengths. Inertial oscillations ensue if $\Sigma >1$ - otherwise disturbances must in general be treated as an initial value problem. We reflect upon these results and offer a speculation.

Most of our knowledge on the nonlinear development of the magneto-rotational instability (MRI) relies on the results of numerical simulations employing the shearing box (SB) approximation. A number of difficulties arising from this approach have recently been pointed out in the literature. We thoroughly examine the effects of the assumptions made and numerical techniques employed in SB simulations. This is done in order to clarify and gain better understanding of those difficulties as well as of a number of additional serious problems, raised here for the first time, and of their impact on the results. Analytical derivations and estimates as well as comparative analysis to methods used in the numerical study of turbulence are used. Numerical experiments are performed to support some of our claims and conjectures. The following problems, arising from the (virtually exclusive) use of the SB simulations as a tool for the understanding and quantification of the nonlinear MRI development in disks, are analyzed and discussed: (i) inconsistencies in the application of the SB approximation itself; (ii) the limited spatial scale of the SB; (iii) the lack of convergence of most ideal MHD simulations; (iv) side-effects of the SB symmetry and the non-trivial nature of the linear MRI; (v) physical artifacts arising on the too small box scale due to periodic boundary conditions. The computational and theoretical challenge posed by the MHD turbulence problem in accretion disks cannot be met by the SB approximation, as it has been used to date. A new strategy to confront this challenge is proposed, based on techniques widely used in numerical studies of turbulent flows - developing (e.g., with the help of local numerical studies) a sub-grid turbulence model and implementing it in global calculations.

We study the saturation near threshold of the axisymmetric magnetorotational instability (MRI) of a viscous, resistive, incompressible fluid in a thin-gap Taylor-Couette configuration. A vertical magnetic field, Keplerian shear and no-slip, conducting radial boundary conditions are adopted. The weakly non-linear theory leads to a real Ginzburg-Landau equation for the disturbance amplitude, like in our previous idealized analysis. For small magnetic Prandtl number (P \ll 1), the saturation amplitude scales as P^{2/3} while the magnitude of angular momentum transport scales as P^{4/3}. The difference with the previous scalings (~P^{1/2} and P respectively) is attributed to the emergence of radial boundary layers. Away from those, steady-state non-linear saturation is achieved through a modest reduction in the destabilizing shear. These results will be useful to understand MRI laboratory experiments and associated numerical simulations.

A 2D model representing the dynamical interaction of dust and gas in a planetary channel is explored. The two components are treated as interpenetrating fluids in which the gas is treated as a Boussinesq fluid while the dust is treated as pressureless. The only coupling between both fluid states is kinematic drag. The channel gas experiences a temperature gradient in the spanwise direction and it is adverse the constant force of gravity. The latter effects only the gas and not the dust component which is considered to free float in the fluid. The channel is also considered on an f-plane so that the background vorticity gradient can cause any emerging vortex structure to drift like a Rossby wave. A linear theory analysis is explored and a nonlinear amplitude theory is developed for disturbances of this arrangement. It is found that the presence of the dust can help generate and shape emerging convection patterns and dynamics in the gas so long as the state of the gas exceeds a suitably defined Rayleigh number appropriate for describing drag effects. In the linear stage the dust particles collect quickly onto sites in the gas where the vorticity is minimal, i.e. where the disturbance vorticity is anticylonic which is consistent with previous studies. The nonlinear theory shows that, in turn, the local enhancement of dust concentration in the gas effects the vigor of the emerging convective roll by modifying the effective local Rayleigh number of the fluid. It is also found that without the f-plane approximation built into the model the dynamics there is an algebraic runaway caused by unrestrained growth in the dust concentration. The background vorticity gradient forces the convective roll to drift like a Rossby wave and this causes the dust concentration enhancements to not runaway.

We show by means of a perturbative weakly nonlinear analysis that the axisymmetric magnetorotational instability (MRI) of a viscous, resistive, incompressible rotating shear flow in a thin channel gives rise to a real Ginzburg-Landau equation for the disturbance amplitude. For small magnetic Prandtl number (${\cal P}_{\rm m}$), the saturation amplitude is $\propto \sqrt{{\cal P}_{\rm m}}$ and the resulting momentum transport scales as ${\cal R}^{-1}$, where $\cal R$ is the {\em hydrodynamic} Reynolds number. Simplifying assumptions, such as linear shear base flow, mathematically expedient boundary conditions and continuous spectrum of the vertical linear modes, are used to facilitate this analysis. The asymptotic results are shown to comply with numerical calculations using a spectral code. They suggest that the transport due to the nonlinearly developed MRI may be very small in experimental setups with ${\cal P}_{\rm m} \ll 1$.

The linear normal-mode stratorotational instability (SRI) is analytically reexamined in the inviscid limit where the length scales of horizontal disturbances are large compared their vertical and radial counterparts. Boundary conditions different than channel walls are also considered. This quasi-hydrostatic, semi-geostrophic (QHSG) approximation allows one to examine the effect of a vertically varying Brunt-Vaisaila frequency, $N^2$. It is found that the normal-mode instability persists when $N^2$ increases quadratically with respect to the disc vertical coordinate. However we also find that the SRI seems to exist in this inviscid QHSG extreme only for channel wall conditions: when one or both of the reflecting walls are removed there is no instability in the asymptotic limit explored here. It is also found that only exponential-type SRI modes (as defined by Dubrulle et al. 2005) exist under these conditions. These equations also admit non-normal mode behaviour. Fixed Lagrangian pressure conditions on both radial boundaries predicts there to be no normal mode behaviour in the QHSG limit. The mathematical relationship between the results obtained here and that of the classic Eady (1949) problem for baroclinic instability is drawn. We conjecture as to the mathematical/physical nature of the SRI. The general linear problem, analyzed without approximation in the context of the Boussinesq equations, admits a potential vorticity-like quantity that is advectively conserved by the shear. Its existence means that a continuous spectrum \emph{is a generic feature of this system}. It also implies that in places where the Brunt-Vaisaila frequency becomes dominant the linearized flow may two-dimensionalize by advectively conserving its vertical vorticity.

The purpose of this paper is to explore how accretion discs manifest the phenomenon of transient growth on a global scale. We investigate analytically the time response of a thin accretion disc to particular axisymmetric perturbations. To facilitate an analytical treatment we replace the energy equation with a general polytropic assumption. The asymptotic expansion of Klu\'zniak & Kita (2000), which extended the method of Regev (1983) to a full steady polytropic disc (with $n=3/2$), is further developed and implemented for both the steady (for any polytropic index) and time-dependent problems. The spatial form and temporal behaviour of selected dynamical disturbances are studied in detail. We identify the perturbation space which leads to transient growth and provide analytical solutions which manifest this expected transient growth behaviour. Three terms (physical causes) responsible for the appearance of transient growth are identified. Two depend explicitly on the viscosity while the third one is relevant also for inviscid discs. The main conclusion we draw is that the phenomenon of transient growth exists in discs on a global scale.

The linear stability of accretion disks is revisited. The governing equations are expanded asymptotically and solved to first order in the expansion parameter $\epsilon$ defined by the ratio of the disk's vertical thickness to its radial extent. Algebraically growing solutions are found for global perturbations on the radial accretion flow of thin inviscid compressible Keplerian disks. The algebraic temporal behavior is exhibited in the vertical velocities and the thermodynamic variables and has the form $t\sin\Omega_0 t$ locally in the disk where $\Omega_0$ is the Keplerian rotation rate. The physical implications and relations to the Solberg-Hoiland stability criteria are discussed.

We present here both analytical and numerical results of hydrodynamic stability investigations of rotationally supported circumstellar flows using the shearing box formalism. Asymptotic scaling arguments justifying the shearing box approximation are systematically derived, showing that there exist two limits which we call small shearing box (SSB) and large shearing box (LSB). The physical meaning of these two limits and their relationship to model equations implemented by previous investigators are discussed briefly. Two dimensional (2D) dynamics of the SSB are explored and shown to contain transiently growing (TG) linear modes, whose nature is first discussed within the context of linear theory. The fully nonlinear regime in 2D is investigated numerically for very high Reynolds (Re) numbers. Solutions exhibiting sustained dynamics are found and manifest episodic but recurrent TG behavior and these are associated with the formation and long-term survival of coherent vortices. The life-time of this spatio-temporal complexity depends on the Re number and the strength and nature of the initial disturbance. We show results for a case in which the dynamical activity persists for the entire duration of the simulation (hundreds of box orbits). Some combinations of the Re number and the initial perturbation spectrum (and strength) show the transiently growing phenomenon to ultimately fade away or, in some instances, to be absent altogether. Because the SSB approximation used here is equivalent to a 2D incompressible flow, the dynamics can not depend on the Coriolis force. Therefore, three dimensional (3D) simulations are needed in order to decide if this force indeed suppresses nonlinear hydrodynamical instability in rotationally supported disks in the shearing box.

We study the weakly nonlinear evolution of acoustic instability of a plane- parallel polytrope with thermal dissipation in the form of Newton's law of cooling. The most unstable horizontal wavenumbers form a band around zero and this permits the development of a nonlinear pattern theory leading to a complex Ginzburg-Landau equation (CGLE). Numerical solutions for a subcritical, quintic CGLE produce vertically oscillating, localized structures that resemble the oscillons observed in recent experiments of vibrated granular material.