> LNK%` 8bjbj .B̟̟FDDDD`$($&&&&&&$h/JJPPP$P$PPPPTDP0P0PPd<^PLJJF
$ $ USING NEURAL NETWORKS TO SOLVE THE INVERSE PROBLEM OF EXTRACTING SALINITY FROM THE UNESCO EQUATION OF STATE FOR SEA WATER AND TO SPEED UP THE EVALUATION OF THIS EQUATION
Vladimir M. Krasnopolsky, Laurence C. Breaker, Dmitri Chalikov, and Desiraju B. Rao
National Centers for Environmental Prediction, Washington, D.C. 20233
We consider two related problems which arise in oceanic modeling and data assimilation. (1) In most ocean models, the UNESCO International Equation of State for Seawater (UES) [1] is evaluated at each point of a threedimensional grid for each time step. For highresolution models, the solution of this equation consumes a significant part of the overall computation time. (2) In data assimilation systems, the assimilation of temperature into ocean circulation models which employ the full equation of state, without making corresponding adjustments to the salinity can lead to problems in the assimilation process. To perform such an adjustment, we must invert the UES to obtain salinity as a function of temperature, density and depth (or pressure). Unfortunately, it is not a simple matter to extract salinity, given temperature and density, from the UES since this represents what is essentially an inverse problem. Alternately, we could use a simpler and more readily invertible version of the equation of state for seawater such as that given by [2], but then we sacrifice accuracy, particularly for extreme values. Here we propose solutions for both problems using neural networks (NNs)  a technique which wellsuited for nonlinear modeling and inversion.
The UES gives the following expression for the density r (kg/m3) as described in [1],
INCLUDEPICTURE "http://polar.ncep.noaa.gov/marine.meteorology/marine.winds/NNs/eq1.gif" \* MERGEFORMATINET (1)
where r is the density of seawater in kg/m3, T is the temperature in C, S is the salinity in psu, P is the pressure, and K(T,S,P) is the secant bulk modulus [1].
The empirically based UES equation (1) is given over a threedimensional domain for D = { 2 < T < 40C, 0 < S < 40 psu, and 0 < P < 10000 decibars}. This domain represents all possible combinations of T, S, and Z which are encountered globally. Mathematically, the functions (T,S,0) and K(T,S,P) are represented by multidimensional polynomials and, as a result, the density (1) is the ratio of two, threedimensional polynomials which contain more than 40 parameters.
Due to the limitations of the UES described above, we have developed an alternate approach for obtaining simpler and faster local parameterizations for (T,S,Z), using neural networks which are excellent devices for approximating multidimensional functions [3]. The neural network (NN) approach has also allowed us to invert the density to obtain salinity as a function of temperature, density and depth, S(T,r,Z).
We have developed regional NN representations for density and salinity for a threedimensional subdomain, DS, where DS = { 2 < T < 35C, 5 < S < 38 psu, and 0 < Z < 5700 m}. In this threedimensional volume, 4,000 points (Ti, Si, Zi) were generated. The UES was used to estimate the density of seawater, ri, for each point. This simulated data set (ri; Ti, Si, Zi) was used as a proxy for experimental data in order to train the NNs to extract density and salinity. Two expressions were obtained:
r = rNN (T,S,Z) (2a)
S = SNN (T,,Z) (2b)
where rNN and SNN are the NN representations which can be expressed as follows and f = r or S,
INCLUDEPICTURE "http://polar.ncep.noaa.gov/marine.meteorology/marine.winds/NNs/eq2.gif" \* MERGEFORMATINET (3)
where the matrix Wji, the vectors Bj, wj, and b collectively represent the weights and biases of the NN, and a and b are scaling parameters. All of these coefficients are different for than they are for S, and r in each case they are determined during the process of training the NN.
The NN for estimating is 23 times faster than USE. Creating density lookup tables for each model level may make the calculations of the sea water density even faster. However, for inversion of UES for salinity, the NN (2b) provides a fast and robust solution which there is no better alternative to the best of our knowledge. To evaluate the errors in using the NN approach to estimate salinity, we used the UES to estimate the density of seawater, .0NP4D#$01DE18
(*NRnpr~hg6H*]hg6]jhgU hgH*hgOJQJhghg5\RD M .B:1F1Z1n111 $$Ifa$gdggdg$a$gdg8djlptjlprv~^`bdhjtv$&24"*****.*0*6*8*>*@*B*f*h*l***++++++$+&+T+V+X+Z+\+d+f+Uj$ hgUjhgUhgOJQJhg6H*]hg6]hgSri, at 16,000 points (Ti, Si, Zi). Then the NN for SNN (2b) was applied to calculate a new salinity, si, using the corresponding values (Ti, ri, Zi). Then the differences (Si  si) were utilized to estimate the accuracy of the NNderived salinities (first line in Table). To further evaluate the quality of the NNderived salinities, the UES was applied again, this time to the triad (Ti, si, Zi) to recalculate the density of seawater, r'i. If the NNobtained values for salinity were perfect, then the density, r'i, would be equal to ri. The differences between these two values, (ri  r'i), were then used to further estimate the accuracy of the salinitytrained NN to retrieve salinities in terms of the density (second line in Table). Table shows that the NN estimates of salinity (5b) have an RMS error of 0.1 psu. In terms of the related error in density, this accuracy corresponds to an RMS error of 0.08 kg m3.
UnitsMin errorMax errorMean errorRMS errorpsu0.330.850.000.10Kg m30.270.710.000.08Two particular applications of NNs have been considered here: (1) accelerating numerical estimation of a complicated multivariate function by modeling this function using a simple NN, and (2) inversion of a multivariate function using NNs. There are many similar problems in meteorology and oceanography which also may benefit through the application of NNs.
REFERENCES
[1] Fofonoff, N.P., and R.C. Millard Jr., "Algorithms for computation of fundamental properties of seawater", UNESCO technical paper in marine science 44, UNESCO, 1983
[2] Mamayev, O.I., 1975: Temperature  Salinity Analysis of the World Ocean Waters. Elsevier Scientific, Amsterdam
[3] Chen, T., and H. Chen, "Approximation capability to functions of several variables, nonlinear functionals, and operators by radial basis function neural networks," Neural Networks, vol. 6, pp.904910, 1995
jlnp....0.2.4........*1.1:1D1F1X1Z1l1n1111111111228hgCJaJhg5\ hgH*hg6H*]hgOJQJhg'1111111E<<<<< $$Ifa$kdT$$Ifr!0634abytg1111112E<<<<< $$Ifa$kd($$Ifr!0634abytg2244L64788E@@@@@>gdgkd$$Ifr!0634abytg,1h/ =!"#$%$ DdJL
S(Aeq1bT@u՝>FPT`DnXT@u՝>FPTPNG
IHDRbWPLTEƦFخ 000PPP```ppp@ @ @@@@@@``@```@@@@@@@@@@@@ @ @@ @ @ @@@@@@@@@@@@`@`@@`@`@`@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @ @@@@@@``@```@@@@@@ @ @@@@@@``@```@@@@@@ @ @@@@@@``@```@@@@@נ༴ tRNSރYbKGDHgIFg8^DcmPPJCmp0712HsIDATx^۶ {UTHHV2Y]Iyş@T$EE` ,j ZBvVcːfv_o,~?>cn/ROV4'?
oX:*0>bBiՋ3,ZZ;Pxtbir1lhfs"E9Wuc4*;P5 t@ҩsS
LEWK\z/Ui}fSQtTJUk\,JՋF0BiPkX:b 7wRJ.ZKoJX\r՝XL%iC'~*54Xn^Cg;K
m3 XrW~Myz]+X=?KZIz,5ɂxJ;%`hTR)9t0/'Q+r(k9~
Ϝ~PT[7x,`Sp@?#n w^:МQ7=BC9X<==d"ETD2 Tcdq?>e߱kX1kkgPKŮP
¨IENDB`0
Dd L
S(Aeq2b qӞ/ӒIl h nd qӞ/ӒIPNG
IHDRFoҺPLTEƦFخ 000PPP```ppp@ @ @@@@@@``@```@@@@@@@@@@@@ @ @@ @ @ @@@@@@@@@@@@`@`@@`@`@`@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @ @@@@@@``@```@@@@@@ @ @@@@@@``@```@@@@@@ @ @@@@@@``@```@@@@@נ༴ tRNSރYbKGDHgIFg8^DcmPPJCmp0712HsIDATx^\r+9[nGRe)~sx=As{rzOfW00cDTD=T"KMr`]575Sɇ>NC6MBDWLc1j]Ct9rj7
i]zqYFO+)
ꥐ:m ;?"qg{}G%5^M8L^
dz9Ktj:zRpXm?'iy="lo;?+ۣ9G *5%^l Ӈ
ߨ%!ݛqւ%
Ӽ7^/{sԬ{IVl]5`C{W4A{RohmSۃ6l8hYBi{ )dlz`,K'ih1G:c\
Q=4Tk!kY}&֦rl>O2O4m'{8uzg]T)9"ZϢ'@+9X{P8^a!Xe/cjDNL>~)`jȁŤ :g70D=וM4{u@nQf.vH7A<%Y6x^=4+ ?:Wa[kr=$X
K'=7ݚ,X'P<
r9y6{m"F͌^#
wŤ<}DDV'?nKlD:˚'nuS:S9.fA=6X=}%"KB:r)MջApMO+KXx7+fXlygrokpa'us=_bndmW7҈s&&H9KHr;QGy7,[{xVڛm·fɁR~~fs
^xds6,ˀ=nƅƅώIpl%A[Z9d)=zUV3JvW}[C'zIENDB`$$If!vh555555{#v5#v#v#v#v{:V065/34ytg$$If!vh555555{#v5#v#v#v#v{:V065/34ytg$$If!vh555555{#v5#v#v#v#v{:V065/34ytg@@@NormalCJ_HaJmH sH tH DA@DDefault Paragraph FontRi@RTable Normal4
l4a(k@(No ListB^`BgNormal (Web)dd[$\$FBDM@! 8*
stzHTqDH000000000000000000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 000000DM@! 8*
stzHTqDHX0X0X0X0X0X0X0X0X0X0X0X0X0X0X0X0X0X0X0X0X0X0X0X0X0X0X0X0X0X0X0X0X0X0X0@0@0@0@0@0
81128 8sFCC AtGBLC"DQEL?"FG~HId
9>==CffH7=CCBHHooH=*urn:schemasmicrosoftcom:office:smarttags PlaceType=*urn:schemasmicrosoftcom:office:smarttags PlaceName>*urn:schemasmicrosoftcom:office:smarttags
PostalCode9*urn:schemasmicrosoftcom:office:smarttagsState8*urn:schemasmicrosoftcom:office:smarttagsCity9 *urn:schemasmicrosoftcom:office:smarttagsplaceP ~7:NV ij
<>efhj\^`b %&
V]?BX`H!&3
7
/ x

DH333333333333DEHHI29QgtzH@EEpCcEE11111F`@```$@`````@@`*UnknownGz Times New Roman5Symbol3&z Arial"qh]FbF$24::2IX)?g2USING NEURAL NETWORKS TO SOLVE THE INVERSE PROBLEM OF EXTRACTING SALINITY FROM THE UNESCO EQUATION OF STATE FOR SEA WATER AND TO SPEED UP THE EVALUATION OF THIS EQUATION renichols renicholsOh+'0LXlx
USING NEURAL NETWORKS TO SOLVE THE INVERSE PROBLEM OF EXTRACTING SALINITY FROM THE UNESCO EQUATION OF STATE FOR SEA WATER AND TO SPEED UP THE EVALUATION OF THIS EQUATIONrenicholsNormal.dotrenichols1Microsoft Office Word@^в@B_T@T՜.+,0hp
DOC/NOAA/NWS/NCEP:'USING NEURAL NETWORKS TO SOLVE THE INVERSE PROBLEM OF EXTRACTING SALINITY FROM THE UNESCO EQUATION OF STATE FOR SEA WATER AND TO SPEED UP THE EVALUATION OF THIS EQUATIONTitle
!#$%&'()*+,./0123456789:<=>?@ABDEFGHIJMRoot Entry F@TOData
"1TableWordDocument.BSummaryInformation(;DocumentSummaryInformation8CCompObjq
FMicrosoft Office Word Document
MSWordDocWord.Document.89q