Backup of Projects(No. 19)

 
 
 

波およびサブメソスケールダイナミクスを考慮した黒潮海域における海洋構造の再現---超高解像度領域海洋モデルとスペクトル波浪モデルの連成並列計算による地球環境研究のパラダイムシフト †

 

　黒潮は北太平洋亜熱帯循環の西端に形成される西岸境界流であり，その流量は毎秒約3000万m3（全世界の総河川流量の30倍に相当）にも達する．北緯15度付近を西進する北赤道海流はフィリピン沖で南進するミンダナオ海流と北進する黒潮に分岐し，黒潮は陸棚に沿って琉球列島西側を北上，トカラ海峡で方向を転じ，日本南岸へ至る．日本沿岸では，四国・紀伊半島沖から伊豆海嶺の間において10数年に一度大きく離岸する「大蛇行」によって特徴付けられるように，頻繁にその流路を大きく変動させる．黒潮の主流は表層500m程度までに集中しており，傾圧性の強い流れであるため，蓄積された位置エネルギーが運動エネルギーへと開放される傾圧不安定によって大蛇行が生じるものと考えられる．流路パターンは概ねその平均流量によって規定されており，低流量で直進，高流量で蛇行，中流量では多重平衡状態となり，その両者の間をカオス的に行き来する．特に中流量時に遠州灘沖で観測される黒潮流路変動はメソスケール渦（準地衡流バランスによる中規模渦）によると考えられている．このように様々な時空間スケールで変動する黒潮流路であるが，流路変動のメカニズムは未だ解明されていない．特に急峻な沿岸・海底地形との相互作用や，波浪などによる小スケール（サブメソスケール）での運動量の水平・鉛直混合過程の影響など，スケール間の相互作用の解明が急務となっている．

　一方，海洋流動に対する波の影響には，Stokes driftとコリオリ力の相互作用，それに伴う海洋表層Ekmanバランスの改変，海洋渦とStokes driftとの相関，砕波による波から流れへの運動量輸送，Langmuir循環による表層混合の促進などがある．反対に波動場は，流れによってドップラーシフトと屈折変形を受けるため，平均流や海洋渦を変化させるフィードバック作用を持っている．このように，波は地球規模での海洋のダイナミクスや海洋乱流のエネルギーバランスを考える上で潜在的に重要な役割を果たしているにも関わらず，McWilliams et al.（2004）によって波形勾配に関する弱非線形近似などの仮定の下に位相平均されたprimitive方程式が整理されるまではこれらを記述する術がなかったため，一部のごく単純化された理論や要素実験を除くと，統一的に検討されることはなかった．これに対し我々は，McWilliams理論を強非線形波のregimeまで拡張し，上記のような波の影響を全て考慮した一般的な理論体系を構築し，海洋モデルROMSに組み込むことに成功した（Uchiyama et al., 2010）．同じ北太平洋亜熱帯循環の東岸境界流であるCalifornia海流に対する予備的な実験（Uchiyama et al., 2010, Ocean Sciences Meeting）によれば，流れに対する波の影響は顕著であり，平均流速を最大で30％強化し，海洋－大気界面プロセスを左右する混合層厚さを20％程度変化させるなどの事実が明らかになっている．流量も気象変動もより大きい西岸の黒潮海域では，さらに強い影響を受けることが予想される．

　以上のような背景および予備的な検討成果を踏まえ，本研究では，最先端の領域海洋モデルROMSおよびスペクトル波浪モデルを用い，マルチスケール多段ネスティングにより，グローバル海洋データからスタートして西太平洋全域（解像度10km），東アジア沿岸（解像度3km），および日本沿岸（平解像度1km）へ4段階でダウンスケールすることにより，グローバルからサブメソスケールまでの現象を再現可能な超高解像度シミュレーションを行い，海盆規模（basin scale）の流れと波動，複雑地形，内部潮汐，メソ・サブメソスケール渦や乱流などを全て考慮した形で黒潮の統計的特徴を再現・抽出する．そして波浪を考慮した場合としない場合の双子実験を行い，黒潮流路変動に代表される海洋流動と波浪との双方向的な相互作用について定量的に評価することを試みるものである．鍵となるのは超高解像度モデルによるサブメソスケール現象の再現と，波－流れ相互作用を考慮する点に集約され，最新の理論と計算技術を駆使することにより黒潮流路動態研究に新たなブレークスルーをもたらすことが大いに期待される．

 

surface vorticity in the Kuroshio xtentionsurface vorticity in the Kuroshio Extention

 
Figure: （左）高解像度モデルROMS-L2（解像度1km）の計算領域．（右）黒潮続流域における海洋表層の温度（上）と無次元渦度（下）の瞬間像．右の図は黒枠をそれぞれ拡大している．流軸−1度程度の領域でサブメソスケール渦が多数発生していることが分かる．
 
 
↑

瀬戸内海全域を対象とした高解像度海洋モデルの開発に基づく外洋影響を考慮した広域流動の動態把握と内湾環境保全戦略 †

 

　瀬戸内海では，東京湾や伊勢湾と同様に，高度経済成長期以降の環境悪化に伴い海域に流入する窒素やリンの総量規制を行ってきたが，近年の海域水質は横ばい，あるいは底層酸素量などの項目ではむしろ悪化傾向にあると報告されている．瀬戸内海は種々の有用水産生物が生息している生産力の高い海域であるが，水質低下とシンクロして水産資源の減少が懸念されており，漁獲量は1987年頃から減少傾向にあり，2000年以降はピーク時（1982年）の半分程度で推移している．その主たる原因は海域の富栄養化や藻場干潟の消失などによる生息環境の悪化であるとされ，水産資源の回復と適切な管理が望まれる．そのためには海域環境を適切に評価し，最適な環境保全戦略を策定する必要があり，アセスメント・海況予報ツールとしての瀬戸内海の物質輸送と流れを精緻に再現する数値流動モデルを開発することが急務となっている．

　瀬戸内海は多くの海峡と複雑な地形，大小様々な島嶼を有するため，これらを正確に表現することが流動モデリング成功への鍵となる．また，内湾の物質拡散や乱流エネルギー収支において主要な役割を果たすサブメソスケール渦は1～数km程度の空間規模であるため，このような沿岸域の海洋構造を正確に再現するためには，最低でも500m程度以下の水平解像度で評価しなければならない．一方で，瀬戸内海の平均的な流動は黒潮流路変動の影響を強く受けていること（駒井ら，2009），湾内を通過する流量は約10年周期の変動し，多くの年では豊後水道から流入して紀伊水道から流出すること（藤原ら，2006），その流量は瀬戸内海に流入する全一級河川の総流量の10倍にも達することなどが分かりつつある（中山ら，2009）．太平洋全体の風成循環の一部である黒潮流路もまた，10年程度の周期で間欠的に直進から大蛇行へレジームシフトすることがよく知られており，結局，瀬戸内海の流動は，内因的な潮汐残差流や密度流に加え，黒潮などを通じてグローバルな環境変動の影響を強く受けて形成されている．

　本研究ではまず，瀬戸内海環境アセスメントのための中核ツールとして，四国沖太平洋側から瀬戸内海にかけての広範な海域に対応した新世代の超高解像度3次元海洋流動モデルを開発する．次いで，2004年の黒潮大蛇行時を含む最近15年間分の再解析を行い，瀬戸内海の広域流動の中長期的な動態把握，瀬戸内海を構成する湾・灘（大阪湾，播磨灘，広島湾等々）間の海水交換と平均滞留時間，稚魚類・動植物プランクトンの行動モデルを組み込んだラグランジュ粒子の3次元追跡を行い，瀬戸内海内部の物質・生物分散と，紀伊水道，豊後水道，関門海峡を通じた外洋との海水交換機構を解明することを主たる目的とする．

　具体的には，データ同化された西太平洋海域渦解像海洋モデルJCOPE2（Miyazawaら，2009）から開始し，ROMS（Regional Oceanic Modeling System，Shchepetkin & McWilliams，2005）を用いて合計3段階でダウンスケールする．つまり，JCOPE2（解像度約10km）→ROMS-L1（同2km） →ROMS-L2（同500m）と順次入れ子状に領域を収斂させ，外洋からの地球規模環境シグナルを取り込みつつ，微細な地形特性や内因的なサブメソスケール渦などの乱流をも詳細に再現することの可能な枠組みを構築する．瀬戸内海流動に対して本質的な重要性を持つ河川からの淡水流入も当然ながら考慮しつつ，15年間の再解析を実施する．また，2年目以降には，ラグランジュ粒子の時間逆行計算と水質－生態系モデルとのカップリングを行い，浅海生物の起源および黒潮を考慮した瀬戸内海生物ネットワーク形成について調査，流れの相関を評価することにより，戦略的な内湾環境保全計画の策定に直接貢献することをプロジェクト全体のゴールとして設定している．

 

vort_L1.jpgvort_L2.jpg

 
Figure: 瀬戸内海における表層無次元渦度分布．（左）低解像度モデル（ROMS-L1，解像度2 km），（右）高解像度モデル（ROMS-L2，解像度600 m）．
 
 
↑

Vortex force理論に基づく新しい位相平均プリミティブ方程式による海浜流場の3次元構造の解明---海岸・海洋工学の新展開 †

 

　沿岸域流動の時空間的構造の実態把握と定量評価は，地球物理学的には海洋循環の境界条件を与え，土木工学的には，a) 津波や高潮など沿岸災害の評価・予測，b) 海岸侵食対策，c) 海洋構造物の性能設計，d) 生態系や環境の保全技術の確立等に対して第一義的に重要である． Michael Longuet-Higginsによって1950年代後半から1960年代に提案されたradiation stressは，風波によって駆動された沿岸の流れをグロスに記述しうるほぼ唯一の理論であった．弱非線形波，WKB近似などを駆使し簡便な数学的表記に帰着させることにより，工学分野では一定の成功を収めていたが，一方では3次元への拡張性が極めて低く，密度成層や地球自転の効果が欠如しているなど，著しく一般性の欠けた概念であった．これに対して我々は，海洋表層で生じるLangmuir循環を記述するために1970年代にCraikとLeibovichによって導出された，Navier－Stokes方程式非線形項のHelmholtz分解に基づくvortex force理論の拡張性の高さに着目し，マルチスケール漸近展開を用いて波動に関して位相平均されたプリミティブ方程式を導出し，領域海洋循環モデルROMSに組み込んだ（Uchiyamaら，2009；2010）．この理論とモデルは，原理的に海洋物理学に求められるほぼ全ての現象を統一的に包含しているため一般性が非常に高く，精度的にもradiation stressをはるかに凌駕し，なおかつこれまで別個に扱われていた沿岸流動と地球規模の海洋流動とをシームレスに接合することが可能な画期的なツールであり，その応用を通じて海岸工学，海洋物理学分野にパラダイムシフトをもたらすことが期待される．

　本研究では，上記の枠組みを応用し，海岸工学分野における典型的な流体力学的問題である海浜流系統の発達過程を精緻に調査し，今まで全く把握されていなかった海浜流の3次元像を明らかにすることを主たる目的とする． 3次元海浜流場では，1）砕波による岸向き表層流とそれを補償する沖向き低層流に代表される鉛直循環流の形成，2）沿岸方向圧力勾配に伴う1）の変形と流れの3次元化，3） 1）の機構に伴う鉛直シア不安定の発生，4）渦度のストレッチ作用などによるエネルギーカスケードの促進， 5） vortex forceによる岸沖方向のstreak構造の発生，などが新しい現象として発見されるものと予想している．これらはいずれも本研究の枠組みがあって初めて再現・解析可能なものであり，本研究の海岸・海洋工学分野に対する学術的インパクトは計り知れない．

 
3D rip current
 
Figure: rip channel周辺に形成された非定常3次元離岸流の瞬間構造．カラー：barotropic渦度，黒矢印：表層流速，赤矢印：底層流速．
 
 
↑

福島第一原発を放出源とした放射性セシウムの沿岸海洋分散に関する研究 †

 

　2011年3月11日に発生した東日本大震災に伴う巨大津波により，東京電力福島第一原発の原子炉3機で水蒸気爆発が生じ，3月12日から現在に至るまで，炉心や燃料棒の冷却のために海水や淡水が営々と注入されている．注水された水は放射能汚染水となって施設内に貯まり続け，再循環システム稼働以前には重大な漏水事故が惹起された．判明しているだけで4月16日に2号機から4,700 TBq，4月4〜10日に緊急計画放出により 150 GBq，5月10〜11日に3号機から 20 TBq（いずれも推定値）の放射性物質が海洋へ放出されたと報告されている．

　電中研・JAMSTEC/文科省・仏 SIROCCO による先行研究は，いずれも放出後の放射性物質の海洋での分散は速やかに行われ，その大半は黒潮続流により北太平洋中緯度域方向へ輸送されることを示唆していた．しかしながら，沿岸域の流れは海岸海底地形，河川流出，波浪等の影響を受けるため，準地衡流的な沖合流動によるものとは異なる分散パターンを生じる可能性がある．特に半減期 30 年のセシウム137（137Cs）は海岸周辺の底質に吸着することにより沿岸域を再循環し，中長期にわたって海洋汚染を引き起こす可能性も否定できない．本研究では，領域海洋循環モデル ROMS（Shchepetkin & McWilliams, 2005）を用いたダウンスケーリング高解像度数値実験により沿岸域での放射性物質の分散特性を定量化することを最終的な目的とする．これまでのところ，2段階のネスティングにより水平解像度 1 km までのダウンスケーリングに成功，その結果を用い，再現性の確認と沿岸域での放射性物質（137Cs）濃度出現パターンを解析している．

 
Fukushima Dispersal
 
Figure: 福島第1原発を放出源とした放射性セシウムの海洋分散．放出後13日目の海洋表層の流速分布とセシウム濃度の水深積分値．（左）河川を考慮したケース，（右）考慮しないケース．
 
 
↑

南カリフォルニア湾の沿岸域流動・物質分散に対するサブメソスケール乱流および波浪の影響に関する研究 †

 

　米国西海岸・南カリフォルニア湾の中央部付近に位置するサンタモニカ湾（SMB）は，後背地に人口密集都市ロサンゼルスを抱え，人間活動の影響を色濃く受ける半開放性海域である．SMB は沖合を南進するカリフォルニア海流と，岸近くの中層部を北進するカリフォルニア反流から構成されるカリフォルニア海流システム（CCS）の縁辺部に位置し，春期の北西風による沿岸湧昇や，陸域からの負荷の影響を受ける生産性の高い海域である．この SMB 中央部の沖合約 8km（5 mile）地点の海底には Hyperion Treatment Plant（HTP）からの 2 次処理水排水口があり，平均毎秒約 15 トンの処理水が恒常的に放出されている．ロサンゼルス市では老朽化に伴う排水口の更新を予定しており，排水口位置を沿岸約 1.6km（1 mile）の地点に移設することを第 2 案として検討している．本研究では，領域海洋循環モデル（ROMS；Shchepetkin andMcWiilams, 2005, Ocean Modell.）をベースに，外洋からの遠隔シグナルや，地形性・フロント不安定によるサブメソスケール渦などを取り込むべく，多段ネスティングを用いたダウンスケーリング実験を行い，高解像度（水平解像度約 75m） SMB モデルにより処理水の分散プロセスを解析し，排水口位置の違いによる海洋構造の変化や分散パターンの相違などについて検討した．

 
Santa Monica Bay Hyperion sewage effluent
 
Figure: 米国Santa Monica湾における下水処理水をトレーサーとしたサブメソスケール海洋分散．（上）排水口設置位置は1マイル沖（diversion），（下）5マイル沖（ノーマル）．
 
 
↑

ハドソン川における3次元非定常密度成層流の乱流特性とサーマルプリュームの時空間的応答に関する研究 †

 

　米国ニューヨーク州とニュージャージー州の境界を南下するハドソン川は，河床勾配が極めて緩いためニュージャージー湾の影響を強く受ける感潮域が最上流まで続く，典型的な大陸河川である．本研究では，ハドソン川中流の約80kmの区間に対して地球自転と密度成層の影響を考慮した3次元静水圧モデリングを行い，Indian Point周辺における非定常3次元流れおよびそれに伴う乱流構造を抽出するととみに，Indian Point左岸側に位置する火力発電所による河川水の取水・温排水に伴うサーマルプリュームの非定常分散過程に関する解析を実施している（神戸大）．また，LES/Fluentシステムによる小領域の定常非静水圧モデリングを併せて行い，POD解析により典型的な物理モードの抽出を試みている（立命館大学）．

 
Hudson Estuary Plume
 
Figure: ハドソン川中流域表層における流速，相対渦度，水温，発電所をソースとした仮想無次元トレーサー濃度（対数スケール）の瞬間像．
 
 
↑

紀伊半島沿岸における台風通過に伴う海洋構造変化に関する研究 †

 
  • 担当：西井達也
  • 共同研究相手：京都大学防災研究所

Coming soon.

 
 
↑

琉球諸島周辺海域における黒潮系水塊の波及効果に関する研究 †

 
  • 担当：上平雄基
  • 共同研究相手：沖縄科学技術大学院大学（OIST）

　琉球諸島西岸沖約200kmの東シナ海大陸棚に沿って北上する黒潮暖流は浮遊幼生，栄養塩等の輸送に大きな役割を果たしているため，この海域の生態系保全を考える上では黒潮の波及効果を正確に把握することが重要となる．黒潮は琉球諸島という明確な幾何的境界条件の影響を受けて独特の乱流場が形成される可能性が予想される．よって，本海域を対象とした精緻な海洋モデリングを行い，観測データ等と相互補完しながら海域の流動環境に対する理解を深化させることが急務である．
　本研究では，まず，領域海洋循環モデルROMS（RegionalOcean Modeling System）を用いた2段階のネスティングにより水平解像度約10km➡3km➡1kmとしてダウンスケーリングを行い，黒潮の波及効果を解析し得るサブメソスケール解像海洋モデリングシステムを開発する．出力モデル結果と現地観測データや衛星データとの比較を行い，黒潮流路，メソスケール変動渦運動エネルギー，黒潮流路上の鉛直断内の密度構造などに関する再現性を確認後，琉球諸島の中間海域におけるサブメソスケール乱流の構造を解析し，黒潮系水塊の波及効果を定量化することを試みた．

 
 
↑

大阪湾・播磨灘海域における海洋物質分散プロセスに関する高解像度モデリング †

 
  • 担当：神吉亮佑
  • 共同研究相手：建設工学研究所，神戸市下水道局，海洋研究開発機構（JAMSTEC）

　大阪湾・播磨灘海域の水質は，その巨大な後背地人口のため，河川水や下水処理水の影響を特に受ける.しかし，下水処理水の詳細な分散過程に関しては，未解明な部分が多く残されている．
　大阪湾・播磨灘は狭隘な海峡に囲まれた閉鎖性海域であることに加え，埋め立て地などによる複雑な海岸線を有するため，形成される流動場も複雑な構造を持っている．一方，栗山（2013）によれば，瀬戸内海全域の水収支は四国沖を東進する黒潮流路の季節変動の影響を強く受けつつも，年間を通じて豊後水道から紀伊水道方向への時計回りの流れが卓越することが分かっている．また，淀川や大和川などの一級河川からの浮力の影響により，特徴的なエスチュアリ循環が形成されることが知られている．したがって，潮汐，海上風や河川などに代表される外力条件を正確に考慮することに加え，瀬戸内海全体の流れと複雑な海岸線を同時に表現することが，本海域の流動・物質分散モデリングの成功の鍵となる．そのためには，瀬戸内海全域モデルからのダウンスケーリングによる高解像度モデルが必要となる．
　以上のことから，本研究では，①大阪湾・播磨灘を対象として多段ネスティングを用いたダウンスケーリング海洋モデルを開発し，②外洋影響を考慮した瀬戸内海全体の流れと、複雑で詳細な地形情報を取り入れた高解像度モデリングを実施して流れ場を再現すること，さらに，③密度プリュームとして移流分散される垂水処理場からの処理水の広域分散過程をシミュレートし，大阪湾・播磨灘における処理水の挙動や湾内への滞留プロセスの解析，を行った．

 

tpas at surfacetpas at bottom

 
Figure: （左）計算開始16日後の表層での処理水分布（右）計算開始16日後の底層での処理水分布．大阪湾での処理水の分散過程において，表層と底層で違いが見られる．
 
 
↑

海浜地形変化に及ぼす波-流れ相互作用の影響に関する研究 †

 
  • 担当：宮崎 大
  • 共同研究相手：

　海岸侵食は，国土保全，環境保全，防災的な観点から土木工学的に非常に重要な問題である．海岸環境の精緻なアセスメントに向けて，海浜地形変化予測技術の高精度化が不可欠である．海浜地形は，波，流れ，地盤系に生じる不安定現象に支配され，離岸流に伴うrip channel, cusp地形など周期的に変化する場合が多く知られている．既往の研究として，Kaida and Uchiyama (2012)は，Yu and Slinn(2003), Weirら(2011)の研究を発展させ海岸侵食に密接に関わる海浜流のうち，離岸流の発達・抑制プロセスを詳細に解析し，波―流れ相互作用のうち，流れから波へのフィードバック機構(Current Effects on Wave),以後CEWと呼ぶが，CEWによって離岸流の発達が抑制され，岸近くにとどめられるメカニズムを力学的に明らかにした．また，Garnierら(2008)はrip channelの自励的な発達機構について検討しており，沿岸方向に一様なバー型海浜地形に対し，波―流れ相互作用を簡易的に取り入れた水理モデルに漂砂モデルを組み合わせることにより，周期的なrip channelの形成を再現したが，CEWのうち波の屈折効果の考慮が不十分であった．したがって，CEWを正確に考慮した場合に生じる離岸流の抑制効果がrip channelの形成過程に及ぼす影響については未解明なままである．
　そこで本研究での目的は，包括的な数値実験によりCEWの有無による海岸地形変化の差異を定量化し，そのメカニズムについて考察することである．そのためにまず，Uchiyama et al.(2009)の海浜流＋波浪変形モデルに漂砂モデルを組み込み，流体―地盤系の不安定によって生じるrip channelの形成・発達過程を再現する．次に，移動床数値実験によりCEWの有無によるrip channelの発達過程の差を定量的に評価する．さらに，十分に発達したrip channel地形を用いた固定床数値計算を行い，CEWの有無による流れ，波，漂砂量の変化を解析しCEWがrip channel形状に及ぼすメカニズムを明らかにした．

 
 
↑

波ー流れ相互作用による海浜流系へ及ぼす影響に関する研究 †

 
  • 担当：甲斐田秀樹
  • 共同研究相手：
  • 参考URL 離岸流の発達に及ぼす流れから波へのフィードバック機構の影響，流れから波へのフィードバック機構による離岸流の発達抑制メカニズム．
    離岸流は，砕波帯や沿岸域における物質輸送や海底地形変化等に重大な影響を及ぼす因子であると同時に，しばしば沖向きの速い流れを伴うことから，海浜利用者の安全性にも強く関与している．そのため，離岸流の定量的な発生・発達過程の予測は工学的に極めて重要な課題であり，詳細な実験，観測等に加えて，モデルによる精緻な解析を行うことが不可欠である．
    下図に示すように，流れから波へのフィードバック機構（Current Effects on Waves，以下CEWと略称する）を考慮しない場合は，離岸流が沖へと過剰に発達してしまい，十分な精度で離岸流を再現できないが，CEWを考慮するとその発達が適切に抑制されることが知られている（Haasら，1998；Yu･Slinn, 2003．）．CEWによる離岸流の抑制効果について，YS03は波のエネルギー平衡方程式中のradiation stressによる仕事の変化が，Weirら（2011）は様々なCEWの効果のうち，流れによる波の屈折作用（ドップラーシフトを通じた波数変化）が効果的であることをそれぞれ示した．しかしながら，CEWによる離岸流場の変調に対する力学構造の変化についての解析は手つかずのままであり，未だ詳細なメカニズムの解明には至っていない．そこで本研究では，離岸流によるCEWを介した波浪変形機構と，それに伴う海浜流場の運動量収支構造の変化に関する詳細な解析を行い，CEWによる離岸流の発達抑制メカニズムについて運動学的・力学的に明解に示した．
    また，波の入射角の違いに基づく，離岸流場・離岸流場から沿岸流場への遷移過程・沿岸流場についての力学的解析も行っており，その中でCEWが流れ場へ及ぼす影響の大きさが明らかになりつつある．
     
    rip current
    Figure: 本研究で得られた離岸流の発達パターンの一例．上図：CEWなしのケース，下図：CEWありのケース．カラーは渦度（1/s），コンターはLagrange流速による流線をそれぞれ表している．離岸流の流路であるrip channelはy=128mに，沿岸砂州の頂部はx=80mに位置している．
     
     
↑

超高解像度三次元海洋モデルを用いた砕波帯ー陸棚相互作用の解析 †

 
  • 担当：甲斐田秀樹
  • 共同研究相手：カリフォルニア大学ロサンゼルス校（UCLA）
    Coming soon.
Updated on 2021-07-16 (Fri) 14:07:08