于慧珍,馬 艷*,郭麗娜,郭飛燕,孟 赫,黃 容(.青島市氣象局,山東 青島 66003;.山東省青島生態環境監測中心,山東 青島 66003)
自《大氣污染防治行動計劃》實施以來,我國74 個重點城市2016年PM2.5質量濃度均值與2013年相比下降了35%,而臭氧濃度與污染發生頻次有所上升[1-2].研究發現臭氧污染對人體有很大的危害[3],也會對植物和生態系統造成損害,導致農林植物減產、經濟效益下降等[4-6].
臭氧污染對不同的天氣形勢非常敏感,在不同的天氣形勢下臭氧污染的來源以及影響臭氧的關鍵天氣要素不同[7-8],引入環流分型的臭氧預報模型的預報改進效果明顯[9-10].天氣形勢及其演變決定了污染物的時空分布,而氣象因子直接影響污染過程的持續時間和嚴重程度[11-12].大氣流場不僅直接輸送臭氧,其輸送攜帶的前體物也會影響局地的臭氧污染[13-15].在不同的大氣環流場下,氣象因子如輻射、溫度和相對濕度等則是驅動局地臭氧污染的關鍵因素[16].
國內對臭氧的研究區域主要集中在京津冀、長江三角洲、珠江三角洲地區,對其他地區的研究還比較少.不同地區的臭氧時空分布特征不同,影響臭氧污染的天氣因素也不同[17-18].研究發現華北臭氧污染容易發生在大陸高壓和低壓前部等強暖性天氣背景場下[19-20],而華東地區臭氧污染的典型環流形勢為副高控制和副高西北側控制[8].重慶市主城區高臭氧濃度對應每日最高溫度區間為35℃以上以及相對濕度區間70%以下[21].洪瑩瑩等[22]分析指出在珠江三角洲氣象條件變化導致臭氧濃度上升的貢獻率為29.8%,而排放的變化引起臭氧濃度下降的貢獻率為7.1%.
從全國臭氧濃度分布看,山東區域的臭氧污染比較嚴重,而且有日益增強的趨勢[23-24],山東省內不同城市間的臭氧分布也存在明顯的局地性特征[25].膠東經濟圈位于山東半島沿海地區,包括青島、濰坊、威海、煙臺、日照五市,是我國對日韓開放最前沿,21世紀海上絲綢之路與新亞歐大陸橋經濟走廊交匯的關鍵區域,沿黃省份和上合組織國家主要出海口.據統計,2020年,五市的常住人口3243 萬人,地區生產總值3.1 萬億元,分別占全省的31.9%、42.5%.前人對青島市臭氧分布特征[26]和臭氧污染VOCs 化學特征和來源進行了分析[27-28].對影響此區域臭氧污染發生的環流形勢以及關鍵天氣因子還不清楚.膠東經濟圈在不同季節受到低緯度和高緯度天氣系統的影響差異較大,且三面鄰海,不同風向時受海陸影響差異大,因此有必要對此區域影響臭氧的天氣因子進行研究.本文通過對大氣環流進行客觀分型,研究不同天氣形勢下膠東經濟圈的臭氧污染特征,找出容易發生臭氧污染的天氣形勢,并對比不同天氣形勢下影響臭氧污染的關鍵天氣要素以及氣流傳輸路徑,為膠東經濟圈臭氧污染防治工作提供參考.
臭氧實況數據為青島、煙臺、威海、濰坊和日照5 個站點的2014~2021年每天的日最大8h平均.根據環境空氣質量指數(AQI)技術規定,換算為空氣質量分指數(IAQI),當IAQI≤100 沒有臭氧污染,101~150 為臭氧輕度污染,151~200 為臭氧中度污染,201~300 為臭氧重度污染,>300 為臭氧嚴重污染.在本研究中,當5 個站點都沒有臭氧污染時定義為非污染日,有一個及以上站點出現污染時為污染日.
天氣形勢分析采用水平分辨率為0.25°×0.25°,時間分辨率為1h 的歐洲中期天氣預報中心第五代大氣再分析資料(ERA5),包括位勢高度、溫度、濕度、風、降水等.后向軌跡追蹤采用美國國家環境預報中心(NCEP)提供的全球資料同化系統(GDAS)數據,水平分辨率為1°×1°,時間分辨率為6h,包括不同層次的溫度、濕度、風等物理量.
采用的環流分型方法為旋轉T 模態主成分分析(PCT)[29-30],環流分型基于ERA5 逐日08:00(北京時)的500hPa 位勢高度場和海平面氣壓場,范圍為30°~45°N,110°~130°E.
利用斯皮爾曼等級相關分析方法來確定不同等級的臭氧污染與天氣要素的相關性,并進行顯著性檢驗.根據前人對不同城市影響臭氧的天氣因子的研究結果[17,31],選擇影響臭氧光化學反應和傳輸擴散的太陽輻射(Radiation)、降水量(Rain)、總云量(Cloud),不同高度層(850hPa、925hPa、地面)的風(緯向風U,經向風V,風速)、溫度(T)、相對濕度(Rh)、氣壓(Pres),925hPa 和2m 的溫差(T_(925hPa-2m)),邊界層高度(Blh),分別計算它們與臭氧的相關系數,找出關鍵天氣因子.
后向軌跡模式采用美國國家海洋大氣研究中心空氣資源實驗室開發的 HYSPLIT5.2.1 版本.HYSPLIT 模式可用于氣塊軌跡、傳輸、擴散、化學轉化和沉積的模擬.最常用的模型應用之一是后向軌跡分析,以確定氣團的起源和建立源-受體關系[32-33].研究顯示青島地區臭氧典型日變化峰值出現在午后13:00 時或14:00 時[26],本研究中計算每日14:00 時開始的24h 后向氣流軌跡,軌跡起始高度500m(500m 高度的風場能夠反映邊界層的平均流場特征)[34].軌跡分組采用的是HYSPLIT 模式中的聚類分析方法,根據軌跡的空間相似度對其進行分類.
2014~2021年,膠東經濟圈5 個城市中臭氧污染日最多的是濰坊428d,其次是日照198d,威海183d,煙臺167d,最少的是青島148d.各等級臭氧污染日最多的均為濰坊,重度污染只出現在濰坊.輕度污染日最少的是青島129d,中度污染日最少的是威海11d.
根據本文對膠東經濟圈臭氧污染日的定義,2014~2021年總計污染571d,占總日數的20%,其中輕度污染488d,中度污染77d,重度污染6d.污染發生在3~11月,中度污染發生在4~9月,重度污染發生在5~7月,其中污染最多出現在6月129d,其次是5月118d.臭氧污染連續出現的概率較高,持續2d 以上的污染有460d,占總污染日的81%,持續2d以上的中度污染有46d,持續2d 以上的重度污染有3d.污染最長持續時間為10d,中度以上污染最長持續時間4d,重度污染最長持續時間3d.
利用PCT 方法基于500hPa 位勢高度場和海平面氣壓場將天氣形勢分成9 種類型,type1 高空為西風槽,地面對應北高南低的形勢;type2 高空為西風槽底的平直西風,地面為鞍型場;type3、type4 和type7高空為槽后西北氣流,地面為高壓,三個類型的主要區別是西風槽和高壓的位置;type5 高空為西風槽,地面為東北低壓槽;type6 高空為平直西風,地面為高壓底部;type8 和type9 高空為弱高壓脊,地面為高壓后部,兩個類型的主要區別是高壓的位置.不同類型出現的日數和臭氧污染日數如圖1.臭氧污染主要出現在第1、4、5、9 四類天氣形勢下,這四類天氣型出現的日數占總分型日數的53%,這四類中的臭氧污染日數占總污染日的97%,全部的中度以上污染日都出現在這4 類天氣型下.

圖1 不同天氣類型出現日數以及不同等級臭氧污染日數分布Fig.1 The number of days under different weather patterns,days of different levels of ozone pollution under each weather pattern
對比type1 的污染日和非污染日的天氣形勢,污染日半島位于西風槽后(圖2a),非污染日半島位于西風槽前(圖略).污染日南方低壓較弱,受弱反氣旋控制,半島為弱南風或弱北風(圖2e).非污染日南方低壓較強,受南方低壓外圍氣流的影響,半島為偏東風.這一天氣型下的臭氧污染出現在春夏秋三季,其中8~9月最多(圖3).這一天氣型的臭氧污染日占總污染日的28%,污染持續2d 以上的有72d,占這一型污染日的45%,最長持續5d(表1).

表1 不同天氣型的臭氧污染日比例、持續污染概率和污染最長持續時間Table 1 Proportion of ozone pollution day,probability of persistent pollution,and longest duration of pollution under different weather patterns

圖2 不同天氣型的臭氧污染日合成平均的環流形勢Fig.2 Synthetic circulation patterns on ozone pollution days under different weather types

圖3 不同天氣型臭氧污染日數的月分布Fig.3 Monthly distribution of ozone pollution days under different weather types
對比type4 污染日和非污染日的天氣形勢,高空均為西風槽,地面為高壓,污染日西風槽和地面高壓的位置(圖2b,f)比非污染日偏東.污染日有兩個高壓中心分別位于內蒙古和山東半島以南,半島為偏西風或偏北風(圖2f),非污染日的高壓位于內蒙古地區,半島為西北風.這一天氣型下的臭氧污染出現在春秋兩季,其中4~5月最多(圖3).臭氧污染日占比5%,污染持續2d 以上的有4d,占這一型污染日的13%,最長持續2d(表1).
對比type5 污染日和非污染日天氣形勢,主要差別是西風槽的位置和東北低壓、海上高壓的強度.污染日半島位于西風槽后(圖2c),非污染日半島位于西風槽前.污染日地面東北低壓的強度比非污染日強,海上高壓強度比非污染日弱,污染日為西南風(圖2g),非污染日半島為南風.這一型出現在春夏秋3個季節,其中5~7月最多(圖3).臭氧污染日占比58%,污染持續2d以上的有243d,占這一型污染日的73%,最長持續7d(表1).
type9 污染日受西風槽后西北氣流影響(圖2d),非污染日高空受弱脊控制.污染日和非污染日半島均位于海上高壓后部,為偏南風(圖2h).這一天氣型下臭氧污染出現在春秋兩季,其中3~4月最多(圖3).臭氧污染日占比6%,污染持續2d以上的有13d,占這一型污染日的41%,最長持續3d(表1).
綜上所述,同一天氣型下是否發生臭氧污染主要與半島相對于西風槽的位置、西風槽的強度、地面高低壓系統的位置和強度有關.4 類天氣型下臭氧污染都發生在高空槽后西北氣流控制下,主要區別是地面天氣系統,type1 為弱反氣旋環流控制,type4為大陸高壓控制,type5 為東北低壓槽影響,type9 為海上高壓后部.為了方便表述,后文中type1、type4、type5 和type9 分別稱為弱反氣旋型、大陸高壓型、東北低壓槽型和海上高壓后部型.
根據影響臭氧光化學反應和污染擴散的可能天氣要素[17,34],選擇代表不同高度的溫度、相對濕度、風場、穩定度、降水、云量的天氣要素進行相關性計算.如表2所示,不同天氣型下,925hPa 的溫度(T_925hPa)與臭氧呈顯著正相關,即925hPa 溫度越高,越有利于臭氧污染的發生.從水平分布看,污染日925hPa 存在暖舌自內陸伸向半島地區,污染強時對應的暖舌也更強(圖略).溫度與太陽輻射相互關聯,當太陽輻射強時,氣溫升高,為光化學反應產生臭氧提供有利條件,有利于臭氧污染的發生.其他要素與臭氧的相關關系,不同天氣型存在明顯的差異.

表2 不同天氣型下,不同天氣要素區域平均的日平均值與臭氧的相關系數Table 2 Correlation coefficient between area-averaged daily mean values of different weather elements and ozone under different weather patterns
弱反氣旋和東北低壓槽天氣型下,臭氧與不同天氣要素的相關關系一致,與濕度的相關比與溫度的相關高(表2).與光化學反應相關的要素中,臭氧與太陽輻射(Radiation)、925hPa 溫度(T_925hPa)和2m最高溫度(Tmax_2m)正相關,與降水(Rain)、云量(Cloud)、850hPa 相對濕度(Rh_850hPa)、925hPa 相對濕度(Rh_925hPa)和2m 的相對濕度(Rh_2m)負相關,即太陽輻射越強、溫度越高、云量越少、濕度越小,越有利于臭氧污染的發生.太陽輻射、溫度、云量和相對濕度是相互關聯的,相對濕度低時大氣中水汽凝結成云和產生降水的可能性小,到達地表的太陽輻射高,溫度升高,為光化學反應產生臭氧提供有利條件,有利于臭氧污染的發生.與擴散傳輸條件相關的要素中,臭氧與低層緯向風(U_850hPa,U_925hPa,U_10m)、925hPa 和2m 溫差(T_(925hPa-2m))正相關,與850hPa 和925hPa 經向風(V_850hPa,V_925hPa)負相關,即低層偏西風分量越大、偏南分量越小,低層大氣越穩定,越有利于臭氧污染的發生.大氣穩定時,有利于臭氧的積累.不同風向對臭氧的影響與天氣形勢場有關,弱反氣旋天氣型下非污染日受低層南側低壓影響,半島為偏東南風,污染日受弱反氣旋控制,半島為偏西北風;東北低壓槽天氣型下的非污染日低層為西南氣流,污染日偏南分量減小,為西西南氣流.雖然弱反氣旋和東北低壓槽天氣型下各氣象要素與臭氧相關性基本一致,但是影響臭氧的關鍵要素存在一些差別,弱反氣旋天氣型下相關性最高的要素是太陽輻射、925hPa 相對濕度、云量、925hPa 緯向風和10m 緯向風,包括與光化學反應條件和擴散傳輸條件相關的要素,東北低壓槽天氣型下相關性最高的要素是太陽輻射、云量、三層的相對濕度(850hPa、925hPa 和2m),均是光化學反應條件相關的要素.
大陸高壓和海上高壓后部天氣型下,臭氧與溫度的相關比與濕度的相關高(表2).大陸高壓天氣型下臭氧與光化學反應相關的要素850hPa 和925hPa溫度(T_850hPa、T_925hPa)呈正相關,即溫度越高,越有利于臭氧污染的形成.臭氧和相對濕度相關的物理量不顯著相關,這是因為大陸高壓型為槽后西北氣流控制,整體比較干,云量少,污染日和非污染日的差異不明顯.與擴散傳輸條件相關的經向風(V_850hPa,V_925hPa,V_10m)正相關,即偏北風分量越大,越不有利于臭氧污染發生.與925hPa 和2m 溫差(T_(925hPa-2m))正相關,與邊界層高度(Blh)負相關,即低層大氣越穩定、邊界層高度越低,越有利于臭氧污染的發生.海上高壓后部天氣型下臭氧與光化學反應相關的要素太陽輻射(Radiation),925hPa溫度(T_925hPa),2m 最高溫度(Tmax_2m)正相關,與云量和925hPa 相對濕度(Rh_925hPa)負相關,即太陽輻射越強、溫度越高、濕度越低,越有利于臭氧污染的形成.海上高壓后部型下的臭氧與擴散傳輸條件相關的要素不顯著相關.
分別計算白天平均(08:00~20:00 時)和夜間平均(20:00~08:00 時)的各天氣要素與臭氧的相關,發現白天的太陽輻射、溫度、濕度、云量、降水與臭氧的相關性明顯高于夜間,進一步說明這些量是與光化學反應有關,主要是白天的光化學反應條件影響臭氧污染發生.而白天的風和穩定度與臭氧的相關系數與夜間差異不明顯,且不同類別結果不同,說明白天和夜間的擴散傳輸都會影響臭氧污染發生.
850 hPa、925hPa、地面2m 三層的溫度(T)和相對濕度(Rh)與臭氧的相關性一致,其中925hPa 的相關性最高,一方面可能是因為對流層低層臭氧濃度極值發生在925hPa 附近[35-36],另一方面地面氣象要素受海陸分布影響較大,對于個別沿海站點,臭氧與地面溫濕不顯著相關(單個站點臭氧與氣象要素的相關結果略),因此臭氧與區域平均的地面溫濕的相關小于925hPa.與降水的相關性低于與云量的相關,可能是因為降水數據是離散的,站點數據的代表性較差.除了表2 中顯示的要素,其他要素與臭氧的相關結果顯示,臭氧與850hPa、925hPa 和地面的氣壓值均不顯著相關.同一高度層緯向風(U)或經向風(V)與臭氧的相關比風速的相關大,除了弱反氣旋型外其他三種天氣型下的臭氧與850hPa、925hPa 和10m 風速都不相關,說明臭氧污染主要受風向的影響.
為了更直觀的對比低層風對于臭氧污染的影響,并找出影響膠東經濟圈臭氧污染的氣流來源,對不同天氣型的污染日和非污染日分別進行后向氣流軌跡追蹤,然后對軌跡進行聚類,對比同一天氣型的污染日和非污染日,以及不同天氣型污染日的氣流軌跡特征,結果如下:弱反氣旋型的污染日氣流主要由江蘇北部近海(31%)、山東西北部(20%)和半島北部(27%)輸送.非污染日氣流主要由東南海上輸送(51%)和半島北部近海(38%)輸送(圖4a).污染日和非污染日氣流軌跡主要區別是來源于內陸或海上,這和臭氧與低層緯向風、濕度的相關性一致.

圖4 不同天氣型在污染日(黑色)和非污染日(灰色)的24h 后向軌跡聚類Fig.4 Cluster analysis of the 24h backward trajectory for pollution day(black)and no-pollution day(grey)
大陸高壓型的污染日氣流主要由山東中部(39%)、天津(23%)和安徽北部(23%)輸送.非污染日氣流主要由河北北部、蒙古及以北(65%)輸送和山東中部(35%)輸送(圖4b).污染日和非污染日氣流軌跡的主要區別是氣流源地的緯度,污染日氣流源地主要在40°N 以南,非污染日源地主要在40°N 以北,這和臭氧與低層經向風、溫度的相關性一致.
東北低壓槽型的污染日氣流主要由山東西部(19%)、安徽北部(22%)、江蘇南部(21%)和青島近海(30%)輸送.非污染日氣流由較遠的東海(33%)、南部海上(25%)、安徽西部(22%)和河北東北部(20%)輸送(圖4c).污染日和非污染日氣流軌跡主要差別是來源于偏西內陸或東南海上,這和臭氧與低層緯向風、濕度的相關性一致.
海上高壓后部型的污染日氣流主要由江蘇北部(34%)、安徽中部(22%)、江蘇南部沿海(22%)和半島東部(22%)輸送.非污染日氣流主要由東南海上輸送(35%),南部近海(31%)和安徽北部(34%)輸送(圖4d).污染日和非污染日氣流軌跡的主要區別是來源于內陸或海上.
污染日的氣流主要來源于內陸和南部近海.其中內陸來源主要是本省和臨近省市,山東省內影響最大,平均占比32%,江蘇省占19%,安徽省占17%,京津冀占8%.山東近海和江蘇近海占16%.
2014年5月26~31日膠東經濟圈連續6d 發生臭氧污染,其中29~31日為臭氧重度污染.膠東經濟圈以及周邊總計9 個城市的臭氧觀測結果顯示,5月25日膠東經濟圈無污染,西南內陸的棗莊出現輕度污染,26~28日單站開始陸續出現污染,29~31日出現大范圍污染,其中濰坊和東營發生重度污染,6月1日污染過程結束(圖5a).這一次臭氧污染過程最早開始且污染最嚴重的地區位于膠東經濟圈的西北和西南,東南沿海地區開始晚且污染輕.下面對這一連續臭氧污染過程進行天氣形勢和關鍵氣象要素的分析,研究連續臭氧污染發生的天氣機制.

圖5 臭氧污染指數和氣象因子隨時間的變化Fig.5 Time series of ozone pollution index and weather elements
按照本文的天氣分型,5月25~31日為type5(污染日為東北低壓槽型),6月1日為type1(污染日為弱反氣旋型).對于半島地區,5月25日位于高空500hPa 西風槽前,隨著西風帶系統的東移,26~30日受西風槽后的西北氣流影響,31日受高壓脊控制,6月1日再次受到西風槽前西南氣流影響(圖6a~d).這一次臭氧污染發生在高空槽后西北氣流或者高壓脊控制下,在臭氧污染發生前和結束時受西風槽前西南氣流影響,與2.2小節中污染日和非污染日500hPa 形勢對比結果一致.從地面形勢看,5月25~30日半島地區受內陸西北或西南氣流影響,5月31日開始受到偏南氣流影響,6月1日受東南海上氣流影響(圖6e~h),這與2.4 小節中污染日和非污染日氣流來源對比結果一致.

圖6 一次臭氧污染過程的環流形勢Fig.6 Circulation patterns during a ozone pollution process
從光化學反應條件看,5月25日~6月1日未發生降水,是否產生臭氧污染主要差異是太陽輻射、低層溫度和濕度(圖5b).5月25日,受西風槽前輻合上升的影響,膠東經濟圈的925hPa 相對濕度達到80%,云量達到7 成,太陽輻射、低層925hPa 和地面2m的溫度比26~31日的低,未發生臭氧污染.26~30日,受西風槽后西北氣流影響,925hPa 相對濕度明顯減小,日平均值小于30%,云量低于4 成,太陽輻射增強,低層925hPa 和地面2m 溫度升高,為光化學反應產生臭氧提供有利條件.31日,受高壓脊控制,雖然云量增加至9 成,但是低云量仍小于1 成(圖略),太陽輻射和低層溫度仍較高,有利于臭氧污染的維持.6月1日,受西風槽前西南氣流影響,低層925hPa 相對濕度增加至50%,太陽輻射減弱,低層溫度下降,光化學反應條件變差.從擴散傳輸條件看,5月26~31日925hPa與2m 溫度差大于5月25日和6月1日(圖5b),說明臭氧污染發生時低層空氣穩定度較大,有利于臭氧污染的積累.5月25~30日,低層(925hPa/地面)受西北或西南氣流影響,有利于內陸地區的臭氧向沿海地區輸送,污染面積增大.31日開始受弱的偏南氣流影響,6月1日東南風加強,受海上氣流影響,臭氧污染過程結束.
上述分析結果說明,這一次連續臭氧污染過程發生在高空西風槽后或高壓脊控制下,太陽輻射強、濕度低、云量少有利于低層溫度升高為光化學反應產生臭氧提供有利條件,低層受內陸西北或西南氣流影響,臭氧污染范圍擴大.當高空受西風槽前西南氣流影響時,濕度和云量增加使得太陽輻射減弱、低層溫度降低不利于光化學反應產生臭氧,低層受到海上東南氣流影響,臭氧污染過程結束.
3.1 2014~2021年膠東經濟圈臭氧污染日總計571d,占總日數的20%.臭氧污染在3~11月都有可能發生,其中5~6月最多.發生2d 以上持續污染的概率較大,占總污染日的81%.
3.2 產生臭氧污染的天氣型主要為4 類,高空均為西風槽后西北氣流,地面分別對應弱反氣旋型(type1),大陸高壓型(type4),東北低壓槽型(type5),海上高壓后部型(type9),臭氧污染日占總污染日的比例分別為28%,5%,58%和6%,弱反氣旋型和東北低壓槽型出現在春夏秋三季,大陸高壓型和海上高壓后部型出現在春秋兩季.
3.3 4類天氣型中臭氧都與925hPa 溫度顯著正相關.不同天氣型下,與臭氧相關的關鍵天氣要素不同.弱反氣旋型和東北低壓槽型中臭氧與濕度的相關比溫度高,大陸高壓型和海上高壓后部型中臭氧與溫度的相關比濕度高.弱反氣旋型的關鍵因子為太陽輻射,云量,925hPa 相對濕度,925hPa 和10m 緯向風.東北低壓槽型的關鍵因子是太陽輻射,云量,850hPa、925hPa 和2m 的相對濕度.大陸高壓型的關鍵因子為850hPa 和925hPa 溫度,925hPa 和10m 經向風,大氣穩定度.海上高壓后部型的關鍵因子為太陽輻射、925hPa 溫度、2m 最高溫度、云量和925hPa相對濕度.
3.4 從氣流后向軌跡看,大陸高壓型下是否發生臭氧污染主要與氣流來源的緯度有關,當氣流來源于40°N 以南時,容易發生污染.其他三個天氣型主要與氣流來源于內陸還是海上有關,當氣流來源于內陸時,容易發生污染.
3.5 膠東經濟圈臭氧污染日的氣流主要來源于內陸和南部近海.其中山東省內最多,占比32%,其次是江蘇省占19%,安徽省占17%,山東近海和江蘇近海占16%.
3.6 一次連續臭氧污染過程分析顯示,當連續多日受高空西風槽后西北氣流或高壓脊控制,低層受內陸西北或西南氣流影響時,有利于膠東經濟圈臭氧污染的發生和持續;當受到高空西風槽前西南氣流和低層東南海上氣流影響時,臭氧污染過程結束.