[MIS13-03] Isolating oceanic contributions to ocean bottom pressure observations
by using numerical ocean models to detect transient seafloor displacements
Keywords:Crustal deformation, Ocean Bottom Pressure , Numerical Ocean Modeling
In this research, We focused on isolating non-tidal oceanic contributions by using several numerical ocean models. Ocean modeling has hardly ever been used to calculate OBP, however comparison has begun with expanding of OBP observation (e.g. Muramoto et al., 2018; Fredrickson et al., 2019). We show the applied models in this research as follows; #1 HYCOM (Cummings et al., 2013), #2 GLORYS (Ferry et al., 2010), #3 JCOPE2M (Miyazawa et al.,2017), #4 ECCO2 (Menemenlis et al., 2011), #5 Single-layer Ocean Model (SOM) (Inazu et al., 2012). Each model except #5 has tens of vertical layers to explain pycnocline. Horizontal resolution of model #1, #2, #3, #4 is 1/12° (about 10 km), though 1/4° is employed for model #4, also output from #3 is available only in Northwest Pacific. Model skills of each models were compared by amplitude after subtracting modeled ocean variations from observed OBP time series. We employed RMS noise level (3σ) and power spectrum of the observed OBP for evaluating the model skills. We compiled OBP datasets (10~12months) from several subduction zones around Pacific rim (Japan Trench (JT) 19 sites, Nankai Trough (NT) 51sites, Hikurangi Trough (HT) 7sites, Cascadia subduction zone (CSZ) 27sites, Chile trench (CT) 5sites).
We report that there were different characteristics of noise correction ability within 5 models for each region. Especially in HT, we can see the high predictability of SOM driven by wind stress and sea level pressure in short periodic band (shorter than 20 days) than other models driven by wind stress and heat/freshwater flux (e.g. 3σ comparison; before model correction 6.24cm; after SOM correction 3.57cm; after HYCOM correction 4.65cm ). This result implied that OBP variations of HT is significantly affected by sea level pressure, actually HT is located under the effect of low-pressure zone around Antarctica which causes strong synoptic scale atmospheric disturbance. Based on this result, We calculated wind stress and sea level pressure driven OBP variations separately using the SOM framework. After confirming the linearity of wind stress and sea level pressure to OBP variation, we combined multi-layer model output and sea level pressure driven OBP from SOM, then compared with observations. As a result, We got better noise correction ability from all multi-layer models combined with sea level pressure component in HT (e.g. GLORYS wind 4.88cm; wind + prs 4.02 cm ). In addition to the HT case, comparison in NT also showed large contributions of sea level pressure to OBP variations because all models became to be improved by combining sea level pressure component. From this result, We concluded sea level pressure contribution as well as wind stress should be considered in order to predict consistent OBP variations and contribute to detect centimeter-scale vertical seafloor displacements.