Research Opportunities
I am always looking for bright, enthusiastic, self-motivated graduate students to work with me in the Earth and Planetary Sciences at McGill University in Montreal, Canada. An undergraduate background in one of the following areas is preferred: geophysics, climate physics, general physics, mathematics, computer science, statistics and earth sciences with numerical modeling experience. If you are interested in working with me, please get in touch by emailing me at natalya.gomez@mcgill.ca. (Note that I may take quite a while to respond, as I tend to start reaching out to interested candidates all at once.) Topics for graduate research could be related, but not limited to:
|
|
Research in the News
Here are links to some articles for the general public related to my research:
- The Influence of Earth Deformation and Gravity on Future Antarctic Ice Loss
- Ice Sheets in Peril? Gravity to the Rescue
- Popsicle Earth
- Antarctic Ice Loss and Sea Level Rise
- Collapse of the West Antarctic Ice Sheet (1), (2), (3)
Sea-Level Change

Understanding Planet Earth gives an introductory overview of topics in Earth and Planetary Sciences for non-specialists. One might think that if an ice sheet melts, sea level will rise evenly everywhere, as if the oceans were filling up like a bathtub. In reality, gravitational, Earth deformational and rotational effects following the retreat and advance of ice sheets and mountain glaciers lead to markedly non-uniform patterns of sea-level change. These patterns are often referred to as "sea-level fingerprints" because they are unique to the geometry and timescale of the changes in grounded ice distribution.
As an example, the figure on the right shows predictions of the sea-level change immediately following the rapid collapse of the West Antarctic Ice Sheet (left) and marine sectors of the East Antarctic Ice Sheet (right). The plots are normalized by the global average value of sea-level rise, so that they can be scaled up or down for any amount of ice sheet retreat. If you melted enough ice from the West Antarctic to raise sea level by 1 meters assuming that the melt water were spread evenly over the oceans, the US East Coast would actually see 1.35 meters of sea-level rise. Similarly, if you melted enough ice for 5 meters of global average sea-level rise, the US East Coast would see 6.5 meters of sea-level rise (top left plot). In general, the patterns are characterized by a greater-than-average sea level rise far away from the location of ice loss and a large sea-level fall in the vicinity of a retreating ice sheet (the opposite applies for an advancing ice sheet). The bottom plots show polar projections of the same fields in the top plots, zooming in on the region of ice loss.
Sea-level modeling can give insight into a wide variety of problems in climate and solid Earth science, including (but not limited to): understanding the response of the planet to ice-age and modern climate changes; predicting the contributions from ice sheet retreat to future sea-level rise; and processing and interpreting measurements of sea-level change, glacial isostatic adjustment and solid deformation.
RELATED PUBLICATIONS:
Hay, C., Lau, H., Gomez, N., Austermann, J., Powell, E., Mitrovica, J. X., Latychev, K., and Wiens, D. (2017). Sea-level fingerprints in a region
of complex Earth structure: The case of WAIS. Journal of Climate, 30, 1881-1892. doi: 10.1175/JCLI-D-16-0388.1
Gomez, N., Gregoire, L.J., Mitrovica, J.X., Payne, A.J. (2015). Laurentide-Cordilleran Ice Sheet saddle collapse as a contribution to
meltwater pulse 1A, Geophys. Res. Lett., 42, doi: 10.1002/2015GL063960.
Hay, C., Mitrovica, J.X., Gomez, N., Creveling, J.R., Austermann, J. and Kopp, R.. (2014) The sea-Level Fingerprint of Ice-Sheet Collapse
During Interglacial Periods. Quat. Sci. Rev. 87: 60–69. doi: 10.1016/j.quascirev.2013.12.022.
Mitrovica, J.X., Gomez, N., Morrow, E., Hay, C., Latychev, K. and Tamisiea, M.E.. (2011) On the robustness of predictions of sea level
fingerprints. Geophys. J. Int. 187, 729-742. doi: 10.1111/j.1365-246X.2011.05090.x.
Gomez, N., Mitrovica, J.X., Tamisiea, M. and Clark, P.. (2009) A New Projection of Sea-Level Change in Response to a Collapse of the
Antarctic Ice Sheet. Geophysical Journal International. 180, 623-634. doi: 10.1111/j.1365-246X.2009.04419.x.
Mitrovica, J.X., Gomez, N., Clark, P. (2009) The Sea-Level Fingerprint of West Antarctic Collapse. Science. 323, 753. doi: 10.1126/science.1166510.
As an example, the figure on the right shows predictions of the sea-level change immediately following the rapid collapse of the West Antarctic Ice Sheet (left) and marine sectors of the East Antarctic Ice Sheet (right). The plots are normalized by the global average value of sea-level rise, so that they can be scaled up or down for any amount of ice sheet retreat. If you melted enough ice from the West Antarctic to raise sea level by 1 meters assuming that the melt water were spread evenly over the oceans, the US East Coast would actually see 1.35 meters of sea-level rise. Similarly, if you melted enough ice for 5 meters of global average sea-level rise, the US East Coast would see 6.5 meters of sea-level rise (top left plot). In general, the patterns are characterized by a greater-than-average sea level rise far away from the location of ice loss and a large sea-level fall in the vicinity of a retreating ice sheet (the opposite applies for an advancing ice sheet). The bottom plots show polar projections of the same fields in the top plots, zooming in on the region of ice loss.
Sea-level modeling can give insight into a wide variety of problems in climate and solid Earth science, including (but not limited to): understanding the response of the planet to ice-age and modern climate changes; predicting the contributions from ice sheet retreat to future sea-level rise; and processing and interpreting measurements of sea-level change, glacial isostatic adjustment and solid deformation.
RELATED PUBLICATIONS:
Hay, C., Lau, H., Gomez, N., Austermann, J., Powell, E., Mitrovica, J. X., Latychev, K., and Wiens, D. (2017). Sea-level fingerprints in a region
of complex Earth structure: The case of WAIS. Journal of Climate, 30, 1881-1892. doi: 10.1175/JCLI-D-16-0388.1
Gomez, N., Gregoire, L.J., Mitrovica, J.X., Payne, A.J. (2015). Laurentide-Cordilleran Ice Sheet saddle collapse as a contribution to
meltwater pulse 1A, Geophys. Res. Lett., 42, doi: 10.1002/2015GL063960.
Hay, C., Mitrovica, J.X., Gomez, N., Creveling, J.R., Austermann, J. and Kopp, R.. (2014) The sea-Level Fingerprint of Ice-Sheet Collapse
During Interglacial Periods. Quat. Sci. Rev. 87: 60–69. doi: 10.1016/j.quascirev.2013.12.022.
Mitrovica, J.X., Gomez, N., Morrow, E., Hay, C., Latychev, K. and Tamisiea, M.E.. (2011) On the robustness of predictions of sea level
fingerprints. Geophys. J. Int. 187, 729-742. doi: 10.1111/j.1365-246X.2011.05090.x.
Gomez, N., Mitrovica, J.X., Tamisiea, M. and Clark, P.. (2009) A New Projection of Sea-Level Change in Response to a Collapse of the
Antarctic Ice Sheet. Geophysical Journal International. 180, 623-634. doi: 10.1111/j.1365-246X.2009.04419.x.
Mitrovica, J.X., Gomez, N., Clark, P. (2009) The Sea-Level Fingerprint of West Antarctic Collapse. Science. 323, 753. doi: 10.1126/science.1166510.
Ice Sheet - Sea Level Interactions
Marine ice sheets like the West Antarctic ice sheet rest on bedrock that lies below the sea surface elevation and interact with the surrounding ocean. The ice loss across the grounding line of a marine ice sheet depends strongly on the depth of water there, such that if the grounding line retreats a little bit into deeper water, ice loss across the grounding line will increase leading to further, runaway retreat (see top frames of videos A and B below, showing a simplified, 1-D ice sheet-shelf model retreating down reversed beds with different reversed slopes). Past models of marine ice sheets have considered the impact of global average sea-level rise following ice-sheet retreat, but as we saw above, sea-level actually falls at the edge of a retreating ice sheet. Collaborators and I have shown that these gravitationally self-consistent sea-level changes have a stabilizing influence on marine ice sheets, acting to slow down and in some cases halt grounding line migration. The bottom frames below show simulations in which a sea-level coupling has been incorporated into the ice sheet model, demonstrating the slowdown (A) and halting (B) of the grounding line as compared to the upper frames in which no sea-level changes have been incorporated into the ice-sheet model.
(A) Simulation on a steeper reversed bedslope of -0.7 m/km in which the sea-level fall at the grounding line acts to compensate for some of the deepening of the bed, slowing down the ice sheet retreat.
|
(B) Simulation on a shallower reversed bedslope of -0.3 m/km in which the sea-level fall at the grounding line is large enough to fully compensate for the deepening of the bed and stop the ice-sheet retreat.
|
In recent work, we coupled a more realistic, 3-D ice sheet-shelf model to a global sea-level model that incorporates viscoelastic deformation of the solid Earth and applied the coupled model to simulate the evolution of the Antarctic Ice Sheet over the last deglaciation. These results demonstrated the stabilizing influence of the sea-level coupling in a more realistic setting, confirming the conclusions of our earlier work with simplified ice-sheet models. We can use the coupled ice sheet - sea level model as a tool to take an integrated view of the evolution of ice sheets, sea level and the solid Earth.
RELATED PUBLICATIONS:
Gomez, N., Latychev, K. and Pollard, D.. (2018) A coupled ice sheet-sea level model incorporating 3D Earth structure: Variations in
Antarctica during the last deglacial retreat. J. Climate, 31, 4041-4054. doi: 10.1175/JCLI-D-17-0352.1.
Pollard, D., Gomez, N., DeConto, R. and Han, H.. (2018) Estimating modern elevations of Pliocene shorelines using a coupled ice sheet-
Earth-sea level model. J. Geophys. Res. - Earth Surface. doi: 10.1029/2018JF004745.
Pollard, D., Gomez, N. and DeConto, R.. (2017). Variations of the Antarctic Ice Sheet in a coupled ice- sheet – Earth - sea level model:
sensitivity to viscoelastic Earth structure. J. Geophys. Res. Earth Surface., 122. doi: 10.1002/2017JF004371.
Gomez, N., Pollard, D., and Holland, D.. (2015) Sea-level feedback lowers projections of future Antarctic ice-sheet mass loss. Nature
Communications 6, 8798. doi: 10.1038/ncomms9798.
Gomez, N., Pollard, D., and Mitrovica, J.X.. (2013) A 3-D coupled ice sheet - sea level model applied to Antarctica through the last
40 ky. Earth Planet. Sci. Lett. 384, 88-99. doi: 10.1016/j.epsl.2013.09.042.
Gomez, N., Pollard, D., Mitrovica, J.X., Huybers, P. and Clark, P.U.. (2012) Evolution of a Coupled Marine Ice Sheet - Sea Level Model. J.
Geophys. Res. 117, 9 pp. doi: 10.1029/2011JF002128.
Gomez, N., Mitrovica, J.X., Huybers, P. and Clark, P.U.. (2010) Sea Level Change As a Stabilizing Influence on Marine Ice Sheets. Nature
Geoscience. 3, 850-853. doi: 10.1038/ngeo1012.
Gomez, N., Latychev, K. and Pollard, D.. (2018) A coupled ice sheet-sea level model incorporating 3D Earth structure: Variations in
Antarctica during the last deglacial retreat. J. Climate, 31, 4041-4054. doi: 10.1175/JCLI-D-17-0352.1.
Pollard, D., Gomez, N., DeConto, R. and Han, H.. (2018) Estimating modern elevations of Pliocene shorelines using a coupled ice sheet-
Earth-sea level model. J. Geophys. Res. - Earth Surface. doi: 10.1029/2018JF004745.
Pollard, D., Gomez, N. and DeConto, R.. (2017). Variations of the Antarctic Ice Sheet in a coupled ice- sheet – Earth - sea level model:
sensitivity to viscoelastic Earth structure. J. Geophys. Res. Earth Surface., 122. doi: 10.1002/2017JF004371.
Gomez, N., Pollard, D., and Holland, D.. (2015) Sea-level feedback lowers projections of future Antarctic ice-sheet mass loss. Nature
Communications 6, 8798. doi: 10.1038/ncomms9798.
Gomez, N., Pollard, D., and Mitrovica, J.X.. (2013) A 3-D coupled ice sheet - sea level model applied to Antarctica through the last
40 ky. Earth Planet. Sci. Lett. 384, 88-99. doi: 10.1016/j.epsl.2013.09.042.
Gomez, N., Pollard, D., Mitrovica, J.X., Huybers, P. and Clark, P.U.. (2012) Evolution of a Coupled Marine Ice Sheet - Sea Level Model. J.
Geophys. Res. 117, 9 pp. doi: 10.1029/2011JF002128.
Gomez, N., Mitrovica, J.X., Huybers, P. and Clark, P.U.. (2010) Sea Level Change As a Stabilizing Influence on Marine Ice Sheets. Nature
Geoscience. 3, 850-853. doi: 10.1038/ngeo1012.
Implications of GIA and Climate Change For Arctic Canada
The location of future shorelines in Arctic Canada will impact the livelihood of native communities as well as their claim to land and resources. In the flat terrain of the Hudson Bay region (see the photo below, taken near Fort Albany, Ontario), future shoreline migration can be expected to occur due to both glacial isostatic adjustment associated with the last Ice Age and sea-level changes associated with ongoing melting of the polar ice sheets and mountain glaciers in response to future climate warming. Greenland ice loss would result in a small gradient of sea level fall across the bay, whereas Antarctic ice loss would cause a relatively uniform sea level rise. Changing sea levels in turn affect ocean tides, which are highly sensitive to changes in water depths. Tides in Hudson Bay are among the world's largest, thus the region is particularly interesting to study the interactions between sea level and tide changes. For her M.Sc., Anna-Mireilla Hayden investigated the impacts of future ice sheet retreat on sea level and tides in the region.
The publications below describe work I was involved in to model future shoreline migration in the Hudson Bay region, explore its impact on native coastal communities, and educate young native Canadians about these issues.
The publications below describe work I was involved in to model future shoreline migration in the Hudson Bay region, explore its impact on native coastal communities, and educate young native Canadians about these issues.
RELATED PUBLICATIONS:
Tsuji, L.J.S., Daradich, A., Gomez, N., Hay, C., and Mitrovica, J.X.. (2016). Sea-Level Change in the Western James Bay Region of Sub-Arctic
Ontario: Emergent Land and Implications for Treaty No. 9. Arctic, 69, 1, 99-107. doi: 10.14430/arctic4542
Karagatzides, J.D., Kozlovic, D.R., De Iuliis, G., Liberda, E.N., General, Z., Liedtke, J., McCarthy, D.D., Gomez, N., Metatawabin, D. &
Tsuji, L.. (2011) Youth Environmental Science Outreach in the Mushkegowuk Territory of Subarctic Ontario, Canada, Applied
Environmental Education & Communication, 10, 201-210. doi: 10.1080/1533015X.2011.669684
Tsuji, L., Gomez, N., Mitrovica, J. (2009) Post-Glacial Adjustment and Global Warming in Sub-Arctic Canada: Implications for Islands
of the James Bay Region. Arctic. 64, 24 pages. doi: 10.14430/arctic176
Tsuji, L.J.S., Daradich, A., Gomez, N., Hay, C., and Mitrovica, J.X.. (2016). Sea-Level Change in the Western James Bay Region of Sub-Arctic
Ontario: Emergent Land and Implications for Treaty No. 9. Arctic, 69, 1, 99-107. doi: 10.14430/arctic4542
Karagatzides, J.D., Kozlovic, D.R., De Iuliis, G., Liberda, E.N., General, Z., Liedtke, J., McCarthy, D.D., Gomez, N., Metatawabin, D. &
Tsuji, L.. (2011) Youth Environmental Science Outreach in the Mushkegowuk Territory of Subarctic Ontario, Canada, Applied
Environmental Education & Communication, 10, 201-210. doi: 10.1080/1533015X.2011.669684
Tsuji, L., Gomez, N., Mitrovica, J. (2009) Post-Glacial Adjustment and Global Warming in Sub-Arctic Canada: Implications for Islands
of the James Bay Region. Arctic. 64, 24 pages. doi: 10.14430/arctic176
Energy and Environmental Issues

As a member of the Harvard Graduate Consortium on Energy and the Environment, I engaged in courses, seminars and research on the issues, consequences and strategies related to the future of energy and the environment. I am particularly interested in taking a broad, interdisciplinary perspective on topics related to climate change and sea-level rise, but I have also done research in other areas.
For example, I worked with Meredith Langstaff, another graduate student in the Consortium, to investigate the increase of renewable energy penetration in Aruba. In particular, we looked into ways to incentivize electricity load management technologies such as electric vehicles and electric thermal storage to accommodate intermittency issues associated with wind and solar energy. We had the opportunity to travel to Aruba to meet with many of the key players in the Energy and Government sectors and wrote a report for the Aruban Government on our work. An article from the Carbon War Room discusses Aruba’s progress towards renewable energy and sustainability.
NAVIGATE TO TOP OF PAGE
For example, I worked with Meredith Langstaff, another graduate student in the Consortium, to investigate the increase of renewable energy penetration in Aruba. In particular, we looked into ways to incentivize electricity load management technologies such as electric vehicles and electric thermal storage to accommodate intermittency issues associated with wind and solar energy. We had the opportunity to travel to Aruba to meet with many of the key players in the Energy and Government sectors and wrote a report for the Aruban Government on our work. An article from the Carbon War Room discusses Aruba’s progress towards renewable energy and sustainability.
NAVIGATE TO TOP OF PAGE