Abs ac Ca bonyl sul ide (OCS) is used o quan i y he ca bon cap u e po en ial o he biosphe e because
o i s di ec co ela ion wi h CO2 up ake du ing pho osyn hesis. Howe e , o cons ain he u ban biosphe e
signal, i is necessa y o e alua e po en ial an h opogenic sou ces. We conduc ed wo sampling campaigns in
he Me opoli an A ea o Ba celona (AMB), Spain, du ing May ( ull COVID lockdown) and Oc obe 2020
o measu e he spa ial dis ibu ion and a iabili y o OCS in ou u ban land uses as ollows: buil , u ban
o es , u ban pa k, and pe i-u ban ag icul u e. The OCS backg ound le els de e mined a Tibidabo (442m
asl) we e app oxima ely 484±20 pp and 407±8 pp o May and Oc obe 2020, espec i ely, and ag eed
wi h o he seasonal su eys conduc ed in Eu ope du ing ha same pe iod. A e aged emissions we e in he
ange o +Δ12±40 pp o he ci y and we e +Δ9.4±40.9 pp o u ban +Δ22.1±48 pp o u ban g een
+Δ20.7±42.9 pp o ag icul u al and −Δ4.8±19.6 pp o o es . The u ban alues anged om neu al o
abo e backg ound, sugges ing nea by an h opogenic and ma ine emissions such as +Δ150 pp in Mon juic,
which is downwind o Ba celona's ha bo . Du ing he c op-g owing season in May, he ag icul u al a eas
consis en ly showed alues below he backg ound (up o −Δ76 pp in Ga à, up ake) a 7:00 UTC when he land
b eezes we e dominan , while la e in he mo ning, when he sea b eeze a e de eloped, he plan sink is masked
by he anspo o ma ine emissions. U ban o es s loca ed no h o Tibidabo showed OCS alues up o −Δ70
pp , sugges ing signi ican up ake by u ban o es s. We conclude ha de e mining he u ban biosphe e signal
using OCS as a ace is mo e complex han expec ed because he ma ine and an h opogenic emissions om he
po s ongly impac he spa ial- empo al dis ibu ion o OCS.
Plain Language Summa y Ca bonyl sul ide (OCS) is a gas used o quan i y he ca bon cap u e
po en ial o he biosphe e. Howe e , ci ies ha e po en ial an h opogenic OCS sou ces ha can impede a clea
cons ain o he u ban biosphe e signal. We conduc ed wo sampling campaigns in he Me opoli an A ea o
Ba celona (AMB), Spain, du ing May ( ull COVID lockdown) and Oc obe 2020 o measu e he a iabili y o
OCS in ou u ban land uses: buil , u ban o es , u ban pa k, and pe i-u ban ag icul u e. The OCS backg ound
le els de e mined a Tibidabo (we e app oxima ely 484±20 pp and 407±8 pp o May and Oc obe 2020,
espec i ely). The u ban buil si es we e ypically abo e backg ound, sugges ing nea by an h opogenic and
ma ine emissions. Du ing he c op-g owing season, he ag icul u al a eas consis en ly showed alues below
he backg ound sugges ing plan OCS up ake when he land b eezes we e dominan , while la e in he mo ning,
when he sea b eezes de eloped, he in luence o ma ine emission esul ed in alues abo e he backg ound.
U ban o es s showed signi ican OCS up ake. We conclude ha de e mining he u ban biosphe e signal using
OCS as a ace is mo e complex han expec ed because he ma ine and an h opogenic emissions om he po
s ongly impac he spa ial- empo al dis ibu ion o OCS.
ESTRUCH ETAL.
© 2023. The Au ho s.
This is an open access a icle unde
he e ms o he C ea i e Commons
A ibu ion-NonComme cial-NoDe i s
License, which pe mi s use and
dis ibu ion in any medium, p o ided he
o iginal wo k is p ope ly ci ed, he use is
non-comme cial and no modi ica ions o
adap a ions a e made.
Explo ing he In luence o Land Use on he U ban Ca bonyl
Sul ide Budge : A Case S udy o he Me opoli an A ea o
Ba celona
Ca me Es uch1,2 , Sau eu Bel iso3 , Alba Badia1 , Ve onica Vidal1,4 , Roge Cu coll1,5,
Mi eia Udina6 , Claudia G ossi5 , Josep-An on Mo guí1,7, Rica d Segu a1, Se gi Ven u a1,
Yolanda Sola6 , and Ga a Villalba1,8
1Ins i u de Ciència i Tecnologia Ambien als (ICTA), Uni e si a Au ònoma de Ba celona, Ba celona, Spain, 2Eu eca , Cen e
Tecnològic de Ca alunya, Clima e Change Resea ch Depa men , Ampos a, Spain, 3Labo a oi e des Sciences du Clima e de
l’En i onnemen , LSCE, Pa is-Saclay, F ance, 4Depa men o Compu e A chi ec u e and Ope a i e Sys ems, Uni e si a
Au ònoma de Ba celona (UAB), Ba celona, Spain, 5Ins i u e o Ene gy Techniques, Uni e si a Poli ècnica de Ca alunya,
Ba celona, Spain, 6Depa amen de Física Aplicada–Me eo ologia, Uni e si a de Ba celona, Ba celona, Spain, 7Depa amen
de Biologia E olu i a, Ecologia i Ciències Ambien als, Uni e si a de Ba celona, Ba celona, Spain, 8Depa men o
Chemical, Biological and En i onmen al Enginee ing, Uni e si a Au ònoma de Ba celona (UAB), Ba celona, Spain
Key Poin s:
• The u ban buil si es we e ypically
abo e backg ound, sugges ing nea by
an h opogenic and ma ine emissions
• Ag icul u al a eas consis en ly
showed alues below he backg ound
sugges ing plan Ca bonyl sul ide
(OCS) up ake when he land b eezes
we e dominan
• De e mining u ban biosphe e signal
using OCS is complex because ma ine
and an h opogenic emissions impac
i s spa ial- empo al dis ibu ion
Suppo ing In o ma ion:
Suppo ing In o ma ion may be ound in
he online e sion o his a icle.
Co espondence o:
G. Villalba,
[email p o ec ed]
Ci a ion:
Es uch, C., Bel iso, S., Badia, A.,
Vidal, V., Cu coll, R., Udina, M.,
e al. (2023). Explo ing he in luence
o land use on he u ban ca bonyl
sul ide budge : A case s udy o he
me opoli an a ea o Ba celona. Jou nal
o Geophysical Resea ch: A mosphe es,
128, e2023JD039497. h ps://doi.
o g/10.1029/2023JD039497
Recei ed 23 JUN 2023
Accep ed 26 NOV 2023
Au ho Con ibu ions:
Concep ualiza ion: Ca me Es uch,
Sau eu Bel iso, Ga a Villalba
Da a cu a ion: Josep-An on Mo guí,
Ga a Villalba
Funding acquisi ion: Ga a Villalba
In es iga ion: Ca me Es uch, Sau eu
Bel iso, Alba Badia, Ve onica Vidal,
Roge Cu coll, Mi eia Udina, Claudia
G ossi, Josep-An on Mo guí, Rica d
Segu a, Se gi Ven u a, Yolanda Sola,
Ga a Villalba
10.1029/2023JD039497
RESEARCH ARTICLE
1 o 21
Jou nal o Geophysical Resea ch: A mosphe es
ESTRUCH ETAL.
10.1029/2023JD039497
2 o 21
1. In oduc ion
Ca bonyl sul ide (he ea e OCS) is an a mosphe ic ace gas used o s udy he CO2 seques a ion capaci y o
ege a ion (Remaud e al.,2023; Whelan e al.,2018) because i is des oyed h ough pho osyn hesis in he
same way as CO2 wi hou p oducing any emissions h ough espi a ion (P o oschill-K ebs & Kesselmeie ,1992).
Acco ding o he la es es ima ions, e es ial ege a ion and soils akes up om 530 o 670GgSy
−1 (Remaud
e al.,2023) he la ges na u al sink o a mosphe ic OCS. Su ace oceans a e he mos impo an na u al sou ce
o OCS con ibu ing di ec ly and indi ec ly, ha is, om he a mosphe ic oxida ion o ca bon disul ide and
dime hylsul ide, be ween 270 and 320 (Lenna z e al.,2017). The global budge o OCS is pa ially balanced by
an h opogenic emissions, which a e es ima ed o be 400GgSy
−1 (Campbell e al.,2015). Te es ial sinks and
na u al and an h opogenic sou ces a e gene ally sepa a ed in space and ime, and exhibi dis inc seasonal and
spa ial a ia ions (Remaud e al.,2023). Hence, land use changes a e expec ed o ha e a signi ican in luence on
he global budge o a mosphe ic OCS.
Al hough many s udies ha e e alua ed OCS budge s a egional and global scales, li le is known abou OCS
sou ces and sinks a he u ban scale. Unde s anding land use in luence on he budge o a mosphe ic OCS a he
u ban scale is much mo e di icul o es ablish because sou ces and sinks can coexis o e small dis ances. Fo
example, Commane e al.(2013) indi ec ly assessed OCS sou ces along he highly u banized coas line o Cali-
o nia whe e hey de ec ed biogenic OCS emissions om we lands in Sac amen o Bay, an h opogenic emissions
nea a s o age acili y and e ine y in he po o Los Angeles, and OCS loss du ing he day ime due o pho osyn-
he ic up ake in ege a ed a eas. Besis e al.(2021) also conduc ed a s udy ha iden i ied po en ial sou ces o OCS
a he ci y le el and obse ed he spa ial he e ogenei y in luencing OCS mixing a ios o Thessaloniki, G eece,
and iden i ied h ee loca ions whe e OCS concen a ions in inne -ci y a eas we e h ee imes highe han hose a
he sea on and whe e he was e managemen sys em was a sou ce o OCS. Acco ding o Zumkeh e al.(2018),
an h opogenic emissions o OCS due o aluminum, esiden ial and indus ial coal, TiO2, sol en s and i es in he
Me opoli an A ea o Ba celona con ibu e o 65% (178MgSy
−1) o he o al emissions in he Ca alonia egion
(276MgSy
−1, see Figu e S1 in Suppo ing In o ma ionS1). Howe e , a s udy by Bel iso e al.(2023) demon-
s a ed ha he Zumkeh e al.(2018) in en o y o F ance la gely o e es ima ed OCS an h opogenic emissions
by app oxima ely an o de o magni ude. A p e ious s udy by he same g oup used di ec measu emen s om a
sampling owe and compu ed back ajec o ies wi h he FLEXPART anspo model o ind ha he Pa is subu -
ban a ea was ac ing mainly as an OCS sink (Bel iso e al.,2022). Ano he u ban s udy o Innsb uck, Aus ia,
measu ed he hou ly exchange a es o OCS wi h an eddy co a iance owe and showed ha u ban OCS emissions
a e signi ican ly highe a noon du ing weekdays han du ing weekends, which sugges s ha oad a ic could be
a ne sou ce o OCS emissions.
The use o OCS as a pho osyn hesis ace has ecen ly ecei ed a en ion as a po en ial me hod o de e mine
he con ibu ion o he u ban biosphe e owa d educing he ca bon oo p in (Mallik e al., 2016; Villalba
e al.,2021). Mallik e al.(2016) epo ed he i s eco ds o annual OCS mixing a ios o e an Indian ci y,
which we e subs an ially highe han oposphe ic alues, and hose alues changed seasonally wi h an obse ed
86% d op in OCS in No embe when ai masses changed om oceanic o con inen al. Villalba e al.(2021)
also obse ed changes in OCS ela ed o ai masses coming om ma ine and con inen al sou ces. This wo k
highligh ed he di icul y o cha ac e izing he biosphe e signal when he inpu masses a e con inen al due o he
he e ogenei y, hey p esen in OCS con en compa ed o s able ma ine condi ions.
This s udy aims o con ibu e o he unde s anding o he dynamics o OCS budge s in u ban a eas by explo ing
o wha ex en di e en u ban land uses can ha e an impac on a mosphe ic OCS mixing a ios in e ms o he
pa icula u ban clima e, geog aphy, and opog aphy a ec ing he anspo and mixing o OCS sou ces and
sinks. Mo i a ed by he p esen lack o in si u OCS con inuous measu emen s o p ope ly in e su ace OCS
(Remaud e al.,2022), we conduc ed an OCS measu emen campaign in he Me opoli an A ea o Ba celona
(AMB). We selec ed ou di e en land uses as sampling si es as ollows: o es , ag icul u al land, u ban pa ks
and impe ious a eas wi h he in en ion o unde s anding he ole o u ban land use on he OCS budge in a
coas al ci y whe e we would expec a mixing a io g adien be ween he sea (whe e OCS is p oduced) and he land
(OCS cap u e in ege a ed a eas). The selec ion o sampling a eas was based on hei po en ial as a sink o OCS
based on pho osyn he ic a es (Campbell e al.,2017; Villalba e al.,2021; Yang e al.,2018) o a local sou ce
based on an h opogenic in en o ies (Lee & B imblecombe,2016; Yan e al.,2019; Zumkeh e al.,2018). The
OCS sampling campaigns we e unded and conduc ed wi hin he p ojec In eg a ed Sys em Analysis o U ban
Me hodology: Ca me Es uch, Ga a
Villalba
P ojec Adminis a ion: Ga a Villalba
Resou ces: Ga a Villalba
Supe ision: Sau eu Bel iso, Ga a
Villalba
Visualiza ion: Ca me Es uch, Ve onica
Vidal, Roge Cu coll, Ga a Villalba
W i ing – o iginal d a : Ca me Es uch,
Ga a Villalba
W i ing – e iew & edi ing: Ca me
Es uch, Sau eu Bel iso, Alba Badia,
Roge Cu coll, Mi eia Udina, Claudia
G ossi, Ga a Villalba
21698996, 2023, 24, Downloaded om h ps://agupubs.onlinelib a y.wiley.com/doi/10.1029/2023JD039497 by Readcube (Lab i a Inc.), Wiley Online Lib a y on [10/01/2024]. See he Te ms and Condi ions (h ps://onlinelib a y.wiley.com/ e ms-and-condi ions) on Wiley Online Lib a y o ules o use; OA a icles a e go e ned by he applicable C ea i e Commons License
Jou nal o Geophysical Resea ch: A mosphe es
ESTRUCH ETAL.
10.1029/2023JD039497
3 o 21
Vege a ion and Ag icul u e (URBAG), which was inanced by he Eu opean Resea ch Council wi h he o e all
aim o unde s anding he in luence o g een in as uc u es on he u ban CO2 budge . As a coas al ci y wi h
in ense sea b eezes du ing he day, he AMB is an ideal case s udy o de e mine he con ibu ion o he u ban
biosphe e and o he an h opogenic sou ces o he o e all u ban OCS budge . The campaigns we e conduc ed
du ing he sp ing and ea ly au umn o 2020 and almos coincided wi h he OCS sp ing maximum and au umn
minimum (Mon zka e al.,2007). Fu he mo e, he May campaign ook place du ing he ull lockdown due o he
COVID pandemic, which helped cons ain he biosphe e signal gi en he educ ion in an h opogenic ac i i ies
ha a e po en ial sou ces o OCS such as ehicle exhaus (Zumkeh e al.,2017). We analyzed he obse ed OCS
concen a ions in e ms o land-use in luence while unde s anding he ole o local wea he condi ions, synop ic
e en s, and he plane a y bounda y laye by using a mosphe ic anspo models, clima e obse a ions and a mos-
phe ic ace s. Mo e speci ically, we a emp ed o answe he ollowing ques ion: Wha is he land-use in luence
on he u ban OCS budge ?
2. Me hods
2.1. S udy Si e
The wo expe imen al campaigns ook place in he Me opoli an A ea o Ba celona (AMB), loca ed in he NE
co ne o Spain, as shown in Figu e1. The AMB, wi h a o al su ace a ea o app oxima ely 640km
2 and
o e 3.3 million inhabi an s (16,160 people/km
2), is he se en h la ges u ban a ea and eigh h mos densely
Figu e 1. (a) A map o Spain showing he Me opoli an A ea o Ba celona (AMB, ed ec angle). O og aphy o he AMB
(b) and land use (c), whe e iplica es o colo ed do s ep esen he upwind, a ge , and downwind sampling loca ions (d).
Backg ound a mosphe ic OCS le els we e in es iga ed a he upwind sampling poin s a he si es o Tibidabo (highes al i ude
in he PBL), Ga à and Poble Nou (closes o he sho e). These si es a e depic ed wi h colo ed do s wi h whi e ci cles in b and
c. The ed diamond ep esen s he loca ion o ICTA and he blue one whe e he adiosondes we e launched and he adon
measu emen s we e pe o med (c).
21698996, 2023, 24, Downloaded om h ps://agupubs.onlinelib a y.wiley.com/doi/10.1029/2023JD039497 by Readcube (Lab i a Inc.), Wiley Online Lib a y on [10/01/2024]. See he Te ms and Condi ions (h ps://onlinelib a y.wiley.com/ e ms-and-condi ions) on Wiley Online Lib a y o ules o use; OA a icles a e go e ned by he applicable C ea i e Commons License
Jou nal o Geophysical Resea ch: A mosphe es
ESTRUCH ETAL.
10.1029/2023JD039497
4 o 21
popula ed egion in he Eu opean Union. I is also one o he mos indus ialized a eas in he wes e n EU (Que ol
e al.,2004), specializing in ood p oduc ion, chemicals, ehicle manu ac u ing, and gas p oduc ion and dis i-
bu ion (Ins i u d’Es adís ica de Ca alunya,2023). The AMB is bounded by he Besós Ri e in he no heas , he
Llob ega Ri e in he wes , he Collse ola moun ain ange o he no h, and he Medi e anean Sea o he sou h,
as shown in Figu e1. The AMB has a Medi e anean d y subhumid clima e wi h an annual p ecipi a ion o
app oxima ely 600mm, which is mainly concen a ed du ing sp ing and au umn, and a mean annual empe a u e
o 16.5°C. The o og aphy o he ci y ex ends om sea le el o 516m abo e sea le el (a.s.l.) a he moun ain ange
loca ed no hwes o he main u ban a eas.
Mo e han hal o he AMB's su ace a ea is co e ed by ege a ion in a ious o ms, including u ban o es (43%),
ag icul u al ields (8%) and u ban pa ks (4%), while 41% is co e ed by impe ious and o he in as uc u es (4%).
Land uses a e shaped by he ollowing physical cons ain s: c ops a e loca ed in he coas al plains o med by he
del a o he Llob ega Ri e , while u ban o es s composed o pine and holm oaks a e loca ed in he Collse ola
moun ain ange (Figu e1). U ban g een spaces a e sca e ed h oughou he me opolis and popula ed wi h o na-
men al species anging om acacias o plan ain ees and o namen al sh ubs and g asses.
2.2. Field Campaigns
Ai samples we e aken a he eigh si es shown wi h colo ed do s in Figu e1c, which we e selec ed o ep e-
sen he he e ogenei y o land use h oughou he AMB and a e summa ized in o ou land-use ca ego ies (see
Table1): ag icul u al a ea (loca ed in he Ga à and P a a eas), u ban o es (Tibidabo and Collse ola), u ban
g een (Mon juic and Guina dó) and u ban (Sag ada Familia and Poble Nou). We ook ai samples a each si e
(he ea e e e ed o as he a ge si e) and upwind and downwind o each a ge si e in an a emp o ollow he
ai mass low du ing i s ajec o y o e he a ge si e. The upwind and downwind sampling si es we e de e -
mined based on his o ical wind pa e ns and ajec o ies ypical o each a ge loca ion unde simila a mosphe ic
condi ions expec ed du ing he campaign and a e a ailable om he Me eo ological Se ice o Ca alonia (www.
me eo.ca ) and shown in Tex S1 in Suppo ing In o ma ionS1. No e ha in May and Oc obe 2020, obse ed
Campaign Si e (ele a ion) Land use ype Da e CWTs a SWC b PBLH ange (m) OCS 7:00 UTC (pp ) OCS 9:30 UTC (pp )
May 2020 Tibidabo (175–442m) Fo es 19/05 C N 19–441 478±20 494±9
Ga à (0–2m) Ag icul u al 21/05 U N 11–249 482.1±NA 495.5±10.7
Poble Nou (0–4m) U ban 21/05 U N 44–2,217 461±20 NA
Mon juic (42–98m) U ban g een 20/05 C N 327–883 485±14 490±5
P a (0–2m) Ag icul u al 25/05 E E 737–1,182 502±NA 508±2
Sag ada (37–55m) U ban 26/05 E E 771–1,900 460±6 NA
Collse ola (107–442m) Fo es 27/05 E E 395–1,079 441±39 497±3
Guina dó (115–169m) U ban g een 28/05 A SE 590–1,524 505±8 612±17
Oc obe 2020 Ga à (0–2m) Ag icul u al 23/10 U NW 21–5,019 401±9 403±3
Poble Nou (0–4m) U ban 15/10 CN N 287–571 405±23 426±45
Guina dó (115–169m) U ban g een 19/10 ASW SW 69–934 389±20 392±16
Sag ada (37–55m) U ban 20/10 SW SW 29–1,085 444±19 460±93
Mon juic (42–98m) U ban g een 21/10 CSW SW 92–734 425±NA 485±87
P a (0–2m) Ag icul u al 22/10 SW SW 103–495 468±24 490±37
No e. An unab idged e sion o his able wi h upwind and downwind OCS including a de ailed da a o backg ound measu emen s concen a ions (see Figu e4) is
a ailable in Tex S1 in Suppo ing In o ma ionS1.
aCWT: Ci cula ion wea he ype. The di e en ypes o CWTs a e cyclonic (C), an icyclonic (A), pu e ad ec i e (N, S, E, W, NE, SE, NW, SW), hyb id cyclone-ad ec i e
(CN, CNE, CE, CSE, CS, CSW, CW, CNW), hyb id an icyclone-ad ec i e (AN, ANE, AE, ASE, AS, ASW, AW, ANW, and unde ined (U)).
bSWC: Synop ic wind
componen summa izing he p essu e pa e n, on s, wind di ec ion and speed.
Table 1
Ci cula ion Wea he Type (CWT), Synop ic Wind Componen (SWC) and PBLH Ex ac ed F om WRF Simula ions Exp essed in Me e s Abo e G ound Le el (a.g.l),
A e aged OCS Concen a ions o All Measu emen s Taken Each Day Du ing he May and Oc obe Campaigns, and Exp essed in pp (Mean Values F om he Th ee
Si es, Mean Values±S anda d De ia ions)
21698996, 2023, 24, Downloaded om h ps://agupubs.onlinelib a y.wiley.com/doi/10.1029/2023JD039497 by Readcube (Lab i a Inc.), Wiley Online Lib a y on [10/01/2024]. See he Te ms and Condi ions (h ps://onlinelib a y.wiley.com/ e ms-and-condi ions) on Wiley Online Lib a y o ules o use; OA a icles a e go e ned by he applicable C ea i e Commons License
Jou nal o Geophysical Resea ch: A mosphe es
ESTRUCH ETAL.
10.1029/2023JD039497
5 o 21
wind pa e ns some imes di e ed om his o ical da a. The loca ions o he upwind, a ge , and downwind si es
a e shown o each poin using he same colo in Figu e1c. The spa ial dis ance among he h ee sampling poin s
o each si e anges om 500 o 4,000m. Sampling imes we e chosen o cap u e inland winds du ing he ea ly
mo ning (i.e., 4:00 and 7:00 UTC am) and he ypical sea b eezes ha de elop du ing he mo ning and ea ly a e -
noon (i.e., 9:30 and 11:00 UTC am). In each campaign, each place was measu ed one day, a he s ipula ed hou s.
The i s campaign ook place be ween 18 May and 28 May 2020, which was du ing a pe iod o in ense anspo
es ic ions due o he COVID-19 pandemic, and esul ed in highly signi ican ai pollu ion educ ions (Badia
e al.,2021). An h opogenic emissions d as ically dec eased, especially o he a ia ion (95%), oad anspo (up
o 85%), and mo e mode a ely acco ding o shipping (21%) and indus y (11%) sec o s (Gue a a e al.,2021).
The second campaign ook place om he 13 o he 27 o Oc obe 2020, when he an h opogenic ac i i y e u ned
o no mal.
The empe a u e anged om 20°C o 24°C du ing May and 14°C o 20°C du ing Oc obe . Samples we e collec ed
a 7:00 UTC and 9:30 UTC, and o some speci ic days chosen a he beginning and a he end o he campaign,
we pe o med addi ional sampling a 4:00 UTC and 11:00 UTC o he Tibidabo, Ga à and Poble Nou upwind
loca ions o de e mine he backg ound concen a ion o OCS. The wea he was s able du ing bo h campaigns
wi h one day o sligh p ecipi a ion du ing each sampling campaign, which accumula ed 2.6 and 1.6mm du ing
May 29 and Oc obe 26, espec i ely. The ai mass in luence was mos ly cyclonic and ad ec i e wi h some days
cha ac e ized by an icyclonic ci cula ion. Synop ic winds we e mos ly no h and eas du ing he May campaign
and p edominan ly sou hwes du ing Oc obe as summa ized in Table S2 in Suppo ing In o ma ionS1.
We eco ded empe a u e, ai p essu e, wind di ec ion and wind speed a each sampling si e using a po able
anemome e (Skywa ch Xplo e 4). The ela i e humidi y was ob ained om he Me eo ological Se ice o Ca a-
lonia o he i s campaign (Se ei Me eo ològic de Ca alunya,2023), and addi ionally, i was egis e ed on si e
using a digi al hyg ome e (TFA Dos mann®, p ecision ±4%) du ing he Oc obe campaign.
2.3. OCS Measu emen s and Analysis
Samples we e collec ed in cylind ical bo osilica e glass lasks (1L olume, No mag Labo und P ozess echnik
GmbH, Ge many) wi h Kel-F PCTFE al es i ed a bo h ends. This lask ma e ial showed he lowes pe mea-
ion o gases compa ed o o he sealing ma e ials (S u m e al.,2004). To ensu e he quali y o he measu emen
echnique and epea abili y o he expe imen , wo lasks connec ed o each o he in a po able de ice we e
simul aneously illed in se ies (see Tex S2 in Suppo ing In o ma ionS1, U ban si es, o de ice pic u es). The
inle was ins alled a 3.5m abo e g ound le el (m a.g.l.) and ai , a e passing h ough a magnesium pe chlo a e
desiccan ube o emo e mois u e and a il e o emo e pa icles, was pumped o he lasks wi h he help o a
diaph agm pump (KNF N84.4 ANDC). Be o e aking he sample, ai was lushed h ough he lasks o a leas
10min. A ha ime, we checked ha he low was s eady and be ween 2.0 and 3.5Lmin
−1. The lask sampling
ook abou o 10min o be comple ed. A e lushing, he exi al e was closed, and he lasks we e p essu ized
wi h sample ai o app oxima ely 1.5ba . Samples om he a ge , upwind and downwind si es we e aken a he
exac same momen o simula e an Eule ian app oxima ion in which he sampling poin emains ixed (Leelőssy
e al.,2016).
Flasks we e sen o he Labo a oi e des Sciences du Clima e de l’En i onnemen , LSCE, Pa is-Saclay, F ance,
o OCS analyses. Flasks we e s o ed o a maximum o 2mon hs be o e analysis. OCS was analyzed ollowing
he analy ical me hod desc ibed in Bel iso e al.(2016, 2022). In his p ocedu e, a mosphe ic OCS (500ml)
is apped using an au oma ed sampling sys em (En ech P7100 p econcen a o ) and hen analyzed by gas
ch oma og aphy wi h pulsed lame pho ome y de ec ion (PFPD, Va ian Model 3,800, Bel iso e al.,2016). Cali-
b a ion was pe o med e e y week using a calib a ion gas supplied by Ai P oduc s (Bel iso e al.,2016). The
di e ence be ween he assigned alue o an NOAA long- e m ai s anda d (NOAA-2004 scale) and he a e age
o mul iple analyses o ha same s anda d was smalle han 1.5%. The d i be ween calib a ions was assessed on
a daily basis using a sho - e m a ge gas (Bel iso e al.,2016).
2.4. De e mining he Plane a y Bounda y Laye Heigh
I is necessa y o de e mine he plane a y bounda y laye heigh (PBLH) o unde s and i s in luence on he
accumula ion and dilu ion o OCS concen a ion a each si e o e he u ban a ea. PBLH measu emen s we e
21698996, 2023, 24, Downloaded om h ps://agupubs.onlinelib a y.wiley.com/doi/10.1029/2023JD039497 by Readcube (Lab i a Inc.), Wiley Online Lib a y on [10/01/2024]. See he Te ms and Condi ions (h ps://onlinelib a y.wiley.com/ e ms-and-condi ions) on Wiley Online Lib a y o ules o use; OA a icles a e go e ned by he applicable C ea i e Commons License
Jou nal o Geophysical Resea ch: A mosphe es
ESTRUCH ETAL.
10.1029/2023JD039497
6 o 21
no possible a he eigh a ge si es, bu a leas he PBLH was calcula ed om measu emen s aken on he
campus o he Uni e si y o Ba celona (Figu e1, label: Radiosonde/Radon) using wo di e en echniques as
ollows: (a) based on obse ed measu emen s om a adiosonde (RSD) s a ion ins alled on he oo o he Physics
Depa men o he Uni e si y o Ba celona (41°23′03″ N 2°07′01″E), which is pa o he Global Me eo ological
Ne wo k (GCOS,2016); and (b) ia a lase ceilome e (CL-31, Vaisala Inc., Finland) loca ed a he same place
wi h a measu emen ange up o 7.6km, 10-m e ical esolu ion and a empo al esolu ion o 16s. Fo he la e ,
he PBLH is es ima ed om he ceilome e da a using he Vaisala Bounda y-Laye View so wa e (BL-VIEW)
Enhanced G adien me hod (VAISALA,2020)—see de ails in Ga cia-Dalmau e al.(2021)— ollowed by selec-
ion algo i hms acco ding o he me hodology o Lo e ane e al.(2016).
To de e mine he PBLH based on adiosonde measu emen s, we ollowed a obus nume ical p ocedu e p oposed
by Liu and Liang(2010) using midday adiosonde (RSD) obse a ions om he selec ed pe iods o s udy and
conside ing 1100 UTC as he launching ime. The PBLH was compu ed ollowing he ai pa cel me hod o a
con ec i e domina ed bounda y laye (Holzwo h,1964; Seibe e al.,2000; S ull,1988) as he heigh whe e ai
pa cel e olu ion ollowing a d y adiaba ic e olu ion om su ace in e sec s wi h he en i onmen empe a u e
p o ile. The calcula ed PBLH om RSD was isually co ec ed o sol e he me hod e o s and mos ly ela ed o
he p esence o se e al empe a u e in e sion laye s. A de ailed explana ion o he me hod and how we applied i
o ou case s udy is p esen ed in Suppo ing In o ma ionS1 (Tex S3).
In his s udy, we also analyzed adon, which is a adioac i e and noble gas known o be a good a mosphe ic ace
(G ossi e al.,2012), o be e unde s and he e olu ion o he PBLH and o co obo a e he PBLH calcula ed om
he ceilome e obse a ions using he me hod p oposed by G i i hs e al.(2013). Hou ly a mosphe ic concen a-
ions o Rn we e measu ed du ing Oc obe 2020 using an A mosphe ic Radon MONi o (ARMON) a he Ba ce-
lona School o Indus ial Enginee ing s a ion (41°23′03″ N 2°07′01″E, 10m a.g.l) a he same si e ha adiosonde
and ceilome e da a we e eco ded and iden i ied wi h he label Radiosonde/Radon in Figu e1c. The moni o ing
was designed and buil by he Uni e si a Poli ècnica de Ca alunya as desc ibed by G ossi e al.(2012). Due o
COVID-19 es ic ions, no adon da a we e eco ded o May 2020 because he equipmen had o be shu down.
The es ima ions o PBLH using adiosonde and ceilome e we e also compa ed wi h PBLH simula ed using he
Wea he Resea ch and Fo ecas ing model (WRF). WRF simula ions we e pe o med o ob ain in o ma ion on
he main anspo p ocesses and he e ical mixing condi ions h oughou he AMB as desc ibed in Sec ion2.6.
2.5. WRF Simula ions o he U ban A mosphe e
The WRF 4.3 model (Skama ock e al.,2021) coupled wi h he mul ilaye u ban canopy scheme BEP-BEM
(Salamanca e al.,2010) was used in his s udy o simula e he u ban a mosphe e o e he AMB o be e unde -
s and he in luence o local wea he condi ions and PBL de elopmen on OCS concen a ions. The WRF con ig-
u a ion used o his s udy consis s o 3 wo-way nes ed domains wi h ho izon al esolu ions o 9, 3, and 1km
co e ing he Ibe ian Peninsula, he Ca alan e i o y, and he AMB, espec i ely, wi h 45 e ical laye s up o
100hPa (see Figu e S2 in Suppo ing In o ma ionS1). The physics pa ame e iza ion o he plane a y bounda y
laye is he Boulac scheme, which has been shown o be he mos sui able o esol e he u bulence in ensi y in
con ec i e si ua ions o he AMB (Segu a e al.,2021). We included 11 u ban local clima e zone (LCZ) clas-
si ica ions (S ewa & Oke,2012) o he inne mos domain and each had a speci ic alue o he mal, adia i e
and geome ic pa ame e s o he buildings and g ound, which a e used by he BEP-BEM o compu e he hea
and momen um luxes in u ban a eas (Gilabe e al.,2021; Ribei o e al.,2021; Yu & S einbe ge ,2011). We
conduc ed simula ions o May and Oc obe 2020 using me eo ological ini ial and bounda y condi ions om
ERA5 (Cope nicus,2023). Table S3 in Suppo ing In o ma ionS1 u he desc ibes he model cha ac e is ics
and expe imen al con igu a ions. The model simula ions we e alida ed using empe a u e, wind speed and
wind di ec ion obse a ions om 14 me eo ological s a ions o he Me eoca ne wo k (h ps://www.me eo.ca /
obse acions/xema).
2.6. FLEXPART Back T ajec o ies
FLEXPART back ajec o ies we e calcula ed o de e mine he p ecedence o ai sampled a he a ious si es
o es ablish whe he he ai mass in luence was local o egional. FLEXPART (FLEXible PARTicle dispe sion
model) is a Lag angian anspo and dispe sion model ha is used o calcula ing he long- ange and mesoscale
21698996, 2023, 24, Downloaded om h ps://agupubs.onlinelib a y.wiley.com/doi/10.1029/2023JD039497 by Readcube (Lab i a Inc.), Wiley Online Lib a y on [10/01/2024]. See he Te ms and Condi ions (h ps://onlinelib a y.wiley.com/ e ms-and-condi ions) on Wiley Online Lib a y o ules o use; OA a icles a e go e ned by he applicable C ea i e Commons License
Jou nal o Geophysical Resea ch: A mosphe es
ESTRUCH ETAL.
10.1029/2023JD039497
7 o 21
dispe sion o ai pollu an s (S ohl e al.,2005). The FLEXPART e sion
adap ed o he WRF mesoscale model (B ioude e al.,2013) was used o
calcula e back ajec o ies o ai ace s om each o he sampling si es a
he sampling imes. The back ajec o ies we e compu ed om he ou pu o
he WRF model a he nes ed domain o 1km wi h 16 e ical le els om
0 o 2,000m a.g.l. Each back ajec o y was calcula ed o 24h a e he
sampling ime. The numbe o pa icles eleased o he a mosphe e in he
simula ion was 2,000. F om he 3D back ajec o ies, only he laye s om
0 o 100m abo e g ound le el we e used o he s udy o oo p in s, as i is
whe e he exchange among he land/ocean - a mosphe e akes place.
3. Resul s
3.1. OCS Sample Repea abili y and Accu acy in U ban En i onmen s
We pe o med a p ecision analysis o he GC measu emen s o lask-ai
samples acco ding o he me hods de eloped by Bel iso e al. (2022) o
assess he epea abili y and accu acy o he OCS measu emen s in an u ban
en i onmen . Eigh pai s o lasks in o al we e selec ed andomly om he
May and Oc obe campaigns, and each lask was analyzed wice consecu-
i ely o OCS. The a e age di e ence be ween duplica ed measu emen s,
ha is, he in a lask di e ence, was 9.8±8.7 pp (median=8.2 pp , in e -
qua ile ange (IQR)=3.6–12.7 pp , n=16, 90 h pe cen ile=16.5 pp , 6.3%
o ou lie s), as shown in Figu e2a. I was sligh ly highe han ha epo ed
by Bel iso o samples collec ed in a u al a ea o e whea and apeseed
c ops (Bel iso e al.,2022), i.e., 9.3±7.2 pp (median=7.5 pp , in e qua ile
ange (IQR)=4.5–12.5 pp , n=58, 5.2% o ou lie s, bu in ha case, he
90 h pe cen ile was 19.5 pp ). Figu e2c shows ha among he 54 lasks (each
analyzed wice consecu i ely), ou displayed mixing a ios highe han he
90 h pe cen ile (i.e., 22.9 pp ), while he median and IQR o he in a lask
di e ence (7.7 pp and 3.3–12.1 pp , espec i ely) appea o be highly
consis en wi h he da a displayed in Figu e2a. Consis ency is also ound in
e ms o he p opo ion o abe an de e mina ions o OCS mixing a ios (i.e.,
6.3% (Figu es2a) and 7.4% (Figu e2c)). I was concluded ha he p ecision
o GC measu emen s o lask-ai samples emained unchanged o e he wo
campaigns and was consis en wi h ha epo ed by Bel iso e al.(2022).
The a e age di e ence be ween lasks o he same pai (i.e., he in e lask
di e ence) was also analyzed in a duplica e manne and was 10.1±14.0 pp
(median=4.8 pp , IQR=0.05–18.6 pp , n=8, no ou lie , Figu e2b). Hence, he in e lask di e ence du ing he
AMB s udy displayed highe a iabili y han ha o Bel iso e al.(2022), ha is, 5.5±4.5 pp (median=4.9 pp ,
IQR=2.2–8.0 pp , n=29, a 6.9% p opo ion o ou lie s was epo ed in ha s udy). The in e lask di e ence was
also assessed om a single analysis o he con en o each lask. The median, IQR and 90 h pe cen ile we e 12.1
pp , 4.3–22.5 pp and 45.0 pp , espec i ely, wi h n=95 (Figu e2d). We epo ed he exis ence o six ou lie s
(i.e., a p opo ion o 6.3%) exhibi ing in e lask di e ences highe han 45.0 pp and up o 360 pp . The in e lask
di e ence calcula ed om pai s o lasks each analyzed only once o OCS (Figu e2d) displays highe a iabili y
han ha om pai s o lasks analyzed wice consecu i ely (Figu e2b).
Consequen ly, al hough he ou ou lie s displayed in Figu e2c should ha e been disca ded, we a bi a ily ein-
co po a ed hose exhibi ing in a lask di e ences in he ange o 22.9 pp (i.e., he 90 h pe cen ile) and 36.8 pp .
Acco ding o Bel iso e al.(2022), da a om pai s o lasks exhibi ing in e lask di e ences highe han 14 pp
should be disca ded. No e ha his h eshold is app oxima ely wice ha se by Mon zka e al.(2007). Al hough
he in e lask di e ence calcula ed om pai s o lasks each analyzed only once o OCS could o igina e om
abe an da a (Figu e2d), nobody e e documen ed he high equency dis ibu ion o OCS mixing a ios nex
o po en ial u ban sou ces o ha gas. In o he wo ds, al hough lasks a e supposed o be simul aneously illed
in pai s, he p ocess o p essu iza ion o wo lasks connec ed in se ies could gene a e concen a ion g adien s
Figu e 2. Analysis p ecision o OCS measu emen s om lask-ai samples
collec ed in he AMB showing he 10 h, 25 h, median, 75 h and 90 h
pe cen iles. Panels (a, c) show Δ in e lask di e ences in OCS concen a ion
be ween wo consecu i e analyses o he con en o each lask (median o lask
1 minus ha o lask 2), whe eas (b, d) show he Δ in e lask di e ences in
OCS concen a ion be ween lasks o he same pai . The diamonds co espond
o he ou lie s (>90 h pe cen ile), which a e o e 23 pp in panel (c) and
o e 45 pp in (d). Ou lie s a e labeled acco ding o he loca ion and ime o
sample collec ion. The codes o each sampling si e a e si e, loca ion and ime
UTC (si es: Tib: Tibidabo, P a: P a , Pob: Poble Nou, Mon: Mon juic, Gui:
Guina dó; loca ion: U: Upwind, D: Donwind, T: Ta ge ). No e ha his o ical
wind da a was used o ha pu pose. The label colo s ep esen he campaigns
and co espond o ed o May and blue o Oc obe .
21698996, 2023, 24, Downloaded om h ps://agupubs.onlinelib a y.wiley.com/doi/10.1029/2023JD039497 by Readcube (Lab i a Inc.), Wiley Online Lib a y on [10/01/2024]. See he Te ms and Condi ions (h ps://onlinelib a y.wiley.com/ e ms-and-condi ions) on Wiley Online Lib a y o ules o use; OA a icles a e go e ned by he applicable C ea i e Commons License
Jou nal o Geophysical Resea ch: A mosphe es
ESTRUCH ETAL.
10.1029/2023JD039497
8 o 21
e lec ing he he e ogeneous dis ibu ion o OCS nex o po en ial u ban sou ces o ha gas, in ou case samples
ook 10min o ill ou in se ies. In ligh o his possibili y, we decided o a bi a ily alida e he six ou lie s
epo ed in Figu e2d and in es iga e whe he in e lask di e ences highe han 45 pp would be si e dependen
o would esul om abe an da a.
3.2. E alua ion o WRF Simula ions
We an WRF BEP-BEM simula ions in May 2020 and Oc obe 2020 o ob ain high- esolu ion me eo ological
da a, a mosphe ic s abili y, and PBLH du ing he sampling pe iods. S a is ical analyses o he oo -mean-squa e
e o (RMSE), mean bias (MB) and co ela ion ac o (R) we e used o e alua e pe o mances o he model
compa ed wi h g ound-based me eo ological obse a ions om Me eoca (see Sec ion2.6). The nea -su ace
ela i e humidi y (RH), empe a u e (T) and wind speed (WS) s a is ical analyses a e displayed in Figu es S3
and S4 in Suppo ing In o ma ionS1 o campaigns in May and Oc obe , espec i ely. The e is an o e all good
pe o mance o he modeled RH, T and WS in compa ison wi h he obse a ions wi h an a e age e o o ±14%
o RH, ±1.7 o 1.9°C o T and ±1.65ms −1 o WS. This s a is ical analysis is a ailable in Table S4 in Suppo ing
In o ma ionS1.
Du ing he May campaign, he e we e p edominan no he ly winds om he no h om he 18 h un il he 21s ,
and he ci cula ion wea he ypes we e hyb id-cyclonic o cyclonic egimes, while he second hal (25 h h ough
he 28 h) was cha ac e ized by winds om he eas wi h a Medi e anean Sea ai mass in luence and ad en i e
egimes, wi h he excep ion o he 28 h, which was an icyclonic wi h a sou heas synop ic wind componen and
local ma i ime ai mass in luence. Du ing he Oc obe campaign, wes e ly and sou he ly winds cha ac e ized by
ad en i e ci cula ion wea he ypes p e ailed o mos o he pe iod wi h a cyclonic no h in luence and egional
con inen al ai mass in he i s 2days swi ching o a ma i ime in luence o he es o he campaign. Table1
summa izes hese wea he pa e ns o each day o he wo campaigns.
3.3. Plane a y Bounda y Laye Heigh Calcula ions and E alua ion
We had adiosonde and ceilome e da a o calcula e he PBLH in one loca ion o he Me opoli an A ea o
Ba celona, which we compa ed o he WRF-BEP BEM simula ed PBLH and es ablished ha he WRF model
pe o med easonably well in e ms o bounda y laye de elopmen . Fu he mo e, WRF-BEP BEM was alida ed
in e ms o empe a u e, wind di ec ion and speed, and ela i e humidi y. Thus, we we e able o employ WRF
me eo ological ou pu such as wind di ec ion and speed o explain he obse ed OCS concen a ions, mainly
in e m o hei accumula ion in o low o high PBLH and possible anspo pa e ns. Fo he same eason, he
a mosphe ic s abili y classes, based on he adon ace me hod, may gi e mo e in o ma ion on he a mosphe ic
episodes analyzed in o he s udy.
The PBLH om he WRF BEP-BEM simula ions was compa ed agains ceilome e es ima ions and adiosonde
obse a ions a 11:00 UTC o each day o he campaigns (see Figu e3). All da a ime se ies ollow he expec ed
PBLH diu nal pa e n wi h a apid g ow h o he PBLH du ing he mo ning, a maximum ha was eached a
midday, and a dec ease du ing he a e noon. The PBLH calcula ed by he WRF model is gene ally highe han
ha es ima ed by he ceilome e and adiosonde. Du ing he i s pe iod o May and i s days o Oc obe when
he e a e con inen al egional ai masses wi h p e ailing la ge-scale no he n winds, we no iced an impo an
di e ence be ween he ceilome e , adiosonde, and WRF alues in some cases up o 150m di e ence a 11
UTC. In his si ua ion, he la ge-scale winds domina e he bounda y-laye scales, which limi he capaci y o he
WRF model o es ima e he PBLH and leads o high unce ain y in he adiosonde obse a ions (S ull,1988).
The in e se co ela ion be ween a mosphe ic adon concen a ions and he ceilome e -based PBLH is uni o mly
main ained du ing he en i e Oc obe campaign as you can see in Figu e S5 in Suppo ing In o ma ionS1 how
PBLH a e a is maximum a 00:00 UTC while a mosphe ic adon concen a ion a e a i s lowes , and he opposi e
beha io is seen a ound 12:UTC (see Figu es S5 and S6 in Suppo ing In o ma ionS1), which u he indica es
an o e es ima ion o he PBLH by WRF du ing s ong wind e en s in ea ly Oc obe . The highes adon ac i i y
concen a ions (be ween 0.01 and 0.05 Bqm
−3) a e obse ed du ing nigh ime (21–2h ) when he PBLH eaches
i s minimum, whe eas he lowes alues (be ween 1 and 2Bqm
−3) a e obse ed du ing diu nal hou s (11–16h )
when
222Rn is dilu ed due o a well-de eloped PBL.
The o e es ima ed PBLH modeled alues in May a e consis en wi h an o e p edic ion o he 2m empe a u e
seen o a s a ion loca ed 1.2km om he adiosonde launching loca ion. In con as , be e ag eemen be ween
21698996, 2023, 24, Downloaded om h ps://agupubs.onlinelib a y.wiley.com/doi/10.1029/2023JD039497 by Readcube (Lab i a Inc.), Wiley Online Lib a y on [10/01/2024]. See he Te ms and Condi ions (h ps://onlinelib a y.wiley.com/ e ms-and-condi ions) on Wiley Online Lib a y o ules o use; OA a icles a e go e ned by he applicable C ea i e Commons License
Jou nal o Geophysical Resea ch: A mosphe es
ESTRUCH ETAL.
10.1029/2023JD039497
9 o 21
he modeled and obse ed PBLH (bo h adiosonde and ceilome e ) is seen du ing he second hal o he May
campaign om he 25 h o 28 h o May (Figu e3a), when he eas e ly winds lead o a mo e s able a mosphe e.
Oc obe showed be e ag eemen o e all among he PBLH alues ob ained om he ceilome e , RSD and WRF
model due o highe a mosphe ic s abili y. The sou h-eas e ly weak winds gene ally lead o s able condi ions and
local scale dominance whe e he PBLH can be well es ima ed.
3.4. WRF-FLEXPART Foo p in s
FLEXPART simula ions showed ha he sampling campaign co e ed di e en synop ic si ua ions and di e en
ai mass o igins. Tex S4 in Suppo ing In o ma ionS1 shows ha du ing he May campaign, wo main di ec ions
we e obse ed as ollows: one coming om he no hwes ollowing he Llob ega Ri e basin (days 18, 19, and
21 o May) and ano he om he no heas wi h mo e sea in luence (days 20 and 25 o May) and some o he s
wi h mo e land in luence (such as days 26, 27, and 28 o May). Gene ally, he simula ion a 09:30 UTC shows
mo e sea o igin o he ai masses a i ing a sampling poin s han he simula ions a 07:00 UTC in acco dance
wi h he sea‒land b eeze egimes.
Du ing he Oc obe campaign, he ai masses mainly came om he no hwes ollowing he Llob ega Ri e
basin (days 15, 19, and 27 o Oc obe ), al hough a si ua ion wi h s ong winds coming om he sou h was
obse ed be ween he 20 h and 21s o Oc obe . A s onge b eeze egime is obse ed du ing his mon h wi h
gene al changes in he o igin o ai masses be ween he simula ions a 07:00 and 09:30 UTC.
Figu e 3. The ime se ies o modeled PBLH (blue line) oge he wi h he Ceilome e ( ed line) and he Radiosonde (blue
shade) o (a) May and (b) Oc obe campaigns. Shaded blue a eas indica e he PBLH using wo di e en de ini ions o i s
calcula ion. Ve ical lines co espond o he 4 a.m., 7 a.m., 9:30 a.m., and 11 a.m. UTC sampling imes on a daily basis.
21698996, 2023, 24, Downloaded om h ps://agupubs.onlinelib a y.wiley.com/doi/10.1029/2023JD039497 by Readcube (Lab i a Inc.), Wiley Online Lib a y on [10/01/2024]. See he Te ms and Condi ions (h ps://onlinelib a y.wiley.com/ e ms-and-condi ions) on Wiley Online Lib a y o ules o use; OA a icles a e go e ned by he applicable C ea i e Commons License
Jou nal o Geophysical Resea ch: A mosphe es
ESTRUCH ETAL.
10.1029/2023JD039497
16 o 21
he no h whe e he ag icul u al ields a e loca ed. This indica es ha he OCS d awdown migh no be en i ely a
local e ec o ag icul u e bu could also be due o a mo e egional e ec o OCS aken up by soil and ege a ion
du ing he nigh all o e Ca alonia and hen ad ec ed by he land b eeze o e he AMB. Howe e , by 9:30 UTC,
OCS concen a ions inc eased (up o an ΔOCS o +116 o he downwind si e) when ai masses we e coming
om he u ban and po a eas and despi e PBLH de elopmen , which eached 1,380m ha mo ning. Du ing he
Oc obe campaign, sou hwes e ly winds p edomina ed bo h a 7:00 and 9:30 UTC, which b ough in ai masses
om he po and Ga à u ban a ea wi h high indus ial ac i i y esul ing in ΔOCS be ween +82 and +40 pp a
7:00 UTC. A 9:30 UTC, ΔOCS alues dec eased sligh ly and had less a iabili y anging om +49 o +55 pp ,
and his indica ed a mo e homogenous ai mix u e due o a mo e de eloped PBL, which inc eased om 21–84
o 380–519m o ha loca ion. Du ing Oc obe , ai anspo had much mo e in luence on he OCS budge han
he ag icul u al land use as i did du ing he May campaign when he ields we e in ull pho osyn he ic ac i i y.
The opposi e beha io was obse ed a he El P a ag icul u al si e, whe e a deep PBLH (1,200m) and egional
ai mass in luence occu ed in May, while a mo e s able si ua ion and lowe PBLH (691m) domina ed in Oc obe .
In gene al, he alues we e simila o he backg ound OCS o May wi h ΔOCS alues o +25, +37, and +22 pp
a 7:00 UTC and e en highe enhancemen s a 9:30 (+17, +44 and +173 pp , Figu e7b). FLEXPART back ajec-
o ies in May showed ai masses coming om he eas wi h longe esiden ial imes loca ed in he sea bu s ill
ecei ed in luence om he ai po and he po o Ba celona. Du ing he Oc obe campaign, El P a had highe
OCS enhancemen s a bo h sampling imes (ΔOCS o +70, +33, and +80 pp o 7:00 UTC; +7, +122, and +77
pp o 9:30 UTC), which, gi en he low PBLH, sugges s a local OCS sou ce. The FLEXPART showed ha ai
masses came om he eas indica ing ma ine in luence and/o om he po o Ba celona, whe e ac i i ies, such
as ai na iga ion and ma i ime a ic, had eco e ed 100% a e he lockdown. The lowe OCS concen a ions
in May compa ed o Oc obe highligh he signi ican sink o ag icul u e du ing he g owing season as was seen
o he Ga à loca ion.
Figu e 7. OCS enhancemen (Δ OCS) in pp calcula ed as he measu e o OCS minus backg ound OCS (484 pp o May
and 407 pp o Oc obe ) o each campaign a he wo u ban pa ks as ollows: (a) Ga à and (b) El P a . The black a ow
shows he wind di ec ion and speed measu ed on si e, and he absence o an a ow indica es ha he e was no wind blowing
a he momen o measu emen . A small map embedded in he lowe igh -hand co ne o each panel shows he esidence ime
on a loga i hmic scale o he ai masses a i ing a he sampling si e as de e mined wi h he FLEXPART model.
21698996, 2023, 24, Downloaded om h ps://agupubs.onlinelib a y.wiley.com/doi/10.1029/2023JD039497 by Readcube (Lab i a Inc.), Wiley Online Lib a y on [10/01/2024]. See he Te ms and Condi ions (h ps://onlinelib a y.wiley.com/ e ms-and-condi ions) on Wiley Online Lib a y o ules o use; OA a icles a e go e ned by he applicable C ea i e Commons License
Jou nal o Geophysical Resea ch: A mosphe es
ESTRUCH ETAL.
10.1029/2023JD039497
17 o 21
To summa ize, he luc ua ions in OCS mixing a ios in u ban pa ks seemed o be mo e hea ily in luenced by
he di ec ion and in ensi y o he wind. In u ban o es s, he e was a clea educ ion o OCS compa ed o u ban
pa ks and buil en i onmen . Du ing he May campaign wi h he signi ican a ic educ ion due o he COVID
es ic ions, we saw a educ ion o OCS concen a ions compa ed o Oc obe - leading us o he conclusion ha
an h opogenic sou ces in u ban a eas a e signi ican . The pe i-u ban ag icul u al a eas did no show lowe alues
o OCS as was expec ed, mos likely due o he high in luence o he ai po and he seapo ha we e close by.
In gene al, i was di icul o de e mine i he e was any signi ican OCS d awdown by ege a ion, excep in
u ban o es ed a eas, while OCS mixing a ios well abo e backg ound le els e idenced signi ican u ban sou ces
o OCS.
4.2. Conside a ions o OCS S udies in U ban Ecosys ems
An h opogenic emissions a e he p ima y sou ce o OCS on a plane a y scale, ollowed by oceanic lux (Remaud
e al.,2023), hey a e signi ican a he local and u ban scales. In ou s udy, he peaks o OCS concen a ions we e
mos ly due o ai masses coming om he po o he ci y, as was he case o Mon juic in Oc obe o Ga à du ing
he May campaign, and eached alues o 665 pp o Mon juic and 600 pp o Ga à (up o ΔOCS o +258 and
+193 pp , espec i ely).
The e a e a ious po en ial sou ces o OCS in u ban a eas, such as he e osion o au omobile i es agains
asphal , ehicle exhaus umes, na u al gas combus ion and indus ial p ocesses including oil e ine ies and
aluminum p oduc ion (Yan e al.,2019). In ac , in hei in en o y, Yan e al. quan i ied he o al emissions in
2015 o China alone as 174GgSy
−1 wi h he highes con ibu ions coming om indus ial OCS emission
p ocesses (120GgSy
−1). In ha same s udy, he o al ehicle exhaus was esponsible o 3GgSy
−1, which
is app oxima ely hal o he yea ly global es ima e o 6GgSy
−1 gi en by Lee and B imblecombe(2016). The
la e s udy also indica ed ha he hea y uel oil used in ships is a signi ican sou ce o he OCS budge and
p oduces 30GgSy
−1. The AMB quali ies o all hese po en ial sou ces and i has a hea ily ansi ed ma ine
po ha also has a combined cycle gas u bine powe plan ha uns on na u al gas as well as an in ense use
o p i a e ehicles wi hin he ci y.
Figu e 7. (Con inued)
21698996, 2023, 24, Downloaded om h ps://agupubs.onlinelib a y.wiley.com/doi/10.1029/2023JD039497 by Readcube (Lab i a Inc.), Wiley Online Lib a y on [10/01/2024]. See he Te ms and Condi ions (h ps://onlinelib a y.wiley.com/ e ms-and-condi ions) on Wiley Online Lib a y o ules o use; OA a icles a e go e ned by he applicable C ea i e Commons License
Jou nal o Geophysical Resea ch: A mosphe es
ESTRUCH ETAL.
10.1029/2023JD039497
18 o 21
Addi ionally, a b oad a ie y o sul ides a e p oduced by was ewa e whe e ola ile o ganic compounds (VOC)
such as me hane hiol (MeSH), dime hyl sul ide (DMS), ca bon disul ide (CS2), dime hyl disul ide(DMDS), dime-
hyl isul ide (DMTS) and OCS a e gene a ed in anoxic en i onmen s o sewe s (Lee & B imblecombe,2016),
and i has been no ed ha he mic obial ac i i y in hese anoxic sedimen s is closely linked o he p oduc ion
o OCS (Ki z e al.,2019). Sul ide emissions om u ban was ewa e ne wo ks a e equen ly de ec ed due o
hei bad odo as is he case o Ba celona (Eijo-Río e al.,2015) and Thessaloniki (Besis e al.,2021). In ou
s udy, he sampling si es o Sag ada Familia in Oc obe (ΔOCS o +150 pp ) and Poble Nou (ΔOCS o +50
pp ) concu wi h loca ions o a combined sewe sys em p e iously epo ed o ha e high peaks o H2S, which
indica ed a po en ial OCS sou ce om was ewa e . Un o una ely, cons aining his signal om o he po en ial
u ban sou ces was ou o he scope o he p esen s udy. Howe e , he obse a ions p o ided by ou measu emen
campaigns could se e o es ima e an h opogenic OCS sou ces when compa ed o simula ions o OCS mole
ac ions combined wi h emission in en o ies and plan up ake luxes as was ecen ly done o he Ne he lands
by Zanche a e al.(2023).
In e ms o ege a ed land uses, we ound one case o a signi ican OCS d awdown in ag icul u al si es du ing he
g owing season o May. In addi ion o pho osyn he ic up ake (Campbell e al.,2008,2017; Mon zka e al.,2007),
wa e logged soils ha e been epo ed o ac as sinks o OCS when plan ed wi h c ops (Yi e al.,2008). The
ag icul u al ields o he Llob ega basin a e mos ly i iga ed by inunda ion, which could u he in ensi y he
OCS up ake and explain he OCS d awdown a 7:00 UTC when pho osyn he ic ac i i y was high. We also ound
alues indica ing up ake du ing he g owing season in Collse ola Fo es , which is in ag eemen wi h o he s udies
e alua ing o es OCS up ake (Bel iso e al.,2016; Campbell e al.,2008,2017).
Finally, de e mining he in luence o land use on he o e all u ban OCS budge equi es a eliable backg ound
measu emen . In gene al, backg ound le els o OCS a e associa ed wi h ai masses coming om he sea because
he s abili y p o ided by oceans esul s in cons an OCS concen a ions, a leas a a seasonal scale (Mon zka
e al.,2007), as was he case o San F ancisco (Villalba e al.,2021). Howe e , his is no applicable o he AMB
because he con ec ion egime is qui e peculia . The AMB has a unique land‒sea ai mass exchange in which
sea b eezes unde speci ic condi ions do no b ing new ai masses bu a he ecycle ai masses ha o igina ed
o e he u ban a ea, howe e , a deepe esea ch should be done in his ega d, since ma ine con ibu ion canno
be neglec ed. This ype o p ocess has been obse ed by Jaén e al.(2021), who desc ibed how ozone plumes
s ay o e he sho e due o ad ec i e p ocesses du ing hea episodes such as hea wa es. The he mal componen
o con ec i e p ocesses c ea ing episodes o ai eci cula ion on he Spanish eas e n coas has also been docu-
men ed by Millán e al.(2000), whe e hey ou line how he sea b eeze can a o he o ma ion o laye s o s a -
i ied ai ese oi s in heigh s whe e compounds accumula e ha can la e e u n o he land. This pa icula sea
b eeze ci cula ion cons ic s local en ila ion and exace ba es he accumula ion o an h opogenic gases emi ed
in he AMB.
5. Conclusions
We pe o med wo measu emen campaigns o OCS mixing a ios in he Me opoli an A ea o Ba celona (AMB)
o explo e o wha ex en di e en u ban land uses can ha e an impac on he OCS mixing a ios measu ed on
lask samples ega ding a pa icula u ban clima e, geog aphy, and opog aphy o AMB. Ai samples we e aken
a loca ions ep esen ing ou di e en land uses as ollows: u ban pa ks, o es s, ag icul u al ields and impe i-
ous u ban landscapes.
Based on he in e lask compa ison desc ibed in 2.3, we de e mined ha he OCS measu emen unce ain y can
be as high as +53 pp , which e lec s ha nea -su ace measu emen s in u ban he e ogeneous en i onmen s esul
in highe unce ain ies han in homogenous u al en i onmen s. We also compa ed a e age OCS concen a ions
±σ wi h alues obse ed om ICOS- LSCE in Gi -su -Y e e (GIF), which we e aken a he same momen as
ou campaign (Figu e S3 in Suppo ing In o ma ionS1). We obse ed how mean OCS concen a ions and σ in
he AMB we e sys ema ically highe han in GIF, which sugges s ha an h opogenic emissions could a ec he
inc ease in unce ain y.
The abili y o cap u e a pho osyn he ic signal appea s weak in ege a ed a eas wi hin he ci y, such as pa ks, and
is elega ed o o es s and pe i-u ban ag icul u al a eas. U ban pa ks showed alues signi ican ly highe han he
es ablished backg ound alue, bo h in May (i.e., up o +144 pp ) and in Oc obe (i.e., up o +150 pp ), which
21698996, 2023, 24, Downloaded om h ps://agupubs.onlinelib a y.wiley.com/doi/10.1029/2023JD039497 by Readcube (Lab i a Inc.), Wiley Online Lib a y on [10/01/2024]. See he Te ms and Condi ions (h ps://onlinelib a y.wiley.com/ e ms-and-condi ions) on Wiley Online Lib a y o ules o use; OA a icles a e go e ned by he applicable C ea i e Commons License
Jou nal o Geophysical Resea ch: A mosphe es
ESTRUCH ETAL.
10.1029/2023JD039497
19 o 21
shows he in luence o an h opogenic sou ces, possibly om he po and he ci y cen e , which masks he possi-
bili y o cap u ing he signal om he ege a ion. The u ban o es sugges ed OCS up ake alues up o −33 pp ,
which indica es a signi ican pho osyn hesis sink in dense ege a ion masses e en wi hin ci y limi s.
The u ban a ea showed a a ie y o signals om almos neu al o abo e he backg ound alue, which sugges s
nea by emission sou ces. In ou case, we ound he explana ion o some p ocess alues ha may be aligned wi h
u ban a ic bu abo e all wi h he indus ial ac i i y o he po and ai po (e.g., high OCS alues in Mon juic
when he wind came di ec ly om he po ). We can also hink o o he biogenic p ocesses ha p oduce OCS
p ocesses linked o human ac i i y on u ban land such as he was e sys em. These p ocesses should be s udied
in mo e de ail in u u e wo ks. We ha e no been able o cap u e an e ec o u ban ege a ion on OCS up ake;
ins ead, he pe i-u ban a ea o ag icul u e and o es s may ha e mo e impo an oles in he OCS dynamics in
ci ies.
Las , ag icul u al a eas showed a clea d awdown du ing he g owing season o he loca ion o Ga à, whe e we
we e able o disca d egional in luences and in ense ai mixing. OCS enhancemen s in Oc obe e lec he loss o
pho osyn he ic up ake and he e u n o business as usual a ic and nea by ai po and po ac i i ies.
O e all, he a iabili y o OCS in he u ban ecosys em was highe han we expec ed. The URBAG campaigns o
OCS measu emen s showed e idence ha he use o OCS as a ace gas o de e mine he con ibu ion o he u ban
biosphe e o he CO2 budge is a complica ed ask and ha he ege a ion signal is di icul o de ec because local
an h opogenic ac i i ies ha e a signi ican in luence. The u ban ecosys em is complex in e ms o land up ake,
ma ine emissions, and an h opogenic emissions om he ci y po , and s ongly impac s he spa ial- empo al
dis ibu ion o OCS in he AMB and impedes a clea cons ain o he biosphe e signal.
Da a A ailabili y S a emen
The da a acqui ed du ing he URBAG measu emen campaigns used o he analysis o land in luence in his
s udy a e a ailable a ZENODO ia h ps://doi.o g/10.5281/zenodo.8072833 wi h C ea i e Commons A ibu ion
4.0 In e na ional license (Villalba, G.,2023).
Re e ences
Badia, A., Langemeye , J., Codina, X., Gilabe , J., Guile a, N., Vidal, V., e al. (2021). A ake-home message om COVID-19 on u ban ai
pollu ion educ ion h ough mobili y limi a ions and elewo king. npj U ban Sus ainabili y, 1, 35. h ps://doi.o g/10.1038/s42949-021-00037-7
Bel iso, S., Abadie, C., Mon agne, D., Hadja , D., T opée, D., Viale es, L., e al. (2022). Ca bonyl sul ide (COS) emissions in wo ag oecosys-
ems in cen al F ance. PLoS One, 17(12), 1–15. h ps://doi.o g/10.1371/jou nal.pone.0278584
Bel iso, S., Pison, I., Pe i , J. E., Be che , A., Remaud, M., Simon, L., e al. (2023). The Z-2018 emissions in en o y o COS in Eu ope: A semi-
quan i a i e mul i-da a-s eams e alua ion. A mosphe ic En i onmen , 300, 119689. h ps://doi.o g/10.1016/j.a mosen .2023.119689
Bel iso, S., Rei e , I. M., Loube , B., G os, V., La hiè e, J., Mon agne, D., e al. (2016). A op-down app oach o su ace ca bonyl sul ide
exchange by a Medi e anean oak o es ecosys em in sou he n F ance. A mosphe ic Chemis y and Physics, 16(23), 14909–14923. h ps://doi.
o g/10.5194/acp-16-14909-2016
Besis, A., Geo giadou, E., & Sama a, C. (2021). Odo -ac i e ola ile o ganic compounds along he sea on o Thessaloniki, G eece. Implica ions
o sou ces o nuisance odo . Science o he To al En i onmen , 799, 149388. h ps://doi.o g/10.1016/j.sci o en .2021.149388
B ioude, J., A nold, D., S ohl, A., Cassiani, M., Mo on, D., Seibe , P., e al. (2013). The Lag angian pa icle dispe sion model FLEXPART-WRF
e sion 3.1. Geoscien i ic Model De elopmen , 6, 1889–1904. h ps://doi.o g/10.5194/gmd-6-1889-2013
Campbell, J. E., Ca michael, G. R., Chai, T., Mena-Ca asco, M., Tang, Y., Blake, D. R., e al. (2008). Pho osyn he ic con ol o a mosphe ic
ca bonyl sul ide du ing he g owing season. Science, 322(5904), 1085–1088. h ps://doi.o g/10.1126/science.1164015
Campbell, J. E., Whelan, M. E., Be y, J. A., Hil on, T. W., Zumkeh , A., S ineciphe , J., e al. (2017). Plan up ake o a mosphe ic ca bonyl sul ide
in coas edwood o es s. Jou nal o Geophysical Resea ch: Biogeosciences, 122(12), 3391–3404. h ps://doi.o g/10.1002/2016JG003703
Campbell, J. E., Whelan, M. E., Seib , U., Smi h, S. J., Be y, J. A., & Hil on, T. W. (2015). A mosphe ic ca bonyl sul ide sou ces om an h opogenic
ac i i y: Implica ions o ca bon cycle cons ain s. Geophysical Resea ch Le e s, 42(8), 3004–3010. h ps://doi.o g/10.1002/2015GL063445
Commane, R., He ndon, S. C., Zahnise , M. S., Le ne , B. M., McManus, J. B., Munge , J. W., e al. (2013). Ca bonyl sul ide in he plane-
a y bounda y laye : Coas al and con inen al in luences. Jou nal o Geophysical Resea ch: A mosphe es, 118(14), 8001–8009. h ps://doi.
o g/10.1002/jg d.50581
Cope nicus (2023). ERA5 hou ly da a on single le els om 1940 o p esen . Re ie ed om h ps://cds.clima e.cope nicus.eu/cdsapp#!/da ase /
eanalysis-e a5-single-le els? ab=o e iew
Eijo-Río, E., Pe i -Boix, A., Villalba, G., Suá ez-Ojeda, M. E., Ma in, D., Amo es, M. J., e al. (2015). Municipal sewe ne wo ks as sou ces o
ni ous oxide, me hane and hyd ogen sulphide emissions: A e iew and case s udies. Jou nal o En i onmen al Chemical Enginee ing, 3(3),
2084–2094. h ps://doi.o g/10.1016/j.jece.2015.07.006
Ga cia-Dalmau, M., Udina, M., Bech, J., Sola, Y., Mon olio, J., & Jaén, C. (2021). Pollu an concen a ion changes du ing COVID-19 lockdown
in ba celona and su ounding egion: Modi ica ion o diu nal cycles and limi ed ole o me eo ological condi ions in p ess.
GCOS. (2016). The global obse ing sys em o clima e implemen a ion needs. Wo ld Me eo ological O ganiza ion 200.
Gilabe , J., Ven u a, S., Segu a, R., Ma illi, A., Badia, A., Llasa , C., e al. (2021). Aba ing hea wa es in a coas al Medi e anean ci y: Wha can
cool oo s and ege a ion con ibu e? U ban Clima e, 37, 100863. h ps://doi.o g/10.1016/J.UCLIM.2021.100863
Acknowledgmen s
This wo k has been made possible hanks
o he inancial suppo o he Eu opean
Resea ch Council (ERC) Consolida o
p ojec : In eg a ed Sys em Analysis
o U ban Vege a ion and Ag icul u e
(818002-URBAG), he Spanish Minis y
o Science, Inno a ion and Uni e si-
ies, h ough he “Ma ia de Maez u”
p og amme o Uni s o Excellence
(CEX2019-000940-M), and he unding
and ecogni ion awa ded o he esea ch
g oup Sos enip a (2021 SGR 00734) by
he Depa men o Resea ch and Uni e -
si ies o he Gene ali a de Ca alunya.
This wo k also has been g an ed by he
Spanish Minis y o Science and Inno a-
ion MCIN AEI/10.13039/501100011033
unde con ac PID2020-113614RB-C21,
by he Spanish Go e nmen unde G an
PRE2018-085425 and by he Ca a-
lan Go e nmen unde con ac 2021
SGR 00574. The au ho s acknowledge
he compu e esou ces a PICASSO
and he echnical suppo p o ided
by he Uni e sidad de Málaga (G an
RES-AECT-2020-2-0004).
21698996, 2023, 24, Downloaded om h ps://agupubs.onlinelib a y.wiley.com/doi/10.1029/2023JD039497 by Readcube (Lab i a Inc.), Wiley Online Lib a y on [10/01/2024]. See he Te ms and Condi ions (h ps://onlinelib a y.wiley.com/ e ms-and-condi ions) on Wiley Online Lib a y o ules o use; OA a icles a e go e ned by he applicable C ea i e Commons License
Jou nal o Geophysical Resea ch: A mosphe es
ESTRUCH ETAL.
10.1029/2023JD039497
20 o 21
G i i hs, A. D., Pa kes, S. D., Chambe s, S. D., McCabe, M. F., & Williams, A. G. (2013). Imp o ed mixing heigh moni o ing h ough a
combina ion o lida and adon measu emen s. A mosphe ic Measu emen Techniques, 6(2), 207–218. h ps://doi.o g/10.5194/am -6-207-2013
G ossi, C., A nold, D., Adame, J. A., López-Co o, I., Bolí a , J. P., De La Mo ena, B. A., & Va gas, A. (2012). A mosphe ic
222Rn concen a-
ion and sou ce e m a El A enosillo 100 m me eo ological owe in sou hwes Spain. Radia ion Measu emen s, 47(2), 149–162. h ps://doi.
o g/10.1016/j. admeas.2011.11.006
Gue a a, M., Jo ba, O., So e , A., Pe e in, H., Bowdalo, D., Se adell, K., e al. (2021). Time- esol ed emission educ ions o a mosphe ic chem-
is y modelling in Eu ope du ing he COVID-19 lockdowns. A mosphe ic Chemis y and Physics, 21(2), 773–797. h ps://doi.o g/10.5194/
acp-21-773-2021
Holzwo h, G. C. (1964). Es ima es o mean maximum mixing dep hs in he con iguous Uni ed S a es. Mon hly Wea he Re iew, 92(5), 235–242.
h ps://doi.o g/10.1175/1520-0493(1964)092<0235:eommmd>2.3.co;2
Ins i u d’Es adís ica de Ca alunya (2023). Idesca . Re ie ed om h ps://www.idesca .ca /
Jaén, C., Udina, M., & Bech, J. (2021). Analysis o wo hea wa e d i en ozone episodes in Ba celona and su ounding egion: Me eo ological
and pho ochemical modeling. A mosphe ic En i onmen , 246, 118037. h ps://doi.o g/10.1016/j.a mosen .2020.118037
Ki z, F., Gómez-B andón, M., Ede , B., E emadi, M., Spielmann, F. M., Hamme le, A., e al. (2019). Soil ca bonyl sul ide exchange in ela ion
o mic obial communi y composi ion: Insigh s om a managed g assland soil amendmen expe imen . Soil Biology and Biochemis y, 135,
28–37. h ps://doi.o g/10.1016/j.soilbio.2019.04.005
Lee, C. L., & B imblecombe, P. (2016). An h opogenic con ibu ions o global ca bonyl sul ide, ca bon disul ide and o ganosul ides luxes.
Ea h-Science Re iews, 160, 1–18. h ps://doi.o g/10.1016/J.EARSCIREV.2016.06.005
Leelőssy, Á., Mona, T., Mészá os, R., Lagzi, I., & Ha asi, Á. (2016). Eule ian and Lag angian app oaches o modelling o ai quali y. h ps://
doi.o g/10.1007/978-3-319-40157-7_5
Lenna z, S. T., Ma andino, C. A., Von Hobe, M., Co es, P., Quack, B., Simo, R., e al. (2017). Di ec oceanic emissions unlikely o accoun
o he missing sou ce o a mosphe ic ca bonyl sul ide. A mosphe ic Chemis y and Physics, 17(1), 385–402. h ps://doi.o g/10.5194/
acp-17-385-2017
Liu, S., & Liang, X.-Z. (2010). Obse ed diu nal cycle clima ology o plane a y bounda y laye heigh . Jou nal o Clima e, 23(21), 5790–5809.
h ps://doi.o g/10.1175/2010JCLI3552.1
Lo e ane , C., & Pi inge , M. (2016). Mixing-heigh ime se ies om ope a ional ceilome e ae osol-laye heigh s. Bounda y-Laye Me eo ology,
161(2), 265–287. h ps://doi.o g/10.1007/s10546-016-0169-2
Mallik, C., Chand a, N., Venka a amani, S., & Lal, S. (2016). Va iabili y o a mosphe ic ca bonyl sul ide a a semi-a id u ban si e in wes e n
India. Science o he To al En i onmen , 551–552, 551–552. h ps://doi.o g/10.1016/j.sci o en .2016.02.014
Millán, M. M., Man illa, E., Sal ado , R., Ca a alá, A., Sanz, M. J., Alonso, L., e al. (2000). Ozone cycles in he wes e n Medi e anean basin:
In e p e a ion o moni o ing da a in complex coas al e ain. Jou nal o Applied Me eo ology, 39(4), 487–508. h ps://doi.o g/10.1175/1520-0
450(2000)039<0487:OCITWM>2.0.CO;2
Mon zka, S. A., Cal e , P., Hall, B. D., Elkins, J. W., Conway, T. J., Tans, P.P., & Sweeney, C. S. (2007). On he global dis ibu ion, seasonali y,
and budge o a mosphe ic ca bonyl sul ide (COS) and some simila i ies o CO2. Jou nal o Geophysical Resea ch, 112(D9), 1–15. h ps://
doi.o g/10.1029/2006JD007665
P o oschill-K ebs, G., & Kesselmeie , J. (1992). Enzyma ic pa hways o he consump ion o ca bonyl sulphide (COS) by highe plan s. Bo anica
Ac a, 105(3), 206–212. h ps://doi.o g/10.1111/j.1438-8677.1992. b00288.x
Que ol, X., Alas uey, A., Ruiz, C. R., A iñano, B., Hansson, H. C., Ha ison, R. M., e al. (2004). Specia ion and o igin o PM10 and PM2.5 in
selec ed Eu opean ci ies. A mosphe ic En i onmen , 38, 6547–6555. h ps://doi.o g/10.1016/j.a mosen .2004.08.037
Remaud, M., Che allie , F., Maignan, F., Bel iso, S., Be che , A., Pa ou e, A., e al. (2022). Plan g oss p ima y p oduc ion, plan espi a ion and
ca bonyl sul ide emissions o e he globe in e ed by a mosphe ic in e se modelling. A mosphe ic Chemis y and Physics, 22(4), 2525–2552.
h ps://doi.o g/10.5194/acp-22-2525-2022
Remaud, M., Ma, J., K ol, M., Abadie, C., Ca w igh , M. P., Pa a, P., e al. (2023). In e compa ison o a mosphe ic ca bonyl sul ide (T ansCom-
COS; Pa One): E alua ing he impac o anspo and emissions on oposphe ic a iabili y using g ound-based and ai c a da a. Jou nal o
Geophysical Resea ch: A mosphe es, 128(6), 1–23. h ps://doi.o g/10.1029/2022JD037817
Ribei o, I., Ma illi, A., Falls, M., Zona o, A., & Villalba, G. (2021). Highly esol ed WRF-BEP/BEM simula ions o e Ba celona u ban a ea wi h
LCZ. A mosphe ic Resea ch, 248, 105220. h ps://doi.o g/10.1016/J.ATMOSRES.2020.105220
Salamanca, F., K po, A., Ma illi, A., & Clappie , A. (2010). A new building ene gy model coupled wi h an u ban canopy pa ame e iza ion
o u ban clima e simula ions-pa I. o mula ion, e i ica ion, and sensi i i y analysis o he model. Theo e ical and Applied Clima ology,
99(3–4), 331–344. h ps://doi.o g/10.1007/s00704-009-0142-9
Segu a, R., Badia, A., Ven u a, S., Gilabe , J., Ma illi, A., & Villalba, G. (2021). Sensi i i y s udy o PBL schemes and soil ini ializa ion using
he WRF-BEP-BEM model o e a Medi e anean coas al ci y. U ban Clima e, 39, 100982. h ps://doi.o g/10.1016/j.uclim.2021.100982
Seibe , P., Bey ich, F., G yning, S. E., Jo e, S., Rasmussen, A., & Te cie , P. (2000). Re iew and in e compa ison o ope a ional me hods o he
de e mina ion o he mixing heigh . A mosphe ic En i onmen , 34(7), 1001–1027. h ps://doi.o g/10.1016/S1352-2310(99)00349-0
Se ei Me eo ològic de Ca alunya. (2023). El emps a Ca alunya. Re ie ed om h ps://www.me eo.ca /
Skama ock, W. C., Klemp, J. B., Dudhia, J., Gill, D. O., Liu, Z., Be ne , J., e al. (2021). A desc ip ion o he ad anced esea ch WRF model
Ve sion 4.3. h ps://doi.o g/10.5065/1d h-6p97
S ewa , I. D., & Oke, T. R. (2012). Local clima e zones o u ban empe a u e s udies. Bulle in o he Ame ican Me eo ological Socie y, 93(12),
1879–1900. h ps://doi.o g/10.1175/BAMS-D-11-00019.1
S ohl, A., Fo s e , C., F ank, A., Seibe , P., & Wo awa, G. (2005). Technical no e: The Lag angian pa icle dispe sion model FLEXPART e sion
6.2. A mosphe ic Chemis y and Physics, 5(9), 2461–2474. h ps://doi.o g/10.5194/acp-5-2461-2005
S ull, R. B. (1988). An In oduc ion o bounda y laye me eo ology. Kluwe Academic Publishe s.
S u m, P., Leuenbe ge , M., Si ignano, C., Neube , R. E. M., Meije , H. A. J., Langen elds, R., e al. (2004). Pe mea ion o a mosphe ic gases
h ough polyme O- ings used in lasks o ai sampling. Jou nal o Geophysical Resea ch: A mosphe es, 109(D4), D04309. h ps://doi.
o g/10.1029/2003jd004073
URBAG. (2023). Full da ase o ca bonyl sul ide URBAG measu emen campaign May and Oc obe 2020. Re ie ed om h ps://zenodo.o g/
eco ds/8072833
VAISALA (2020). In es iga ion o bounda y laye s uc u es wi h ceilome e . VAISALA. Re ie ed om h ps://www. aisala.com/en/p oduc s/
so wa e/bl- iew
Villalba, G. (2023). Full da ase o ca bonyl sul ide URBAG measu emen campaign May and Oc obe 2020 [Da ase ]. Zenodo. h ps://doi.
o g/10.5281/zenodo.8072833
21698996, 2023, 24, Downloaded om h ps://agupubs.onlinelib a y.wiley.com/doi/10.1029/2023JD039497 by Readcube (Lab i a Inc.), Wiley Online Lib a y on [10/01/2024]. See he Te ms and Condi ions (h ps://onlinelib a y.wiley.com/ e ms-and-condi ions) on Wiley Online Lib a y o ules o use; OA a icles a e go e ned by he applicable C ea i e Commons License
Jou nal o Geophysical Resea ch: A mosphe es
ESTRUCH ETAL.
10.1029/2023JD039497
21 o 21
Villalba, G., Whelan, M., Mon zka, S. A., Came on-Smi h, P.J., Fische , M., Zumkeh , A., e al. (2021). Explo ing he po en ial o using
ca bonyl sul ide o ack he u ban biosphe e signal. Jou nal o Geophysical Resea ch: A mosphe es, 126(13), e2020JD034106. h ps://doi.
o g/10.1029/2020JD034106
Whelan, M. E., Lenna z, S. T., Gimeno, T. E., Weh , R., Wohl ah , G., Wang, Y., e al. (2018). Re iews and syn heses: Ca bonyl sul ide as a
mul i-scale ace o ca bon and wa e cycles. Biogeosciences, 15(12), 3625–3657. h ps://doi.o g/10.5194/bg-15-3625-2018
Yan, Y., Li, R., Peng, L., Yang, C., Liu, C., Cao, J., e al. (2019). Emission in en o y o ca bonyl sul ide (COS) om p ima y an h opogenic
sou ces in China. En i onmen al Pollu ion, 247, 745–751. h ps://doi.o g/10.1016/j.en pol.2019.01.096
Yang, F., Qubaja, R., Ta a ino , F., Ro enbe g, E., & Yaki , D. (2018). Assessing canopy pe o mance using ca bonyl sul ide measu emen s.
Global Change Biology, 24(8), 3486–3498. h ps://doi.o g/10.1111/gcb.14145
Yi, Z., Wang, X., Sheng, G., & Fu, J. (2008). Exchange o ca bonyl sul ide (OCS) and dime hyl sul ide (DMS) be ween ice paddy ields and
he a mosphe e in sub opical China. Ag icul u e, Ecosys ems & En i onmen , 123(1–3), 116–124. h ps://doi.o g/10.1016/j.agee.2007.05.011
Yu, J., & S einbe ge , Y. (2011). Ve ical dis ibu ion o mic obial communi y unc ionali y unde he canopies o Zygophyllum dumosum and
Hammada scopa ia in he Nege Dese , Is ael. Mic obial Ecology, 62(1), 218–227. h ps://doi.o g/10.1007/s00248-011-9846-3
Zanche a, A., Kooijmans, L. M. J., an Heu en, S., Sci o, A., Schee en, H. A., Mamma ella, I., e al. (2023). Sou ces and sinks o ca bonyl sul ide
in e ed om owe and mobile a mosphe ic obse a ions in The Ne he lands. Biogeosciences, 20(16), 3539–3553. h ps://doi.o g/10.5194/
bg-20-3539-2023
Zumkeh , A., Hil on, T. W., Whelan, M., Smi h, S., & Campbell, J. E. (2017). G idded an h opogenic emissions in en o y and a mosphe ic ans-
po o ca bonyl sul ide in he U.S. Jou nal o Geophysical Resea ch, 122(4), 2169–2178. h ps://doi.o g/10.1002/2016JD025550
Zumkeh , A., Hil on, T. W., Whelan, M., Smi h, S., Kuai, L., Wo den, J., & Campbell, J. E. (2018). Global g idded an h opogenic emissions
in en o y o ca bonyl sul ide. A mosphe ic En i onmen , 183, 11–19. h ps://doi.o g/10.1016/j.a mosen .2018.03.063
Re e ences F om he Suppo ing In o ma ion
Hong, S.-Y., Kim, J.-H., Lim, J., & Dudhia, J. (2006). The WRF single momen mic ophysics scheme (WSM). Jou nal o he Ko ean Me eo o-
logical Socie y, 42, 129–151.
Seidel, D. J., Ao, C. O., & Li, K. (2010). Es ima ing clima ological plane a y bounda y laye heigh s om adiosonde obse a ions: Compa ison
o me hods and unce ain y analysis. Jou nal o Geophysical Resea ch, 115(D16), D16113. h ps://doi.o g/10.1029/2009JD013680
21698996, 2023, 24, Downloaded om h ps://agupubs.onlinelib a y.wiley.com/doi/10.1029/2023JD039497 by Readcube (Lab i a Inc.), Wiley Online Lib a y on [10/01/2024]. See he Te ms and Condi ions (h ps://onlinelib a y.wiley.com/ e ms-and-condi ions) on Wiley Online Lib a y o ules o use; OA a icles a e go e ned by he applicable C ea i e Commons License