Studying The Ripple Effects Of Shrinking Arctic Sea Ice


Scientists on the research vessel Akademik Fedorov spent a week or so setting up a network of scientific monitoring equipment up to about 25 miles from the MOSAiC ship.

Ravenna Koenig/NPR

Arctic sea ice is one of the most dramatic indicators of the changing climate. Ice cover on the Arctic Ocean is in some months about half what it was decades ago, and its thickness has shrunk, by some estimates 40%.

Changes in the ice may also mean a host of other changes, in the Arctic system and around the globe. To better understand this, scientists have frozen an icebreaker alongside an Arctic ice floe that they will observe for a whole year.

The project is called MOSAiC, for Multidisciplinary drifting Observatory for the Study of Arctic Climate. And the primary questions they’re trying to answer: what are the causes of diminishing Arctic ice, and what are the consequences?

At just about 5 degrees from the North Pole, ocean physicist Tim Stanton from the Naval Postgraduate School stands next to a hole in the ice, surrounded by boxes of tools and equipment.

“I’ve got to just get the ‘hair dryer,'” he says, eyeing two electrical connectors for a science buoy that need to be warmed up in the 18 degrees Fahrenheit temperature.

A hair dryer? He clarifies: “Well, it’s an electrical what-do-you-call-it… heat gun,” he says. “It will frizz your hair, that’s for sure!”

Stanton is in the middle of a grueling eight hour process to install the buoy about 15 miles from the spot where the MOSAiC ship, the German icebreaker Polarstern, is moored.

It’s part of a network of equipment that’s being distributed around the Polarstern and will operate independently throughout the next year. It will provide additional data to what’s being collected at the central research camp on the ice next to the ship.

Ocean physicist Tim Stanton with the buoy system he’s installing in the ice. The aim is to get a better sense of the ocean factors that may be driving Arctic ice melt.

Ravenna Koenig/NPR

The buoy is a big banana yellow device, with a whole bunch of scientific bells and whistles that hang below it in the water.

“The flux package mounts on here,” says Stanton, pointing to a cylindrical instrument with sensors on it that will run up and down a metal rail hanging vertically in the water. “And that’s what measures the transport of heat, salt and momentum in the water column.”

Stanton wants to collect data on those attributes of the ocean because he thinks it may help explain why sea ice is disappearing as fast as it is.

“At first glance it must be obvious, right? You add heat, you melt ice,” he says. “But it is so complicated.”

As more sea ice melts in the summertime, it’s contributing fresher water to the top of the ocean. The saltier ocean water, which sits lower because it’s more dense, can create a barrier that prevents the fresher water from going down.

If that top water is trapped near the surface all summer, Stanton thinks it can absorb a lot more heat from the sun, and lead to even more ice melting.

“You can get these fresh warm layers that, when a little bit of wind comes along, does a little bit of mixing, really melts the heck out of the ice,” he says.

Tim Stanton installing a science buoy with the help of student Rosalie McKay. The buoy will measure heat, salt and momentum in the upper layer of the ocean over the course of a year.

Ravenna Koenig/NPR

While Stanton is asking questions about things that are going on below the ice, other scientists are looking at things going on above it.

Jessie Creamean of Colorado State University, for example, is out on the ice testing a device that collects and counts tiny particles in the atmosphere called aerosols.

“Alright little aerosol sampler, do well today,” Creamean says, closing a pelican case about the size of a carry-on piece of luggage. She’s tested it before in Colorado, but today’s experiment is to see how well it does in the cold.

People may be most familiar with aerosols created by pressurized cans like hairspray, but that’s just one kind. Aerosols can also come from natural sources like dust, pollen, fungi, or sea salt, and they’re actually the seeds that clouds need to form and grow.

In the Arctic, scientists think that microbes in the ocean, like bacteria or algae, can generate aerosols. And Creamean hypothesizes that less ice on the Arctic Ocean could mean more aerosols getting blown from the water into the atmosphere, and seeding more clouds.

Scientist Jessie Creamean moves a portable aerosol sampler out onto the ice to test it in the cold conditions.

Ravenna Koenig /NPR

The mechanism for that could be twofold: through more sunlight getting to the ocean as sea ice decreases, and potentially causing more growth of microbes, and also through the increased contact between the ocean and atmosphere.

MOSAiC scientists are interested in clouds because they’re important for regulating temperature, similar to a thermostat. Depending on the season, whether clouds are over water or ice, and the properties of the clouds, they can wind up cooling or warming the earth below.

Scientist Jessie Creamean and her portable aerosol sampler in the lab on the research vessel Akademik Fedorov.

Ravenna Koenig/NPR

“That affects how much heat can basically help melt the sea ice, or it can actually reflect sunlight from the sea ice,” says Creamean. “So it has a big role in controlling how much sea ice we have here.”

Creamean and Stanton are among hundreds of scientists from different disciplines trying to better understand this changing region.

“We’re looking at the interactions in the system,” says Matthew Shupe, an atmospheric scientist with the University of Colorado and the National Oceanic and Atmospheric Administration, and one of the coordinators for the expedition.

“How the atmosphere interacts with the sea ice, how the ocean interacts with the sea ice, the ecosystem, the biogeochemical processes,” he says.

The overarching goal of collecting all this data is to improve the way the Arctic is represented in climate models. Those are the computer simulations scientists use to estimate things like how much the earth could warm in the next 50 years.

The better you reflect how reality works in simulation, the better a prediction you’ll get. But because so little is known about how the Arctic Ocean system works, Shupe says predictions for how the Arctic will respond to climate change vary significantly.

The primary questions MOSAiC is asking: what are the causes and consequences of diminishing Arctic sea ice?

Ravenna Koenig/NPR

“The Arctic is a place where the models agree the least,” he says. “So that tells us that we’re missing something.”

Projecting changes in the Arctic — such as when the Arctic Ocean will see its first ice-free summer — is obviously important for the local ecosystem, for Arctic communities, and for anyone interested in doing commercial activity in the region.

But this research will also help scientists figure out how changes in the Arctic will impact other places on earth. For example, it may contribute to scientists’ understanding of the possible connections between warming in the Arctic and extreme weather events at mid-latitudes.

“We need to understand the physics, and ultimately improve our models that can help answer those questions for us,” says Shupe.

It will also help scientists anticipate the speed at which the Greenland ice sheet could melt, raising global sea level, and improve projections for how much global temperature will rise in the coming years.

By drifting across the Arctic Ocean for the next year and observing how all the smaller pieces of the Arctic system fit together, scientists hope they can bring these big picture questions into clearer focus.

Across Arctic Canada, sea ice levels are again at record lows

Across Arctic Canada, sea ice levels are again at record lows

“We could be up to a month behind where we should be at this time of year”

“Freeze-up is very, very late this year,” says Gjoa Haven’s Willie Aglukkaq, who took this photo on Oct. 18. “Can still drive boats to mainland, which is very unusual for Gjoa Haven. Never seen the ocean open this late in my lifetime.” (Photo by Willie Aglukkaq)

By Dustin Patar

Last October, when Juuta Sarpinak of Igloolik set off on a seal hunt, he took his snowmobile. For the same trip exactly a year later, he’d need a boat.

Across Nunavut and the rest of the Arctic, the same scenario is playing out.

“It’s the latest I’ve ever seen it freeze up in my lifetime,” said Willie Aglukkaq, who’s lived in Gjoa Haven for nearly 40 years.

He’s right.

As of Oct. 15, sea ice levels hit 5.118 million square kilometres, making them the lowest on record, which began in 1968.

“Depending on the area, we could be up to a month behind where we should be at this time of year,” said Gilles Langis, senior ice forecaster with the meteorological service of Canada.

But it’s not the first time this year ice levels have been at record lows.

April also saw the least amount of sea ice ever recorded.

This was then followed by a record-setting July.

In August, the ice loss slowed.

“We were tied for the second-lowest [amount of ice] with 2007 and 2016, but quite a bit above 2012,” said Walt Meier, a senior research scientist at the National Snow and Ice Data Center.

“But now we’re actually below 2012 levels for this time of year because it’s a slow freeze-up—2012 recovered pretty quickly, but we’ve been very sluggish in growing ice this year.”

As of Oct. 15, sea ice levels dropped below those from 2012, which previously held the record. (Graph courtesy of the NSIDC)

With sea ice growth, the big driver is air temperature.

“If you have an area of open water and the temperatures above it are cooler than the water temperature, it’ll lose heat into the atmosphere and the water will start cooling off to the point where it reaches freezing,” said Langis.

But Arctic air temperatures through October have been unseasonably warm, much like the previous six months.

“It was a pretty extreme summer,” said Meier.

From April onward, every month has ranked in the top three warmest—in terms of Arctic air temperature—on record, with May and August now holding top spots.

According to Meier, because of this, there was an abundance of ice-free water throughout the summer that absorbed a lot of solar energy and warmed up quite a bit.

“That heat is just taking a while to dissipate into the atmosphere and for the ocean to cool enough to grow ice.”

(Photo courtesy of UNEP)

What’s going on this year is representative of something larger.

“Under the influence of global heating caused by human-induced greenhouse gases emissions, we have seen a sharp decrease in the extent of Arctic sea [ice] since 1979,” says Pascal Peduzzi, Director of GRID-Geneva, in a press release published by the UN Environment Programme (UNEP).

Much like the scenario currently playing out across the Arctic this year, declining sea ice has amplified Arctic warming over the last several decades.

According to the press release, “Temperatures increased by around 0.5°C per decade between 1982 and 2017, primarily due to increased absorbed solar radiation accompanying sea ice loss since 1979. This is twice as fast as the global average.”

This season isn’t an exception.

But while the year-to-year trend continues on, ice coverage this year is on the rebound.

“It’s just finally starting to freeze,” said Aglukkaq.

Over the last 10 days, since Oct. 20, sea ice coverage has increased by almost 1.5 million square kilometres, including the water near Gjoa Haven.

If the current pace of ice growth continues, 2019 may soon once again be above 2012 levels.

315 billion-tonne iceberg breaks off Antarctica

Image captionThe EU’s Sentinel-1 satellite system captured these before and after images

The Amery Ice Shelf in Antarctica has just produced its biggest iceberg in more than 50 years.

The calved block covers 1,636 sq km in area – a little smaller than Scotland’s Isle of Skye – and is called D28.

The scale of the berg means it will have to be monitored and tracked because it could in future pose a hazard to shipping.

Not since the early 1960s has Amery calved a bigger iceberg. That was a whopping 9,000 sq km in area.

Amery is the third largest ice shelf in Antarctica, and is a key drainage channel for the east of the continent.

The shelf is essentially the floating extension of a number of glaciers that flow off the land into the sea. Losing bergs to the ocean is how these ice streams maintain equilibrium, balancing the input of snow upstream.

So, scientists knew this calving event was coming. What’s interesting is that much attention in the area had actually been focussed just to the east of the section that’s now broken away.

This is a segment of Amery that has affectionately become known as “Loose Tooth” because of its resemblance in satellite images to the dentition of a small child. Both ice areas had shared the same rift system.

Loose ToothImage copyrightNASA
Image captionLoose Tooth pictured in the early 2000s. D28 is seen forming to the left

But although wobbly, Loose tooth is still attached. It’s D28 that’s been extracted.

“It is the molar compared to a baby tooth,” Prof Helen Fricker from the Scripps Institution of Oceanography told BBC News.

Prof Fricker had predicted back in 2002 that Loose Tooth would come off sometime between 2010 and 2015.

“I am excited to see this calving event after all these years. We knew it would happen eventually, but just to keep us all on our toes, it is not exactly where we expected it to be,” she said.

The Scripps researcher stressed that there was no link between this event and climate change. Satellite data since the 1990s has shown that Amery is roughly in balance with its surroundings, despite experiencing strong surface melt in summer.

“While there is much to be concerned about in Antarctica, there is no cause for alarm yet for this particular ice shelf,” Prof Fricker added.

Amery ice shefImage copyrightRICHARD COLEMAN/UTAS
Image captionAmery experiences a lot of summer surface melt, but the data indicates it is in equilibrium

The Australian Antarctic Division will however be watching Amery closely to see if it reacts at all. The division’s scientists have instrumentation in the region.

It’s possible the loss of such a big berg will change the stress geometry across the front of the ice shelf. This could influence the behaviour of cracks, and even the stability of Loose Tooth.

D28 is calculated to be about 210m thick and contains some 315 billion tonnes of ice.

The name comes from a classification system run by the US National Ice Center, which divides the Antarctic into quadrants.

The D quadrant covers the longitudes 90 degrees East to zero degrees, the Prime Meridian. This is roughly Amery to the Eastern Weddell Sea.

D28 is dwarfed by the mighty A68 berg, which broke away from the Larsen C Ice Shelf in 2017. It currently covers an area more than three times as big.

Nearshore currents and winds will carry D28 westwards. It’s likely to take several years for it to break apart and melt completely.

Running the numbers on an insane scheme to save Antarctic ice


It would take a lot. Like a real lot.

Antarctica's Pine Island Glacier sheds some icebergs. Could we... sort of... put them back?
Enlarge / Antarctica’s Pine Island Glacier sheds some icebergs. Could we… sort of… put them back?

Imagine, if you will, the engineers of the king’s court after Humpty Dumpty’s disastrous fall. As panicked men apparently competed with horses for access to the site of the accident, perhaps the engineers were scoping out scenarios, looking for a better method of reassembling the poor fellow. But presumably none of those plans worked out, given the dark ending to that fairy tale.

A recent study published in Science Advances might be relatable for those fairy tale engineers. Published by Johannes Feldmann, Anders Levermann, and Matthias Mengel at the Potsdam Institute for Climate Impact Research, the study tackles a remarkable question: could we save vulnerable Antarctic glaciers with artificial snow?

Keeping our cool

Antarctica’s ice is divided into two separate ice sheets by a mountain range, with the smaller but much more vulnerable West Antarctic Ice Sheet representing one of the biggest wildcards for future sea level rise. In 2014, a study showed that two of the largest glaciers within that ice sheet—known as the Pine Island Glacier and Thwaites Glacier—had likely crossed a tipping point, guaranteeing a large amount of future ice loss that would continue even if global warming were halted today.

Much of the bedrock beneath the West Antarctic Ice Sheet is actually below sea level, though it’s buried below kilometers of solid ice. This makes for situations where the bed beneath the ice slopes down as you go inland from the coast. That’s inherently unstable, and once a glacier starts retreating downslope, the invading water provides an increasing floating force that reduces the sliding friction that slows the seaward flow of ice.

In the case of the Pine Island and Thwaites Glaciers, it seems that this is exactly what’s happening. Although this process can take centuries to fully play out, this portion of the ice sheet contains enough ice to raise global sea level by more than a meter.

Is there some extraordinary measure that could prevent that loss and preserve these glaciers? It’s the kind of question people will often ask, and scientists (who know the scale of these things) generally ignore as implausible.

But in this case, the researchers decided to go wild. Using a computer model of the ice sheet, they simulated the effects of adding huge amounts of ice near the front of these two glaciers. The idea works like this: Where a glacier meets the sea, it transitions from grounded to floating. Behind this “grounding line,” the glacier sits on the bedrock and sediment beneath; in front it gets thinner and floats as an ice shelf. To preserve the glacier, you need to keep that grounding line from retreating downhill. Thicken the ice on the inland side of the grounding line, and the thickness of ice flowing over the line and into the ice shelf increases—its weight keeps the grounding line pinned in place.

This map shows bedrock elevation beneath the ice sheet, with the white box highlighting the area of the Pine Island and Thwaites Glaciers where snow would be added in this scenario.
Enlarge / This map shows bedrock elevation beneath the ice sheet, with the white box highlighting the area of the Pine Island and Thwaites Glaciers where snow would be added in this scenario.

The researchers played around with different amounts of ice added to the glaciers for different periods of time, ranging from 10 year treatments to 50 years. Spreading it out over a longer period could mean a less preposterous addition of ice each year, but they found that the total amount has to increase if you do it that way. So in the end, the scenario they selected was 7,400 billion tons of ice added over 10 years. That was enough to restabilize these glaciers, preventing their inexorable decline.

Two for one special

To put that into context, removing that much seawater from the ocean would lower global sea level by about 2 millimeters per year. Current total sea level rise is a little over 3 millimeters per year, so it would be like nearly halting sea level rise… by bailing water out of the ocean. We can call that a bonus positive.

This analysis is more about what it would take than what such a scheme would look like, but the basic options are to pump water up and hose it around—hoping it freezes quickly—or to freeze it into snow like the world’s most awkwardly located ski resort.

Here, the researchers transition to listing all the reasons this is impractical and all the negative impacts it could have. For starters, the seawater would have to be desalinated since salt would probably affect the physics and behavior of the ice. Simply pumping that much water up the 640 meters and spreading it over an area nearly the size of West Virginia would require the power of something like 12,000 wind turbines—and that’s without the very substantial energy requirements for desalination and snow-making.

“The practical realization of elevating and distributing the ocean water would mean an unprecedented effort for humankind in one of the harshest environments of the planet,” the researchers write.

The impacts on Antarctic ecosystems could also be huge. Pumping that water out of the sea near the coast would significantly alter the circulation of water, which might even become somewhat self-defeating, as it could bring more warm water up against the ice shelf, increasing melt.

In the Potsdam Institute’s press release, Levermann puts it this way: “The apparent absurdity of the endeavour to let it snow in Antarctica to stop an ice instability reflects the breath-taking dimension of the sea-level problem. Yet as scientists we feel it is our duty to inform society about each and every potential option to counter the problems ahead.”

And to be clear, this is in addition to halting climate change—the scenario the numbers are based on assumes the temperatures don’t keep rising. But as the alternative is eventual inundation of parts of the world’s coastal cities, an argument can be made that the cost could be worth paying. It least it gives us an idea just how hard it would be to put Humpty Dumpty back together again.

When is Arctic ice like a magnet?

Arctic meltponds
Arctic meltponds (Image courtesy: Kenneth Golden)

model that simulates magnets can also reproduce pools of water on Arctic sea ice from just one real-world measurement. Researchers in the US and UK adapted the Ising statistical model of phase transitions in ferromagnets to recreate melt-pond patterns. Characterizing the distribution of meltwater at small scales could improve climate models and predict ice loss better.

As climate changes, warming is expected to be especially rapid at high latitudes. In the Arctic, a lengthening melt season over the last few decades has already reduced the volume and extent of sea ice, with year-on-year ice-loss rates generally exceeding predictions.

One reason for the failure of models to predict this loss accurately might lie in how they account for melt ponds. The formation and growth of these ponds, which occurs at the transition between sea ice and open ocean, includes a positive feedback loop that makes the system especially sensitive.

“Pond evolution largely controls sea-ice albedo, a key parameter in climate modelling and one of the most important — and least understood — processes in determining the role of sea ice in the climate system,” says Kenneth Golden of the University of Utah, US.

Real-life and modelled meltponds

As melt ponds form on the sea-ice surface, highly-reflective ice and snow are replaced by darker pools of water. The meltwater absorbs more solar radiation and warms up more than the ice that it replaces; the extra energy can trigger additional melting in its surroundings. When the melt pond penetrates the full thickness of the ice, warmer ocean water floods in from below, accelerating the process.

Existing descriptions of pond formation in global climate models consider the overall volume of meltwater but not its surface distribution. Yet, because the albedo change caused by melt ponds is a surface process, and because the rate at which ponds conduct heat to their surroundings depends on their perimeter, understanding the rules governing melt pond sizes is crucial for climate modelling.

Sea-ice spin

The standard Ising model comprises a lattice of interacting particles, each of which is assigned a spin value that is either up or down. To capture the detailed geometry of melt ponds, Golden, with colleagues at Northumbria University, UK, the University of Dayton, US, and the University of Utah, US, created a version where the lattice represents the sea-ice surface, and each node a pixel of either ice or liquid water.

Starting with a random input configuration that does not resemble the real-world, each node interacts with its neighbours until the system settles on a local low-energy state. With a lattice spacing of 1 metre – the length scale over which Arctic ice exhibits significant topographic differences – the pattern that emerges from the energy-minimization process closely matches the melt-pond distribution seen in real life. For example, both real and modelled ponds scale in size according to the same power law, and both form more complex, fractal shapes when they grow larger than 100 sq. m.“The approach could ultimately lead to a framework for representing pattern formation occurring at spatial scales smaller than the grid spacing used in global climate models, which currently track meltwater volume without representing its spatial distribution,” says Golden.

Other parameters that describe ferromagnetic behaviour in the original Ising model also have real-world analogues in the version adapted for sea-ice. The global magnetic field, for example, which conventionally governs how likely particle spins are to align as up or down, now corresponds to solar energy input, making each lattice point more or less likely to be water or ice. The strength of coupling between neighbouring particle spins, meanwhile, now describes heat flow between water and ice in adjacent pixels.

Although Golden and colleagues ran their model with zero global field and an infinite coupling strength, changing these parameters after the initial process could perturb a realistic pond arrangement from its metastable state into an alternative low-energy configuration. In this way, the researchers might simulate how sea ice evolves as melt ponds respond to changing environmental conditions.

Golden and colleagues reported their findings in New Journal of Physics.

Marine ice sheet instability amplifies and skews uncertainty in projections of future sea-level rise

Alexander A. RobelHélène Seroussi, and Gerard H. Roe
  1. Edited by Isabel J. Nias, Goddard Space Flight Center, Greenbelt, MD, and accepted by Editorial Board Member Jean Jouzel June 11, 2019 (received for review March 20, 2019)


The potential for collapse of the Antarctic ice sheet remains the largest single source of uncertainty in projections of future sea-level rise. This uncertainty comes from an imperfect understanding of ice sheet processes and the internal variability of climate forcing of ice sheets. Using a mathematical technique from statistical physics and large ensembles of state-of-the-art ice sheet simulations, we show that collapse of ice sheets widens the range of possible scenarios for future sea-level rise. We also find that the collapse of marine ice sheets makes worst-case scenarios of rapid sea-level rise more likely in future projections.


Sea-level rise may accelerate significantly if marine ice sheets become unstable. If such instability occurs, there would be considerable uncertainty in future sea-level rise projections due to imperfectly modeled ice sheet processes and unpredictable climate variability. In this study, we use mathematical and computational approaches to identify the ice sheet processes that drive uncertainty in sea-level projections. Using stochastic perturbation theory from statistical physics as a tool, we show mathematically that the marine ice sheet instability greatly amplifies and skews uncertainty in sea-level projections with worst-case scenarios of rapid sea-level rise being more likely than best-case scenarios of slower sea-level rise. We also perform large ensemble simulations with a state-of-the-art ice sheet model of Thwaites Glacier, a marine-terminating glacier in West Antarctica that is thought to be unstable. These ensemble simulations indicate that the uncertainty solely related to internal climate variability can be a large fraction of the total ice loss expected from Thwaites Glacier. We conclude that internal climate variability alone can be responsible for significant uncertainty in projections of sea-level rise and that large ensembles are a necessary tool for quantifying the upper bounds of this uncertainty.

In marine ice sheets, the grounding line is a critical boundary where ice flowing from the ice sheet interior becomes thin enough to float in ocean water. The grounding line location and ice flux are sensitive functions of the depth of the surrounding ocean (13). When the bedrock beneath the grounding line is reverse sloping (i.e., deepens toward the ice sheet interior), a small retreat of the grounding line onto deeper bed leads to greater ice flux and therefore more retreat. This positive flux feedback leads to the potential for rapid and irreversible retreat wherever the bed is reverse sloping, which has been termed the “marine ice sheet instability” (4). Rapid ice loss from the Antarctic ice sheet through this instability will likely drive sea-level rise beyond the next century (56). However, projections of future sea-level rise are uncertain due to imperfect representation of ice sheet processes in models, unknown future anthropogenic emissions, and the internal variability of future climate forcing of ice sheets. Even with improvements in ice sheet models and climate projections, there will always remain a component of sea-level projection uncertainty that cannot be reduced due to the fundamentally unpredictable internal variability of the climate system which causes ice sheet change. This fundamental lower bound in the uncertainty of projections due to internal climate variability (for ice sheets and other elements of the climate system) has been termed “irreducible uncertainty” (7). The inevitability of uncertainty remaining in sea-level projections necessitates robust modeling of the ice sheet dynamical factors which produce uncertainty in future ice sheet projections, for which there is no existing theoretical framework. In particular, it is critically important to constrain the upper bounds of uncertainty in sea-level rise projections, which have a disproportionate influence on planning for coastal adaptation measures (8).

Ice sheets evolve in response to changes in climate variables (i.e., climate “forcing”), such as snowfall and ocean temperature. The amount and structure of uncertainty in projections of the future ice sheet contribution to sea-level rise can thus be determined through large ensembles of ice sheet model simulations, which are plausible realizations of the future evolution of an ice sheet in response to climate forcing (see Fig. 1Upper for a conceptual illustration of an ensemble). Each ensemble member is distinguished by selecting either a unique set of model parameters or one realization of future variable climate forcing from a distribution of possibilities. At a particular point in time, the statistical properties of the full ensemble represent the probability distribution of ice sheet state, represented by a variable such as extent or volume (conditioned on the probability of the model parameters or climate variability having particular values). The spread (or standard deviation [SD]) of the probability distribution quantifies the amount of uncertainty in the projection. In a probability distribution that is symmetric, the skewness is 0, and the projected change in ice sheet volume is equally likely to fall above or below the mostly likely projection (Fig. 1Lower Left). A negative skewness of the probability distribution indicates that the probability of ice sheet volume turning out to be below the most likely projection is greater than that of ice sheet volume turning out to be above the most likely projection (and vice versa for positive skewness). Put another way, a negative skewness indicates that the probability of worst-case scenarios (i.e., more ice sheet mass loss and corresponding sea-level rise than is expected from the highest-likelihood projection) is greater than the probability of best-case scenarios (of less sea-level rise than the highest-likelihood projection).

Fig. 1.

(Upper) Conceptual illustration of an evolving ensemble of ice sheet model simulations, where each solid line is a single plausible realization of the future ice sheet evolution. Dashed colored lines indicate two time slices, for which probability distributions are provided (Lower Left and Lower Right). The statistical properties of the probability distribution of ensemble simulations change with time from a narrow, symmetric distribution at early times (blue dashed line, Lower Left) to a wide, asymmetric distribution with negative skewness (red dashed line, Lower Right), making the probability of ice sheet volumes below the most likely projection many times greater than that of ice sheet volumes above the most likely projection (arrows show example of ice sheet volumes above and below the most likely ice sheet volume projection).

In this study, we develop a framework for determining how marine ice sheet processes produce uncertainty in projections of future sea level. We do so using 2 complementary approaches: 1) Stochastic perturbation analysis of a simple model of marine ice sheet evolution with one evolving quantity and 2) statistical analysis of a large number of simulations of the future evolution of a West Antarctic glacier using a state-of-the-art numerical ice sheet model with many thousands of evolving quantities. These 2 approaches represent end members of the hierarchy of modern ice sheet models. They provide both a theoretical framework for understanding the sources of uncertainty in projections of the future ice sheet contribution to sea level and an application of this framework to an actual glacier that is thought to be undergoing the marine ice sheet instability.

Stochastic Perturbation Theory

Over the past several decades, much of the observed increase in Antarctic ice loss has been caused by ocean-driven ice sheet melting (9). Our goal in this study is thus to quantify the uncertainty in projected ice sheet state that is caused by uncertainty in ocean-induced ice sheet melt, although the principles involved are extendable to uncertainty in other climate forcing (or uncertainty in glaciological parameters). We start by mathematically and numerically analyzing the uncertainty in ice sheet state simulated by a minimal model of grounding-line migration for a glacier under climate forcing (1011)


[1]where L is the distance of the grounding line from the glacier onset, and hg=ρwρibghg=−ρwρibg is the thickness of ice at the grounding line, which depends on the ratio of seawater (ρwρw) and ice (ρiρi) densities and the local bedrock depth (bgbg). This model, which is derived and discussed in more detail in Robel et al. (11), tracks the total mass balance of the glacier, captured in 3 processes. The first term captures ice entering as snowfall over the glacier surface at accumulation rate P (averaged through glacier geometry: L and hghg). The second term captures ice leaving the glacier due to ice flow through the grounding line (γhβ1gγhgβ−1). The third term captures ocean-induced melt at the grounding line (η). The parameters β and γ are related to the balance of glaciological processes which contribute to setting the velocity of ice at the grounding line (e.g., gravitational driving stress, basal sliding, and ice shelf buttressing). Observations (12), mathematical analyses (1313), and numerical simulations (111415) have all shown that grounding-line velocity is generally a nonlinear function of grounding-line ice thickness (as we assume in Eq. 1), even during periods of transient grounding-line migration, although there is perhaps some variation in the exact value of β that applies in such a situation. Still, this minimal model is meant as a tool to understand the processes which drive uncertainty in simulations of marine ice sheet instability, not a means for making actual predictions of ice sheet change. As we show later on, the conclusions drawn from this minimal model are reproduced in a state-of-the-art ice sheet model which does not make the same simplifying assumptions.

In Eq. 1, the rate at which the grounding line migrates in response to ocean-induced melting (or freezing) is η=η¯+η(t)η=η¯+η′(t), consisting of a time-averaged component (η¯η¯) and a time-variable component (η(t)η′(t)). The time-averaged ocean forcing may be uncertain and so is drawn from a Gaussian distribution with SD σPσP. The time-variable ocean forcing is a first-order autoregressive Gaussian noise process with interannual SD σFσF, and decorrelation timescale τFτF. Other studies have shown using a complex spatially resolved model that ocean-induced grounding-line variability is filtered through the frequency-dependent response of the ice shelf to ocean forcing (16). In this minimal model, we assume a simpler form for η in which all ocean-induced grounding-line migration is the result of melting directly at the grounding line and neglects effects from sub-ice shelf melt and buttressing. While this assumption will likely produce some quantitative difference than that of a model which includes sub-ice shelf melt beyond the grounding line, the more complicated ice-shelf–resolving simulations in the next section show that the qualitative aspects of our mathematical analysis do not appear to be changed by such detailed considerations.

Fig. 2 shows ensembles (with 10,000 simulations each) of simulated grounding-line migration over a bed of constant slope, calculated using Eq. 1 (details in Materials and Methods). The only forcing during the simulation period shown is ocean-induced grounding-line migration (η). The ensemble spread (represented in Fig. 2 by the interquartile range) captures the uncertainty in projected grounding-line position due to uncertainty in ocean forcing. On a forward-sloping bed (green shading/lines), small interannual forcing in the grounding-line position (σF=100σF=100 m/y, η¯=0η¯=0) produces an ensemble spread that remains small, bounded, and symmetric. For the same stochastic forcing on a reverse-sloping bed (orange shading/lines), all ensemble members retreat, as the ensemble spread grows rapidly without bound and becomes skewed (negatively) toward more retreat.

Evolution of a 10,000-member ensemble of minimal model (Eq. 1) simulations of grounding-line retreat over idealized constant-slope bed topography and variability in ocean forcing. (A) Ensemble interquartile range (shading spans 25th percentile to 75th percentile). (B) SD derived from numerically calculated ensemble statistics (solid lines) and analytic predictions from stochastic perturbation theory (circles). (C) Skewness derived from numerically calculated ensemble statistics (solid lines) and analytic predictions from stochastic perturbation theory (circles). Negative skew indicates more retreated grounding line. There is no analytic approximation available for skewness under autocorrelated forcing (see SI Appendix for discussion of stochastic perturbation theory). Pink shading and lines are simulations on a reverse-sloping bed (bx=4×10−3) with constant ocean forcing, selected from a Gaussian distribution. Blue shading and lines are simulations on a reverse-sloping bed (bx=4×10−3) with interdecadal variability in ocean forcing (τF=10 y). Orange shading and lines are simulations on a reverse-sloping bed (bx=3.5×10−3) with interannual variability in ocean forcing (τF=1 y). Green shading and lines are simulations on a forward-sloping bed (bx=−10−3) with interannual variability in ocean forcing (τF=1 y). In all simulations, P=0.35 m/y, γ=7.8×10−9 m−2.75⋅y−1, and β=4.75. In simulations with variable ocean forcing (blue, orange, green), η¯=0 and σF=100 m/y. In simulations with uncertainty in constant ocean forcing (pink), σP=50 m/y.

” data-icon-position=”” data-hide-link-title=”0″>Fig. 2.

Fig. 2.

Evolution of a 10,000-member ensemble of minimal model (Eq. 1) simulations of grounding-line retreat over idealized constant-slope bed topography and variability in ocean forcing. (A) Ensemble interquartile range (shading spans 25th percentile to 75th percentile). (B) SD derived from numerically calculated ensemble statistics (solid lines) and analytic predictions from stochastic perturbation theory (circles). (C) Skewness derived from numerically calculated ensemble statistics (solid lines) and analytic predictions from stochastic perturbation theory (circles). Negative skew indicates more retreated grounding line. There is no analytic approximation available for skewness under autocorrelated forcing (see SI Appendix for discussion of stochastic perturbation theory). Pink shading and lines are simulations on a reverse-sloping bed (bx=4×103bx=4×10−3) with constant ocean forcing, selected from a Gaussian distribution. Blue shading and lines are simulations on a reverse-sloping bed (bx=4×103bx=4×10−3) with interdecadal variability in ocean forcing (τF=10τF=10 y). Orange shading and lines are simulations on a reverse-sloping bed (bx=3.5×103bx=3.5×10−3) with interannual variability in ocean forcing (τF=1τF=1 y). Green shading and lines are simulations on a forward-sloping bed (bx=103bx=−10−3) with interannual variability in ocean forcing (τF=1τF=1 y). In all simulations, P=0.35P=0.35 m/y, γ=7.8×109γ=7.8×10−9 m−2.75⋅y−1, and β=4.75β=4.75. In simulations with variable ocean forcing (blue, orange, green), η¯=0η¯=0 and σF=100σF=100 m/y. In simulations with uncertainty in constant ocean forcing (pink), σP=50σP=50 m/y.

The growth in ensemble spread (i.e., uncertainty) occurs because the marine ice sheet instability amplifies (rather than damps) small perturbations from stochastic ocean forcing. These growing perturbations accumulate over time, leading to a divergence between ensemble members, which each experience a different series of perturbations from ocean forcing. The skewness of the ensemble can be understood physically and from an analysis of Eq. 1. For a retreating ice sheet, the rate of retreat is set by the difference between two fluxes: The accumulation flux and the grounding-line flux. If the grounding-line flux is more sensitive to the position of the grounding line than the accumulation (which it is for sufficiently nonlinear grounding-line flux), the net effect will be that the rate of retreat of more-retreated ensemble members will be greater than the rate of retreat of the less-retreated ensemble members. The result is that the retreat of the ensemble becomes progressively more negatively skewed with time. In ensembles where all simulations are advancing, the most advanced ensemble members accelerate faster, producing a positive skew.

Observations indicate that climate forcing of glaciers in Antarctica (and elsewhere) exhibits strong variability on decadal timescales (1718). In our minimal model, when stochastic forcing has decadal persistence, the ensemble spread and skewness grow considerably faster (blue shading/lines in Fig. 2). Such an amplified glacier response to temporal persistence in forcing agrees with previous model studies of mountain glaciers (19) and periodically forced ice streams (162021).

For temporal ocean variability, stochastic perturbation theory (22) (SI Appendix) provides a theoretical framework for determining the physical processes which control the amplification and skewing of uncertainty in ice sheet projections. Analytic approximations for the spread and skewness of ensembles derived from stochastic perturbation theory (circles in Fig. 2 B and C) match well with numerically calculated ensemble statistics. This theoretical framework shows that ensemble spread grows exponentially with a rate that is proportional to the bed slope and the nonlinearity in grounding-line ice flux (β in Eq. 1). Thus, when the bed is forward sloping (negative), the ensemble spread remains bounded. When the bed is reverse sloping (positive), the marine ice sheet instability causes the ensemble spread to grow exponentially without bound. The ensemble variance is also proportional to the decorrelation timescale of the forcing, implying greater uncertainty in projections when climate forcing is persistent on longer timescales.

As alluded to above, the skewness of the ensemble is caused by the changing rate of grounding-line migration over a reverse-sloping bed. It can be shown analytically (see SI Appendix for details) that when the grounding-line ice flux is sufficiently nonlinear with respect to ice thickness (β>3β>3 in Eq. 1), then ensembles of a retreating grounding line will tend to be skewed toward more retreat. Conversely, when grounding-line ice flux is linear or weakly nonlinear (β<3β<3), then ensembles will tend to be skewed toward less retreat. In a wide range of realistic settings, we expect ensembles to skew toward more retreat during the retreating phase of the marine ice sheet instability because of the high-degree nonlinearity of grounding-line ice flux (β=4.75β=4.75 in ref. 1β=5β=5 in ref. 2β=4β=4 in ref. 3). In other words, the fact that the probability distribution is skewed in the direction of more sea-level rise is a fundamental consequence of the strong nonlinearity inherent in grounding-line dynamics.

Uncertainty in the time-averaged ocean forcing (η¯η¯; pink shading/lines in Fig. 2) produces even more ensemble spread (for relatively less uncertainty, σP=50σP=50 m/y) and further indicates the importance of the marine ice sheet instability for amplifying and skewing uncertainty in ice sheet projections. Uncertainty in the time-averaged climate forcing can be thought of as a limiting case of the response to an initial impulse with an infinite decorrelation time (τFτF→∞). However, for such a nonstochastic case, there is no formal limit from stochastic perturbation theory in which the ensemble spread can be predicted.

Large Ensembles of Thwaites Glacier Instability.

To demonstrate that the intuition gained from the theoretical framework developed in the previous section applies to realistic glacier models, we simulate ensembles of the future retreat of Thwaites Glacier in West Antarctica using the Ice Sheet System Model (ISSM), a state-of-the-art finite-element model of ice sheet flow (23). Thwaites Glacier rests on a reverse-sloping bed and is currently retreating rapidly, which is argued to be the result of the marine ice sheet instability (2425). In ISSM, as in other models, ocean-induced ice sheet melting is parameterized with a depth-dependent melt rate, with maximum melt rate, MmaxMmax, prescribed at some depth (25). We treat MmaxMmax as a first-order autoregressive noise process that varies monthly with a prescribed decorrelation timescale (τFτF), a mean of 80 m/y, no long-term trend, and no variation in space. Many observations and models indicate that subice shelf melt rates at glaciers in the Amundsen Sea, and elsewhere in Antarctica, exhibit strong variability on decadal (and longer) timescales (172627). A short run of a regional ocean model simulation for the Amundsen Sea region (SI Appendix) produces variability in MmaxMmax with interannual SD of 1.4 m/y. This estimate of variability is likely an underestimate since the ocean simulation was run for only 15 y (the time period over which reanalysis forcing is available) and did not include coupled ocean–atmosphere feedbacks. Thus, for our baseline ensemble of Thwaites Glacier, our conservative estimate for the statistics of ocean-induced melt variability is τF=10τF=10 y and σM=1.4σM=1.4 m/y. In reality, we also expect that MmaxMmax varies in space, due to (for example) the Coriolis effect on ocean circulation in the subshelf cavity, which would quantitatively (but not necessarily qualitatively) affect our results.

Fig. 3A shows the evolution of the probability distribution of a 500-member ensemble of ISSM-simulated ice volume at Thwaites Glacier in response to decadal variability in subice shelf melt rate. All ensemble members initialized with the modern state of Thwaites Glacier eventually reach complete deglaciation (Fig. 3B), in agreement with previous studies (2425). The rate of grounding-line migration (Fig. 3C) experiences significant variability over the course of the retreat due to the stochastic ocean forcing and the presence of forward-sloping “speed bumps” in the bed topography, both of which can slow the rate of retreat or even cause advance for short durations (28). Even with the relatively conservative assumption that there is no variability in surface mass balance and a small amplitude of subshelf melt variability (representing <2%<2% of the time-averaged subshelf melt), the spread in the ensemble spans ∼20 cm of uncertainty in projected sea-level rise during periods of fast retreat (i.e., the green probability distribution functions [PDFs] in Fig. 3A), with a probability distribution skewed in the direction of lower ice volume (greater contribution to sea level). This uncertainty is over 25% of the entire sea-level rise due to deglaciation of Thwaites Glacier and 50% of the median sea-level rise achieved during those periods of fast retreat. This spread between simulations amounts to instantaneous differences of hundreds of kilometers in the grounding-line position (Fig. 3D). Following this growth of uncertainty during the centuries of most rapid retreat, the ensemble then contracts and skews in the opposite direction as individual ensemble members achieve complete deglaciation of Thwaites Glacier, due to the limited model domain used in our simulations. In simulations of the entire Antarctic Ice Sheet in which the marine ice sheet instability spreads to other glaciers (56), we would expect even faster amplification of uncertainty as multiple glaciers become involved in deglaciation.

Fig. 3.

Evolution of a 500-member ensemble of ISSM simulations of Thwaites Glacier evolution over 500 y (where year 0 in model time is the modern glacier state) in response to decadal variability and constant average in maximum subice shelf melt rate. (A) Evolution of ensemble PDF over time, plotted every 25 y, with probability on the y axis and Thwaites Glacier ice volume (in cm sea-level equivalent [SLE]) on the x axis. (B) Black lines are simulated ice volume contained in Thwaites Glacier catchment in cm SLE for all ensemble members. (C) Black dots are evolving grounding-line migration rates for all ensemble members (based on the centroid of the 2D grounding line). (D) Snapshots (red, orange, and pink lines) of grounding-line positions at year 635 in model time, from 5th percentile, 50th percentile, and 95th percentile ice volume ensemble members.

In Fig. 4, we compare Thwaites Glacier ensemble statistics, given a comparable amplitude (σM=1.4σM=1.4 m/y) of variability in MmaxMmax, but differing degrees of temporal persistence. As predicted by theory, ensemble spread (Fig. 4A) increases with longer persistence in forcing variability (proportional to τF−−√τFSI Appendix). During the century of fastest retreat, multidecadal ocean variability (yellow line; τF=30τF=30 y) produces skewed uncertainty that is nearly 50% of the total ice loss from Thwaites glacier or ∼40 cm of uncertainty in projected sea-level rise. Some studies have suggested that Antarctic glaciers may be subject to such multidecadal variability in forcing through low-frequency coupled modes of the ocean–atmosphere system (26) or sporadic detachment of very large tabular icebergs (29).

Fig. 4.

Evolution of uncertainty and skewness of uncertainty in four Thwaites Glacier ensembles (500 simulations each). Three of the ensembles have variability in ocean forcing specified using a first-order autoregressive model: Including variability at interannual (τF=1.1τF=1.1 y, blue line), interdecadal (τF=10τF=10 y, red line; Fig. 3), and multidecadal (τF=30τF=30 y, yellow line) timescales. One ensemble has no temporal variability, but the constant maximum subshelf melt rate is uncertain and so drawn from a Gaussian distribution (with σM=5σM=5m/y and the same mean as other ensembles, purple line). (A) “Fractional projection uncertainty” given by the ratio of ensemble spread (measured by ±2σ±2σ of ensemble) to total ice loss at the end of simulation: 4σvol/μvloss4σvol/μvloss. (B) Ensemble skewness, with negative skewness representing a distribution skewed toward lower total ice volume (more ice loss and more sea-level rise).

Poorly constrained subice shelf properties [such as roughness (30)] and the small scale of the turbulent ice–ocean boundary layer make it difficult to accurately simulate even the time-averaged subice shelf melt rate given some change in global climate. Consequently, we also consider uncertainty in the time-averaged subshelf basal melt rate (which may also result from uncertainties in future anthropogenic emissions) by keeping MmaxMmax constant in time, but varying it between ensemble members (drawing from a Gaussian distribution with SD of 5 m/y). As in the minimal model, this ensemble (purple line in Fig. 4) has a very strong amplification of skewed uncertainty due to the accumulation of differences in subshelf basal melt rate among ensemble members over the course of the instability. For several centuries, the spread in this ensemble amounts to nearly the entire signal of ice loss from Thwaites Glacier (i.e., some ensemble members have retreated completely while others have lost almost no ice at all), skewed in the direction of more ice loss throughout most of the early period of the simulation.

Discussion and Conclusions

Studies of the future evolution of the Antarctic Ice Sheet have estimated the uncertainty in future sea-level rise due to poorly constrained model parameters (563134). Other studies (35) have investigated the role of internal climate variability in the Greenland Ice Sheet contribution to sea-level rise, but with simulations that were too short to capture much of the marine ice sheet instability that may occur in the future. No study has provided a theoretical framework explaining the role of ice sheet dynamics in setting the amount and structure of uncertainty in sea-level rise projections. We provide such a theoretical framework in this study and find that ice sheet instabilities are amplifiers of uncertainty, which is a common mathematical property of unstable nonlinear systems (22). Although there are processes not considered here that might stabilize (36) or further destabilize (6) an ice sheet, our analysis shows that we should expect more rapid instability (of any kind) to cause more rapid uncertainty growth. Indeed, the theoretical framework developed in this study applies to a sufficiently broad set of assumptions regarding ice sheet dynamics, such that we expect that any type of ice sheet instability, regardless of the processes involved, will experience rapid growth in the ice sheet projection uncertainty during periods of most rapid instability. Integrating the contribution of marine ice sheet instability over many glaciers also integrates the uncertainty of each glacier’s future evolution, potentially leading to considerable uncertainty in sea-level projections, as has been seen in ensemble studies of total ice sheet contribution to future sea-level rise (56323435). We have shown that model ensembles can be used to quantify a range of possible scenarios for future sea-level rise, including potentially catastrophic scenarios of rapid sea-level rise. However, large model ensembles can be prohibitively expensive when extended to the entire Antarctic Ice Sheet. To fully capture the complete range of possible Antarctic futures, we will need efficient methods for uncertainty quantification (3237) and model order reduction that captures the complexities of ice sheet dynamics (3134). Such sophisticated methods will ensure that we can make the most useful sea-level projections beyond 2100 for those stakeholders who depend on them.

‘World’s most dangerous glacier’ could cause catastrophic sea level rise, study warns

A glacier in West Antarctica, known as “the world’s most dangerous,” could completely melt away and cause a rapid and “catastrophic” sea-level rise, a new study warns.

The study, published in the scientific journal PNAS, notes that the Thwaites Glacier is at a proverbial “tipping point” that could cause a neverending flow of ice into the world’s oceans.

“If you trigger this instability, you don’t need to continue to force the ice sheet by cranking up temperatures. It will keep going by itself, and that’s the worry,” said the study’s lead author and Georgia Tech professor Alex Robel, in a statement. “Climate variations will still be important after that tipping point because they will determine how fast the ice will move.”

The Thwaites Glacier is seen above.

The Thwaites Glacier is seen above. (NASA/OIB/Jeremy Harbeck)


NASA JPL scientist Helene Seroussi, who worked on the study along with Robel, said that the glacier could lose all of its ice over the next 150 years. “That would make for a sea level rise of about half a meter (1.64 feet),” Seroussi added in the statement.

According to the National Oceanic and Atmospheric Administration, sea levels “continue to rise at a rate of about one-eighth of an inch per year.”

The Thwaites Glacier is “the largest single source of uncertainty in projections of future sea-level rise,” according to the study’s abstract. If it were to collapse, it would make “worst-case scenarios of rapid sea-level rise more likely in future projections,” the abstract added.

“If you trigger this instability, you don’t need to continue to force the ice sheet by cranking up temperatures. It will keep going by itself, and that’s the worry.”

— Georgia Tech professor Alex Robel

The Antarctic ice sheet has more than 50 times the amount of ice than the mountain glaciers in the world combined, and eight times as much ice in the Greenland ice sheet, Robel added in the statement.

If and when the glacier becomes unstable, the after-effects would be considered “catastrophic.”

Thwaites Glacier acts like a giant cork that holds back the West Antarctic Ice Sheet.

Thwaites Glacier acts like a giant cork that holds back the West Antarctic Ice Sheet. (NASA/James Yungel)

“Once [the] ice is past the grounding line and only over water, it’s contributing to sea level because buoyancy is holding it up more than it was before,” Robel said. “Ice flows out into the floating ice shelf and melts or breaks off as icebergs.”

The grounding line is the line between where the ice sheet rests on the seafloor and where it extends over the water.

As is common for most sea-level studies, the time scale for the study was in centuries and simulations showed that the ice loss for the Thwaites Glacier started after 600 years. However, the researchers warn that if ocean temperatures continue to rise, the instability in the glacier could occur much faster than many expect.


“It could happen in the next 200 to 600 years. It depends on the bedrock topography under the ice, and we don’t know it in great detail yet,” Seroussi added.

Earlier this year, NASA scientists found a massive hole two-thirds the size of Manhattan under the Thwaites Glacier, a fact that left them disturbed.

The huge cavity – which is approximately 1,000 feet tall, about as tall as New York City’s Chrysler Building – is growing at the bottom of the glacier and is large enough to have once contained 14 billion tons of ice, according to NASA. Most of that ice has melted over the past three years.

Researchers say world’s second-largest emperor penguin colony has been wiped out

Researchers say what was once the world’s second-largest colony of emperor penguins has “now all but disappeared” after changes in sea-ice conditions made their typical breeding grounds highly unstable.

A group of researchers from the British Antarctic Survey (BAS) published their findings in the Antarctic Science journal on Thursday. The team said in a statement that they studied “very high resolution satellite imagery to reveal the unusual findings.”

According to their research, satellite imagery showed that the emperor penguin colony at Halley Bay in Antarctica had drastically decreased over the past three years on account of breeding failures caused by severe changes in local environmental conditions.

“For the last 60 years the sea-ice conditions in the Halley Bay site have been stable and reliable,” the team said. “But in 2016, after a period of abnormally stormy weather, the sea-ice broke up in October, well before any emperor chicks would have fledged.”

The group said the conditions were repeated the following two years, leading to “the death of almost all the chicks at the site each season.”

“The colony at Halley Bay colony has now all but disappeared, whilst the nearby Dawson Lambton colony has markedly increased in size, indicating that many of the adult emperors have moved there, seeking better breeding grounds as environmental conditions have changed,” the researchers said.

Peter Fretwell, the lead author of the report and a remote sensing specialist at BAS, said the team has been studying the population of penguins at the Halley Bay colony and other nearby colonies for years using the high resolution satellite imagery.

“These images have clearly shown the catastrophic breeding failure at this site over the last three years,” Fretwell said. “Our specialized satellite image analysis can detect individuals and penguin huddles, so we can estimate the population based on the known density of the groups to give reliable estimate of colony size.”

Phil Trathan, a penguin expert with BAS who co-authored the report, said “it is impossible to say whether the changes in sea-ice conditions at Halley Bay are specifically related to climate change, but such a complete failure to breed successfully is unprecedented at this site.”

“Even taking into account levels of ecological uncertainty, published models suggest that emperor penguins numbers are set to fall dramatically,” he said, adding that the penguins are likely to lose between 50 percent and 70 percent “of their numbers before the end of this century as sea-ice conditions change as a result of climate change.”

The researchers said they plan to continue to study the colony’s response to the changing sea-ice conditions to help other scientists gain “vital information about how this iconic species might cope with future environmental change.”

Arctic breaks records with 5-year heat streak

Arctic warming continues at about double the rate of the rest of the planet, NOAA says

A polar bear climbs out of the water to walk on the ice in the Franklin Strait in the Canadian Arctic Archipelago in 2017. A new report says Arctic warming continues at about double the rate of the rest of the planet. (David Goldman/Associated Press)

The Arctic had its second-hottest year on record in 2018, part of a warming trend that may be dramatically changing earth’s weather patterns, according to a report released on Tuesday by the U.S. National Oceanographic and Atmospheric Administration.

“Arctic air temperatures for the past five years have exceeded all previous records since 1900,” according to the annual NOAA study, the 2018 Arctic Report Card, which said the year was second only to 2016 in overall warmth in the region.

It marks the latest in a series of warnings about climate change from U.S. government bodies, even as President Donald Trump has voiced skepticism about the phenomenon and has pushed a pro-fossil fuels agenda.



NOAA issues 2018 , region experienced 2nd warmest air temps ever recorded. More details: 

215 people are talking about this

The study said the Arctic warming continues at about double the rate of the rest of the planet, and that the trend appears to be altering the shape and strength of the jet stream air current that influences weather in the Northern Hemisphere.

“Growing atmospheric warmth in the Arctic results in a sluggish and unusually wavy jet-stream that coincided with abnormal weather events,” it said, noting that the changing patterns have often brought unusually frigid temperatures to areas south of the Arctic Circle.

Some examples are “a swarm of severe winter storms in the eastern United States in 2018, and the extreme cold outbreak in Europe in March 2018 known as ‘the Beast from the East.'”

Environmentalists have long warned of rapid warming in the Arctic, saying it threatens imperilled species like polar bears, and is a harbinger of the broader impacts of climate change on the planet.

Trillions in infrastructure damage

Scientists have warned that the region could suffer trillions of dollars worth of climate change-related damage to infrastructure in the coming decades.

But the melting of Arctic ice has piqued the interests of polar nations like the United States, Canada and Russia by opening new shipping routes and expanding access to a region believed to be rich in petroleum and minerals.

The United States and Russia have both expressed an interest in boosting Arctic drilling, and Russia has bolstered its military presence in the north.

The 2018 Arctic Report Card found the Arctic region had the second-lowest overall sea-ice coverage on record. The map shows the age of sea ice in the Arctic ice pack in March 1985 (left) and March 2018 (right). Ice that is less than a year old is darkest blue. Ice that has survived at least 4 full years is white. (Mark Tschudi/University of Colorado/NOAA/

The NOAA report comes weeks after more than a dozen U.S. government agencies released a study concluding that climate change is driven by human consumption of fossil fuels and will cost the U.S. economy hundreds of billions of dollars by the end of the century.

Trump, who has been rolling back Obama-era environmental and climate protections to maximize production of domestic fossil fuels, said of the update to the National Climate Assessment: “I don’t believe it.”

Trump last year announced his intention to withdraw the United States from the 2015 Paris Deal agreed by nearly 200 nations to combat climate change, arguing the accord would kill jobs and provide little tangible environmental benefit.

Antarctic ice melts to January record low


The amount of sea ice around Antarctica has melted to a record low for January, scientists announced this week.

As of January 1, there was 2.11 million square miles of sea ice around the continent, the smallest January area since records began in 1978, according to the National Snow and Ice Data Center.

Sea ice is frozen ocean water that melts each summer, then refreezes each winter. Antarctic sea ice is typically at its smallest in late February or early March, toward the end of summer in the Southern Hemisphere.

“Antarctic sea ice extent is astonishingly low this year, not just near the Ross Ice Shelf, but around most of the continent,” Cecilia Bitz, a polar scientist at the University of Washington, told Grist on Thursday.

Specifically, the area of sea ice around Antarctica on Jan. 1 was 11,600 square miles below the previous record low for that date, set in 2017. It was 726,000 square miles below average – an area roughly twice the size of the state of Texas.

With six to eight weeks remaining in the melt season, the ice center said it remains to be seen whether the Antarctic sea ice will set its all-time record-low minimum.

More: Ancient Antarctic ice sheet collapse could happen again, triggering a new global flood

More: Greenland’s ice sheet melt has ‘gone into overdrive’ and is now ‘off the charts’

Sea ice loss – especially in the Arctic and less so in the Antarctic – is one of the clearest signals of global warming, the National Climate Assessment reported last year.

In addition to human-caused warming of the Earth’s atmosphere and oceans, multiple factors – including the geography of Antarctica, the region’s winds, as well as air and ocean temperatures – affect the ice around Antarctica.

“Although it is too soon for us to isolate what caused the rapid December decline and recent record low extents, it is likely that unusual atmospheric conditions and high sea surface temperatures – important factors in the 2016-2017 record lows – are playing a role,” according to a statement from the National Snow and Ice Data Center.

The center added that even though sea ice occurs primarily in the polar regions, it influences our global climate and weather patterns around the world.