Parameter Identification for Salinity in a Quasilinear Thermodynamic System of Sea Ice

This study is intended to provide a parameter identification method to determine salinity of sea ice by temperature and salinity observations. A quasilinear thermodynamic system of sea ice with unknown salinity is described and its property is proved. Then, a parameter identificationmodel is established and the existence of its optimal solution is discussed.The salinity profile is calculated by the temperature and salinity data, which were measured at Nella Fjord around Zhongshan Station, Antarctica, during the polar night time by the 22nd Chinese Antarctic Research Expedition. Another simulation for temperature profiles during different measurement periods is operated. Results show that better simulations of the salinity and temperature distribution are possible with the estimated parameters than Eicken’s (Eicken 1992) and THESCI’s (Lv et al. 2009).This method will help people understand the salinity evolution of sea ice more thoroughly.


Introduction
Sea ice accounts (in areal extent) for nearly two-thirds of the earth's ice cover [1].It is both an important actor and an indicator of climate change [2].Its salinity is the main difference from water ice.The input of brine from sea ice is important in the ocean interaction, through modification of the density of the surface layer and the induction of thermohaline convection [3].Salinity may also affect largescale sea ice and ocean characteristics [4].Therefore, the physical properties of sea ice are thus strongly dependent on salinity.
Salinity has long fascinated sea ice researchers for the complexity that it produces through depression of freezing and melting temperatures [5].Malmgren [6] discussed the evolution of the vertical salinity profile and recognized its effect on the thermal properties of sea ice.In evaluation of salinity data, most researchers resort to studying mean values integrated over the thickness of an entire floe.Most of this work has concentrated on the Arctic and Antarctic [7][8][9][10].Eicken's ice-core based classification of ice salinity profiles features typical profiles that differ significantly from their Arctic counterparts [8].Although the seasons were not evenly sampled, Eicken's classification provides a general idea of the salinity profile in Antarctic sea ice.
The above researchers' methods to determine the ice salinity are according to salinity observations.The salinity measurement is by hand, and the field data are spare and unsatisfactory due to the difficulties associated with fieldwork, especially during the polar winter; thus it could not help people to thoroughly understand the physical evolution of sea ice salinity just by salinity data.Parameter identification method is effective in solving this problem.Identification refers to the determination from observed data of unknown parameters in the system model such that the predicted response of the model is close, in some well-defined sense, to the process observations [11].In recent years, there are many researchers devoted to the identification problems.In [12][13][14], parameter identification problems were performed via chaotic ant swarm.A semigroup approach to the mathematical analysis of the inverse coefficient problems of identifying the unknown coefficient in the quasilinear parabolic equation with Dirichlet boundary conditions was presented [15].References [16,17] provided the mathematical theory to calculate 2 Mathematical Problems in Engineering the ice densities, specific heats, and thermal conductivities.In [18], an ice salinity profile was determined by the parameter identification method abbreviated by THESCI.
In this paper, a parameter identification method to determine salinity using salinity and temperature observations is presented.Different to [18], a quasilinear thermodynamic system without a source term is chosen.The existence and uniqueness of the solution for the system are proved.The coefficients describing the salinity of sea ice are taken as the identified parameters, the deviation sum of the temperature and the salinity is defined as the performance criterion, and a parameter identification model is constructed; the existence of the identified parameters is presented.There are many powerful optimal algorithms to solve this problem, such as chaotic particle swarm optimization and chaotic ant swarm.To get the optimal solution, an optimization algorithm is constructed on the basis of Genetic Algorithm, Hooke-Jeeves Algorithm, and semi-implicit difference scheme.The numerical results will be compared with Eicken's [8] and THESCI's [18].The main contribution of this paper is that we not only provide a new method to determine the salinity of sea ice but also present its mathematical theory.
The objective of this study is to present a new method to determine the ice salinity and better understand the salinity profiles.Section 2 describes the quasilinear thermodynamic system of sea ice and its property.A parameter identification model is constructed and the existence for the optimal solution is proved in Section 3. Section 4 performs a numerical experiment and compares it to Eicken's [8] and THESCI's [18].Final conclusions follow in Section 5.

Quasilinear Thermodynamic System of Sea Ice
For the purpose of describing the salinity profile more accurately, the thermodynamic system of sea ice not containing a source term at the polar night is presented.Since the gradient variation in the vertical direction is far greater than the one in the horizontal direction, we only consider the heat flux in the vertical direction.Let  representing the depth of the system be the vertical coordinate taken as positive downward and let the origin be at some point on the ice surface, and the unit is in meters; the area we consider is from  1 meter to  2 meter depth, and 0 <  1 <  2 .Let  denote time,   (> 0) the final time, and the unit is in seconds.Let (, ) be the temperature at depth  and time , and the unit is in Kelvins.Let Ω = ( 1 ,  2 ) ⊂ R,  = (0,   ),  = Ω × .For convenience, set I  = {1, 2, . . ., }, and  is any positive integer.According to Fourier's law of heat conduction, the thermodynamic system is described by the following heat conduction equations denoted by QTS: where (, , ) is the thermal conductivity of sea ice and (, , ) is the heat capacity of sea ice, which are expressed by the following formulae on the basis of previous work [19].For more applications, the ice salinity is taken as mean salinities determined for ice normalized-depth which is the ratio of the ice depth and the ice thickness.So the salinity is calculated by the following cubic polynomial function: where From Remark 3, we make the following assumption.
Remark 5. From the expression of (, , ), we get that (, , ) > 0, which squares with the facts of science.
Remark 6.Because (, , ) is the thermal conductivity of sea ice, we can get that (, , ) > 0. In fact, if (, , ) = 0, there is no heat transmitted through sea ice; if (, , ) < 0, which shows that heat is transmitted from low temperature to high temperature.Therefore (, , ) ≤ 0 is not in accordance with the laws of science.
From Remarks 5 and 6, we make the following assumption.

Parameter Identification Model
For the purpose of making the salinity and temperature of sea ice be close to the actual process observations, a parameter identification model is established as PIM: where  mea (, ) and  mea (, ) are the field data of the temperature and salinity and  max and  max denote the maximums among the field data.Here the temperature unit is taken as Kelvin, and the salinity's ppt.Since units are different, the temperature and salinity in () are normalized.
Our goal is to estimate the salinity coefficients such that the temperature and salinity values computed from QTS most closely reproduce these observations.Next section will give the algorithm to solve the identified model.

Algorithm.
In this section, an algorithm is constructed to solve the above identification model.
To solve the model described in Section 2, the ice temperature  in formula () should be calculated.So the first thing is to calculate the thermodynamic system QTS.Here we choose the semi-implicit finite difference scheme abbreviated to SIFDS, which is an approach with unconditional convergence.To get the optimal solution of the identified model, the Genetic Algorithm is used in conjunction with Hooke-Jeeves Algorithm, where Genetic Algorithm and Hooke-Jeeves Algorithm are abbreviated to GA and HJA, respectively.GA is a global search one, which can find the global minimum point in theory.HJA is a local search one.The combination of GA's good global search and HJA's accurate local search is a good manner to solve the model rapidly and effectively.The convergence is obvious.We will present the four-dimensional closed domain  for the identified vector  and not give some value.The initial value for  will be generated randomly in the domain.The stop criterion is a minor error  for () or a large iterative time .If the value of () is less than  or the iterative time is greater than , the program will stop.Next are the algorithm steps.Step 1. Give , , , and generate  0 ∈  randomly; set  = 0.

An Experiment
In this section, an experiment to estimate the parameters describing the ice salinity is performed.The field data are described in Section 4.1.Numerical results and discussions are performed in Section 4.2.
where  and  are temperature numbers of measured spatial spots and temporal spots.
4.1.Field Data.The field data were measured by the 22nd Chinese Antarctic Research Expedition.The experimental site was located in the coastal area off Zhongshan Station (Figure 1). Figure 1 is seen in [21].In the region, polar night occurs from late May to late July.Five survey sections were set up to monitor landfast ice (Figure 1).As shown in Figure 1, Section 3 was located north of the station, where the field data taken in this paper were measured.Landfast ice off Zhongshan usually breaks off in early summer, often due to coincident high wind speeds, monthly peak tides, and/or the effect of penetrating ocean swell [21].A thermistor string was deployed at northern end of Section 3 on 27 March at site B (Figure 1), where the water depth is 35 m.At the time of installation, the ice thickness was 0.270 m, with 0.015 m snow.The thermistor string was installed in an auger hole, which refroze naturally, and the portion of the hole above water level was filled with water of about 3 ppt.The temperature sensors in the first 1.26 m beneath the upper ice surface were spaced at 0.06 m intervals and at 0.12 m intervals below that.Each sensor extended by 0.02 m from the rod.The coherence of these sensors was controlled by selecting sensors whose reading was within ±0.005 K of each other during calibration.From 31 March through 12 December 2006, ice temperatures were measured at 30 min interval and stored on board a DUT-5000 storage system.The data system was designed to hold 20-day data and the battery pack to provide powder for more than 20 days.The measured time for temperature data taken here is from May 27 to June 15, 2006, which is the polar night time.To get more accurate results, we take  1 = 0.36 m and  2 = 0.66 m.In the simulated temperature profiles, we take  1 = 0.12 m and  2 = 0.78 m.
During 2006 ten full length ice cores were collected using a 0.13 m diameter corer within an area of 10 m by 10 m near the center of Section 3. The cores were cut lengthwise to obtain 0.05-0.15m thick vertical sections.Sea ice salinity was measured using a portable conductivity meter DDBT-350 with an accuracy of ±0.5 ppt [21].Here we take the salinity data of 9 ice cores on April 25, May 21, June 18, July 27, August 29, October 2, November 1, November 18, and November 30.
In March 2006, when ice allowed safe travel only adjacent to the shore, sea ice thickness was measured every 5 days at 3 to 5 holes, drilled with a 0.05 m diameter stainless steel auger.The accuracy of the drill hole method was ±0.005 m [21].From late March to late December 2006 the ice thickness was measured every 2-5 days using the drill hole method at a distance of 1 m around the initial measurement sites.The fitting thickness functions for the period are ℎ() = 0.7822 + 8.0 × 10 −8 .The pure ice density, heat capacity, and thermal conductivity are related to temperature and pressure.The parameters values in QTS are given in Table 1.

Results and Discussions.
Using the observed temperature data from June 21 to July 15 and the salinity observations, we get the results of identified parameters shown in Table 2. Figures 2 and 3 present the salinity distributions.Figure 2 gives the mean salinity distributions at the normalized-depth and Figure 4 shows nine salinity profiles at the ice depth on April 25, May 21, June 18, July 27, August 29, October 2, November 1, November 18, and November 30.The salinity unit is in parts per thousand (ppt), that is, the mass of salt in grams per kilogram of sea ice, where the red curves are salinity profiles for PIM, the blue for Eicken, and the green for THESCI.The black dots are observations.The gray lines in Figure 2 are the ranges for the measurement errors.Figure 4 shows the absolute errors for the ice salinity profiles on the nine days.Figure 5 presents the absolute errors and Figure 6 gives the defined temperature error configurations.To examine the validity, we make another simulation for temperature profiles from May 27 to June 20.The defined salinity and temperature errors are presented in Tables 2 and 3. Figure 7 presents the absolute errors and Figure 8 gives the defined temperature errors for the simulated temperature profiles.
In Figure 2, the mean salinity distribution for PIM is shaped like a question mark with low, high, low, and high salinities as seen from the top.This is uniform with the observations and [21].Compared to the salinity observations in Figure 3, there are 28 points for PIM in the range of measurement errors, 23 for Eicken, and 18 for THESCI.This means that, among the three results, PIM is the best.
In Figure 3, as seen from the top, the shape observed in the first subgraph is shaped like a question mark with low, high, low, and high salinities.But it does not happen for the three numerical results.This is because the three are shaped with high, low, low, and high salinities.Compared with Eicken and THESCI, the top salinity for PIM is much closer to the observation.Except for the first subgraph in Figure 4, vertical salinity distributions for PIM show typical "question-mark-shaped" profiles.
In Figure 4, the smallest deviation to the observations for PIM is the top point on November 1 and the largest is the bottom point on November 18. From Figure 5, we can conclude that 40 absolute errors for PIM are the smallest among the three results, 32 for Eicken, and 18 for THESCI.
Since we take the ice depth 0.36 m and 0.66 m as the upper and lower boundaries for the thermodynamic system QTS in Figure 5 and 0.12 m and 0.78 m in Figure 7, the absolute temperature errors in the first and last subgraphs are zero.The error fluctuations to zero for PIM at the remainder six depths in Figure 5 and at the remainder eight depths in Figure 7 are shown.Most absolute errors for PIM are smaller.
In Figures 6 and 8, most errors of the defined temperature for PIM are smaller than those for Eicken.
From Tables 2 and 3, we can get that the errors for PIM are less than those for Eicken and THESCI.
From these tables and figures, we can conclude that the numerical results using PIM correspond to the actual characteristics of the ice salinity and temperature distributions and approach the field data better than those using Eicken's and THESCI's parameters.The parameter identification method is effective when field data are spare and unsatisfactory due to the difficulties associated with fieldwork, especially during the polar winter.Thus the mean salinity function can be applied and our method is valid, and the method can help in interpreting field data and can be used to overcome data gaps.

Conclusions and Further Work
In this paper, we aim to provide a new method to estimate the salinity of sea ice neglecting the influence of snow and ocean on ice.We have gotten the parameter values describing the salinity and operated another numerical simulation during different measurement periods.Compared with the results of Eicken's and THESCI's, the numerical results of the new method indicate that not only the method is valid, but also our model with the optimal parameters can be used, for example, to overcome data gaps.
However, currently the impact of snow on sea ice, of the oceanic heat flux, and of radiation penetrating through snow on top of the sea ice is not included in our mathematical framework.Additionally, we do not account for the impact of brine moving within sea ice on the temperature field.These limitations will be the subject of the future work.

Figure 1 :
Figure 1: Locations of landfast ice survey areas in the vicinity of Zhongshan Station (solid circle) in 2006.

Figure 2 :
Figure 2: The mean salinity profile of sea ice.

Figure 5 :Figure 6 :
Figure 5: The absolute error configurations of the ice temperature at different depths from June 21 to July 5, 2006, for PIM and Eicken.

Figure 7 :(
Figure 7: The absolute error configurations of the ice temperature at different depths from May 27 to June 20, 2006, for PIM and Eicken.

Figure 8 :
Figure 8: The temperature error configurations from May 27 to June 20, 2006, for PIM and Eicken.

Table 2 :
Salinity parameter values and salinity errors of three methods.