MPE Mathematical Problems in Engineering 1563-5147 1024-123X Hindawi Publishing Corporation 10.1155/2014/540379 540379 Research Article Parameter Identification for Salinity in a Quasilinear Thermodynamic System of Sea Ice http://orcid.org/0000-0003-2824-5563 Lv Wei 1 Li Xiaojiao 1 Feng Enmin 2 Peng Haipeng 1 Department of Mathematics Shanghai University Shanghai 200444 China shu.edu.cn 2 School of Mathematical Sciences Dalian University of Technology Dalian 116024 China dutmba.com 2014 1072014 2014 19 03 2014 20 06 2014 14 7 2014 2014 Copyright © 2014 Wei Lv et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

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 identification model 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.

1. Introduction

Sea ice accounts (in areal extent) for nearly two-thirds of the earth’s ice cover . It is both an important actor and an indicator of climate change . 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 . Salinity may also affect large-scale sea ice and ocean characteristics . 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 . Malmgren  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 . Eicken’s ice-core based classification of ice salinity profiles features typical profiles that differ significantly from their Arctic counterparts . 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 . In recent years, there are many researchers devoted to the identification problems. In , 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 . References [16, 17] provided the mathematical theory to calculate the ice densities, specific heats, and thermal conductivities. In , 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 , 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  and THESCI’s . 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  and THESCI’s . Final conclusions follow in Section 5.

2. 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 z 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 l1 meter to l2 meter depth, and 0<l1<l2. Let t denote time, tf(>0) the final time, and the unit is in seconds. Let T(z,t) be the temperature at depth z and time t, and the unit is in Kelvins. Let Ω=(l1,l2)R, I=(0,tf), Q=Ω×I. For convenience, set In={1,2,,n}, and n 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: (1)C(z,t,T)Tt=z(k(z,t,T)Tz),inQ,T(z,0)=T0(z),inΩ¯,T(l1,t)=TU(t),inI-,T(l2,t)=TL(t),inI-, where k(z,t,T) is the thermal conductivity of sea ice and C(z,t,T) is the heat capacity of sea ice, which are expressed by the following formulae on the basis of previous work . 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: (2)C(z,t,T)=ρpicpi+γS(z,t)(T-273.15)2,inQ,k(z,t,T)=kpi+βS(z,t)T-273.15,inQ,S(z,t)=S(z^)=q1+q2z^+q3z^2+q4z^3,inQ,z^=zh(t),inQ, where ρpi, cpi, and kpi denote the density, specific heat, and thermal conductivity of pure ice, respectively; γ and β are known positive constants; z^ is the normalized- depth; h(t) is the ice thickness; S(z,t) is the mean salinity of sea ice; and q1, q2, q3, and q4 are the unknown parameters describing the ice salinity, which will be estimated.

Assumption 1.

There exist qLi>0, qUi>0, and qLiqUi, for any (z,t)Q¯, such that qLiqiqUi, iI4.

Assumption 2.

Consider TL(t)C1(I-;R+), TU(t)C1(I-;R+), T0(z)C1(Ω¯;R+), h(t)C1(I-;R+).

Remark 3.

From Assumptions 1 and 2, we can get S(z,t)C1(Q¯;R+). And for any (z,t)Q¯, there exist hM>0, SMZ>0, SMT>0, SL>0, SU>0, SLSU, such that |h|hM, |St|SMT, |Sz|SMZ, and SLS(z,t)SU.

From Remark 3, we make the following assumption.

Assumption 4.

Consider βSU<273.15kpi.

Remark 5.

From the expression of C(z,t,T), we get that C(z,t,T)>0, which squares with the facts of science.

Remark 6.

Because k(z,t,T) is the thermal conductivity of sea ice, we can get that k(z,t,T)>0. In fact, if k(z,t,T)=0, there is no heat transmitted through sea ice; if k(z,t,T)<0, which shows that heat is transmitted from low temperature to high temperature. Therefore k(z,t,T)0 is not in accordance with the laws of science.

From Remarks 5 and 6, we make the following assumption.

Assumption 7.

Take Ψ(z,0)=T0(z), Ψ(l1,t)=TU(t), Ψ(l2,t)=TL(t), Ψ(z,t)H2+α,1+α/2(Q¯), and satisfy (3)[Ψt(z,0)-a11(z,0,Ψ(z,0),Ψz(z,0))|z=l1,l2×Ψzz(z,0)+a(z,0,Ψ(z,0),Ψz(z,0))]|z=l1,l2=0, where (4)a11(z,0,Ψ(z,0),Ψz(z,0))=k(T)C(T),a(z,0,Ψ(z,0),Ψz(z,0))=-kzTzC(T).

Remark 8.

From Assumption 1 and Remark 3, we define (5)Qad={i=14q=(q1,q2,q3,q4)0<qLjqjqUj,jI4,SLi=14qiz^i-1SU}; then Qad is a compact convex subset of R4.

From the above analysis, we can get the following results.

Theorem 9.

Suppose Assumptions 17 hold; then there exists a unique solution T(z,t)H2+α,1+α/2(Q¯) satisfying QTS.

Proof.

Take T/t=Tt, T/z=Tz, 2T/z2=Tzz, Ψ(z,0)=T0(z), Ψ(l1,t)=TU(t), Ψ(l2,t)=TL(t); then QTS can be written as (6)Tt-k(z,t,T)C(z,t,T)Tzz-kz(z,t,T)C(z,t,T)Tz=0,inQ,(7)Ψ(z,0)=T0(z),inΩ¯,(8)Ψ(l1,t)=TU(t),inI-,(9)Ψ(l2,t)=TL(t),inI-, where (10)k(z,t,T)C(z,t,T)=kpi(T-273.15)2+βS(T-273.15)ρpicpi(T-273.15)2+γS,inQ,kz(z,t,T)C(z,t,T)=β(Sz(T-273.15)-STz)ρpicpi(T-273.15)2+γS,inQ,Sz(z,t)=q2h(t)+2q3zh2(t)+3q4z2h3(t),inQ. Let (11)a11(z,t,T,Tz)=k(z,t,T)C(z,t,T)=kpi(T-273.15)2+βS(T-273.15)ρpicpi(T-273.15)2+γS,a(z,t,T,Tz)=-kz(z,t,T)TzC(z,t,T)=-βTz(Sz(T-273.15)-STz)ρpicpi(T-273.15)2+γS.

The expression of a11(z,t,T,Tz) is none of Tz, so a11(z,t,T,Tz) is written as a11(z,t,T). Then (6) can be written as (12)Tt-a11(z,t,T)Tzz+a(z,t,T,Tz)=0,inQ. From Remarks 5 and 6, for (z,t)Q¯,T, we have (13)a11(z,t,T)=k(T)C(T)>0. From the expression of a(z,t,T,Tz), for b10, we have (14)Ta(z,t,T,0)=0-b1T2. For (z,t)Q¯,|T|M=273.15-βS/kpi, pTz, we can get that a11(z,t,T) and a(z,t,T,p) are continuous.

And (15)a11(z,t,T)p=0,a11(z,t,T)T=-βSρpicpi(T-273.15)2-2kpiSγ(T-273.15)-βγS2[ρpicpi(T-273.15)2+γS]2,a11(z,t,T)z=(T-273.15)2Sz[βρpicpi(T-273.15)-γkpi][ρpicpi(T-273.15)2+γS]2. It is obvious that a11/p, a11/T, and a11/z are continuous, so a11(z,t,T) is differential with respect to z, T, and p for (z,t)Q¯, |T|M,p.

Take (16)ν=kpi(M-273.15)2+βSL|M-273.15|ρpicpi(M+273.15)2+γSU>0,μ=kpi(M+273.15)2+βSU(M+273.15)ρpicpi(M-273.15)2+γSL>0; then for (z,t)Q¯,|T|M,p, we have that νa11(z,t,T(z,t))μ.

For (z,t)Q¯,|T|M,p, we have (17)|a11/z|1+|p|=|ρpicpikpi(T-273.15)4-ρpicpiβ(Sz-S)(T-273.15)3|(1+|p|)[ρpicpi(T-273.15)2+γS]2+γkpiSz(T-273.15)2(1+|p|)[ρpicpi(T-273.15)2+γS]2ρpicpikpi(M+273.15)4+ρpicpiβ(SMZ+SU)(M+273.15)3[ρpicpi(M-273.15)2+γSL]2M1,|a11T|=|βSρpicpi(T-273.15)2-2kpiSγ(T-273.15)-βγS2|[ρpicpi(T-273.15)2+γS]2βSUρpicpi(M+273.15)2+2kpiSUγ(M+273.15)+βγSU2[ρpicpi(M-273.15)2+γSL]2M2,|a|(1+|p|)2=|βSp2-βSz(T-273.15)p|(1+|p|)2[ρpicpi(T-273.15)2+γS]βmax{SU,SMZ(M+273.15)}ρpicpi(M-273.15)2+γSLM3. So (18)|a11p|(1+|p|)3+|a11T|(1+|p|)2+|a11z|(1+|p|)+|a(z,t,T,p)|=(|a11/z|1+|p|+|a11T|+|a(z,t,T,p)|(1+|p|)2)(1+|p|)2(M1+M2+M3)(1+|p|)2.

For any (z,t)Q¯, we can get (19)|a11t|=|St(T-273.15)2[βρpicpi(T-273.15)-γkpi][ρpicpi(T-273.15)2+γS]2|SMT(M+273.15)2[βρpicpi(M+273.15)+γkpi][ρpicpi(M-273.15)2+γSL]2. So a11/t is bound. And a11/z and a11/T are also bound; we can get that a11(z,t,T) is Lipschitz continuous on z,t,T.

For (z,t),(z,t)Q¯; T,T, |T|M,|T|M; p,p, and |p|Mp,|p|Mp, set Λ=(l1,l2)×(0,tf)×(-M,M), and we have that (20)a11(z,t,T)Hα,α/2(Λ¯)=sup(z,t,T),(z,t,T)Λ¯,|z-z|2l2-2l1|a11(z,t,T)-a11(z,t,T)||z-z|α+sup(z,t,T),(z,t,T)Λ¯,|t-t|2tf|a11(z,t,T)-a11(z,t,T)||t-t|α/2+sup(z,t,T),(z,t,T)Λ¯,|T-T|2M|a11(z,t,T)-a11(z,t,T)||T-T|α+maxΛ¯|a11(z,t,T)|=sup(z,t,T),(z,t,T)Λ¯,|z-z|2l2-2l1|a11(z,t,T)-a11(z,t,T)||z-z|·|z-z|1-α+sup(z,t,T),(z,t,T)Λ¯,|t-t|2tf|a11(z,t,T)-a11(z,t,T)||t-t|·|t-t|1-α/2+sup(z,t,T),(z,t,T)Λ¯,|T-T|2M|a11(z,t,T)-a11(z,t,T)||T-T|·|T-T|1-α+maxΛ¯|a11(z,t,T)|<. Therefore, a11(z,t,T) is Hölder continuous in the variable t with the exponent α/2 and in z, T, and p with α.

Obviously, we can also get that a(z,t,T,p) is Hölder continuous in the variable t with the exponent α/2 and in the variables z, T, and p with the exponent α.

From the above analysis, Assumption 7 and Theorem 5.2 , we can get that there exists a unique solution T(z,t)H2+α,1+α/2(Q¯) satisfying QTS, which completes our proof.

3. 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: (21)minJ(q)=T(z,t;q)-Tmea(z,t)TmaxL2(Q)2+S(z^;q)-Smea(z^)SmaxL2[0,1]2s.t.T(z,t;q)H2+α,1+α/2(Q¯),qQad, where Tmea(z,t) and Smea(z,t) are the field data of the temperature and salinity and Tmax and Smax 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 J(q) 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.

3.1. Existence of Optimal Parameter

The next theorem presents the existence of the identified parameters for the identification model.

Theorem 10.

Suppose Assumptions 17 hold; then there exists one optimal solution q*Qad satisfying PIM.

Proof.

For any parameter qQad, let {qn}Qad be a sequence, such that qn-qQad0 as n. Let Tn=T(qn) and T=T(q) be the solutions of QTS. Then we have that (22)Tn-TH2+α,1+α/2(Q¯)=sup(z,t),(z,t)Q¯,|z-z|2l2-2l1|(z,t)Dz2(Tn-T)(z,t)aaaaaaaaaaaaaaaaaxaaaa-Dz2(Tn-T)(z,t)|aaaaaaaaaaaaaaaaaaa×(|z-z|α)-1+sup(z,t,T),(z,t,T)Q¯,|z-z|2l2-2l1|(z,t)Dt(Tn-T)(z,t)aaaaaaaaaaaaaaaaaxaaaaaa-Dt(Tn-T)(z,t)|aaaaaaaaaaaaaaaaaaa×(|z-z|α)-1+sup(z,t,T),(z,t,T)Q¯,|t-t|2tf|(z,t)Dz(Tn-T)(z,t)aaaaaaaaaaaaaxaaaaaaaa-Dz(Tn-T)(z,t)|aaaaaaaaaaaaaaaaaaa×(|t-t|(1+α)/2)-1+sup(z,t,T),(z,t,T)Q¯,|t-t|2tf|(z,t)Dt(Tn-T)(z,t)aaaaaaaaaaaaxaaaaaaaaa-Dt(Tn-T)(z,t)|aaaaaaaaaaaaaaaaaaa×(|t-t|α/2)-1+sup(z,t,T),(z,t,T)Q¯,|t-t|2tf|(z,t)Dz2(Tn-T)(z,t)aaaaaaaaaaaaaxaaaaaaaa-Dz2(Tn-T)(z,t)|aaaaaaaaaaaaaaaaaaa×(|t-t|α/2)-1+maxQ¯|Tn(z,t)-T(z,t)|+maxQ¯|Dz(Tn-T)(z,t)|+maxQ¯|Dz2(Tn-T)(z,t)|+maxQ¯|Dt(Tn-T)(z,t)|0(n).

So we can get that the mapping qT(q) is strongly continuous in H2+α,1+α/2(Q¯).

Next we get the existence of optimal solution.

For any qQad, set X1(T;q)=T(z,t;q)-Tmea(z,t), X2(z^;q)=S(z^;q)-Smea(z^); then (23)J(q)=X1(T;q)Tmea(z,t)L2(Q¯)2+X2(z^;q)Smea(z,t)L2[0,1]20. Obviously, X1(T;q) and X2(z^;q) are continuous on Qad; thus J(q) is continuous on Qad. Since Qad is a nonempty compact convex set, and the mapping qT(q) is strongly continuous in H2+α,1+α/2(Q¯), there exists q*Qad, for all qQad, such that J(q*)J(q); that is q*Qad is the optimal solution of PI.

Next section will give the algorithm to solve the identified model.

3.2. 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 T in formula J(q) 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 PS for the identified vector q and not give some value. The initial value for q will be generated randomly in the domain. The stop criterion is a minor error ε for J(q) or a large iterative time IT. If the value of J(q) is less than ε or the iterative time is greater than IT, the program will stop. Next are the algorithm steps.

Step 1.

Give PS, ε, IT, and generate q0PS randomly; set l=0.

Step 2.

Solve QTS for T(ql) by SIFDS, and calculate J(ql) by GA and HJA.

Step 3.

If J(ql)ε or lIT, set q*=ql, and output q*; else, go to Step 2 with l=l+1.

4. 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.

The salinity and temperature errors are defined as (24)SE(q)=190(s=190(S(z^s;q)-Smea(z^s))2),TE1(q)=1i-·j-j=1j-i=1i-(T(zi,tj;q)-Tmea(zi,tj))2,TE2(q)=j=1j-i=1i-|T(zi,tj;q)-Tmea(zi,tj)|j=1j-i=1i-|Tmea(zi,tj)|,TEZ(z;q)=1j-(j=1j-(T(zi,tj;q)-Tmea(zi,tj))2),aaaaaaaaaaaaasaaaaaaaaaaai=1,,i-,TET(t;q)=1i-(i=1i-(T(zi,tj;q)-Tmea(zi,tj))2),aaaaaaaaaaaasaaaaaaaaaaaaj=1,,j-, where i- and j- 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 . 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.

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

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 . 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 l1=0.36 m and l2=0.66 m. In the simulated temperature profiles, we take l1=0.12 m and l2=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.15 m thick vertical sections. Sea ice salinity was measured using a portable conductivity meter DDBT-350 with an accuracy of ±0.5 ppt . 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 . 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 h(t)=0.7822+8.0×10-8t. 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.

List of model parameters and constants.

Parameters Symbols Value Unit
Density of pure ice ρ pi 916 Kg m−3
Specific heat of pure ice c pi 2093 J kg−1 K−1
Thermal conductivity of pure ice k pi 2.03 W m−1 K−1
γ 1.715 × 10 7 J Km−3 ppt−1
β 0.1172 W m−1 ppt−1
Observed minimum salinity S L 2.2 Ppt
Observed maximum salinity S U 13 ppt

4.2. 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.

Salinity parameter values and salinity errors of three methods.

q 1 q 2 q 3 q 4 SE
PIM 2.1667 32.5360 −81.5137 53.6941 1.3671
Eicken 1.89 29.19 −77.49 56.38 1.7808
THESCI 3.9998 19.6095 −58.6298 48.0169 2.8175

Temperature errors for PIM and Eicken.

TE1 TE2 TE1 TE2
(6/21–7/5) (6/21–7/5) (5/27–6/20) (5/27–6/20)
PIM 0.1055 0.0317% 0.0986 0.0267%
Eicken 0.1082 0.0329% 0.1053 0.0284%

The mean salinity profile of sea ice.

The ice salinity profiles on April 25, May 21, June 18, July 27, August 29, October 2, November 1, November 18, and November 30.

The absolute salinity profiles on April 25, May 21, June 18, July 27, August 29, October 2, November 1, November 18, and November 30.

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

The temperature error configurations from June 21 to July 5, 2006, for PIM and Eicken.

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

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

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 . 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.

5. 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.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

This work was supported by the National Natural Science Foundation of China under Grant nos. 11101262, 11171050, and U1232112, Dalian University of Technology Special Fund under Grant no. DUTTX2011106, Key Disciplines of Shanghai Municipality under Grant no. S30104, and a grant of “The First-Class Discipline of Universities in Shanghai.” The authors are grateful to all the 2006 winter-over crews at Zhongshan Station in the 22nd CHINARE for their field work support.

Maykut G. A. Untersteiner N. Some results from a time- dependent thermodynamic model of sea ice Journal of Geophysical Research 1971 76 6 1550 1575 10.1029/JC076i006p01550 2-s2.0-0015019401 Serreze M. C. Holland M. M. Stroeve J. Perspectives on the Arctic's shrinking sea-ice cover Science 2007 315 5818 1533 1536 10.1126/science.1139426 2-s2.0-33947383369 Carmack E. C. Untersteiner N. Nijhoff M. Circulation and mixing in ice-covered waters The Geophysics of Sea Ice 1986 146 Dordrecht, The Netherland Springer 641 712 Nato ASI Series B Vancoppenolle M. Fichefet T. Goosse H. Simulating the mass balance and salinity of Arctic and Antarctic sea ice. 2. Importance of sea ice salinity variations Ocean Modelling 2009 27 1-2 54 69 10.1016/j.ocemod.2008.11.003 2-s2.0-59249083593 Hunke E. C. Notz D. Turner A. K. Vancoppenolle M. The multiphase physics of sea ice: a review for model developers The Cryosphere 2011 5 4 989 1009 10.5194/tc-5-989-2011 2-s2.0-81255144693 Malmgren F. On the properties of sea ice, The Norwegian North Polar Expedition with the “Maud” 1918–1925 John Griega Bokst 1927 1 5 1 67 Schwarzacher W. Pack-ice studies in the Arctic Ocean Journal of Geophysical Research 1954 64 2357 2367 Eicken H. Salinity profiles of Antarctic sea ice: field data and model results Journal of Geophysical Research 1992 97 C10 15545 15557 2-s2.0-0027070224 Vancoppenolle M. Fichefet T. Bitz C. M. Modeling the salinity profile of undeformed Arctic sea ice Geophysical Research Letters 2006 33 21 L21501 10.1029/2006GL028342 2-s2.0-34248513498 Gough A. J. Mahoney A. R. Langhorne P. J. Williams M. J. M. Haskell T. G. Sea ice salinity and structure: a winter time series of salinity and its distribution Journal of Geophysical Research C: Oceans 2012 117 3 C03008 10.1029/2011JC007527 2-s2.0-84863393859 Omatu S. Seinfeld J. H. Distributed Parameter Systems 1989 Oxford, UK Clarendon Press Li L. Yang Y. Peng H. Wang X. Parameters identification of chaotic systems via chaotic ant swarm Chaos, Solitons and Fractals 2006 28 5 1204 1211 10.1016/j.chaos.2005.04.110 ZBL1121.90426 2-s2.0-29144460641 Peng H. Li L. Yang Y. Liu F. Parameter estimation of dynamical systems via a chaotic ant swarm Physical Review E: Statistical, Nonlinear, and Soft Matter Physics 2010 81 1 016207 10.1103/PhysRevE.81.016207 2-s2.0-75349093230 Li L. Yang Y. Peng H. Fuzzy system identification via chaotic ant swarm Chaos, Solitons and Fractals 2009 41 1 401 409 10.1016/j.chaos.2008.01.011 ZBL1198.90419 2-s2.0-67349261317 Ozbilge E. Demir A. Identification of the unknown coefficient in a quasi-linear parabolic equation by a semigroup approach Journal of Inequalities and Applications 2013 2013, article 212 7 10.1186/1029-242X-2013-212 MR3063984 2-s2.0-84879641936 Lv W. Feng E. Li Z. A coupled thermodynamic system of sea ice and its parameter identification Applied Mathematical Modelling. Simulation and Computation for Engineering and Environmental Systems 2008 32 7 1198 1207 10.1016/j.apm.2007.03.006 MR2409376 ZBL1187.80024 2-s2.0-40649086673 Lv W. Feng E. Li Z. Properties and optimality conditions of a three-dimension non-smooth thermodynamic system of sea ice Applied Mathematical Modelling 2009 33 5 2324 2333 10.1016/j.apm.2008.07.002 MR2494289 2-s2.0-58249136518 Lv W. Feng E. Lei R. Parameter identification for a nonlinear non-smooth thermodynamic system of sea ice International Journal of Thermal Sciences 2009 48 1 195 203 10.1016/j.ijthermalsci.2008.03.003 2-s2.0-54849436877 Ebert E. E. Curry J. A. An intermediate one-dimensional thermodynamic sea ice model for investigating ice-atmosphere interactions Journal of Geophysical Research 1993 98 6 10085 10109 2-s2.0-0027881535 Ladyzenskaja O. A. Solonnikov V. A. Ural'ceva N. N. Linear and Quasilinear Equations of Parabolic Type 1968 Providence, RI, USA American Mathematical Society MR0241822 Lei R. Li Z. Cheng B. Zhang Z. Heil P. Annual cycle of landfast sea ice in Prydz Bay, east Antarctica Journal of Geophysical Research C: Oceans 2010 115 2 C02006 10.1029/2008JC005223 2-s2.0-77249165592