A 198.8 m deep borehole was drilled through ice to subglacial bedrock in the northwestern marginal part of Princess Elizabeth Land, ~12 km south of Zhongshan Station, in January–February 2019. Three years later, in February 2022, the borehole temperature profile was measured, and the geothermal heat flow (GHF) was estimated using a 1-D time-dependent energy-balance equation. For a depth corresponding to the base of the ice sheet, the GHF was calculated as 72.6 ± 2.3 mW m−2 and temperature −4.53 ± 0.27°C. The regional averages estimated for this area based, generally, on tectonic setting vary from 55 to 66 mW m−2. A higher GHF is interpreted to originate mostly from the occurrence of metamorphic complexes intruded by heat-producing elements in the subglacial bedrock below the drill site.

Very few studies have emphasized the effects of high-pressure sintering on snow density evolution, even though snow as a type of engineering material is widely used in construction engineering in cold regions for snow pavement, snow runway and polar infrastructure. This study presents new experimental results of snow densification under high pressures of up to 100 MPa for a temperature range from −3.5 to −17.3°C and uniaxial compression at the temperature of −10°C and constant strain rates from 5 × 10−4 to 10−1 s−1. Results reveal that density evolution of snow to ice under high-pressure sintering can be achieved in a wide temperature range within a duration as short as 5 min. The compressive strength of snow-sintered ice was ~1.2–2.2 times as large as that of water-frozen ice reported by previous work. The orthogonal experiment showed that pressure is a more significant factor affecting the final density in comparison with sintering temperature and time. The increased rates of ice fabrication, low limitations on temperature and reliable sintered snow strength indicate that snow-ice engineering, such as airport construction in Greenland and Antarctica, can be improved by high-pressure sintering of snow to overcome the harsh environment.

Hot-point drills have been widely used for drilling boreholes in glaciers, ice caps and ice sheets. A hot-point drill melts ice through the thermal head at its bottom end. Penetration occurs through a close-contact melting (CCM) process, in which the ice is melted, and the meltwater is squeezed out by the exerted force applied on the thermal head. During the drilling, a thin water film is formed to separate the thermal head from the surrounding ice. For the hot-point drill, the rate of penetration (ROP) is influenced by several variables, such as thermal head shape, buoyancy corrected force (BCF), thermal head power (or temperature) and ice temperature. In this study, we developed a model to describe the CCM process, where a constant power or temperature on the working surface of a thermal head is assumed. The model was developed using COMSOL Multiphysics 5.3a software to evaluate the effects of different variables on the CCM process. It was discovered that the effect of thermal head shape and the cone angle of conical thermal head on ROP is less significant, whereas the increase in the BCF and the power (or temperature) of the thermal head can continuously enhance the ROP.

In glaciology, snow–firn temperature at 10 m is considered a representation of the mean annual air temperature at the surface (MAAT) of the studied site. Although MAAT is an important parameter in ice-sheet investigations, it has not been widely measured in Antarctica. To measure the 10 m snow–firn temperature in Antarctica, a shallow hot-point drill system is designed. In this simple and lightweight system, a hot-point drill can melt boreholes with a diameter of 34 mm in the snow–firn to a depth of 30 m and a temperature sensors string can measure the borehole temperature precisely. In the 2018/19 field season, 16 boreholes along the Zhongshan–Dome A traverse were drilled, and the borehole temperature was measured. Although certain problems existed pertaining to the hot-point drill, a total depth of ~244 m was successfully drilled at an average penetration rate of ~10 m h−1. After borehole drilling, ~12–15 h were generally required for the borehole to achieve thermal equilibrium with the surroundings. Preliminary results demonstrated that the 10 m snow–firn temperature along the traverse route was affected by the increasing altitude and latitude, and it decreased gradually with an increase in the distance from Zhongshan station.

A series of new synthetic armored cables were developed and tested to ensure that they were suitable for use with the RECoverable Autonomous Sonde (RECAS), which is a newly designed freezing-in thermal ice probe. The final version of the cable consists of two concentric conductors that can be used as the power and signal lines. Two polyfluoroalkoxy jackets are used for electrical insulation (one for insulation between conductors, and the other for insulation of the outer conductor). The outer insulation layer is coated by polyurethane jacket to seal the connections between the cable and electrical units. The 0.65 mm thick strength member is made from aramid fibers woven together. To hold these aramid fibers in place, a sheathing layer was produced from a polyamide fabric cover net. The outer diameter of the final version of the cable is ~6.1 mm. The permissible bending radius is as low as 17–20 mm. The maximal breaking force under straight tension is ~12.2 kN. The cable weight is only ~0.061 kg m−1. The mechanical and electrical properties and environmental suitability of the cable were determined through laboratory testing and joint testing with the probe.

The Antarctic subglacial drilling rig (ASDR) is designed to recover 105 mm-diameter ice cores up to 1400 m depth and 41.5 mm-diameter bedrock cores up to 2 m in length. In order to ensure safe and convenient drilling, drilling auxiliaries are designed to support fieldwork and servicing. These auxiliaries are subdivided into several systems for power supply, drill tripping in the borehole, ice core and chip processing, and drill servicing and maintenance. The required equipment also includes two generators, a drilling winch with a cable, logging winch with a cable, control desk, pipe handler with a fixed clamp, chip chamber vibrator, centrifuge, emergency devices and fitting and electrical tools. Additionally, several environmental protective measures such as a new liquid-tight casing with a thermal casing shoe and a bailing device for recovering drilling fluid from the borehole were designed. Most of the auxiliaries were tested during the summer of 2018–2019 near Zhongshan Station, East Antarctica while drilling to the bedrock to a depth of 198 m.

A new, modified version of the cable-suspended Ice and Bedrock Electromechanical Drill (IBED) was designed for drilling in firn, ice, debris-rich ice and rock. The upper part of the drill is almost the same for all drill variants and comprises four sections: cable termination, a slip-ring section, an antitorque system and an electronic pressure chamber. The lower part of the IBED comprises an auger core barrel, reamers, a core barrel for ice/debris-ice drilling and a conventional geological single-tube core barrel or custom-made double-tube core barrel. First, the short and full-scale field versions of the IBED were tested at an outdoor testing stand and a testing facility with a 12.5 m-deep ice well. Then, in the 2018–2019 summer season, the IBED was tested in the field at a site ~12 km south of Zhongshan Station, East Antarctica, and a ~6 cm bedrock core was recovered from a 198 m-deep borehole. A total of 18 d was required to penetrate the ice sheet. The retrieved core samples of blue ice, basal ice and bedrock provided valuable information regarding the Earth's paleo-environment.

In many cases, the efficiency and safety of a drilling project depend on the reliability of the electrical and electronic control system, as the process progresses without visual access of the operator. The electrical and electronic system provides and regulates the power supply for the drill, collects and monitors the drill data during the whole operating process, and sends and receives the control instructions and feedback signals. The entire system is composed of the surface, borehole and software subsystems. The surface subsystem serves for operating the drilling process, transmitting the drilling and environmental data, and supplying power for the drill motor and downhole control system. The borehole subsystem is generally intended for borehole data acquisition, drill motor control, power regulation and communication. The software subsystem is designed for human–computer interaction, data processing and storage, and programming of signal acquisition and transmission of data. The control system of Antarctic subglacial drilling rig was tested during the 2018–2019 summer season near Zhongshan Station, East Antarctica, in the course of drilling to the bedrock at a depth of 198 m. It exhibited a steady and efficient performance without significant system failures.

The British Antarctic Survey Rapid Access Isotope Drill is an innovative new class of electromechanical ice drill, which has recently been used to drill the deepest dry hole drilled by an electromechanical auger drill. The record-breaking depth of 461.58 m was drilled in just over 104 hours at Little Dome C. The drill collects ice chippings, for water stable isotope analysis, rather than an ice core. By not collecting a core the winch can be geared for speed rather than core breaking and is lightweight. Furthermore, emptying of the chippings is performed by simply reversing the drill motor on the surface reducing the overall drilling time significantly. The borehole is then available for instrumentation. We describe the drill in its current state including modifications carried out since it was last deployed. Test seasons and the lessons learned from each are outlined. Finally, future developments for this class of drill are discussed.

In this paper, we propose a compact scheme to numerically study the coupled stochastic nonlinear Schrödinger equations. We prove that the compact scheme preserves the discrete stochastic multi-symplectic conservation law, discrete charge conservation law and discrete energy evolution law almost surely. Numerical experiments confirm well the theoretical analysis results. Furthermore, we present a detailed numerical investigation of the optical phenomena based on the compact scheme. By numerical experiments for various amplitudes of noise, we find that the noise accelerates the oscillation of the soliton and leads to the decay of the solution amplitudes with respect to time. In particular, if the noise is relatively strong, the soliton will be totally destroyed. Meanwhile, we observe that the phase shift is sensibly modified by the noise. Moreover, the numerical results present inelastic interaction which is different from the deterministic case.

We study the construction of symplectic Runge-Kutta methods for stochastic Hamiltonian systems (SHS). Three types of systems, SHS with multiplicative noise, special separable Hamiltonians and multiple additive noise, respectively, are considered in this paper. Stochastic Runge-Kutta (SRK) methods for these systems are investigated, and the corresponding conditions for SRK methods to preserve the symplectic property are given. Based on the weak/strong order and symplectic conditions, some effective schemes are derived. In particular, using the algebraic computation, we obtained two classes of high weak order symplectic Runge-Kutta methods for SHS with a single multiplicative noise, and two classes of high strong order symplectic Runge-Kutta methods for SHS with multiple multiplicative and additive noise, respectively. The numerical case studies confirm that the symplectic methods are efficient computational tools for long-term simulations.

In this paper we propose stochastic multi-symplectic conservation law for stochastic Hamiltonian partial differential equations, and develop a stochastic multi-symplectic method for numerically solving a kind of stochastic nonlinear Schrödinger equations. It is shown that the stochastic multi-symplectic method preserves the multi-symplectic structure, the discrete charge conservation law, and deduces the recurrence relation of the discrete energy. Numerical experiments are performed to verify the good behaviors of the stochastic multi-symplectic method in cases of both solitary wave and collision.

The local one-dimensional multisymplectic scheme (LOD-MS) is developed for the three-dimensional (3D) Gross-Pitaevskii (GP) equation in Bose-Einstein condensates. The idea is originated from the advantages of multisymplectic integrators and from the cheap computational cost of the local one-dimensional (LOD) method. The 3D GP equation is split into three linear LOD Schrödinger equations and an exactly solvable nonlinear Hamiltonian ODE. The three linear LOD Schrödinger equations are multisymplectic which can be approximated by multisymplectic integrator (MI). The conservative properties of the proposed scheme are investigated. It is mass-preserving. Surprisingly, the scheme preserves the discrete local energy conservation laws and global energy conservation law if the wave function is variable separable. This is impossible for conventional MIs in nonlinear Hamiltonian context. The numerical results show that the LOD-MS can simulate the original problems very well. They are consistent with the numerical analysis.