Analysis and prediction of changes in the temperature of the pure freshwater ice column in the Antarctic and the Arctic

A.A. Fedotov,V.V. Kaniber,P.V. Khrapov

DOI: https://doi.org/10.48550/arXiv.2106.08985

2021-06-17

Abstract:This paper investigates the initial boundary value problem for a non-stationary one-dimensional heat equation that simulates the temperature distribution in freshwater ice near the Earth's poles. The mathematical model has been constructed taking into account solid-liquid phase transitions. Data from meteorological stations were used to determine the model parameters, with the help of which the necessary physical and thermophysical characteristics of the computational domain were obtained. For the numerical solution of the problem, the finite volume method (FVM) was used. For the obtained periodic regime, the temperature versus depth dependences for each month were plotted, and the depth of the active layer, as well as the depth of zero annual amplitudes were found for each meteorological station. A forecast of the ice temperature regime for 2100 was modeled for three Representative Concentration Pathway (RCP) scenarios of global warming: moderate RCP2.6, corresponding to the current emissions of RCP7 and adopted at the Paris Agreement in 2015 RCP1.9. By analyzing the impact of an additional 0,5 degrees Celsius of warming on other areas, a reduction in the full range of risks to humanity and the planet as a whole becomes evident with the proper efforts of the global community. Thus, the conducted modeling has confirmed the need to reduce the rate of global warming.

Geophysics,Atmospheric and Oceanic Physics

What problem does this paper attempt to address?

The problem that this paper attempts to solve is the analysis and prediction of the temperature changes of pure freshwater ice columns in the Antarctic and Arctic regions in the context of global warming. Specifically, the research simulates the temperature distribution of freshwater ice near the Earth's poles by constructing a mathematical model of an unsteady - state one - dimensional heat conduction equation considering solid - liquid phase change. The research uses meteorological station data to determine model parameters and conducts numerical solutions through the Finite Volume Method (FVM) to analyze the changes in the ice layer temperature field and determine the time required to reach an unsteady - state periodic state. In addition, based on different global warming scenarios (RCP2.6, RCP7 and RCP1.9), the paper also predicts the ice temperature regime in 2100 to evaluate the impact of different climate paths on the ice layer temperature.

### Research Background
1. **Impact of Global Warming**: Global warming causes glaciers and ice sheets to melt, which in turn leads to a series of environmental problems such as sea - level rise and the ocean absorbing more heat.
2. **Importance of the Polar Regions**: The polar regions are the main areas on the Earth where ice and cold air accumulate, so they are particularly sensitive to climate change responses, and special attention needs to be paid to the temperature changes of their ice layers.

### Mathematical Model
The research constructs an unsteady - state one - dimensional heat conduction equation describing the temperature distribution of the ice layer:
\[ \rho c \frac{\partial u}{\partial t} - \frac{\partial}{\partial z} \left( \lambda \frac{\partial u}{\partial z} \right) + Q \delta(u - u^*) = 0 \]
where:
- \( \rho \): density
- \( c \): specific heat capacity
- \( \lambda \): thermal conductivity coefficient
- \( u(z, t) \): medium temperature
- \( u^* \): phase - change temperature
- \( Q \): latent heat of phase change
- \( \delta(u - u^*) \): Dirac delta function

### Boundary Conditions
1. **Upper Boundary**: At \( z = 0 \), convective heat transfer occurs, and the boundary condition is:
\[ J = h (\theta(t) - u(0, t)) \]
where:
- \( J \): heat flux density at the boundary
- \( h \): convective heat transfer coefficient
- \( \theta(t) \): atmospheric temperature

2. **Lower Boundary**: At \( z = L \), there is a no - heat - flow condition:
\[ J_b = 0 \]

### Physical and Geographical Conditions
The research uses long - term monthly average air temperature, wind speed and snow depth data from multiple meteorological stations to determine the thermal physical property parameters within the calculation domain. These meteorological stations include the Amundsen - Scott Station and the Esperanza Base in Antarctica, and the Kap Morris Jesup Station and the Ernst Krenkel Observatory in the Arctic.

### Prediction Results
1. **RCP2.6 Scenario**: Global warming is controlled below 2°C, and the ice layer temperature rises, but does not reach the melting point.
2. **RCP7 Scenario**: Global warming reaches 4°C, and the ice layer temperature rises significantly, and an active layer (that is, the ice layer begins to melt) appears in some areas.
3. **RCP1.9 Scenario**: Global warming is controlled below 1.5°C, the ice layer temperature changes slightly, and the risk of melting is greatly reduced.

### Conclusion
The research results show that even under the mild RCP2.6 scenario, the increase in ice layer temperature will lead to a significant reduction in the ice sheet area. And the more severe RCP7 scenario will bring more adverse consequences. Therefore, controlling the speed and magnitude of global warming is crucial for slowing down ice layer melting.