True-color renderings of planets orbiting stars of various masses and temperatures, ranging from 2500 K to 8000 K, with the planets orbiting the cooler stars rotating with one side always facing the star. These planets were simulated with ExoPlaSim.
Most of my graduate research focused on the cold edge of the habitable zone. Traditionally, we talk about the habitable zone in terms of distance: the "inner edge" and the "outer edge". But there are a great many factors that affect habitability and climate beyond just distance from the star--CO2, atmospheric mass, water inventory, land masses, etc all play a role. So it makes more sense to talk about the cold edge, which describes the cold limit of habitability across all those parameters.
Relative vorticity (top) and wind speed (bottom) of a tidally-locked planet at different model layers in a high-resolution (75 km) ExoPlaSim simulation, captured in high-cadence output. The substellar point is in the middle.
In order to enable a lot of my research, I have had to make several modifications to the PlaSim climate model. I use PlaSim for the foundation of my work because it is extraordinarily fast (1 year in 45 seconds on 16 CPUs at 32x64x10 resolution), but PlaSim was designed for Earth, and cannot natively model a wide range of habitable exoplanets accurately. Years of modifications and additions to the model have resulted in the ExoPlaSim climate model, of which I am the primary maintainer, which extends PlaSim for planets orbiting stars unlike our Sun, and for different atmospheric pressures. ExoPlaSim also features a Python API which streamlines installation and model configuration, enabling rapid and efficient modelling with a very shallow learning curve. This should make this model accessible to a range of users, including observers with no modelling experience, while permitting experienced theorists and modellers to produce ensembles of models orders of magnitude larger than those possible with other GCMs. ExoPlaSim does not yet have all the physics necessary to model every physical process we will find on Earth-like exoplanets in the habitable zone, but its robustness, reliability, flexibility, and speed make it a good choice for parameter space studies and pathfinding experiments, where the expense of a larger model may not be justified without first testing the waters with ExoPlaSim.
3D renderings of the water content [kg/kg] in the atmosphere of a tidally-locked planet (the dayside is where the big purple region is; the nightside is the other half) and the associated cloud field, generated using output from the LMD Generic model.
Tidally-locked planets represent perhaps one of the greatest challenges in planetary science: they are at once our most promising targets in the search for life and the most attractive targets for follow-up observations, and at the same time also the most different from the planets in our solar system. About 70% of the stars in our galaxy are small, dim, low-mass stars, and planets in their habitable zones are almost all going to be tidally-locked, with one side always facing the star. Due to their proximity to their stars, short orbits, and the small size of those stars, they are the easiest terrestrial planets to observe. But when the planet has a temperate day-side and an inhospitable night-side, and when the very dynamics of the atmosphere operate in a different regime, it becomes difficult to generalize what we know about the habitability of Earth-like climates to tidally-locked planets. I use numerical models at various levels of complexity to probe how basic parameters affect the position of the cold edge of the habitable zone. Models I use include a 1D EBM, a modified version of PlaSim called ExoPlaSim, MITgcm, and the LMD Generic model. ExoPlaSim is the key to this work, because it both runs very quickly (1 year of climate in 45 seconds on 8 CPUs) and contains enough physics to give us useful information about a planet's climate. I use the other models to verify whether my results are robust, or model-dependent.
It's all fine and well if you can identify the conditions that give habitable climates, but what if those conditions are themselves unstable? Earth has several geological feedbacks which connect the atmospheric composition to the planet's interior, and depending on factors like continental configuration and volcanic outgassing rates, CO2 levels can adjust by large amounts over millions of years. Some of my work has focused on examining the consequences of this feedback when outgassing, the source of CO2, is weak. Carbon dioxide sequestration could then continue unabated, pulling the planet into a colder state. Some of our research suggests that this could lead to planets that cycle between being completely frozen and completely temperate, and some of our other research suggests that sometimes planets can get stuck in their frozen states. Not to worry, though--that same research also suggests that when that happens, there are regions of land that remain temperate, even though the rest of the planet is frozen. I use a modified version of PlaSim to compute how CO2 changes over time on different planets and identify whether or not these planets have stable climates.
We are at the cusp of being able to observe and characterize Earth-sized exoplanets. New instruments and bigger telescopes are allowing us to peer further, see smaller, and collect more detail. But our experience doing that with larger planets tells us that this will be an extraordinarily difficult endeavor--particularly because terrestrial climates are extraordinarily complex, with many difficult-to-predict highly-coupled nonlinear feedbacks. We are not yet at the stage where we can take a spectrum of a planet and determine what its climate must be like, and we may not be at a point where we can even try to solve that problem. We can however identify how different features of the climate affect that spectrum. I use the radiative transfer model SBDART together with output from my climate models to compute the reflected light spectra of various climate states, and then look for changes in the spectrum that could be attributed to those changes in climate. As a fun aside, because I compute the full reflected spectrum for each point on the planet, I also have everything I need to compute the color and brightness of each point--allowing me to construct an image of what my model would look like to the human eye, if seen from space.
A true-color rendering of a tidally-locked planet viewed at different phases, simulated with the LMD Generic model, using spectral data computed by SBDART.
Supervisor: Noboru Nakamura
In June 2018, I attended the Rossbypalooza climate modeling summer school at the University of Chicago. Two other students and I worked with Professor Noboru Nakamura on a follow-up to his Science paper on why the jet stream sometimes gets stuck in a holding-pattern of sorts--an event called 'blocking'. His paper proposed using reanalysis data that flow in the jet stream might obey effectively the same equation as traffic flow--that effectively, blocks are atmospheric traffic jams. Over the course of 2 weeks, we wrote and tested a 1D model of the jet stream using this equation to understand under what conditions blocks form, and how that might change as a result of climate change. Professor Nakamura is presenting some of our results at the 2018 AGU Fall Meeting.
Supervisor: Keith Vanderlinde
In the summer of 2016, as my second first-year project at UofT, I worked with Professor Keith Vanderlinde on writing software beamforming algorithms for the CHIME Pathfinder radio telescope. CHIME is a series of half-pipes, with radio antennae along the axis. Traditional telescopes only need one set of receivers at the focal point, because the dish does all the work of gathering the light and interfering it into an image at the receiver. With a telescope like CHIME or the smaller CHIME Pathfinder, a computer needs to compute that using a Fourier transform. By applying different offsets, however, the telescope can form many beams looking at different parts of the sky simultaneously. As a result, CHIME can see a strip of sky spannign almost all the way from the Southern horizon to the Northern horizon--and over the course of a day, as the Earth rotates, it sweeps out the entire sky. CHIME's original mission is to construct a 3D map of neutral hydrogen from redshifts 2.5 to 0.8, but it's also useful for detecting Fast Radio Bursts, or FRBS, which are short, bright, mysterious radio pulses from very distant sources. Because CHIME sees a huge portion of the sky and has to process it with computers, and needs to maintain a high cadence to catch FRBs, the beamforming algorithms have to be screamingly fast--so CHIME uses a huge GPU supercomputer on-site to process several terabytes of data per second. In collaboration with Cherry Ng and a few others, I helped with proof-of-concept and prototype beamforming algorithms, figuring out how to code extremely fast and efficient FFTs for the GPUs powering the CHIME Pathfinder.
The CHIME Pathfinder at the Dominion Radio Astrophysical Observatory near Penticton, BC. CHIME is a series of half-pipes, and uses software to form images of the sky. Photo credit Keith Vanderlinde, Dunlap Institute.
Supervisor: Paul Woodward
In the year between undergrad and grad school (2014-2015), I worked with Professor Paul Woodward at the University of Minnesota on trying to come up with a faster, more efficient way to compute implicit fluid transport in multi-dimensional fluid simulations, using Riemann invariants. We had very promising results in 1D cases, but have yet to find a way to robustly generalize that success to multiple dimensions.
Supervisors: Paul Woodward and William Dai
For my undergraduate astronomy thesis, I worked with Professor Paul Woodward at the University of Minnesota on coming up with a better, more efficient way to compute the Poisson equation for gravity on regular finite-volume cartesian grids. For large, self-gravitating systems like astrophysical disks (ranging from proto-planetary disks to galactic disks), gravity computations can be very expensive and inefficient in highly-parallelized simulations, since the gravitational force on a given position depends on knowledge of the entire system--leading to an N2 problem. Solving a discretized version of the Poisson equation for gravity allows one to iteratively express the problem as simply depending on a local source term and the values of the cell's immediate neighbors. We used a "red-black" relaxation algorithm coupled with manipulations of the domain geometry to reduce the expense and improve the efficiency of this approach for scalable highly-parallelized applications. I continued this work the summer after graduation with Dr. William Dai at the Los Alamos National Laboratory.
Supervisors: Aaron Breneman and Cynthia Cattell
From Summer 2011 until my graduation in Spring 2014, I worked with the Space Physics Lab at the University of Minnesota as an undergraduate research assistant. We studied electromagnetic waves and other phenomena in the solar wind, in particular as they interacted with Earth's magnetosphere. This work involved data from the STEREO spacecraft as well as the Van Allen Space Probes (originally the Radiation Belt Space Probes). I helped to analyze data, write analysis code, designed an early version of the website for our electric field instrument aboard the Van Allen probes, and wrote some end-user software for aquisition and analysis of data from the Van Allen probes.
- GCM Simulations of Unstable Climates in the Habitable Zone
- Habitable Snowballs: Temperate Land Conditions, Liquid Water, and Implications for CO2 Weathering
- Climate diversity in the solar-like habitable zone due to varying background gas pressure
- Climate uncertainties caused by unknown land distribution on habitable M-Earths (2nd author)
- Algorithms for FFT Beamforming Radio Interferometers
- CHIME FRB: An application of FFT beamforming for a radio telescope
- A Machine Learns to Predict the Stability of Tightly Packed Planetary Systems
- STEREO and Wind observations of intense cyclotron harmonic waves at the Earth's bow shock and inside the magnetosheath
- Explaining polarization reversals in STEREO wave data
- ExoPlaSim: Extending the Planet Simulator for Exoplanets
- Fundamental Challenges to Remote Sensing of Exo-Earths
- A Large Repository of 3D Climate Model Outputs for Community Analysis and Postprocessing
- The Long Night: Modeling the Climate of Westeros
- Preliminary Analysis of Planetary Characteristics, Dynamics, and Climates from the Systems Alliance Planetary Survey Catalogue
- Extreme Solar Systems IV - Reykjavik, Iceland, Aug 19-23 2019
- Emerging Researchers in Exoplanet Science V - Cornell University, Ithaca, June 17-19, 2019
- Comparative Climatology of Terrestrial Planets III - Lunar & Planetary Institute, Houston, August 27-30 2018
- Rossbypalooza Summer School - University of Chicago, Chicago, June 11-22 2018
- Emerging Researchers in Exoplanet Science III - Yale University, New Haven, June 12-13 2017
- Emerging Researchers in Exoplanet Science II - Carl Sagan Institute, Cornell University, Ithaca, June 13-14 2016
- GNASH: The anomalous metal-poor stars and convective-reactive nuclear astrophysics - University of Victoria, Victoria, BC, May 25-29, 2015
- NCSA Blue Waters Symposium for Petascale Science and Beyond - Sunriver Resort, Sunriver, Oregon, May 10-13 2015
- American Geophysical Union Fall Meeting - San Francisco, December 9-13 2013