Pub Date : 2024-11-15DOI: 10.1016/j.ocemod.2024.102462
Yongzhi Liu , Minjie Xu , Xianqing Lv
Ecological modeling is an important methodology for studying the spatio-temporal evolution of marine ecosystem. Given the significant role of model parameters as a major source of uncertainty in ecological models, we propose a novel approach by combining the Conditional Nonlinear Optimal Perturbation related to Parameters (CNOP-P) method with the adjoint assimilation method to enhance predictive accuracy. CNOP-P denotes the parameter perturbation that leads to the greatest deviation of the model's development from the reference state. In comparison to other sensitivity analysis methods, this combined approach proves to be more efficient. Considering the nonlinearity of the model structure, the maximum development of the model does not consistently align with the extreme parameter values within the confidence interval. Minor parameter errors can lead to substantial model development, significantly impacting the precision of ecological models. Notably, traditional sensitivity analysis methods such as one-at-a-time (OAT) sensitivity analysis and global sensitivity analysis (GSA) methods fail to capture this characteristic. On the other hand, the GSA methods incur substantial computational costs and tends to overestimate the sensitivity of the most sensitive parameters while underestimating the sensitivity of less sensitive parameters. The combined approach of CNOP-P and adjoint assimilation enables the assimilation of satellite data and the simultaneous optimization of model parameters alongside the CNOP-P calculations. This integration substantially improves both efficiency and precision of the ecological model, thereby improving predictive skill.
{"title":"Improving ecological modeling: Integrating CNOP-P and adjoint assimilation in a coupled ecological model","authors":"Yongzhi Liu , Minjie Xu , Xianqing Lv","doi":"10.1016/j.ocemod.2024.102462","DOIUrl":"10.1016/j.ocemod.2024.102462","url":null,"abstract":"<div><div>Ecological modeling is an important methodology for studying the spatio-temporal evolution of marine ecosystem. Given the significant role of model parameters as a major source of uncertainty in ecological models, we propose a novel approach by combining the Conditional Nonlinear Optimal Perturbation related to Parameters (CNOP-P) method with the adjoint assimilation method to enhance predictive accuracy. CNOP-P denotes the parameter perturbation that leads to the greatest deviation of the model's development from the reference state. In comparison to other sensitivity analysis methods, this combined approach proves to be more efficient. Considering the nonlinearity of the model structure, the maximum development of the model does not consistently align with the extreme parameter values within the confidence interval. Minor parameter errors can lead to substantial model development, significantly impacting the precision of ecological models. Notably, traditional sensitivity analysis methods such as one-at-a-time (OAT) sensitivity analysis and global sensitivity analysis (GSA) methods fail to capture this characteristic. On the other hand, the GSA methods incur substantial computational costs and tends to overestimate the sensitivity of the most sensitive parameters while underestimating the sensitivity of less sensitive parameters. The combined approach of CNOP-P and adjoint assimilation enables the assimilation of satellite data and the simultaneous optimization of model parameters alongside the CNOP-P calculations. This integration substantially improves both efficiency and precision of the ecological model, thereby improving predictive skill.</div></div>","PeriodicalId":19457,"journal":{"name":"Ocean Modelling","volume":"193 ","pages":"Article 102462"},"PeriodicalIF":3.1,"publicationDate":"2024-11-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142702170","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2024-11-14DOI: 10.1016/j.ocemod.2024.102465
Anna Maggiorano, Mark Baird, Clothilde Langlais, Mathieu Mongin, Jennifer Skerratt
Solar radiation propagating through the water column is scattered and absorbed by optically active constituents in the ocean, in particular phytoplankton, coloured-dissolved organic matter (CDOM), suspended inorganic particulate matter (SPIM) and detritus. These wavelength-dependent processes affect the vertical distribution of heating in the water column and its stratification. The continental shelf north-east of Australia, containing the Great Barrier Reef (GBR), is characterised by highly seasonal and intermittent freshwater inputs leading to large sediment and nutrient discharges that strongly impact the water optical properties. While this complex mixture of optically active constituents is known to affect water clarity and the euphotic zone depth in the river plumes, its impact on the ocean circulation and thermal balance is still unclear at the scale of the GBR. In this study, we use a hydrodynamic-optical-biogeochemical ocean model to investigate the feedback between heat absorption by phytoplankton, CDOM and suspended sediments and ocean dynamics in the GBR region. The results show that the attenuation of the vertical heat flux due to phytoplankton, CDOM and SPIM concentrations is stronger on the continental shelf and dominated by the absorption and scattering from suspended sediments. The presence of absorbing constituents in the water column drives a temperature increase at the surface and a decrease below the mixed layer with stronger stratification and greater heat losses to the atmosphere. Inshore, the ocean heat content increases by up to 1% due to optically active constituents. Offshore, absorption by optically active constituents near the surface is compensated by less absorption underneath the mixed layer resulting in a decrease in the ocean heat content of the top 500 m. We find that considering a spatially- and temporally-variable vertical attenuation of heat due to multiple optically-active components improves hydrodynamic model skill. This study highlights the importance of the impact of water clarity and its spatial variability on hydrodynamic processes.
{"title":"Impact of phytoplankton, CDOM, and suspended sediments on the vertical attenuation of light, changing heat content and circulation on a continental shelf: A modelling study of the Great Barrier Reef","authors":"Anna Maggiorano, Mark Baird, Clothilde Langlais, Mathieu Mongin, Jennifer Skerratt","doi":"10.1016/j.ocemod.2024.102465","DOIUrl":"10.1016/j.ocemod.2024.102465","url":null,"abstract":"<div><div>Solar radiation propagating through the water column is scattered and absorbed by optically active constituents in the ocean, in particular phytoplankton, coloured-dissolved organic matter (CDOM), suspended inorganic particulate matter (SPIM) and detritus. These wavelength-dependent processes affect the vertical distribution of heating in the water column and its stratification. The continental shelf north-east of Australia, containing the Great Barrier Reef (GBR), is characterised by highly seasonal and intermittent freshwater inputs leading to large sediment and nutrient discharges that strongly impact the water optical properties. While this complex mixture of optically active constituents is known to affect water clarity and the euphotic zone depth in the river plumes, its impact on the ocean circulation and thermal balance is still unclear at the scale of the GBR. In this study, we use a hydrodynamic-optical-biogeochemical ocean model to investigate the feedback between heat absorption by phytoplankton, CDOM and suspended sediments and ocean dynamics in the GBR region. The results show that the attenuation of the vertical heat flux due to phytoplankton, CDOM and SPIM concentrations is stronger on the continental shelf and dominated by the absorption and scattering from suspended sediments. The presence of absorbing constituents in the water column drives a temperature increase at the surface and a decrease below the mixed layer with stronger stratification and greater heat losses to the atmosphere. Inshore, the ocean heat content increases by up to 1% due to optically active constituents. Offshore, absorption by optically active constituents near the surface is compensated by less absorption underneath the mixed layer resulting in a decrease in the ocean heat content of the top 500 m. We find that considering a spatially- and temporally-variable vertical attenuation of heat due to multiple optically-active components improves hydrodynamic model skill. This study highlights the importance of the impact of water clarity and its spatial variability on hydrodynamic processes.</div></div>","PeriodicalId":19457,"journal":{"name":"Ocean Modelling","volume":"193 ","pages":"Article 102465"},"PeriodicalIF":3.1,"publicationDate":"2024-11-14","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142702168","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2024-11-13DOI: 10.1016/j.ocemod.2024.102458
H.M. Aravind , Tamay M. Özgökmen , Michael R. Allshouse
Lagrangian analyses of oceanic flows provide insight into the various transport pathways in the ocean. This analysis typically relies on a dense set of trajectories that can be computed using high-resolution velocity fields, which are often not available during field experiments. Instruments like drifters and floats are often employed to overcome the limitations imposed by satellite- and radar-based velocity fields, to understand the transport pathways in the ocean. However, the sparsity in available drifter-trajectory data proves prohibitive to obtaining a comprehensive map of the Lagrangian characteristics of the underlying flow. To circumvent these issues, we use Gaussian Process Regression (GPR) to obtain velocity fields from sparse drifter data to generate synthetic trajectories and subsequently estimate two Lagrangian metrics, FTLE and dilation rate. A detailed error analysis is performed for drifter clusters deployed within various dynamical regions in the analytic Bickley jet system. The uncertainties in velocity reconstruction obtained from the GPR method, averaged along particle trajectories, locate Lagrangian confidence regions that are applicable both to synthetic trajectories and the dilation rate field. A sensitivity analysis reveals the role played by factors such as the spatial sampling density and temporal resolution of the drifter data, as well as the effect of position uncertainty as a result of GPS inaccuracy. The method is then applied to the drifter data from the Lagrangian Submesoscale Experiment in 2016 to locate convergent filaments. The results present a marked improvement over direct estimation of area-averaged dilation rates using drifter clusters.
{"title":"Lagrangian analysis of submesoscale flows from sparse data using Gaussian Process Regression for field reconstruction","authors":"H.M. Aravind , Tamay M. Özgökmen , Michael R. Allshouse","doi":"10.1016/j.ocemod.2024.102458","DOIUrl":"10.1016/j.ocemod.2024.102458","url":null,"abstract":"<div><div>Lagrangian analyses of oceanic flows provide insight into the various transport pathways in the ocean. This analysis typically relies on a dense set of trajectories that can be computed using high-resolution velocity fields, which are often not available during field experiments. Instruments like drifters and floats are often employed to overcome the limitations imposed by satellite- and radar-based velocity fields, to understand the transport pathways in the ocean. However, the sparsity in available drifter-trajectory data proves prohibitive to obtaining a comprehensive map of the Lagrangian characteristics of the underlying flow. To circumvent these issues, we use Gaussian Process Regression (GPR) to obtain velocity fields from sparse drifter data to generate synthetic trajectories and subsequently estimate two Lagrangian metrics, FTLE and dilation rate. A detailed error analysis is performed for drifter clusters deployed within various dynamical regions in the analytic Bickley jet system. The uncertainties in velocity reconstruction obtained from the GPR method, averaged along particle trajectories, locate Lagrangian confidence regions that are applicable both to synthetic trajectories and the dilation rate field. A sensitivity analysis reveals the role played by factors such as the spatial sampling density and temporal resolution of the drifter data, as well as the effect of position uncertainty as a result of GPS inaccuracy. The method is then applied to the drifter data from the Lagrangian Submesoscale Experiment in 2016 to locate convergent filaments. The results present a marked improvement over direct estimation of area-averaged dilation rates using drifter clusters.</div></div>","PeriodicalId":19457,"journal":{"name":"Ocean Modelling","volume":"193 ","pages":"Article 102458"},"PeriodicalIF":3.1,"publicationDate":"2024-11-13","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142748208","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2024-11-12DOI: 10.1016/j.ocemod.2024.102459
Mikhail Lavrentiev , Andrey Marchuk , Konstantin Oblaukhov
This paper proposes the use of specialized hardware accelerator based on the Field Programmable Gates Array (FPGA) microchip to compute tsunami wave propagation to assess and manage risks of marine natural disasters, namely, tsunami waves caused by underwater earthquakes. After a sufficiently strong seismic event, many countries and research centres launch extensive computations to estimate the tsunami wave parameters in certain parts of the coast to determine if a declaration of a tsunami alarm is warranted. This requires high computating powers which leads to higher energy costs. The paper demonstrates how an FPGA-based special Calculator (architecture of which has been earlier proposed by the authors), installed on a Personal Computer (PC) could be used to calculate the propagation of a tsunami wave over the entire Pacific Ocean, from the subduction zone offshore Kamchatka Peninsula and Kuril Islands to the coast of Chile. Such calculations offer reliable results within a few minutes and make it possible to obtain the distribution of expected tsunami wave heights along the coast. If the obtained results indicate a danger to the population or possible destruction of infrastructure, it becomes paramount to carry out more detailed calculations to accurately estimate the wave parameters at specific locations along the coast where negative consequences are expected. This requires cluster and/or supercomputer systems, which consume significant energy and hence are expensive. In case the modelling results indicate small values of maximum wave heights at populated coastal areas, population of the near-shore regions can be immediately informed about low amplitude tsunami wave; more detailed studies are not needed. This hence leads to noticeable savings in energy consumption. The paper presents a calculation of the propagation of a tsunami wave across the Pacific Ocean on a personal computer using a FPGA-based hardware acceleration of a computer code execution.
{"title":"Low power computation of transoceanic wave propagation for tsunami hazard mitigation","authors":"Mikhail Lavrentiev , Andrey Marchuk , Konstantin Oblaukhov","doi":"10.1016/j.ocemod.2024.102459","DOIUrl":"10.1016/j.ocemod.2024.102459","url":null,"abstract":"<div><div>This paper proposes the use of specialized hardware accelerator based on the Field Programmable Gates Array (FPGA) microchip to compute tsunami wave propagation to assess and manage risks of marine natural disasters, namely, tsunami waves caused by underwater earthquakes. After a sufficiently strong seismic event, many countries and research centres launch extensive computations to estimate the tsunami wave parameters in certain parts of the coast to determine if a declaration of a tsunami alarm is warranted. This requires high computating powers which leads to higher energy costs. The paper demonstrates how an FPGA-based special Calculator (architecture of which has been earlier proposed by the authors), installed on a Personal Computer (PC) could be used to calculate the propagation of a tsunami wave over the entire Pacific Ocean, from the subduction zone offshore Kamchatka Peninsula and Kuril Islands to the coast of Chile. Such calculations offer reliable results within a few minutes and make it possible to obtain the distribution of expected tsunami wave heights along the coast. If the obtained results indicate a danger to the population or possible destruction of infrastructure, it becomes paramount to carry out more detailed calculations to accurately estimate the wave parameters at specific locations along the coast where negative consequences are expected. This requires cluster and/or supercomputer systems, which consume significant energy and hence are expensive. In case the modelling results indicate small values of maximum wave heights at populated coastal areas, population of the near-shore regions can be immediately informed about low amplitude tsunami wave; more detailed studies are not needed. This hence leads to noticeable savings in energy consumption. The paper presents a calculation of the propagation of a tsunami wave across the Pacific Ocean on a personal computer using a FPGA-based hardware acceleration of a computer code execution.</div></div>","PeriodicalId":19457,"journal":{"name":"Ocean Modelling","volume":"192 ","pages":"Article 102459"},"PeriodicalIF":3.1,"publicationDate":"2024-11-12","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142652228","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
This work examines the use of discrete variance decay of tracers to estimate locally in space and time the numerical mixing caused by different processes during a tracer transport step. Expressions for local discrete variance decay (DVD) rates are directly derived from discrete tracer equations without any assumptions on discrete fluxes of the second moment. They relate the DVD rates to the fluxes of the first moment through the faces of scalar control cell. Mixing associated with advective and diffusive fluxes is thus estimated. The new framework avoids the need for second-moment flux definition when solved directly on finite-volume cell faces but still invokes certain second-moment fluxes when the face DVD rates are partitioned to cells sharing the face. These implied discrete fluxes depend on the partitioning and are non-unique. For third- or higher-order advection schemes, the DVD rates are contaminated by dispersive errors intrinsic to the approach, introducing uncertainty to the locality of any estimates produced by it. Additional temporal averaging or coarse-graining is thus necessary. Through the application of this technique, numerical mixing is found to be correlated with the distribution of eddy kinetic energy. Numerical mixing induced by vertical advection is found to be relatively small and correlated with the distribution of buoyancy fluxes. The explored high-order schemes are found to demonstrate levels of spurious mixing which may locally exceed physical mixing.
这项研究利用示踪剂的离散方差衰减来估算示踪剂传输步骤中不同过程在空间和时间上造成的局部数值混合。局部离散方差衰减(DVD)率的表达式直接来自离散示踪剂方程,无需对第二时刻的离散通量做任何假设。它们将 DVD 率与通过标量控制单元面的第一时刻通量联系起来。这样就可以估算出与平流和扩散通量相关的混合情况。当直接求解有限体积单元面时,新框架避免了第二时刻通量定义的需要,但当面的 DVD 率被分割到共享面的单元时,仍会调用某些第二时刻通量。这些隐含的离散通量取决于分区,并且是非唯一的。对于三阶或更高阶的平流方案,DVD 率会受到该方法固有的分散误差的污染,从而为其产生的任何估计值的定位带来不确定性。因此,有必要进行额外的时间平均或粗粒化。通过应用这种技术,发现数值混合与涡旋动能的分布相关。垂直平流引起的数值混合相对较小,并且与浮力通量的分布相关。研究发现,所探讨的高阶方案显示的虚假混合程度可能会局部超过物理混合程度。
{"title":"Discrete variance decay analysis of spurious mixing","authors":"Tridib Banerjee , Sergey Danilov , Knut Klingbeil , Jean-Michel Campin","doi":"10.1016/j.ocemod.2024.102460","DOIUrl":"10.1016/j.ocemod.2024.102460","url":null,"abstract":"<div><div>This work examines the use of discrete variance decay of tracers to estimate locally in space and time the numerical mixing caused by different processes during a tracer transport step. Expressions for local discrete variance decay (DVD) rates are directly derived from discrete tracer equations without any assumptions on discrete fluxes of the second moment. They relate the DVD rates to the fluxes of the first moment through the faces of scalar control cell. Mixing associated with advective and diffusive fluxes is thus estimated. The new framework avoids the need for second-moment flux definition when solved directly on finite-volume cell faces but still invokes certain second-moment fluxes when the face DVD rates are partitioned to cells sharing the face. These implied discrete fluxes depend on the partitioning and are non-unique. For third- or higher-order advection schemes, the DVD rates are contaminated by dispersive errors intrinsic to the approach, introducing uncertainty to the locality of any estimates produced by it. Additional temporal averaging or coarse-graining is thus necessary. Through the application of this technique, numerical mixing is found to be correlated with the distribution of eddy kinetic energy. Numerical mixing induced by vertical advection is found to be relatively small and correlated with the distribution of buoyancy fluxes. The explored high-order schemes are found to demonstrate levels of spurious mixing which may locally exceed physical mixing.</div></div>","PeriodicalId":19457,"journal":{"name":"Ocean Modelling","volume":"192 ","pages":"Article 102460"},"PeriodicalIF":3.1,"publicationDate":"2024-11-09","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142652230","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2024-11-02DOI: 10.1016/j.ocemod.2024.102461
Jian-Guo Li , Peitao Wang
A model of shallow water equations (SWEs) on a spherical multiple-cell (SMC) grid of 4-level (2.5–5–10–20 km) spatial resolutions is used to simulate tsunami propagation on global ocean surfaces. The unstructured SMC grid retains rectangular cells of the longitude-latitude grid so that efficient finite-difference schemes could be used. It also supports multi-resolutions like mesh refinement to resolve small islands and coastline details while keeping models compact enough to fit into available computers. Two earthquake-induced tsunami cases are simulated and compared with available observations. Results indicate that the modelled tsunami arrival times agree well with observations while tsunami wave heights are underestimated, particularly the observed runups on remote coastal lands. Possible reasons for this underestimation include the smoothing scheme used to suppress numerical oscillations and the missing of initial kinetic energy input from the earthquakes. Another reason is the limitation of the SWEs to describe coastal bores and breaking waves in coastal waters. A possible tsunami scenario induced by a landslide in the Canary Islands is also simulated to assess its potential impact on Atlantic coastal regions. Model results indicate that this kind of tsunami may cause severe damage to local areas but its effects on far fields, like the UK and American coastal regions are small. As the initial landslide disturbance is overly simplified, this study does not give a true representation of a real landslide tsunami but rather a qualitative assessment of its impact on the Atlantic Ocean. More realistic initial condition and improved model representation of coastal processes are needed for further studies of this possible landslide hazard.
{"title":"Global tsunami modelling on a spherical multiple-cell grid","authors":"Jian-Guo Li , Peitao Wang","doi":"10.1016/j.ocemod.2024.102461","DOIUrl":"10.1016/j.ocemod.2024.102461","url":null,"abstract":"<div><div>A model of shallow water equations (SWEs) on a spherical multiple-cell (SMC) grid of 4-level (2.5–5–10–20 km) spatial resolutions is used to simulate tsunami propagation on global ocean surfaces. The unstructured SMC grid retains rectangular cells of the longitude-latitude grid so that efficient finite-difference schemes could be used. It also supports multi-resolutions like mesh refinement to resolve small islands and coastline details while keeping models compact enough to fit into available computers. Two earthquake-induced tsunami cases are simulated and compared with available observations. Results indicate that the modelled tsunami arrival times agree well with observations while tsunami wave heights are underestimated, particularly the observed runups on remote coastal lands. Possible reasons for this underestimation include the smoothing scheme used to suppress numerical oscillations and the missing of initial kinetic energy input from the earthquakes. Another reason is the limitation of the SWEs to describe coastal bores and breaking waves in coastal waters. A possible tsunami scenario induced by a landslide in the Canary Islands is also simulated to assess its potential impact on Atlantic coastal regions. Model results indicate that this kind of tsunami may cause severe damage to local areas but its effects on far fields, like the UK and American coastal regions are small. As the initial landslide disturbance is overly simplified, this study does not give a true representation of a real landslide tsunami but rather a qualitative assessment of its impact on the Atlantic Ocean. More realistic initial condition and improved model representation of coastal processes are needed for further studies of this possible landslide hazard.</div></div>","PeriodicalId":19457,"journal":{"name":"Ocean Modelling","volume":"192 ","pages":"Article 102461"},"PeriodicalIF":3.1,"publicationDate":"2024-11-02","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142652229","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2024-10-31DOI: 10.1016/j.ocemod.2024.102457
Polina Verezemskaya , Sergey Gulev , Bernard Barnier , Jean-Marc Molines , Pedro Colombo , Alexander Gavrikov , Mikhail Kalinin
The dynamics of the North Atlantic subpolar gyre (SPG) largely control the Atlantic Meridional Overturning Circulation (AMOC), thus affecting the long-term climate variability of the global ocean. Many key processes in the North Atlantic SPG, such as Gulf Stream detachment, eddy activity, deep convection and overflows, are frequently incorrectly presented or misrepresented in modern ocean general circulation models. Here, we demonstrate a new regional eddy-resolving model of the subpolar North Atlantic based on the NEMO4 configuration and analyze its ability to adequately represent the dynamics of the subpolar North Atlantic. The model performance is assessed using in situ data (cross-section hydrography, Argo and sea gliders), remote sensing data (satellite sea surface temperature and salinity) and the ANDRO and ARMOR3D products based on optimal interpolation products. In the next step, we analyze the effects of several recent developments in ocean modeling to improve the model solution. The implementation of the self-diffusive momentum advection scheme with additional viscosity allows for the simulation of more realistic Deep Western Boundary Current and Irminger Rings lifecycle characteristics. The surface current feedback parameterization implemented in the forcing function and the modification of the kinetic energy budget in the upper ocean, especially in intensively eddying areas, simulate the behavior predicted by coupled models. While the model is sensitive to the implementation of cool skin–warm layer parameterization, this parameterization does not significantly improve the solution compared to the observational data. Regional implementation of the local-sigma vertical coordinate system in cascading areas allows overflow waters to spread over the ocean bottom more realistically. The remote effect of the implementation of the local-sigma coordinate is traced in the vertical thermohaline structure at the western slopes of the Irminger Sea and Iceland Basin as well as in the Deep Western Boundary Current in the Labrador Sea. Compared to z-partial step setting, local-sigma allows for correct representation of continuous densification of the water in the overturning part of the AMOC and thus North Atlantic Deep Water formation. The potential of the new model configuration for analyzing the subpolar North Atlantic, including convective activity, is discussed.
{"title":"A regional NEMO 4.0 configuration of the subpolar North Atlantic","authors":"Polina Verezemskaya , Sergey Gulev , Bernard Barnier , Jean-Marc Molines , Pedro Colombo , Alexander Gavrikov , Mikhail Kalinin","doi":"10.1016/j.ocemod.2024.102457","DOIUrl":"10.1016/j.ocemod.2024.102457","url":null,"abstract":"<div><div>The dynamics of the North Atlantic subpolar gyre (SPG) largely control the Atlantic Meridional Overturning Circulation (AMOC), thus affecting the long-term climate variability of the global ocean. Many key processes in the North Atlantic SPG, such as Gulf Stream detachment, eddy activity, deep convection and overflows, are frequently incorrectly presented or misrepresented in modern ocean general circulation models. Here, we demonstrate a new regional eddy-resolving model of the subpolar North Atlantic based on the NEMO4 configuration and analyze its ability to adequately represent the dynamics of the subpolar North Atlantic. The model performance is assessed using in situ data (cross-section hydrography, Argo and sea gliders), remote sensing data (satellite sea surface temperature and salinity) and the ANDRO and ARMOR3D products based on optimal interpolation products. In the next step, we analyze the effects of several recent developments in ocean modeling to improve the model solution. The implementation of the self-diffusive momentum advection scheme with additional viscosity allows for the simulation of more realistic Deep Western Boundary Current and Irminger Rings lifecycle characteristics. The surface current feedback parameterization implemented in the forcing function and the modification of the kinetic energy budget in the upper ocean, especially in intensively eddying areas, simulate the behavior predicted by coupled models. While the model is sensitive to the implementation of cool skin–warm layer parameterization, this parameterization does not significantly improve the solution compared to the observational data. Regional implementation of the local-sigma vertical coordinate system in cascading areas allows overflow waters to spread over the ocean bottom more realistically. The remote effect of the implementation of the local-sigma coordinate is traced in the vertical thermohaline structure at the western slopes of the Irminger Sea and Iceland Basin as well as in the Deep Western Boundary Current in the Labrador Sea. Compared to z-partial step setting, local-sigma allows for correct representation of continuous densification of the water in the overturning part of the AMOC and thus North Atlantic Deep Water formation. The potential of the new model configuration for analyzing the subpolar North Atlantic, including convective activity, is discussed.</div></div>","PeriodicalId":19457,"journal":{"name":"Ocean Modelling","volume":"194 ","pages":"Article 102457"},"PeriodicalIF":3.1,"publicationDate":"2024-10-31","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143148258","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2024-10-17DOI: 10.1016/j.ocemod.2024.102448
David E. Gregg , Nigel T. Penna , Christopher Jones , Miguel A. Morales Maqueda
The accuracy of global ocean tide models is assessed in coastal waters of the European North West Shelf to ascertain where higher resolution local (forecast) models are most needed for geophysical and navigational applications, and which global models are most suitable for providing boundary conditions for regional and local tide models. Five recent global ocean tide models (FES2014b, EOT20, TPXO9-atlas-v5, GOT4.10c, and DTU16) are considered, with the models first compared by interpolating them onto common grids and computing the mean absolute deviation at each grid point. Coastline tide gauge and offshore bottom pressure sensor data were collated from several sources to give a total of 279 observation sites for evaluating model accuracy, including observational values from 137 locations that have not previously been released and have therefore not been assimilated into any of the global models tested. The residual errors between each model’s predicted phasor and the corresponding observed phasor were calculated at each observation location, and quantified using the root mean square (RMS) and median absolute residual (MAR) for the eight tidal constituents M2, S2, N2, O1, K1, K2, P1, and Q1. To avoid RMS values being biased by observation point density, a Voronoi-weighted RMS based on the water area of the Voronoi polygon about each observation location was also developed and used. Four zones were defined based on ocean depth to gauge model performance, and model inaccuracy is again demonstrated in near-shore regions. Seven further zones were defined based on geographical areas, which reveals inhomogeneity among the global models. The smallest overall root sum square (RSS) RMS error across all eight constituents arises with FES2014b, although TPXO9-atlas-v5 has the best performance when using the MAR and Voronoi-weighted RMS metrics. Using only the 137 observation sites that have not been assimilated by any model and therefore provide an independent accuracy assessment, FES2014b exhibits the smallest errors at the coastline, with an RSS RMS of 24.46 cm. All models exhibit larger errors with the 137 independent observation sites than with all 279 observation sites, with an average overall increase in RSS RMS error of 12%, and an increase of 30% for coastline tide gauges, highlighting the need for local model development in these areas.
{"title":"Accuracy assessment of recent global ocean tide models in coastal waters of the European North West Shelf","authors":"David E. Gregg , Nigel T. Penna , Christopher Jones , Miguel A. Morales Maqueda","doi":"10.1016/j.ocemod.2024.102448","DOIUrl":"10.1016/j.ocemod.2024.102448","url":null,"abstract":"<div><div>The accuracy of global ocean tide models is assessed in coastal waters of the European North West Shelf to ascertain where higher resolution local (forecast) models are most needed for geophysical and navigational applications, and which global models are most suitable for providing boundary conditions for regional and local tide models. Five recent global ocean tide models (FES2014b, EOT20, TPXO9-atlas-v5, GOT4.10c, and DTU16) are considered, with the models first compared by interpolating them onto common grids and computing the mean absolute deviation at each grid point. Coastline tide gauge and offshore bottom pressure sensor data were collated from several sources to give a total of 279 observation sites for evaluating model accuracy, including observational values from 137 locations that have not previously been released and have therefore not been assimilated into any of the global models tested. The residual errors between each model’s predicted phasor and the corresponding observed phasor were calculated at each observation location, and quantified using the root mean square (RMS) and median absolute residual (MAR) for the eight tidal constituents M2, S2, N2, O1, K1, K2, P1, and Q1. To avoid RMS values being biased by observation point density, a Voronoi-weighted RMS based on the water area of the Voronoi polygon about each observation location was also developed and used. Four zones were defined based on ocean depth to gauge model performance, and model inaccuracy is again demonstrated in near-shore regions. Seven further zones were defined based on geographical areas, which reveals inhomogeneity among the global models. The smallest overall root sum square (RSS) RMS error across all eight constituents arises with FES2014b, although TPXO9-atlas-v5 has the best performance when using the MAR and Voronoi-weighted RMS metrics. Using only the 137 observation sites that have not been assimilated by any model and therefore provide an independent accuracy assessment, FES2014b exhibits the smallest errors at the coastline, with an RSS RMS of 24.46 cm. All models exhibit larger errors with the 137 independent observation sites than with all 279 observation sites, with an average overall increase in RSS RMS error of 12%, and an increase of 30% for coastline tide gauges, highlighting the need for local model development in these areas.</div></div>","PeriodicalId":19457,"journal":{"name":"Ocean Modelling","volume":"192 ","pages":"Article 102448"},"PeriodicalIF":3.1,"publicationDate":"2024-10-17","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142527781","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2024-10-16DOI: 10.1016/j.ocemod.2024.102455
Harikrishnan Sreeshylam , Zhilong Liu , Brian Dzwonkowski , John Lehrter
Accurate estimation of water column temperature is vital for modeling physical and biogeochemical processes. A key process in the thermal dynamics of the upper ocean is the attenuation of solar radiation. In shallow-turbid coastal systems, spatially and temporally varying optical characteristics present challenges for commonly used attenuation parameterization schemes. This study investigates the dependency of temperature with a ROMS model of Mobile Bay, a shallow, turbid estuary, using six different attenuation approaches including three base cases: Conventional approach PS77 based on water type-9; Novel approach SAL relating in situ PAR attenuation to salinity; and Surface trapped irradiance method ST. In addition, these base cases are also tested with surface atmospheric heat flux correction (QC). Simulations were validated against observations from various sources to identify the optimal approach at annual and synoptic scales. While all simulations showed effective temperature performance over an annual cycle, monthly analysis revealed some seasonality, with winter months typically performing better than summer months. The influence of QC notably enhanced temperature performance in both annual and synoptic scales, given that surface heat flux primarily drove temperature changes in this shallow system. The best overall performance was determined to be the ST approach incorporating QC. Conversely, PS77 without QC demonstrated the poorest performance. The SAL model with QC, notably improved performance over PS77 with QC, yet demonstrated comparable yet weaker performance compared to the ST model with QC. The study also implies that neglecting subseasonal validation in long-term regional climate modeling could introduce uncertainty into analyzing events tied to subseasonal temperatures.
准确估算水体温度对物理和生物地球化学过程建模至关重要。太阳辐射衰减是上层海洋热动力学的一个关键过程。在浅湍流沿岸系统中,时空变化的光学特性给常用的衰减参数化方案带来了挑战。本研究采用六种不同的衰减方法(包括三种基本方法),对莫比尔湾(一个浅浊流河口)的 ROMS 模式的温度依赖性进行了研究:基于水体类型-9 的传统方法 PS77;与原位 PAR 衰减和盐度有关的新方法 SAL;以及表面受困辐照度方法 ST。此外,还对这些基本方法进行了地表大气热通量校正(QC)测试。模拟结果与各种来源的观测数据进行了验证,以确定年度和天气尺度上的最佳方法。虽然所有模拟都显示出有效的年周期温度表现,但月度分析显示出一定的季节性,冬季月份的表现通常好于夏季月份。由于地表热通量主要驱动浅层系统的温度变化,因此 QC 的影响明显提高了年尺度和同步尺度的温度性能。总体性能最好的是包含 QC 的 ST 方法。相反,不含 QC 的 PS77 性能最差。采用 QC 的 SAL 模型比采用 QC 的 PS77 有明显改善,但与采用 QC 的 ST 模型相比,其性能相当但较弱。这项研究还表明,在长期区域气候模拟中忽略次季节验证可能会给分析与次季节温度相关的事件带来不确定性。
{"title":"Enhancing model temperature estimations in shallow, turbid, coastal regions: Mobile Bay, Alabama","authors":"Harikrishnan Sreeshylam , Zhilong Liu , Brian Dzwonkowski , John Lehrter","doi":"10.1016/j.ocemod.2024.102455","DOIUrl":"10.1016/j.ocemod.2024.102455","url":null,"abstract":"<div><div>Accurate estimation of water column temperature is vital for modeling physical and biogeochemical processes. A key process in the thermal dynamics of the upper ocean is the attenuation of solar radiation. In shallow-turbid coastal systems, spatially and temporally varying optical characteristics present challenges for commonly used attenuation parameterization schemes. This study investigates the dependency of temperature with a ROMS model of Mobile Bay, a shallow, turbid estuary, using six different attenuation approaches including three base cases: Conventional approach PS77 based on water type-9; Novel approach SAL relating in situ PAR attenuation to salinity; and Surface trapped irradiance method ST. In addition, these base cases are also tested with surface atmospheric heat flux correction (QC). Simulations were validated against observations from various sources to identify the optimal approach at annual and synoptic scales. While all simulations showed effective temperature performance over an annual cycle, monthly analysis revealed some seasonality, with winter months typically performing better than summer months. The influence of QC notably enhanced temperature performance in both annual and synoptic scales, given that surface heat flux primarily drove temperature changes in this shallow system. The best overall performance was determined to be the ST approach incorporating QC. Conversely, PS77 without QC demonstrated the poorest performance. The SAL model with QC, notably improved performance over PS77 with QC, yet demonstrated comparable yet weaker performance compared to the ST model with QC. The study also implies that neglecting subseasonal validation in long-term regional climate modeling could introduce uncertainty into analyzing events tied to subseasonal temperatures.</div></div>","PeriodicalId":19457,"journal":{"name":"Ocean Modelling","volume":"192 ","pages":"Article 102455"},"PeriodicalIF":3.1,"publicationDate":"2024-10-16","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142527780","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2024-10-11DOI: 10.1016/j.ocemod.2024.102451
Linlin Cui , Fei Ye , Y. Joseph Zhang , Haocheng Yu , Zhengui Wang , Saeed Moghimi , Gregory Seroka , Jack Riley , Shachak Pe'eri , Soroosh Mani , Edward Myers , Kyungmin Park , Liujuan Tang , Zizang Yang , Yan-Ming Wang
We demonstrate recent progress made in the simulation of total water level (TWL) at continental scale, using the coastal ocean of US East Coast/Gulf of Mexico coast as an example. A key difference between the continental-scale and small-scale modeling is that the former requires a more accurate vertical datum. Using a geoid-based datum (xGEOID20b), a satellite altimetry product, and a state-of-the-art 3D unstructured-grid model, we significantly improve the accuracy for TWL both near- and off-shore. The average root-mean-square error at all NOAA stations is 14 cm. The non-tidal signals are found to be sensitive to the representation of a large-scale current system near the boundary and extending the domain extent to accommodate this system improves these signals.
{"title":"Total water level prediction at continental scale: Coastal ocean","authors":"Linlin Cui , Fei Ye , Y. Joseph Zhang , Haocheng Yu , Zhengui Wang , Saeed Moghimi , Gregory Seroka , Jack Riley , Shachak Pe'eri , Soroosh Mani , Edward Myers , Kyungmin Park , Liujuan Tang , Zizang Yang , Yan-Ming Wang","doi":"10.1016/j.ocemod.2024.102451","DOIUrl":"10.1016/j.ocemod.2024.102451","url":null,"abstract":"<div><div>We demonstrate recent progress made in the simulation of total water level (TWL) at continental scale, using the coastal ocean of US East Coast/Gulf of Mexico coast as an example. A key difference between the continental-scale and small-scale modeling is that the former requires a more accurate vertical datum. Using a geoid-based datum (xGEOID20b), a satellite altimetry product, and a state-of-the-art 3D unstructured-grid model, we significantly improve the accuracy for TWL both near- and off-shore. The average root-mean-square error at all NOAA stations is 14 cm. The non-tidal signals are found to be sensitive to the representation of a large-scale current system near the boundary and extending the domain extent to accommodate this system improves these signals.</div></div>","PeriodicalId":19457,"journal":{"name":"Ocean Modelling","volume":"192 ","pages":"Article 102451"},"PeriodicalIF":3.1,"publicationDate":"2024-10-11","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142437825","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}