The Great Dark Spot (GDS-89) observed by Voyager 2 was the first of several large-scale vortices observed on Neptune, the most recent of which was observed in 2018 in the northern hemisphere (NDS-2018). Ongoing observations of these features are constraining cloud formation, drift, shape oscillations, and other dynamic properties. In order to effectively model these characteristics, an explicit calculation of methane cloud microphysics is needed. Using an updated version of the Explicit Planetary Isentropic Coordinate General Circulation Model (EPIC GCM) and its active cloud microphysics module to account for the condensation of methane, we investigate the evolution of large scale vortices on Neptune. We model the effect of methane deep abundance and cloud formation on vortex stability and dynamics. In our simulations, the vortex shows a sharp contrast in methane vapor density inside compared to outside the vortex. Methane vapor column density is analogous to optical depth and provides a more consistent tracer to track the vortex, so we use that variable over potential vorticity. We match the meridional drift rate of the GDS and gain an initial insight into the evolution of vortices in the northern hemisphere, such as the NDS-2018.

On 7th August 2019, an impact flash lasting $\sim1$s was observed on Jupiter. The video of this event was analysed to obtain the lightcurve and determine the energy release and initial mass. We find that the impactor released a total energy of $96-151$ kilotons of TNT, corresponding to an initial mass between $190-260$ metric tonnes with a diameter between $4-10$m. We developed a fragmentation model to simulate the atmospheric breakup of the object and reproduce the lightcurve. We model three different materials: cometary, stony and metallic at speeds of $60$, $65 $ and $70$ km/s to determine the material makeup of the impacting object. The slower cases are best fit by a strong, metallic object while the faster cases require a weaker material.

On June 2, 2016 at 10h56m UTC, a $-20.4 \pm 0.2$ magnitude superbolide was observed over Arizona. Fragments were located a few days later and the meteorites were given the name Dishchii'bikoh. We present analysis of this event based on 3 cameras and a multi-spectral sensor observations by the SkySentinel continuous fireball-monitoring camera network, supplemented by a dash cam footage and a fragmentation model. The bolide began its luminous flight at an altitude of $100.2 \pm 0.4$ km at coordinates $\phi = 34.555 \pm 0.002^\circ$N planetographic latitude and $\lambda = 110.459 \pm 0.002^\circ$W longitude, and it had a pre-atmospheric velocity of $17.4 \pm 0.3$ km/s. The calculated orbital parameters indicate that the meteoroid did not belong to any presently known asteroid family. From our calculations, the impacting object had an initial mass of $14.8 \pm 1.7$ metric tonnes with an estimated initial diameter of $2.03 \pm 0.12$ m.

We have investigated the 2009 July impact event on Jupiter using the ZEUS-MP 2 three-dimensional hydrodynamics code. We studied the impact itself and the following plume development. Eight impactors were considered: 0.5 km and 1 km porous (\rho = 1.760 g cm^{-3}) and non-porous (\rho = 2.700 g cm^{-3}) basalt impactors, and 0.5 km and 1 km porous (\rho = 0.600 g cm^{-3}) and non-porous \rho = 0.917 g cm^{-3}) ice impactors. The simulations consisted of these bolides colliding with Jupiter at an incident angle of \theta = 69 degrees from the vertical and with an impact velocity of v = 61.4 km s^{-1}. Our simulations show the development of relatively larger, faster plumes created after impacts involving 1 km diameter bodies. Comparing simulations of the 2009 event with simulations of the Shoemaker-Levy 9 events reveals a difference in plume development, with the higher incident angle of the 2009 impact leading to a shallower terminal depth and a smaller and slower plume. We also studied the amount of dynamical chaos present in the simulations conducted at the 2009 incident angle. Compared to the chaos of the SL9 simulations, where \theta is approximately 45 degrees, we find no significant difference in chaos at the higher 2009 incident angle.

We have studied plume formation after a Jovian comet impact using the ZEUS-MP 2 hydrodynamics code. The three-dimensional models followed objects with 500, 750, and 1000 meter diameters. Our simulations show the development of a fast, upward-moving component of the plume in the wake of the impacting comet that "pinches off" from the bulk of the cometary material ~50 km below the 1 bar pressure level, ~100 km above the depth of greatest mass and energy deposition. The fast-moving component contains about twice the mass of the initial comet, but consists almost entirely (>99.9%) of Jovian atmosphere rather than cometary material. The ejecta rise mainly along the impact trajectory, but an additional vertical velocity component due to buoyancy establishes itself within seconds of impact, leading to an asymmetry in the ejecta with respect to the entry trajectory. The mass of the upward-moving component follows a velocity distribution M(>v) approximately proportional to v^-1.4 (v^-1.6 for the 750 m and 500 m cases) in the velocity range 0.1 < v < 10 km/s.