7#,66(t dx ~*81.:h~~~~~~Salinity Intrusion in the Red River Delta Vu Thanh Ca Department of Civil and Environmental Engineering Saitama University, Urawa, Saitama 338 Japan Abstract Characteristics of salinity intrusion into estuaries of the Red River System, Vietnam are studied based on many year measurement data of salinity concentration at stations along the estuaries and on a numerical model. Computed results of salinity concentration at measurement stations along rivers and tributaries by the numerical model agree satisfactorily with observed data. It was found that in the dry season, the salinity intrusion length is as long as 20 km in the main river and more than 20 km for some tributaries. In the main river and tributaries with high freshwater discharge, the maximum salinity concentration is observed in January while for the tributaries with low freshwater discharge, the maximum salinity concentration is observed in March. 1 Introduction The Lower Red River Delta (Fig. 1) is the most important part of the Bacbo Delta. the main agricultural production supplier for the North of Vietnam. There are two distinguished seasons in this area: the rainy season from May to October and the dry season from November to April of the next year. The difference in freshwater discharge between the two seasons is very large. The largest monthly average freshwater discharge in August is about ten times larger than the smallest one in March (Vu, 1992). In the rainy season, due to large freshwater discharge from upstream, the salinity intrusion problem does not present. However, the dry season. With small freshwater discharge, the salinity intrusion into river estuaries has been a main environmental problem. With relatively large slope of the river bed in this area. the salinity intrusion problem for this river system is not as serious as that in the Mekong River Delta in the South of Vietnam. However, the Lower Red River Delta is the most densely populated area in Vietnam, and consequently salinity intrusion can cause large damage to local economy. Vi (1985) and Ngo (1991) based on the recorded data of salinity concentration at stations along estuaries of the Red River System has drawn some primary remarks on the characteristics of salinity intrusion there. Details of salinity intrusion in each tributary of the river network were not investigated. The aim of this study is to give a detailed view on the characteristics of salinity intrusion in the estuaries of the Red River System. This report is based mainly on two previous publications by Vu (1990) and Vu et al (1991). Using many year recorded data of salinity concentration at stations along the estuaries, monthly-average salinity concentration at each station is computed. The salinity intrusion length in each estuary are also estimated. Details of salinity concentration distributions along the estuaries are studied using a numerical model of the transport and dispersion of salinity. 2 Regime of Salinity Intrusion in the Estuaries of the System Fig. 1 depicts the estuaries of the Red River System together With locations of salinity measurement stations. The data of the salinity concentration at 10 stations have been recorded in the dry season for the period of 33 years (from 1963 up to now). However, in this study, the salinity concentration data for the period of 29 years (from 1963 to 1992) have been used for the evaluation of the regime of salinity intrusion into the estuaries of the Red River System since data for the last four years are not available.  Figs. 2(a-d) depict the monthly average salinity concentrations at the stations along the estuaries of the Red River System in the dry season. In the Day river (Fig. 2a) and the Red River, the main river of the river system, (Fig. 2c), due to relatively large slope of the river bed, and consequently large freshwater discharge, the decrease of salinity concentration with the length from the estuaries is faster than that in the Ninhco and Traly Rivers (2b and 2d). The maximum salinity concentration at stations along the Red River and Day Rivers is observed in January while that in the Ninhco and Traly Branches is observed in March.  In this study, the salinity intrusion length is defined as the length from the river mouth along the river channel to a point where the salinity concentration is l ppt (part per thousand). Fig. 3 depicts a map of contour lines of many-year-average salinity concentration and many--year-average maximum salinity concentration of 4 ppt and 1 ppt. This map is drawn based on the data collected in some special projects (Vi, 1985) and the recorded data at stations. As seen in this figure, the salinity concentration at a place in the river system depends mainly on the distance from the place to the shore-line, not on the distance from the place to the river mouth along the river channel. Thus, even relatively long, almost whole Ninheo river is suffered from salinity intrusion. It has been observed that at the Lieude station, which is about 23 km from the sea along the river channel, the many-year-average salinity in January is 1.976 ppt. As in Fig. 3. the whole area surrounded by Red river and Ninheo river (Haihau and Xuanthuy districts, Namdinh province) is suffered from salinity intrusion. The main agricultural product of this area is rice and the affect of salinity intrusion to rice fields must be prevented.  To protect rice fields from flood and salinity intrusion, dikes have been constructed at both banks along the rivers. However, the construction of dikes can only partly improve the situation. since intake of freshwater from the estuaries is the main source of freshwater for irrigation, and during the dry season. the intake of water from the estuaries is not possible for almost whole region. Also, at places near the coast, the transport of salinity from deep soil layers to the surface in the dry season can pollute the surface soil layer if there is not enough freshwater to wash the salt out. Thus some proper measures should be proposed to reduce effects of salinity intrusion on agricultural activities in the area. 3 Numerical Model 3.1 Governing Equations A numerical model has been developed to study in detail the salinity intrusion problem in the estuaries of the river system. This model is a modification of the flow and salinity intrusion model for a river network (Nguyen, 1985 and Vu, 1990). The governing equations of he model are: Continuity equation  Momentum equation  Salt transport and dispersion equation  Where H is the water level, Q the river discharge. b the river width at the water surface, A is the cross-sectional area of the river. C the Chezy coefficient, g gravitational acceleration, R the average hydraulic radius of the river at the cross-section, q the lateral flow into the river, S the salinity concentration, Dx the apparent longitudinal dispersion coefficient of salinity and G(S) the source or sink of salinity due to lateral flow. Strictly speaking, the system of Equations (1-3) can be applied only for well mixed estuaries. Study results (Vu, 1990) indicated that the estuarine number, defined as the ratio between the volume of seaward freshwater flow in a tidal cycle and volume of tidal prism, for estuaries in the system during dry season varies from 0.099 to 0.1l. Thus, it can be remark that all estuaries can be considered nearly well mixed. And! it is save to apply the system of Equation (1-3) to the river estuaries. 3.2 Numerical Scheme The numerical scheme for the integration of the system of Equations (1-3) is the improvement of that of Nguyen (1985) and Vu (1990). The river network is divided into branches, linked at river junctions. A control volume, staggered grid scheme (Patankar. 1980) is applied for the integration of the system of Equations (1-3) in an individual branch. The branch is divided into segments with finite length x. The water discharge Q and salinity concentration S are computed at the centre of the segment and the water elevation is computed at its ends. A first order upwind scheme is used for the connection terms in Equations (2) and (3). Central difference formulas are used for the salinity dispersion-diffusion term in the Equations (3). Since control volume scheme is used, there is no solution for the water discharge and salinity concentration at the end of the branch. Thus, additional segments with zero thickness (Patankar, 1980) is added to the ends of the branch and the solutions of water discharge and salinity concentration can be determined there. At the river junction, the condition of mass conservation requires that the flux of water towards the junction must equal the ones leaving it; and the solutions for water level an salinity concentration must have unique values there disregarding they are computed from the discretised equations for any of the branches joining the junction. An implicit scheme with the implicit factor theta = 0.667 is used for the time. The above mentioned scheme can produced a tri-diagonal matrix for the unknowns for each branch. By the well known double-sweep algorithm, the unknowns at intermediate points inside the branch can be eliminated to get a system of equations for unknowns of water level, water discharge and salinity at junctions. This system of equations is solved using a standard Gauss elimination scheme. Having, obtained the values of unknowns at junction points, the values of unknowns at nodal points inside a branch can be readily obtained by the Thomas algorithm. Since Equation (2) is nonlinear in term of the unknown water discharge Q, an iterative procedure is used for solving the system of Equations (1) and (2). The solutions of Equations (1) and (2) then can be used for the integration of Equation (3). The boundary conditions for the computation are given water discharge at the upstream boundary points, say, Hanoi station on the Red river and Trieuduong stations on the Luoc river. At Ninhbinh station on the Day river, since there's no available data of water discharge. the boundary conditions at this point and all estuarine boundary points: Nhutan station on the Day river, Phule station on the Ninhco river, Balat station on the Red river and Dinhcu station on the Traly river are given water elevation. The boundary conditions for salinity computations are given salinity concentration at the estuarine boundary points and zero value of salinity at the upstream boundary points. 3.3 Results of Computations and Discussions The numerical model is used to simulate salinity intrusion into the river estuaries in March, 1980. The recorded data of water level, discharge and salinity at stations inside the river system was used to calibrate and verify the model. The Chezy coefficient C and apparent longitudinal dispersion coefficient Dx are determined from the calibration. Figs. 4 depict examples of the comparison between the computed and observed salinity concentration at Lieude station on the Ninhco river. As seen, the model can simulate the real salinity intrusion picture with satisfactory agreement. Results of the computations for other days (not shown) also confirm this remark. Thus, the results of the computations reveal that the assumption of well mixed estuaries can be applied with satisfactory accuracy.  Using the model, details of salinity intrusion into the estuaries of the river system can be investigated. It is found that during the time of flood tide, the contour line of 1 ppt salinity concentration can cover all the Ninhco river while during the time of ebb tide, this line retreats farther than Lieude station. However. during the ebb tide, due to low water level in the rivers, the intake of water for irrigation is difficult. The salinity intrusion length is shortest for the Day river, where the maximum salinity intrusion length is less than 20 km, while for the Red river and Traly river, the salinity intrusion length is about 20 km. From the pattern of daily variation of salinity concentration in Figs. 1 it is remarked that the salinity intrusion has a strong relationship with tide. To study the effect of tide on the salinity intrusion, the numerical model had been used to generate the data at a number of points along the estuaries where the measurement data of salinity are not available. The generated data of salinity concentration were then used for harmonic analysis to obtain the values of harmonic constants of the four main tidal components M2, S2, Ol and Kl for salinity. The computed results (not shown) showed that daily variation of salinity at points along the estuaries can be predicted with reasonable accuracy using harmonic constants of the four tidal components. References Dac N.T. (1985) Mathematical Model for Flow and Salinity Simulation in Deltaic Areas. Report of The Institute of Mechanics, Hanoi, Vietnam. Patankar S.V. (1980) Numerical Heat Transfer and Fluid Flow. Hemisphere Pub. Corp., l91 pp. Ngo T.T. (1991) Salinity Intrusion in the Bac Bo Delta. Proc. 3rd Vietnam National Conference on Marine Science. Vol. 2, Vietnam Natational Centre for Scientific Research (in Vietnamese) 191-197. Vi V.V. (1985) Some Primary Remarks on Salinity Intrusion in the Red River Delta, Proc. 3rd Scientific Conference of The Hydrometeorological Research Institute, Hanoi, Vietnam (in Vietnamese). 177-183. Vu 'T.C. (1990) Characteristics of the Flow and Salinity Intrusion in the Red River Delta. MA Thesis, Asian Institute of Technology, Bangkok, Thailand 101 pp. Vu T.C., V. Suphat and T. Asaeda (1991) Study on Salinity Intrusion in the Red River Delta. Environmental System Research, JSCE, Vol. 22, 213-218. uqTM~!Mp!!  !L``! < G < 0   00  0 0 0 x     `   9 À ᇙ p ь` <31a  Cь` 31#  q` 603  q >00 g9 00>X g 00 0~p@ ?< 00 f88 0pf0 0`Cf f f? f>  f  f c0?      p 8    @000`Lp0x# 0 x0cc' 0 y8 0c c.< 0<8a80<0b c408|~1#00f c5$?|F0q #00c c5$>?Fc< #0ٍ4cc6%69aC>A #6<5$@`1ၘCr1p`@/x`1Cca1p+>1!C @gp`'`1yf>?"`1> ~?yDZ |0>   @!@0!0(`0 pa1>`0!016 `0`a#03<1!<ьa7& 8wc03|1!7& 0ɐ0ccl031ccÁ1 ِ0cx`0319p-~0pc8`01=) 0p2`a838xx % 0pǶaq|  `p`Ϟ?p `0`  ~` ```!0@`#8x0*8 sp0`1 p 1a`8`4">9cf1``4@" LB cg``4@"<pf`1h4 "~Øa>> fl>x3 !Af f<9.AB#c8@*>A& >bA!ca#x`8o ~,@6A I@2EC@2Io0g@6y0bi<8f@ffp@@`{?>90 >8H LCHl l` DdpD,F<400 x`@@@` ` 00  p<?`|>8? p@  800`  @ @ 0 x ~a0 f 130 f 0#0 0|00 0l C0  0d>  p 0F  p  c0`8`@@ @0    0 ٷ8 ُ@   @# @ ̀@"`xx `#`̀`#  ̀`  l@!  d'@! d'@ g$@$  @!?4ٰ@`#<ɰ@!? '4ٰ@! ?%ٰ@!? `çȀ@?`@ @! x~@"p@#`8@#   @$ <p0x # c"!< ! ?# `a8 #p # }&  d> >d6 >l6` }~ 0  y  ?0 ` "  " `1<"` Ǐ4"` lƑ4p"0 dѦ0" $60! $6! < $3  $ 3#`l#0? `l'9 psg !1 p00   ` 0 ` `$@8@$@0 0@&` H& #>I@%  $I4I@&  p$I@#"!2!<``0 #0>|@"8" aa0c@g@~0`n `@`  p`0?'< ?0 # /0$?0$0`< `$0@p#`C% Cy8% Ba͞|x % BAMl`*a AMd`0*?~ 9Ml `&}>Gl0' |>0)#π}&8$ ̀m~8$0m}~ `">Oq  >80`0 `0@`` ``0p 0 ` 080p8p8p` x`"`80"0`"88"`!8 `  8 0`38 0"``@!@` 0!@0 ?8!@"@p @p0?"`!` $<##0``x8`$9%À#<" |8" 88"  `$ 0% " 0?  g ` L0 @ H` @lH @s@<`` `  00 0   ` 00`@p   0 0  00 `   <  | < @ 0  ` 00<  x0 @  0<@ Y6` & 0 0  ?{  ` ws0 %  ~90)c<<+?c<%,?)>D s~!*?D #f7?* Dcg>?$  F ~sg! 0F <$ 0@`ü# f@0@  @0` @p  pg0 xf10  LdO3 Ff60 Ff`   0  @0@ @``! @s#d@19x3gap@!   @@ @@@ !@`x> %@ #%tɐ81  0c*<883BA@*<ɐ= 8!fa@-,oafa`+d`0 afa`-`@ a&a`(@a&a`(@ .afa@(@>3G!f# *`w!c?0+<@p1 |88*  @?|( 1 `>&@1?>c#b??f# !> !fc 00a7  0!1 ``<0 @0 0 B 0F0|F 0 `@   A<6af"0aBf0  l0 0 <0 0 00 0 `?A @A`@ Axc2 `       ]]  \ <l@ '~ t<*~@?,0@/@ .0@@)0@<~~? @vb@@@@8B|0f@~>@~ x"b|p@$ gsf@/ @&&@  / |@64@ 0`/ @3@@ 0<|yM* @#l@ 2͛2~?  <@  . @  Fǁ-@  FƁl-@  w͘la&@  3fF ! @0`@@@#`@ @$0a c0@' q>1`1@' y?!wx3N@' I1a f` "hH@ ( M1! #hH@(` G1a'd A6&hH@@(00 1&d c66l@<1  A1?7 3G@@, @"'`@ ;0<~ (@@ r~~ (@@ 7f (`@  ?2B (0@ 72c $@  ;3bc)@  83~b-@ #<0C`@ , G @)0 L|xϘ>@( x'ف33@( xdG |a#@* ddG a#@+ `FdǴ !#@& `Cdf͍ٿ!@,` Ad@ƈ@ %@@% @f`x*@  </@ ><~m2 @ >B~~Mh3 @@ ?b~M4` @ b6vl=` @B 6gt;@?? @pfFL08`8a0!@ 6~|p 8 @0`9#@ >>< 7>p @?0<3&f`@ 03f  $$`@02 0 $$`@20@1 $$`@810a0a&d`@30?1?0?# @/@@~#f@ ~ >@ "@p p@%@$@%@$@8%?00@' @& @&`@%@`@8$9@@   @# @@#@ @$`@``#@ &@ #9@0( @`0# @ (8@'?@0& @A`&@A(@3 '?@? $ @ '$@ &70@?' @p0! ~@@`( `@0' @@?&6@<@?#$@(<@'<@?'@y0'@`@$8@!@'@7@%@`#`' @7 %@7 #@$ 0$<@,@'@="0@  <@$0@?#@ !p@ " @   @`$ @$ @@$ @@$0@$@8 #@ #@0 #1@`"@">@"3@"@@`@@@  @0`0 @` 90@7 ?@0!@?%@"0`p@!<0p0$$ @0$$`@0&@@0~ #@ @0@&@0%@@)9@@$0@@! @@ '80p@@#@`@&@`@ -@@@@ @@@' 0@@@+ p@@@%? `@ @@@@@#~`@0`@&: @@` @$@@ @&@0@)@@0@@@0@@# p`@@"`` @$0`?0% @` @ @ @@  @'``@@ @`` @%@p`@$`0@$ 8``(p  8@#1  @p  @ ?8  !0  0 @`  0@ p` < 0    0 >|           ? @                 `  `0 p   `  ` 7>  |@  `@ |@  ~`0  &  " f00 ~@ @      #  g @ @ @  y x x  p     @ @ 0`0 `x   0 0             0 x  `@  C$ p8$  q0 8y! ` 8a$ 88#  <*p .  <a0ፀ 84  ??` 6  | )p  $  p+@@| p +@@< 0` Dr s-0 ?0 p3!0??AFp_1?+@ 1 @@78xC1$>8#p;|Gπ8&"pA#p@8̀@ f#`a6p?C H`9;yLEN C#`91@IcD  r#?~@~`Hfo. 0|>f~`HHf'c'<|3~wCrf`l'~0f``~Csgc?@ f0``pb" 3 `"v` 7& `>|` BNv&@ ~}&  @ 0@ 0| s"f R B B"B6b πb|@̈ Ĉ Ĉ Ĉ f f|@ >f_g yȎ | )}! *  / @  0   * `@00"`@0`@0`@`@`` `'$00`&60` / <0 `   p0 ß6` fXX-0  8 `#XX.p   ~`C'&xX.` 0 f`c66l/` 0 p7o/ < ` :0( @ % % & `%?0$?#  ,  |2`| `@2p  `4p   00Aߟ4 0 gdـ48 `l>@؀58 00o@g؀58 0 o@ll؀580`; lLlߘ5p `3ݙgϐ0p  |8+` ) '0"0`$00$0$0$0 #0"0#`0"`00"p00`"00`p" 00`8"0"0 #80#8 0(00'p0'`0(`0``(0 0~0*0 -0  /0 /?00100 2??`?0?`<9@*p`0`- 0009`.00?0?0, 00( ( 0 (0( 0 '  @) 0@(0@(`8`&`0 @(p@`0@)`@% @' 0@) @)@)``@%`@'` @' 0p`@# 00@& p @'0 x @'   @%@' 0@'0@$0``'0`@'00`@&``@&`0(@~'@' @$`' 6`'0@v `'``& 7`&0``$00`$0``$  `$  `$ `"`% 0`$0`$0x` ``"`0`#`0`"@`` `!@`!@x`"``% ``,0`0  >`- 80` `2>B>d 2 ?;>Be< 3 33`BM 3 !!`B2>E 381!`B2fe| 20113f3vd݌ 21!>~> x@ + ``x $|80 $v "~ "  `"8  $`p )  @ +x 0 - 0<< -  2~~` / ɛ3Ff  0 p 92F0 2 y ;2F /693fF / &ٱ1~B  0p  '  "    x p `  l #x0l + 0)  @ 0+ 猏  0) ͛; 0+ L|@ɱ 0,OGɱ 0/lɱ 0.2lAͱ 0- lGo̟ 0& 0` 0  0π      ~cc6b~@0` ?`3  ` @`  g"8a"<a !89#p0!!! `p p 08  0%0$x#1a@ %! 9ǎ8 '#0%7p>0q& cc~p 8"?|.2 ??` 88 ?Î8 |0`.?00p>A 0010?@qIAp0?csp20~p`  QJA2?;ၘ"2s`` |>p [ÀNA ?;`0>~`OOǀMA y 0?!r1??cb`NlcMA0   01 213b30pCz6@foP w?`013 0`0`pCsb`0|@I~`0?0p`0`"c3c 7p1c320`?~`>" ??7  |@ dg>~f`<}π#bD`숀k#bD`Ĉ{!bD`Ĉ;#cd`Ĉ1>"d f f|V.V..  .L``. 8> ||f0c0B0@0p7 a@x00##0a#0a" 0a # C0a 3 ca 3 ca 30|>a #8aa@ 3 |`|  0 @0@0o@0}@0@0@0@@ @@0@@0@@0@@0@@  g0@  8  B F@00!|B f@0|0!@~f@ 1#a``ff 0#c@ Bf 0 #@ B7B  #0  "   "    " "  @c  ~ c ?|  8  @ `    &|0x <@ &81p @ %00@  (00@ xa & A0009 cx~|??p pc ? B<81p8c0s fa0aL>f ` s b 00 B<08~1ǰ| c`0 01`cb! " 0F0B<00 xcOp01c! 4 0F0B<10  !c pQ0v c! < 0F B<00 #c p0b c1b! $  B<080 #c01bc 40D B<10 #B `q f? 0b1c 0D B< 0! #B0 p b1Bc(0 B80! #C0x` Ƹ 03ac 810 B<0!s |xǀy?s8!` (c0 ` .c0 8@`w >0 )8c||ǀ>~a`9 >40p8!xc0f D`a408l!0 `0B `0` !Bpa`0 >41!0`B  >#0b|a `0D >41!1cB` f? 0brAa`0 ?5 !1cB 0 a <BbAa`0Ā ?5A 9!1 b4B4x``ư Fga`1 >4c0!00 s|G}`|x?ǀ|?? =3c1q08x0?p@0 + >|s88@`0 `0 x< x 08    0 00 08 (a00)pap0( `p - `88 p~ : `)9gÌ||~ ;>4 xcb0fDaaxL3 DØa 0 >4& 8a0 `0F ` a`#0BpA 00:4F!a1`B >#0B|>A 0:4!a1cB 0``>0fbAp0?5 a1b0B `a`80BbA 0ŀ?5 a1f49 x``ϰgA 1>4 a01s@xp?x??6ss8p1 # xq00@ @0p8!x0 ? 0    0 0`p`c*`#'>pc &|`` + ``@  p\>B `1s1?ax# `0B<  xc`1aAxL31`#a @C@:~8 0ca1`0!a` #01` !!@AB<!a1 ` ᇃ`<#0F! |!!A@:aa1c0 a@`>0F!`x!!AB<`a1f0! @`0!`!!AB<`9a1f< y xa`x1a0c !0gB< a11c@x?{0s|0n8sc0!@( <pq@@x,@p8;    0 0               00 00000    00 00000 ?0?0 0  @@  0 00000000 0'>0G G$& G$& > ;Ø0 0 000@ @ `0   0 000 `0@00  0 000000000 ` 0        800x000` 00   0p  @ `p  0 0        0  0 @` `@ 0@ @    0p         @0@ @          0P@  8g   @!"    1 0  8!`0  0 < !   ! 0  69"0 ! >0     8 x 0  0@0 8p   xpx |@@@ !@ @ !@ @ !  @@ "?@ #  @ $ @@` $@H $@cH $@2`h $@7`x ! @@ ` @ 0 @   @ ! @   !   !` @ !@@`8 !@@@@ ! @@@ !@@ !@x@@  `  !        < q  A A dA  8   D 0| `0 !  " ` 0  !p 0 ` "q # $  $@ "@p0 % `0 $  $  $@ '  !@@ &D @ %?  "  #A @ "0@  #  $  $@ #@ "@ $ @ &   @ &@ @ "`@ &  !  &  " 0 % ` ! @   @ !   @  @  !  %    %@p@  & @ & 0 &   &`  &   &`  &  @ !@@ `@` # 0  &      @@ "@ @  &@  &0 @ ` @0@@` $`@ @` ( @` '1` $0 ` #`@ (0@ %<0@ ' @` @ ' `0@@ ' <8 @ %  @ "  @ ! @ " 8 @ !@`@0@ !ـ0€@ #  @ &0  (` ! ( @ 8? "p@` ``?3 &g 0 '!   (0` x (  0p *` |x )@2` )@08< *@ *``0p .`  3gf8  /t0x `0 44 x0  4 `88 3 DL@@ 0d@@ 0`x`@  * 0`@p *  ` * F` @ ) # 0 ) 8 8` )  0 (     $@ p )@p @8` (`@0 $ `8` $ %p` (  $900 #0 8 !0`p " "` & 88 &8` +< ` *` *  p )0a@ +`1`  ,@0@80 %`pp0 *`` '@<  @ '@| 0 $@0 (@ @` +0 @  +   (@ 0 00 ,`@00`` +0 +0 + >` +8 '   0 +  @ ,@ )  %@ 0 %0 ` ' 0 # `   @` x   # `~ $ `     p  0   `0@  `@`  @  0 @` @  0 p@     @ `      00 `@  0  `@ `  `@ @@  @`  >  `   `   @     |  ? @ @ @   ` 0@   @ >` @        @   @0  @@ @@ @ @ @ @ @  @` @@` @@ @" @2 @!  !   ` %0@ 01` ) @ x8c;g1 ` )`@@`m" ` * @@   + @ 0@0 0 , @ 0` ,8B c ,B "b0 - b g9&0 (  > 98 @ 0  ` @ 0 0 0 0 00   0@0 @`     @ `      @x @ `<  `f @`<p @@8 p @ @` @   @@    dd  6d|0 @6d |0 @>D  0 @6d0 ` @6<   @@  @@  @` @@ #@@<| (@@n (@0B1c (@ 0@1" ( 0  06 (`0 06 (@ > 06 (@0` 06 (@0p 2`* 1"`0-0 l1c`@`8p.0 l1f;`0.`0pN'=`Ã~|`?s)`` `1fw0!@?`qvg0`1w @ ?`s@`p@p@ p@0` p<`@```~;c8`@@`1v `?81݌w `0Ϝw@`06`06< `070``0 7 ``0 s`?``19c@```?````x``````````````p `?``` p  À xp31 031``pa`?10a~>s30a '30aa 0s0aa0 0 aa0 0a!00 00`3 00 008  0  0  8bNTNDNdxpr  N"N currentpoint ",Times .+b, Symbol( )H( ) t" (+( ))Q(*)x"((<=) q30 dict begin currentpoint 3 -1 roll sub neg 3 1 roll sub 2496 div 928 3 -1 roll exch div scale currentpoint translate 64 43 translate -9 533 moveto /fs 0 def /sf {exch dup /fs exch def dup neg matrix scale makefont setfont} def /f1 {findfont dup /cf exch def sf} def /ns {cf sf} def 384 /Times-Italic f1 (b) show 260 286 moveto /f2 {findfont matrix dup 2 .22 put makefont dup /cf exch def sf} def 384 /Symbol f2 (\266) show 531 286 moveto 384 /Times-Roman f1 (H) show 339 834 moveto 384 /Symbol f2 (\266) show 613 834 moveto 384 /Times-Roman f1 (t) show /thick 0 def /th { dup setlinewidth /thick exch def } def 16 th 237 434 moveto 595 0 rlineto stroke 925 533 moveto 384 /Symbol f1 (+) show 1250 286 moveto 384 /Symbol f2 (\266) show 1458 286 moveto 384 /Times-Italic f1 (Q) show 1295 834 moveto 384 /Symbol f2 (\266) show 1503 834 moveto 384 /Times-Italic f1 (x) show 1227 434 moveto 524 0 rlineto stroke 1875 533 moveto 384 /Symbol f1 (=) show 2202 533 moveto 384 /Times-Italic f1 (q) show end dKMATH? Q bHt+Qx=q-1 " " "dxpr  """ currentpoint ", Symbol .+ ,Times)Q()t"(+( "()x" ( 4Q (=2 (8A"4 .( - * (- ( C * .(K+) gA( d)H(f)x"c(w+( g)|)Q)|)R(A (2 +C (2 +R"9(=) 030 dict begin currentpoint 3 -1 roll sub neg 3 1 roll sub 6656 div 1088 3 -1 roll exch div scale currentpoint translate 64 38 translate 23 355 moveto /fs 0 def /sf {exch dup /fs exch def dup neg matrix scale makefont setfont} def /f1 {findfont dup /cf exch def sf} def /ns {cf sf} def /f2 {findfont matrix dup 2 .22 put makefont dup /cf exch def sf} def 384 /Symbol f2 (\266) show 231 355 moveto 384 /Times-Italic f1 (Q) show 99 903 moveto 384 /Symbol f2 (\266) show 307 903 moveto 384 /Times-Italic f1 (t) show /thick 0 def /th { dup setlinewidth /thick exch def } def 16 th 0 503 moveto 524 0 rlineto stroke 617 602 moveto 384 /Symbol f1 (+) show 1025 355 moveto 384 /Symbol f2 (\266) show 942 903 moveto (\266) show 1150 903 moveto 384 /Times-Italic f1 (x) show 919 503 moveto 435 0 rlineto stroke 1619 355 moveto (Q) show 1906 184 moveto 224 /Times-Roman f1 (2) show 1740 903 moveto 384 /Times-Italic f1 (A) show 1602 503 moveto 474 0 rlineto stroke 1403 355 moveto 384 /Symbol f1 (\346) show 1403 894 moveto (\350) show 1403 675 moveto (\347) show 2088 355 moveto (\366) show 2088 894 moveto (\370) show 2088 675 moveto (\367) show 2353 602 moveto (+) show 2652 602 moveto 384 /Times-Italic f1 (gA) show 3147 355 moveto 384 /Symbol f2 (\266) show 3355 355 moveto 384 /Times-Italic f1 (H) show 3203 903 moveto 384 /Symbol f2 (\266) show 3411 903 moveto 384 /Times-Italic f1 (x) show 3124 503 moveto 547 0 rlineto stroke 3764 602 moveto 384 /Symbol f1 (+) show 4547 355 moveto 384 /Times-Italic f1 (g) show 4733 355 moveto 384 /Times-Roman f1 (|) show 4832 355 moveto 384 /Times-Italic f1 (Q) show 5099 355 moveto 384 /Times-Roman f1 (|) show 5225 355 moveto 384 /Times-Italic f1 (R) show 4115 903 moveto (A) show 4357 732 moveto 224 /Times-Roman f1 (2) show 4867 903 moveto 384 /Times-Italic f1 (C) show 5143 732 moveto 224 /Times-Roman f1 (2) show 5677 903 moveto 384 /Times-Italic f1 (R) show 4066 503 moveto 1860 0 rlineto stroke 6050 602 moveto 384 /Symbol f1 (=) show 6375 602 moveto 384 /Times-Roman f1 (0) show end dMATH Qt+xQ 2 A()+gAHx+g|Q|RA 2 C 2 R=0t   dxpr  " currentpoint ", Symbol .+ ,Times)()AS) )( )t"($+( .)()QS))(7)x"-(O-( \(Z)x"Y (oAD +x ( )S()x" .( h * (  *  .(=) G) ()S))z30 dict begin currentpoint 3 -1 roll sub neg 3 1 roll sub 6432 div 1024 3 -1 roll exch div scale currentpoint translate 64 43 translate 23 318 moveto /fs 0 def /sf {exch dup /fs exch def dup neg matrix scale makefont setfont} def /f1 {findfont dup /cf exch def sf} def /ns {cf sf} def /f2 {findfont matrix dup 2 .22 put makefont dup /cf exch def sf} def 384 /Symbol f2 (\266) show 251 318 moveto 384 /Times-Roman f1 (\() show 416 318 moveto 384 /Times-Italic f1 (AS) show 856 318 moveto 384 /Times-Roman f1 (\)) show 335 866 moveto 384 /Symbol f2 (\266) show 543 866 moveto 384 /Times-Italic f1 (t) show /thick 0 def /th { dup setlinewidth /thick exch def } def 16 th 0 466 moveto 997 0 rlineto stroke 1090 565 moveto 384 /Symbol f1 (+) show 1415 318 moveto 384 /Symbol f2 (\266) show 1643 318 moveto 384 /Times-Roman f1 (\() show 1776 318 moveto 384 /Times-Italic f1 (QS) show 2259 318 moveto 384 /Times-Roman f1 (\)) show 1702 866 moveto 384 /Symbol f2 (\266) show 1910 866 moveto 384 /Times-Italic f1 (x) show 1392 466 moveto 1008 0 rlineto stroke 2492 565 moveto 384 /Symbol f1 (-) show 2899 318 moveto 384 /Symbol f2 (\266) show 2816 866 moveto (\266) show 3024 866 moveto 384 /Times-Italic f1 (x) show 2793 466 moveto 435 0 rlineto stroke 3493 565 moveto (AD) show 3989 661 moveto 224 ns (x) show 4528 318 moveto 384 /Symbol f2 (\266) show 4736 318 moveto 384 /Times-Italic f1 (S) show 4535 866 moveto 384 /Symbol f2 (\266) show 4743 866 moveto 384 /Times-Italic f1 (x) show 4505 466 moveto 448 0 rlineto stroke 3277 355 moveto 384 /Symbol f1 (\346) show 3277 820 moveto (\350) show 3277 601 moveto (\347) show 4965 355 moveto (\366) show 4965 820 moveto (\370) show 4965 601 moveto (\367) show 5261 565 moveto (=) show 5584 565 moveto 384 /Times-Italic f1 (G) show 5867 565 moveto 384 /Times-Roman f1 (\() show 6012 565 moveto 384 /Times-Italic f1 (S) show 6218 565 moveto 384 /Times-Roman f1 (\)) show end dMATHd (AS)t+(QS)x-xAD x Sx()=G(S)efHp:pp  p@p O`@@@@`&`@@`)@`0%@`6, @`71 @@`71 @`71 @`71 @`71 @`71 @`71 @`71 `71 `71 `71 @`71 `71 `71 @@`71 @@`71 @@`71 `71 @`71 @@`71 @`71 @@`71 @@`71 @@`71 @@`6) @@7- @@ 3@@0"3 @1@`- @@`4 D@@`1"@@`3A!  @@`7 $| @@`71| @@`71< @@`71? @@`71~; @@`71@! @@`71 @@`71@? @@`71@= > @@`71@9| @@`71 @q< @@`7 $ p @@`71  ? @@`71 n@@`71 @@`7 $@@`71@@`71?@@`71`@@`71 > L@`71~0L@`71Np@@`71?N0@@`718~0 @`71<?ρ@@`71 @@7-0@@`3(0@0!p40?1P? ?@ 00p@?@`1 @@`/#?ǀ@@`6,0 (@`71 (@>@`710(@@`71(@@`;10hH7@@`<߀10HL7@@`<8108l@s@`< `1|0,@@`71>0x(@`;1P~X8@!`;1xp|@!`:1>8p c`:10P c`:1Pp8!`:1<p|!`;1< `p0#@`;1 `0x@@`:1 `@@`:1C `@`@`:1A`d@`:1C`  `;1@p  @`;1@ @`71`@<@`71`@D|@`<`@1x@4`<x1@>?4@ 6(x>@<3!88>0 d1?@`2 ~p8@`0? @ @`/%Áx@p`7(@ `71  @` `71>8`;~1@`;1`@`71!`@@`710`@`;1``;1`@`;1`@`; 1` _`; 1A` @G`; 1A` @W`; 1 `;1 8|`;1  @`71 `71> @0`:1> @0`;1 ! 0~`;1&! 1`;81?! '`;1!@'8`;x1?1@7@;11`?@:1J @<:1@p`'2(n 4!~<6 = 6 = (@3 8@30 0@1#x8; 1<x8@|;1  8  x;1    |71 0 x  @71  8 @0;1 <8 @0;1 , |@; 1, |@p; 1 ,  @`;1 0,  `0;1, ?  `0;@1 (   `<;1 ( & `71(& `71X, > `<71  >`?718~( @ <`<`718@ >|`?`718@ >\@`718 >\@`?718  \`:18 h  |:18(   a:18 ><  a;18~  c;-8V  a710XT0> 7-0X\0 6'0X<0 6 0 0  7 P 58>P0`6 |:X`P2> 0`p6 80 3@#@:> 0`p9+?: <` 3-?p  <@` 71p> 0@` ;1` p 0``<;1 p   0``~;1 p  ``>:10p   ``710p  0p`<71 p p0|`<;1 p  | ;1 p`  ``;1 p` `` ;1 ` `| ;1 `p ``|`;1`@`8 x``;1 @`80  8`X` ; 1 ` 0`X`71``@ 0`H 71 @`@ 8@H`:1``00@ |@H`:~1p<<@? L`;1x8@ `;1 `@/ `;1 @0x/ % ;1 `0> ~7*0`0h ?68&0@00ho7 0@00  5@  4> | F 5 `  F 2p B 6(|  ?B ;1>0 B`:?1  !C`:18 ?G`:1( >C@`710 C@ 71H C@ 71H  G F710x  C`F71D0  A F71D0hp A 71D0x A`71D0` ÀA`71D` A 71Bx` A`71B88` A`71g8@ A 718@ @ ;1退0|@ @@ 71>@ @@ 71:@0 @` 710:@ @ 71@ @ `71p>@ @| `71P<@ @d 71?X @` @t 71 X@` @t| 6,+@P@ D@|1>@X@d?4!>@P   2 @< 4  <@< 2?? @<`3q0?! @ 4` '0 x @ : @ $ @ `: ~$0 @< :4< @< :46@ :40 @ :4 @`:4  @\`:4aa<  @ `:4qs> @`:4? ? @`71 P @`71 @`71 @<`71@`71π @$`71 @$`71 @$`71 @$`71 @d`71 @D`71@D`71 @`71 @`71<@`71 @`71@`71@`6$`0&)! /&G` @ ||   | <0 08~ 0q0 !0A? 0@ 10`  a0 D CL `L 8 L``;90@8<  ?>    (  0%0 02|;0 @0 `  0|<?10`@8 `$#<000 `D`: +C;8f810 L@D@;1 7?c000f0 F D@D@=9 ;1@@    x x 10@` @ `` ``)  8<`@` ``1  `, `@`À`` /0` @   `@` ` / @@   ``o` ` / @@  `~qa 8 ` /   `c!a1 ` 2%0w `` `08+`8>8` x``0` <6?`?g3c ` 0``0@ <6  0#! ` `| h`a`  ;0 0a`` |x`A`a `;500a ` 0`"`"`<60!0q` 0`"s@8.03@0a `0pyp@@5'0<@@ # ` >@2# ``F ! `??* 2  `0p`H}}@&x8?@>/0` `@  @@ @L ` ?@ ( @`o  @` @` @`!0 @`#0 @`)!0 @`3-00 @`5' 0 @`6?' 00 @`50' 0 @`1 ' 0 ``4 ' 0@`6 ' 0@@71  0@ 5+  03'  0  @4#  0@@3  0 @`5   @`5    @`2  ``3 @`3 ?#0p ``5 ) x0~0 ``7?- 0D @`71;  0G ``71  0G ``71 00 ``71 00? ``71 00` @`71 0?0!` `6, 0'0!` _3& 00!`4" 0O0 ?``3  0 P ``5   P ``5` r W ``5` z03`  ``6   Up ``6` %0Up ``7`?+00 ``93`?? p  ``93`0@0M  ``93` 0L00 ``93` 0H00x ``93` ?00S0 ``93` '0A0S< @`93`0#0@0~< g7,`0#00:`4% !00``5! #0@0 ``50#0@ ``4 "00@9 ``3 6 ƿ@N ``2 &PN ``2 "~ ``30 (0|  ``71? 00 ``710800 ``710|0@0 ``71 >$ 0|  ``71 >$0  ``5+ *@0p L  6`71  @0P L `4$ 0p @0P |4& > 0 <"?``5"0x   ``50 ` < `@40 0` O``1 0`G``5  p O``4 !! G``30'0?  ``7?-~  x``6&0 0 ` ``:100 @ X``:102 00@  ``;<102-|00@8``;106/<00@| ``7106?00@ `710400@X 7* 40@P7%04 0<PO``5 !0 0<P `5 4,04``6 `+``4 0@ "``5  #>``3$8 ``6*8` ``91?0( ``90-( 0@``90-A @``:1 A 0@@``:1 A @@``710 0 @@`@71000 @@`?9@*000 0`9p%00p8`B``8 0P|0@  ``700|0|x``40040<x``3 0x|p(``3 ?0`8``5 0&0?P x``5?+נ00 P`x``71?0  P x``71 00  `P9,``7100 08 ``:1 0 0`8``:1  0```71 0 0 `y `5+ 00; ``6, 0 0 `9< ?7'00A0``7"00A0  ?``6F00~  )` 8B 0 0x``6B A !'``6B0A '``5N0a #;X 0` 60(01 0` 71?8!@1 0``7100!9 0``900!(18=0``:100110À``:10} 1pC0` 8+007 @18A0`:100 @!(0``:10w!(06p(0y!p( ?`8$0a0,  `70 0  `40 q0 `CD 40 q0`GD 70 +~` D 7F0 "`@@ 7B<'A0@@ ` :B-a@ 4 :B1@ < :B1@A0V@@  :f1AB@` :1 Ax>@@ @` 710 A|>@@ 71  @ `:10    s6(0 @  `7$`  `8 `  @_ 5 ` @L 7`\` 6` p< 7?  !8p 7?&px0@ 6' 08@ 71| 71px0< 71 P8$@ 71` @ 71>` @<@ 71x8@ 710V  @ 91V @06)>@0`9%8>  @ 8 8 @  0@@ 4 (a AÀ@~ @@ 5 (! |@^ @ 3 7 @ @ 7"@ @^ @` 7#'< @ P@> @` 9!,   P@ @ :!1   @ @ :!1` <  P@`@ :#1 X@`@ :1` h@`@ :1 F`/h@`@P 71 h@`@ 71 > `|@ @p07)8< ` <@`g 9&8 0`0 8  0_ `@< 6p@p@ 30@@ 6 8@@8 7  @P@ 06%@@009!,>x@@00:!1* @@@ ;1  @@@ 0;1  @@@ 8+> @@@0 8+ ` @@@N 71 `@@@ 4.@b@@ 09-` @`@@?9'` @ :p 0 8 `  *f@ 8` @`@ 9>c @`@0 8 @b@00; $ @f@009)p@@00<4p8' @>@00=4p8 ' @>` 0 =<4 8 ' @>` 0 >4   @`00=4 X 8@6@00=<4 H <@@00=4H  @@0 <0X   @@@9*X  <@@ 07%8  8`0078 0 ! `0068  @>`0058x??0 @@>`0058]  @@@`005$A.0@&`004 )@>>0@.`0071\  0 x@>`0071x 0 @`0071( 0@`@0071( @`@0071( 0 @`@71,  @ @`@07'(   @@07'@(> @@`5(@,: @< @ß05#@, P @< 05@,  `4@,  @@l6`` @<@`6 l?@`@@`x4$Ϙ``@ @ 06O+L``@ @`071 (D``@ @ 8071(@p`@ @`00071x@|@@ @`00710L4@@ @`0071L4@  @`71L@ @`00718L@x@ @`907'8 0@ @9;4) @w07'@ `06!9ß  59` 06` @@68 @@`p5?!8 @<@`p6*0@<@`p0710@4@`710@``071@@`71@@ 71@~ 71@`n`7- @@ `071 @@ 04)@@>07'a6! 5 `8 4 @ 2 L@@ 2" @@ 04)@@ 071`@@`06&@`71@@ 71 @@ 71@@ 71@@ 71@ 71 @@ >4)@@x1&@`-%' `$  #0  0? 0 0  00p ! `p !8 A30BA2~ B2c p"bb "`& B@0&f`0&f  0& bA& 0&0c a& 0"?pF 0"  00 c `  ~ `  `  @  ~ B`? C8 $@ B wpcF %@c8?FD)B@@ÀcaaFFD7 "BF|CC A@ FD?D90"% 8sa#a`bBA dDaD90B+!~ !a1# a`0#B fDF <028C>08C AA#A0`0#B@ Fdf<6 na0AA@#A0`0#B@@0FfÎcA<6 10 A A@#A0`#C@ag |@8.01 !AB@@#A`cc37-01 !AB@@#A`s-#01  AB@a1#a+!01 1 Cc3s!@`(01 A1 fAr` 03 c >@@ 03B>08@ @0  00&0 p @90@Æ 0 08< 0``90` p8 n80`0` :0   00`0A0 80   0`pa0800   00`pa097 `  0`0a0:8| Nx  8?v>`0`00<60` 0f0qt?`810``000<60`0 0a `0` x`00``00<6000 0A p`C0` b00@`a00<6000? 0A8`C0`200` c00<60 0c 0A`C0`0>4090` B`< 0`0 0A`8C0``<8q0```0< 0`00F 0Ab0C0`0a0``1`0<6800F 0AfpѐC0`8q0ap1 <6 <800 tqѐ0 |?~}x 6,` x0ϟ`<>x   ``08@`@0@   @ `@@ ```0?a0s a0@0as a ac10a3 c3311?0a2 c3q  0D0 a2 c20! a2 b>a!8a# c!d0!ac a!| !?>aq 1 <0F p 8= 8@`a@`@@0``@ ` a  11!0 01`1!0 3@31 F#@9BD >0c@ B? `1C@!f? C803`!$  Aq 1<0 FC C W Equation? bHt+Qx=q Equation Qt+xQ 2 A()+gAHx+g|Q|RA 2 C 2 R=0 Equation (AS)t+(QS)x-xAD x Sx()=G(S) Seminar on Environment and Development in Vietnam 6-7 December 1996 Australian National University, Canberra FO File SharingAudio CD 1 PDFABTFLpds 1 >;pFFAudio CD 1 PDFABTFLpds 1 >;pFF,zD\]'(_`uv-W-X6  <xz  @ PR BH B= B8l P P P6*+78j   7 8 [\^_&')*+IJ^_abtuwxabTUjk##$9$:&c&d'['\* % % %!%#%%.%  %  %% %!% %  % % % E* * *6*7++-V-W-Y-Z/ ///222223o3p33445]5^5566 =Oxyz% %k% % % %% % % $ 6@6J%,,63Enop >m)z)2G2R6%% % % % % % % % % % % % @% % % % @% % @% % @% % % z* z  Zyh >   OWX7"""R"S+%:0zzzHH .1@{,, ,-^&'dA A.hnnnnBH-:DW-C 2.0 (AppleTalk) 2ChicagoNew YorkGenevaMonaco Zapf DingbatsPalatinoASCII Zapf ChanceryTimes HelveticaCourierSymbolI Courier ObliqueB Courier BoldBI Courier BoldObliqueBoiDiep NiLienThuyen LocUyenVPSTimesZ UHoaiH 1.1yPhuongThao 1.1 Minh Quanh VPSTimesHArial Courier NewMonotype Corsiva Wingdings > HoangYen 1.1 MT Extra HoangYenH 1.1 PhuongThaoH 1.1  MinhQuanH 1.1 StarTrekTNGTitle StarTrekClassic StarTrekTNGCrille AnhMinhH 1.1 StarTrekClassicMoviesO VI Gio NamW VI SystemX VI MacSystem7VinaPala"Industria Solid"Industria Inline"RIndustria SolidA"SIndustria InlineA+%VI Gio Nam Test:0 Brush ScriptzVtopiaz Anh Minh 1.1z VI UHoai 1.18Embedded_Object_1Embedded_Object_2Embedded_Object_36H`H`_H`D6HHH2G2R2G(,(} 2np >m)z)2G2R66 <x<*<< < <<*U<j* <*62<26