scieee Science in your language
[en] (orig)

Modelling the tidal dispersion of 137 Cs and 239,240 Pu in the English Channel

Abstract

A numerical model to simulate the tidal dispersion of non-conservative radionuclides in the English Channel has been developed. The model solves the shallow water hydrodynamic equations and, simultaneously, the suspended matter equation which includes deposition and resuspension terms. The model considers that radionuclides can be present in three phases: water, suspended matter and bottom sediments. Thus the equations whose solutions give the temporal evolution of activities in the three phases must be solved too. The transfer of radionuclides between the liquid and solid phases has been described in terms of kinetic transfer coe$cients. The model has been applied to simulate the dispersion of 137Cs and 239,240Pu released from a nuclear fuel reprocessing plant at La Hague (France). The model gives a realistic estimation of the activity levels detected in the Channel. The transit time of radionuclides from La Hague to the Dover Strait has also been calculated, using the cross correlation function method, under di!erent weather conditions.

Read accessible full text

Modelling the tidal dispersion of 137 Cs and 239,240 Pu in the English Channel

Author: Periáñez Rodríguez, Raúl
Publisher: Elsevier
Year: 2000
DOI: 10.1016/S0265-931X(99)00125-3
Source: https://idus.us.es/bitstreams/e2e4e4af-e29a-491a-b3ef-a145e1191bba/download
Modelling he idal dispe sion o Cs and
Pu in he English Channel
R. Pe iaH n ez
Depa amen o F&nsica Aplicada, E.U. Ingenie &naT&ecnica Ag &ncola, Uni e sidad de Se illa, C a. U e a km 1,
41013-Se illa, Spain
Abs ac
A nume ical model o simula e he idal dispe sion o non-conse a i e adionuclides in he
English Channel has been de eloped. The model sol es he shallow wa e hyd odynamic
equa ions and, simul aneously, he suspended ma e equa ion which includes deposi ion and
esuspension e ms. The model conside s ha adionuclides can be p esen in h ee phases:
wa e , suspended ma e and bo om sedimen s. Thus he equa ions whose solu ions gi e he
empo al e olu ion o ac i i ies in he h ee phases mus be sol ed oo. The ans e o
adionuclides be ween he liquid and solid phases has been desc ibed in e ms o kine ic
ans e coe$cien s. The model has been applied o simula e he dispe sion o Cs and
Pu eleased om a nuclea uel ep ocessing plan a La Hague (F ance). The model
gi es a ealis ic es ima ion o he ac i i y le els de ec ed in he Channel. The ansi ime o
adionuclides om La Hague o he Do e S ai has also been calcula ed, using he c oss-
co ela ion unc ion me hod, unde di!e en wea he condi ions.
Keywo ds: Model; Tidal; Dispe sion; Radionuclides; English Channel; Cap de la Hague
1. In oduc ion
A nuclea uel ep ocessing plan eleases adionuclides (Sb, Tc, Cs,
Pu) o he sea a Cap de La Hague, loca ed on he F ench coas o he English
Channel (Fig. 1). Some models ha e been de eloped o simula e he dispe sion o
conse a i e adionuclides ( ha emain in solu ion) in pa o he Channel (Salomon,
Gueguenia , O bi & Ba on, 1988) and in he whole channel plus he sou he n No h
Sea (B e on & Salomon, 1995; Salomon, B e on & Gueguenia , 1995; Janin, 1996).
*Tel.: #34-95-4233669; ax: #34-95-4232644.
E-mail add ess: pe [email protected] (R. Pe iaHnez).
Fig. 1. Compu a ional domain. The s a is La Hague nuclea uel ep ocessing plan . The line is he sec ion
along which adionuclide concen a ions a e ob ained (see Sec ion 3.2). Each uni in he x- and y-axis is
5000 m.
Howe e , hese models wo k wi h esidual ci cula ion (a e aged cu en s) and, as
a consequence, hey can only be used o ob ain he long- e m dispe sion o adionucl-
ides ( empo al scale in he o de o yea s).
Recen ly, Pe iaHnez and Regue a (1999) de eloped a nume ical model o simula e he
idal dispe sion o dissol ed adionuclides in he English Channel. The model sol es
he shallow wa e hyd odynamic equa ions o ob ain he ins an aneous wa e el-
oci ies o e he Channel and, simul aneously, he ad ec ion di!usion dispe sion
equa ion is also sol ed. Model p edic ions we e in good ag eemen wi h obse a ions
in he case o Sb and Tc. The model unde es ima ed he concen a ions in he
case o Cs. This can be due o he ac ha his adionuclide is less conse a i e
han Sb and Tc (IAEA, 1985) and hus pa o he eleased Cs is "xed o he
bo om sedimen . Bo om sedimen s may be a sou ce o adionuclides o he wa e
column when concen a ions in wa e a e low enough o allow he deso p ion eac ion
o domina e abso p ion. Du ing he simula ed pe iod he inpu o Cs was smalle
han some mon hs be o e and hus he Cs con ibu ion om he sedimen could be
highe han he inpu om he sou ce. This e!ec was no conside ed in he model and
could explain he unde es ima ion o he Cs concen a ions.
The objec i e o his wo k consis s o ex ending he p e ious model o non-
conse a i e adionuclides. Thus, he shallow wa e hyd odynamic equa ions a e
sol ed and, simul aneously, he suspended ma e equa ion, which includes ad ec ion
and di!usion o suspended pa icles plus he deposi ion and esuspension p ocesses,
and he equa ions ha go e n he dispe sion o adionuclides in he h ee phases
conside ed (wa e , suspended ma e and bo om sedimen s) a e sol ed oo. The
ans e s o adionuclides be ween he liquid and solid phases a e desc ibed in e ms o
kine ic ans e coe$cien s. The model has been applied o simula e he dispe sion o
Cs: he objec i e is o imp o e he esul s o Pe iaHnez and Regue a (1999). I has
also been applied o s udying he dispe sion o Pu. The dispe sion o his
adionuclide, which has a high a$ni y o he solid phases, in he English Channel has
no been s udied be o e wi h nume ical models. Also, he ansi imes o adionuclides
om Cap de La Hague o he S ai o Do e (see Fig. 1) ha e been calcula ed.
The model is p esen ed in he nex sec ion. The ea e , he esul s a e p esen ed and
discussed.
2. The model
2.1. Hyd odynamics
The wo-dimensional shallow wa e hyd odynamic equa ions a e p esen ed in
Pe iaHnez and Regue a (1999) and will no be epea ed he e. They include he non-
linea e ms, a Co iolis e m and he e!ec o wind. The use o a wo-dimensional
model is jus i"ed due o he dominance o ba o opic o e ba oclinic mechanisms in
he shallow and well-mixed wa e s o he Channel (B e on & Salomon, 1995).
Ins ead o using a cons an alue o he bed ic ion coe$cien , as is usual in mos
models, i s alue has been inc eased om 0.0015 in he wes Channel o 0.0875 in he
Do e S ai o ob ain a be e ag eemen be ween he obse ed and compu ed
cu en s. Indeed, P andle (1975) has ound in p e ious modelling wo k ha ic ion
inc eases in he egion a ound Do e and has no ed ha he alue o he ic ion
coe$cien mus be selec ed wi h ca e in his a ea so as o accoun o such an inc ease
in he bo om ic ion.
2.2. Suspended ma e dynamics
As usual in suspended ma e s udies, we will conside ha only pa icles wi h
a diame e (62.5 lm can emain in he wa e column as suspended ma e , since
la ge pa icles will sink apidly o he bo om and, as a consequence, hei ho izon al
mo emen is negligible (Belde son, 1964; Gu bu , Ke shaw & Du ance, 1987; Cla ke,
1995; Pe iaHnez, Ab il & Ga cmHa-LeoHn, 1996a). Indeed, Eisma (1981) has poin ed ou
ha o all p ac ical pu poses muds can be ega ded as synonymous wi h suspended
ma e . The suspended ma e equa ion is
*(Hm)
* #*(uHm)
*x#*( Hm)
*y"*
*xHKV
*m
*x#*
*yHKW
*m
*y
!wQm1!q
qAB #E q+, (1)
whe e mis he suspended ma e concen a ion in g/m,uand a e he componen s o
he wa e eloci y along he xand ydi ec ions, His he ins an aneous wa e dep h and
KVand KWa e he di!usion coe$cien s along he xand ydi ec ions. The las wo
e ms ep esen deposi ion and e osion o he sedimen . The deposi ion e m has been
based on he concep gi en by Teisson (1991), also used by Cla ke (1995): wQis he
se ling eloci y o suspended pa icles, q"(u# is he cu en magni ude and
qAB is a c i ical deposi ion eloci y. Thus, deposi ion occu s only i q(qAB; o he wise
his e m is se o ze o. Eis he e osion cons an (Nicholson & O'Conno , 1986), Mis
some powe o he wa e eloci y, ypically in he ange 2}5 (P andle, 1997) and gi es
he ac ion o small pa icles (diame e (62.5 lm) in he sedimen . This way, he
o mula de i ed by Fukuda and Lick (1980) and La elle, Moj eld and Bake (1984)
has been used o he e osion a e, al hough he ac o has been included o ake
in o accoun ha only pa icles wi h diame e (62.5 lm can be inco po a ed
as suspended ma e in o he wa e column. E osion occu s only i he wa e eloci y
is la ge han a c i ical e osion eloci y, qAC ; o he wise he e osion e m is se
o ze o.
The se ling eloci y o suspended pa icles, a low concen a ions (smalle
han 1000 ppm), inc eases as he suspended ma e concen a ion inc eases. A
s anda d o mula used o ep esen his p ocess is (Pej up, 1988; Meh a, 1989;
Cla ke, 1995):
wQ"am?

, (2)
whe e aand aha e o be ob ained om measu emen s o om model calib a ion.
Finally, a sou ce e m mus be included in he suspended ma e equa ion along he
coas line. This e m ep esen s he inpu o pa icles om uno!o con inen al wa e s
(Pe iaHnez, 1999).
The solu ion o he suspended ma e equa ion p o ides he suspended ma e
concen a ion a each posi ion in he model domain and a each ins an o ime.
2.3. Radionuclide dispe sion
As commen ed abo e, he model conside s ha adionuclides can be p esen in
h ee phases: solu ion, suspended ma e and ac i e bo om sedimen s (pa icles wi h
a diame e (62.5 lm). I is conside ed ha he ans e o adionuclides om wa e o
he solid phase is go e ned by a coe$cien kand he in e se p ocess by a coe$cien
k, which a e deno ed kine ic ans e coe$cien s. The adso p ion o adionuclides
will depend on he su ace o pa icles pe wa e olume uni . This quan i y has been
deno ed as he exchange su ace (Pe iaHnez, Ab il & Ga cmHa-LeoHn, 1996b; Pe iaHnez
& Ma mHnez-Agui e, 1997a; Pe iaHnez, 1999). Thus
k"s(SK#SQ), (3)
whe e SKand SQa e he exchange su aces o suspended ma e and bo om
sedimen s, espec i ely, and sis a pa ame e wi h he dimensions o a eloci y. I is
deno ed as he exchange eloci y (Pe iaHnez e al., 1996b). As a " s app oach,
assuming sphe ical pa icles and a s ep unc ion o he g ain size dis ibu ion o
pa icles, i can be easily ob ained (Pe iaHnez e al., 1996b) ha
SK"3m
oR, (4)
SQ"3¸ 
RH , (5)
whe e ois he suspended pa icle densi y, Ris he mean adius o suspended ma e
and ac i e sedimen pa icles, ¸is he a e aged mixing dep h ( he dis ance o which
he dissol ed phase pene a es he sedimen ) and is a co ec ion ac o ha akes
in o accoun ha no all he mass o he sedimen is in con ac wi h wa e .
The ans e coe$cien kis conside ed cons an . This desc ip ion o he ans e o
adionuclides be ween he dissol ed and he solid phases has been used success ully in
p e ious modelling s udies (Pe iaHnez, Ab il & Ga cmHa-LeoHn, 1996c; Pe iaHnez
& Ma mHnez-Agui e, 1997aPe iaHnez & Ma mHnez-Agui e, 1997b; Pe iaHnez, 1999).
The equa ion ha gi es he empo al e olu ion o adionuclide concen a ion in he
dissol ed phase, CB(Bq/m), is
*(CBH)
* "(ad #di )!kCBH#kCQmH#kAQ¸oQ ;10!jCBH, (6)
whe e (ad #di ) means ad ec i e plus di!usi e anspo o dissol ed adionuclides,
CQ(Bq/g) and AQ(Bq/g) a e adionuclide concen a ions in suspended ma e and
ac i e sedimen pa icles, espec i ely, he sedimen bulk densi y oQis exp essed in
kg/mand jis he adioac i e decay cons an . The ex e nal sou ce o adionuclides
should be added o his equa ion a he poin s whe e i exis s. The ans e coe$cien
kis gi en by
k"s3m
oR#3¸ 
RH . (7)
The equa ion o he empo al e olu ion o speci"c ac i i y in suspended ma e
pa icles is
*(mCQH)
* "(ad #di )#kCBH#kmCQH#(e os!dep)!jmCQH, (8)
whe e kis gi en by he " s e m o Eq. (7) and he e osion and deposi ion e ms a e
gi en by
e os"E q+AQ, (9)
dep"wQmCQ1!q
qAB . (10)
Again, a sou ce e m should be included in Eq. (8) i he e exis s an ex e nal inpu o
adionuclides "xed o solid pa icles.

The equa ion o he empo al e olu ion o speci"c ac i i y in he ac i e ac ion o
he sedimen is
*AQ
* "k
CBH
¸oQ 10 !kAQ#(dep!e os), (11)
whe e kis gi en by he second e m o Eq. (7) and he deposi ion and e osion e ms
a e gi en by
dep"wQmCQ
¸oQ 1!q
qAB 10 , (12)
e os"Eq+AQ
¸oQ
10 . (13)
2.4. Model conxgu a ion
To sol e he equa ions, a spa ial and empo al disc e iza ion is ca ied ou : he
Channel was di ided in o 3750 g id cells ( o ming a ma ix 75;50). The g id ex ends
om 48.33N o 51.03N and om 4.03W o 1.53E; he g id cell size is
*x"*y"5000 m and he ime-s ep is "xed as * "60 s.
The explici "ni e di!e ence scheme desc ibed in Fla he and Heaps (1975) was used
o sol e he hyd odynamic equa ions. Second-o de accu acy "ni e di!e ence schemes
(Kowalick & Mu y, 1993) ha e been applied o sol e he ad ec ion and di!usion
e ms in he suspended ma e and adionuclide dispe sion equa ions. In he case o
ad ec ion, i consis s o a second o de app oxima ion o he upwind scheme. All he
s abili y condi ions (Kowalick & Mu y, 1993) a e sa is"ed by he spa ial and
empo al esolu ions o he model.
The compu a ional domain is p esen ed in Fig. 1. Along he open wes e n bound-
a y, wa e ele a ions we e speci"ed o each ime-s ep om obse a ions (Howa h
& Pugh, 1983). A adia ion condi ion (Kowalick & Mu y, 1993) is applied along he
no heas e n open bounda y (Do e S ai ). Along he coas , he cu en componen
ha is no mal o he bounda y is se o ze o. Only he wo main idal componen s
(M#S) a e conside ed. Good ag eemen wi h obse a ions is ob ained al hough
only wo componen s a e used, as will be seen.
In he case o he dispe sion equa ions, he e is no #ux o adionuclides o
suspended ma e h ough a closed bounda y. Along open bounda ies, he ollowing
bounda y condi ion is applied:
CG" CG , (14)
whe e CG ep esen s he concen a ion o adionuclides o suspended ma e a he
open bounda y and CG  ep esen s he concen a ion jus inside he compu a ional
domain. This condi ion has been p e iously applied in Pe iaHnez (1998a,b,1999) and
Pe iaHnez and Regue a (1999).
A FORTRAN code was de eloped o sol e he equa ions in ol ed in he model and
was implemen ed on a HP SPP 2000 X-Class compu e .
3. Resul s and discussion
3.1. Hyd odynamics and suspended ma e dynamics
The magni ude and di ec ion o he compu ed and obse ed cu en s ha e been
compa ed o 12 poin s in he channel o sp ing, medium and neap ides. Obse ed
and compu ed Mco ange cha s ha e also been compa ed. Resul s a e, in gene al, in
good ag eemen wi h obse a ions. De ails can be seen in Pe iaHnez and Regue a
(1999) and will no be epea ed he e.
The ollowing pa ame e s a e used o sol e he suspended ma e equa ion: di!u-
sion coe$cien s a e "xed as KV"KW"51 m/s and i is assumed ha "0.9 in Eq.
(14) (see Pe iaHnez & Regue a, 1999, o de ails).
One o he main di$cul ies in sedimen anspo modelling is he me hod o
ob aining he e osion and deposi ion h esholds and he e osion cons an . These
pa ame e s a e speci"c o each si e and hey a e o en selec ed a e a model
calib a ion p ocess. Thus hey a e selec ed in such a way ha hey p oduce he bes "
o he model esul s wi h he obse a ions. O cou se, hese pa ame e s mus be
physically ealis ic.
A e some model uns, he e osion cons an has been "xed as
E"2.2;10  gm
s. Uncon en ional uni s a e due o he ac ha Ealso ac s
as a scaling ac o : he e osion a e mus be exp essed in g m s  and i has been
assumed ha M"3.5 (as commen ed abo e, Mis ypically in he ange 2}5).
I has been selec ed, also a e a model calib a ion, ha a"5.7;10  and
a"1.6 in Eq. (2) is wQis exp essed in m/s and min g/m. Since suspended ma e
concen a ions in mos o he Channel a e o he o de o 1 g/m, se ling eloci ies o
he o de o 10  m/s a e ob ained, which is he o de o magni ude ha can be ound
in he li e a u e (P andle, Jago, Jones, Pu die & Tappin, 1993; P andle, 1997)
o he se ling eloci y o cohesi e sedimen s. Indeed, i can be ound in he li e a u e
ha a anges be ween 0.5 and 2 (Pej up, 1988; Eisma, 1993). A alue o 1.6 was
ob ained by Meh a (1989) and was also used in p e ious modelling wo k (Pe iaHnez e
al., 1996a). I can also be seen ha a anges om 10  (Cla ke, 1995) o 10  (Meh a,
1989).
The c i ical deposi ion and e osion eloci ies ha e been aken om li e a u e da a.
A alue o 0.18 m/s has been used o he c i ical deposi ion eloci y (Eisma, 1993;
Pe iaHnez e al., 1996a) and a ypical alue o 0.21 m/s was assumed o he c i ical
e osion eloci y (Hea he shaw, 1981; Pugh, 1987). Indeed, i has al eady been poin ed
ou (Cla ke, 1995) ha due o he small size o cohesi e sedimen s hey end o be slow
alling and hus he e is some lag be ween he end o deposi ion and e osion begin-
ning. Du ing he in e media e ime in e al, he sedimen suspension is main ained by
wa e u bulence wi h no appa en e!ec on e osion o deposi ion. Thus, usually,
qAB(qAC (Cla ke, 1995).
The suspended ma e sou ce- e m om con inen al uno!was selec ed a e
a calib a ion p ocess. Along he F ench coas , i was conside ed ha a mean supply o
0.027;10kg/yea occu s. Along he B i ish coas i was "xed as 0.35;10kg/yea .
O cou se, uno!will no be cons an h ough he yea , bu hese alues may be
conside ed as annual a e ages ha gi e a mean suspended ma e dis ibu ion o e
he Channel ha is, in gene al, in ag eemen wi h he obse ed ea u es o suspended
ma e concen a ions. Fo compa ison, suspended ma e supply o he No h Sea
om di!e en i e s anges om 10(Ems i e ) o 10kg/yea (Thames plus Hum-
be , Rhine plus Meuse i e s), and he supply om coas al e osion (Eas Anglia plus
Holde ness coas s) is 0.7;10kg/yea (Eisma, 1981).
The mean suspended ma e concen a ions compu ed by he model a e p esen ed
in Fig. 2a and a e, in gene al, in ag eemen wi h hose gi en by Eisma and Kal (1987),
shown in Fig. 2b, ob ained om obse a ions. Suspended ma e concen a ions a e
somewha highe along he English coas (some 5 ppm) han along he F ench coas
(some 3 ppm). Also, suspended ma e concen a ions o 5 ppm a e ob ained o he
no h o La Hague. Un o una ely, he e a e no measu emen s in he wes e n pa o
he Channel ( om La Hague o he wes ). Again, he use o a wo-dimensional model
is jus i"ed since Eisma and Kal (1987) ha e ound ha suspended ma e concen a-
ions in he bo om wa e a e a he simila o hose in he su ace wa e . Measu ed
(B ylinski, Dupon & Ben ley, 1984) and compu ed suspended ma e concen a ions
a Do e S ai ha e also been compa ed: measu ed concen a ions 10 and 15 km
om Cap G is-Nez a e 7 and 5 ppm, espec i ely; he co esponding mean compu ed
alues a e 5.5 and 5.2 ppm.
The suspended ma e #ux h ough he Do e S ai in o he No h Sea has been
compu ed as 1.3;10 g/yea . Eisma (1981) es ima ed a #ux o 1.0;10 g/yea ,
al hough he alue ob ained by Van Alphen (1990), 1.7;10 g/yea , is conside ed he
mos accu a e es ima e (Eisma, 1990). Thus he compu ed #ux is in a he good
ag eemen wi h p e ious es ima es.
3.2. Radionuclide dispe sion
The compu ed dis ibu ion o a conse a i e adionuclide (Tc) in he Channel is
p esen ed in Fig. 3 as an example ( om Pe iaHnez & Regue a, 1999). In he case o
a conse a i e adionuclide, s"k"0. I can be seen ha concen a ions a e
highe on he F ench side o he Channel han on he English side. The banded
s uc u e, showing dec easing ac i i ies o! he F ench coas , has been obse ed by
Gueguenia , He mann, Ke shaw, Bailly du Bois and Ba on (1996) and is also
appa en he e. Wa e s o he No man B e on Gul a e a!ec ed by he discha ges om
La Hague due o he exis ence o gy es a ound he Channel Islands (O bi & Salomon,
1988). The Do e S ai ac s as a bo leneck and hus g adien s a e enhanced in his
a ea. This e!ec was also ound by B e on and Salomon (1995).
Some pa ame e s mus be known o simula e he dispe sion o non-conse a i e
adionuclides. The mixing dep h ¸was aken as 0.1 m, ollowing p e ious modelling
wo ks (Pe iaHnez e al., 1996c; Pe iaHnez & Ma mHnez-Agui e, 1997a; Pe iaHnez
& Ma mHnez-Agui e, 1997b; Pe iaHnez, 1999). The densi y o suspended ma e pa -
icles is aken as 2600 kg/m, which is he es ablished alue o soil pa icle densi y
(Ba e , Ga dne & Ga dne , 1972). The mean adius o suspended ma e and ac i e
sedimen pa icles has been selec ed as R"15 lm, since only pa icles wi h a
diame e (62.5 lm can be p esen in he wa e column as suspended ma e , as
Fig. 2. (a) A e age suspended ma e concen a ions (ppm) calcula ed by he model. (b) Suspended ma e
concen a ions (ppm) om obse a ions o Eisma and Kal (1987).
local se ies, ha is de"ned as
h(q)"1C( )C( #q)2
(1C( )2(1C( )2, (15)
whe e his he c oss-co ela ion unc ion, Cand Ca e concen a ions a La Hague
and he cen al Do e S ai , espec i ely, and 12means ime-a e aging. The max-
imum o he c oss-co ela ion unc ion occu s o a gi en alue o q, which is
a s a is ical e alua ion o he ansi ime be ween he wo poin s.
I wind is no included in he model, he compu ed ansi ime due o ides is 4.1
mon hs. Howe e , i he a e age wind, sou hwes 6 m/s (B e on & Salomon, 1995), is
included, he ansi ime is educed o 3.0 mon hs. This alue is in close ag eemen
wi h he p e ious calcula ion o Salomon e al. (1995) wi h a long- e m model. This
dec ease is due o he ac ha a e age wind condi ions p oduce an inc ease in he
esidual d i o he No h Sea ela i e o he idal d i . This compu ed alue is also in
ag eemen wi h he p e ious expe imen al es ima ions o Gueguenia , Bailly du Bois,
Gandon, Salomon, Ba on and Leon (1994), who ound a ansi ime be ween 2 and
8 mon hs. A ansi ime o 3 mon hs implies an a e age eloci y o adionuclides
discha ged a La Hague o 83 km/mon h. O cou se, he ansi ime depends upon
wea he condi ions. Thus, in he case o a cons an eas e ly wind (speed 6 m/s), he
compu ed ansi ime is inc eased o 7.1 mon hs.
O cou se, ansi imes a e a!ec ed by he eal changing and non-homogeneous
wind "elds. Howe e , ou calcula ions allow easonable es ima ions o he o de o
magni ude o he ansi imes. Indeed, 3 and 7.1 mon hs ha e been ob ained o
a e age wind and eas e ly wind condi ions, espec i ely. The es ima e by Gueguenia
e al. (1994), 2}8 mon hs, ep esen s he eal mix u e o di!e en wind condi ions since
i has been deduced om obse a ions. The calcula ions om his modelling s udy a e
in he ange o a ia ion gi en by Gueguenia e al. (1994) and hus hey can be
conside ed as ealis ic alues.
These ansi imes ha e been compu ed o dissol ed adionuclides, conside ing
ha s"k"0. Howe e , i has been obse ed ha i adionuclides ha e a high
a$ni y o he solid phase he ansi ime inc eases signi"can ly. Thus, in he case o
Pu, he ansi ime o Do e is 7.0 mon hs (in a e age wind condi ions) a alue
ha is highe han he ansi ime o a conse a i e adionuclide by a ac o 2.3. The
ansi ime o plu onium had no been calcula ed p e iously.
4. Conclusions
A model o simula e he dispe sion o non-conse a i e adionuclides in he English
Channel has been de eloped. The model sol es he shallow wa e hyd odynamic
equa ions, he suspended ma e equa ion and he equa ions o he empo al e olu-
ion o ac i i y concen a ions in he h ee phases: wa e , suspended ma e and ac i e
bo om sedimen s. The ans e o adionuclides be ween wa e and he solid phases
has been desc ibed in e ms o kine ic ans e coe$cien s.

The suspended ma e sub-model has been es ed by compa ing obse ed and
compu ed suspended ma e dis ibu ions. Also, he suspended ma e #ux o he
No h Sea h ough he Do e S ai has been calcula ed. The compu ed #ux is in
ag eemen wi h p e ious es ima es.
The model has been applied o simula e he dispe sion o Cs and Pu
eleased om he nuclea uel ep ocessing plan a La Hague. In he case o Cs, he
esul s ob ained wi h an ea lie model, in which his adionuclide was conside ed o
emain in solu ion, ha e been imp o ed. This sugges s he impo ance o con-
amina ed sedimen s as a sou ce o adionuclides o he dissol ed phase. The model
gi es a ealis ic es ima ion o Pu ac i i y le els measu ed along he Channel in
wa e , suspended ma e and bo om sedimen s. Also, di!e en dis ibu ion maps ha e
been ob ained o conse a i e and non-conse a i e adionuclides, e ealing he
in#uence o suspended ma e on Cs and Pu dispe sion. Compu ed kB alues o
plu onium a e consis en wi h hose ecommended by IAEA (1985).
T ansi imes o adionuclides om La Hague o he cen al Do e S ai ha e also
been calcula ed using he c oss-co ela ion unc ion me hod. A ansi ime o
3 mon hs has been calcula ed o dissol ed adionuclides wi h a e age wind condi-
ions, in close ag eemen wi h p e ious calcula ions. The ansi ime inc eases o
4 mon hs i only ides a e conside ed. On he o he hand, a alue o 7 mon hs is
ob ained o eas e ly winds. In he case o a adionuclide wi h a high a$ni y o he
solid phases, like Pu, a ansi ime signi"can ly longe han ha o a dissol ed
adionuclide has been ob ained.
Acknowledgemen s
This wo k was pa ially suppo ed by ENRESA, EU con ac FI4PCT960046 and
Spanish CICYT con ac AMB97-1720CE.
Re e ences
Ba e , L. D., Ga dne , W. H., & Ga dne , W. R. (1972). Soil physics. New Yo k: Wiley.
Bailly du Bois, P., Roze , M., Tho al, K., & Salomon, J. C. (1997). Imp o ing knowledge o wa e }mass
ci cula ion in he English Channel using adioac i e ace s. Radiop o ec ion Colloques,32(C2), 63}69.
Belde son, R. H. (1964). Holocene sedimen a ion in he wes e n hal o he I ish Sea. Ma ine Geology,2,
147}163.
Bous , D., Mi chell, P. I., Ga cia, K., Cond en, O., Leon-Vin o, L., & Lecle c, G. (1996). A compa a i e
s udy o he specia ion and beha iou o plu onium in he ma ine en i onmen o wo ep ocessing
plan s. Radiochimica Ac a,74, 203}210.
B e on, M., & Salomon, J. C. (1995). A 2D long e m ad ec ion dispe sion model o he Channel and
Sou he n No h Sea. Pa A: Valida ion h ough compa ison wi h a i"cial adionuclides. Jou nal o
Ma ine Sys ems,6, 495}513.
B ylinski, J. M., Dupon , J., & Ben ley, D. (1984). Condi ions hyd o-biologiques au la ge du Cap G is-Nez
(F ance): P emie s esul a s. Oceanologica Ac a,7, 315}322.
Cla ke, S. (1995). Ad ec i e/di!usion p ocesses in he Fi h o Fo h. Ph.D. hesis. Uni e si y o Wales,
Bango .
Fla he , R. A., & Heaps, N. S. (1975). Tidal compu a ions o Mo ecambe Bay. Geophysical Jou nal o he
Royal As onomical Socie y,42, 489}517.
Eisma, D. (1981). Supply and deposi ion o suspended ma e in The No h Sea. Spec Publs in Ass Sedimen ,
5, 415}428.
Eisma, D. (1990). T anspo and deposi ion o suspended ma e in he No h Sea and he ela ion o coas al
sil a ion, pollu ion and bo om auna dis ibu ion. Aqua ic Sciences,3, 181}216.
Eisma, D. (1993). Suspended ma e in he aqua ic en i onmen . Be lin: Sp inge .
Eisma, D., & Kal , J. (1987). Dispe sal concen a ion and deposi ion o suspended ma e in he No h Sea.
Jou nal o he Geological Socie y, London,144, 161}178.
Fukuda, M. K., & Lick, W. (1980). The en ainmen o cohesi e sedimen s in esh wa e . Jou nal o
Geophysical Resea ch,85(C5), 2813}2824.
Gueguenia , P., Bailly du Bois, P., Gandon, R., Salomon, J. C., Ba on, Y., & Leon, R. (1994). Spa ial and
empo al dis ibu ion (1987}91) o Sb used o ace pa hways and ansi imes o wa e s en e ing he
No h Sea om he English Channel. Es ua ine, Coas al and Shel Science,39,59}74.
Gueguenia , P., He mann, J., Ke shaw, P., Bailly du Bois, P., & Ba on, Y. (1996). A i"cial adioac i i y in
he English Channel and he No h Sea. In P. Gueguenia , P. Ge main, & H. Me i ie , Radionuclides in
he oceans (pp. 121}154). Les Ulis: Les Edi ions de Physique.
Gu bu , P. A., Ke shaw, P. J., & Du ance, J. A. (1987). Modelling he dis ibu ion o soluble and pa icle
adso bed adionuclides in he I ish Sea. In J. C. Gua y, P. Gueguenia , & R. J. Pen ea h, Radionuclides:
A ool o oceanog aphy (pp. 395}407). Else ie : Ams e dam.
Hea he shaw, A. D. (1981). Compa ison o measu ed and p edic ed sedimen anspo a es in idal
cu en s. Ma ine Geology,42,75}104.
He mann, J., Ke shaw, P., Bailly du Bois, P., & Gueguenia , P. (1995). The dis ibu ion o a i"cial
adionuclides in he English Channel, sou he n No h Sea, Skage a and Ka ega , 1990}1993. Jou nal
o Ma ine Sys ems,6, 427}456.
Howa h, M. J., & Pugh, D. T. (1983). Obse a ions o ides o e he con inen al shel o no hwes Eu ope.
In B. Johns, Physical oceanog aphy o coas al and shel seas (pp. 135}185). Ams e dam: Else ie .
IAEA (1985). Sedimen kBand concen a ion ac o s o adionuclides in he ma ine en i onmen . Technical
Repo Se ies 247, Vienna.
Janin, J.M. (1996). Modelisa ion nume ique de la dilu ion des adionuclides en Manche. Pos e p esen ed in
he in e na ional symposium on adionuclides in he oceans. Che bou g, F ance.
Ke shaw, P., & Young, A. (1988). Sca enging o Th in he eas e n I ish Sea. Jou nal o En i onmen al
Radioac i i y,6,1}23.
Kowalick, Z., & Mu y, T. S. (1993). Nume ical modelling o ocean dynamics. Singapo e: Wo ld Scien i"c.
La elle, J. W., Moj eld, H. O., & Bake , E. T. (1984). An in si u e osion a e o a "ne g ain ma ine sedimen .
Jou nal o Geophysical Resea ch,89(C4), 6543}6552.
McKay, W. A., & Pa enden, N. J. (1993). The beha iou o plu onium and ame icium in he sho eline
wa e s o he I ish Sea: A e iew o Ha well s udies in he 1980s. Jou nal o En i onmen al Radioac i i y,
18,99}132.
Meh a, A. J. (1989). On es ua ine cohesi e sedimen suspension beha iou . Jou nal o Geophysical Resea ch,
94(C10), 14303}14314.
Nicholson, J., & O'Conno , B. A. (1986). Cohesi e sedimen anspo model. Jou nal o Hyd aulic
Enginee ing,112, 621}640.
Ny!ele , U. P., Li, Y. H., & San schi, P. H. (1984). A kine ic app oach o desc ibe ace elemen dis ibu ion
be ween pa icles and solu ion in na u al aqua ic sys ems. Geochimica Cosmochimica Ac a,48,
1513}1522.
O bi, A., & Salomon, J. C. (1988). Dynamique de ma ee dans le gol e No mand-B e on. Oceanologica Ac a,
11,55}64.
Pej up, M. (1988). Suspended sedimen anspo ac oss a idal #a . Ma ine Geology,82, 187}198.
Pe iaHnez, R. (1998a). Th ee dimensional modelling o he ide induced dispe sion o adionuclides in he sea.
Jou nal o En i onmen al Radioac i i y,40, 215}237.
Pe iaHnez, R. (1998b). A h ee dimensional pcoo dina e model o simula e he dispe sion o adionuclides in
he ma ine en i onmen . Applica ion o he I ish Sea. Ecological Modelling,114,59}70.
Pe iaHnez, R. (1998c). Modelling he dis ibu ion o adionuclides in deep ocean wa e columns. Applica ion
o H, Cs and Pu. Jou nal o En i onmen al Radioac i i y,38, 173}194.
Pe iaHnez, R. (1999). Th ee dimensional modelling o he idal dispe sion o non conse a i e adionuclides
in he ma ine en i onmen . Applica ion o Pu dispe sion in he eas e n I ish Sea. Jou nal o
Ma ine Sys ems,22,37}51.
Pe iaHnez, R., & Ma mHnez-Agui e, A. (1997a). U and Th concen a ions in an es ua y a!ec ed by phospha e
e ilize p ocessing: expe imen al esul s and a modelling s udy. Jou nal o En i onmen al Radioac i i y,
35, 281}304.
Pe iaHnez, R., & Ma mHnez-Agui e, A. (1997b). A six phases model o simula e he con amina ion by non
conse a i e adionuclides o sedimen s, soils and plan s in a ma sh a ea. Applica ion o he Odiel
Ma sh in sou hwes Spain. Jou nal o En i onmen al Radioac i i y,37,29}54.
Pe iaHnez, R., & Regue a, J. (1999). A nume ical model o simula e he idal dispe sion o adionuclides in he
English Channel. Jou nal o En i onmen al Radioac i i y,43,51}64.
Pe iaHnez, R., Ab il, J. M., & Ga cmHa-LeoHn, M. (1996a). Modelling he suspended ma e dis ibu ion in an
es ua ine sys em. Applica ion o he Odiel Ri e in sou hwes Spain. Ecological Modelling,87, 169}179.
Pe iaHnez, R., Ab il, J. M., & Ga cmHa-LeoHn, M. (1996b). Modelling he dispe sion o non conse a i e
adionuclides in idal wa e s. Pa 1: Concep ual and ma hema ical model. Jou nal o En i onmen al
Radioac i i y,31, 127}141.
Pe iaHnez, R., Ab il, J. M., & Ga cmHa-LeoHn, M. (1996c). Modelling he dispe sion o non conse a i e
adionuclides in idal wa e s. Pa 2: Applica ion o Ra dispe sion in an es ua ine sys em. Jou nal o
En i onmen al Radioac i i y,31, 253}272.
P andle, D. (1975). S o m su ges in he sou he n No h Sea and Ri e Thames. P oceedings o he Royal
Socie y o London A,344, 509}539.
P andle, D. (1997). Tidal cha ac e is ics o suspended sedimen s concen a ions. Jou nal o Hyd aulic
Enginee ing,123, 341}350.
P andle, D., Jago, C. F., Jones, S. E., Pu die, D. A., & Tappin, A. (1993). The in#uence o ho izon al
ci cula ion in he supply and dis ibu ion o ace s. Philosophical T ansac ions o he Royal Socie y,
London A,343, 405}421.
Pugh, D. T. (1987). Tides, su ges and mean sea le el. Chiches e : Wiley.
Salomon, J. C., Gueguenia , P., O bi, A., & Ba on, Y. (1988). A lag angian model o long e m idally
induced anspo and mixing. Ve i"ca ion by a i"cial adionuclide concen a ions. In J. C. Gua y, P.
Gueguenia , & R. J. Pen ea h, Radionuclides: a ool o oceanog aphy (pp. 384}394). Ams e dam:
Else ie .
Salomon, J. C., B e on, M., & Gueguenia , P. (1995). A 2D long e m ad ec ion dispe sion model o he
Channel and Sou he n No h Sea. Pa B: T ansi ime and ans e unc ion om La Hague. Jou nal o
Ma ine Sys ems,6, 515}527.
Teisson, C. (1991). Cohesi e suspended sedimen anspo : easibili y and limi a ions o nume ical model-
ling. Jou nal o Hyd aulic Resea ch,29, 755}769.
Van Alphen, J. (1990). A mud balance o he Belgian-Du ch coas al wa e s be ween 1969 and 1986.
Ne he lands Jou nal o Sea Resea ch,25,10}18.