A method and system for decomposition of the dynamic process of sea surface temperature anomaly of oceanic eddy

Abstract

The application discloses a method and system for decomposing a dynamic process of a marine eddy sea surface temperature anomaly, and the method comprises the following steps: for an eddy region of the sea surface temperature anomaly, a process simulation of a meridional transport is carried out by considering diffusion and damping effects, and a process simulation of vertical motion is carried out by considering the upward and downward flow of the eddy center; a total simulation model of the eddy sea surface temperature anomaly is established according to the process simulation of the meridional transport and the process simulation of the vertical motion, and parameter inversion is carried out to obtain optimal inversion parameters; and the eddy sea surface temperature anomaly caused by the meridional transport and the vertical motion is determined respectively according to the optimal inversion parameters. The application can clearly reflect the physical source of the temperature anomaly by decomposing the sea surface temperature anomaly into the process of the meridional transport and the process of the vertical motion from the perspective of the dynamic process, and has physical interpretability.

Description

Technical Field

[0001] This invention relates to ocean dynamic process modeling and data inversion technology, and in particular to a method and system for decomposing the dynamic processes of ocean eddies and sea surface temperature anomalies. Background Technology

[0002] Ocean eddies are important mesoscale circulation structures widely distributed in the global ocean, typically ranging from tens to hundreds of kilometers in horizontal scale, with lifecycles lasting from weeks to months. They play a crucial role in heat transport, mass distribution, air-sea flux regulation, and global climate change, making them a core topic in ocean dynamics research. Ocean eddies cause significant sea surface temperature anomalies in their interior and periphery by perturbing the surrounding temperature and current fields. Sea surface temperature anomalies are important physical quantities for observing and studying ocean-atmosphere interactions, widely used to analyze the mechanisms by which eddies influence climate systems, cloud cover, precipitation, and atmospheric circulation. Therefore, a deep understanding of the causes and dynamic structures of eddy sea surface temperature anomalies is crucial for revealing the mechanisms of eddy interaction.

[0003] Existing research primarily utilizes satellite remote sensing data and Argo buoy data to observe and analyze the spatial structure of eddy sea surface temperature anomalies. A common approach is to radially average the anomalies and empirically decompose them into monopole and dipole structures. However, these decomposition methods rely on the assumption of eddy axisymmetry and artificially defined radial profiles, lacking rigorous physical and mathematical foundations and failing to fully reflect the coupling of actual physical processes. Especially when multiple physical mechanisms are at play or the background field is non-uniform, this method often exhibits significant bias, making it difficult to accurately identify and quantify the contribution of each process to the eddy sea surface temperature anomaly. Summary of the Invention

[0004] Purpose of the invention: The purpose of this invention is to provide a method and system for decomposing the dynamic process of ocean eddy sea surface temperature anomalies, which is used to analyze and invert the dynamic process of the formation mechanism of ocean eddy sea surface temperature anomalies, and to achieve quantitative characterization and simulation reconstruction of its dynamic mechanism.

[0005] Technical solution: The method for decomposing the dynamic process of sea surface temperature anomalies in ocean eddies according to the present invention includes the following steps:

[0006] For vortex regions with sea surface temperature anomalies, the advection process is simulated considering diffusion and damping effects, and the vertical motion process is simulated considering the rising and sinking flows at the center of the vortex.

[0007] Based on the advection process simulation and vertical motion process simulation, an overall simulation model of vortex sea surface temperature anomaly was established, and parameter inversion was performed to obtain the optimal inversion parameters;

[0008] The optimal inversion parameters are used to determine the vortex sea surface temperature anomalies caused by the advection process and the vertical motion process, respectively.

[0009] Furthermore, at ocean grid points At this location, the vortex of the advection process causes sea surface temperature anomalies. The simulation model is as follows:

[0010]

[0011] in, The diffusion coefficient is... The damping coefficient is... For ocean grid points The vortex background temperature field at that location, and These represent the flow velocities in the east-west and north-south directions, respectively. For the Nabla operator, For time;

[0012] At marine grid points At this location, the vertical motion process vortex causes sea surface temperature anomalies. The simulation model is as follows:

[0013]

[0014] in, Reaction vertical motion intensity The scope of the reaction's influence For ocean grid points The sea level is abnormal.

[0015] Furthermore, the overall simulation model of eddy sea surface temperature anomalies as follows:

[0016]

[0017] in, The integration time is the time for the advection process model. The method is found through parameter inversion. , , , , The optimal parameter values.

[0018] Furthermore, the cost function is solved using numerical optimization methods to obtain... , , , , The optimal parameter values ​​and the cost function are:

[0019]

[0020] in, For ocean grid points Sea surface temperature anomaly observed at the location.

[0021] Furthermore, based on the vortex sea surface temperature anomalies caused by the advection process and the vertical motion process, the relative contributions of the vortex sea surface temperature anomalies caused by the advection process to the vortex sea surface temperature anomalies, and the relative contributions of the vortex sea surface temperature anomalies caused by the vertical motion process to the vortex sea surface temperature anomalies are calculated respectively.

[0022] Furthermore, the relative contribution of the eddy sea surface temperature anomaly caused by the advection process to the eddy sea surface temperature anomaly is as follows:

[0023]

[0024] The relative contribution of the eddy sea surface temperature anomaly caused by vertical motion to the eddy sea surface temperature anomaly is:

[0025]

[0026] in, For ocean grid points At the location, the vortex sea surface temperature anomaly caused by the advection process is determined according to the optimal inversion parameters; For ocean grid points The vertical motion process at the location determined according to the optimal inversion parameters results in vortex sea surface temperature anomalies.

[0027] The dynamic process decomposition system for ocean eddy sea surface temperature anomalies described in this invention comprises:

[0028] Dynamic process decomposition simulation unit: For vortex regions with sea surface temperature anomalies, the advection process is simulated considering diffusion and damping effects, and the vertical motion process is simulated considering the rising and sinking flows at the center of the vortex.

[0029] Overall Simulation and Parameter Inversion Unit: Based on the advection process simulation and vertical motion process simulation, an overall simulation model of the vortex sea surface temperature anomaly is established, and parameter inversion is performed to obtain the optimal inversion parameters;

[0030] Vortex sea surface temperature anomaly reconstruction unit: Based on the optimal inversion parameters, determine the vortex sea surface temperature anomalies caused by the advection process and the vertical motion process, respectively.

[0031] The computer-readable storage medium of the present invention stores a computer program, which, when executed by a processor, implements the method for decomposing the dynamic process of sea surface temperature anomalies in ocean eddies.

[0032] The computer program product of the present invention includes a computer program that, when executed by a processor, implements a method for decomposing the dynamic process of sea surface temperature anomalies in ocean eddies.

[0033] Beneficial effects: Compared with the prior art, the advantages of the present invention are: (1) The present invention decomposes sea surface temperature anomalies into advection transport processes and vertical motion processes from the perspective of dynamic processes, which can clearly reflect the physical source of temperature anomalies and has physical interpretability. (2) The present invention can invert parameters such as damping coefficient, diffusion coefficient, vertical motion intensity, and influence range by parameter inversion, and can invert the optimal parameters that are close to the actual observation data, avoiding the error caused by setting fixed values ​​based on experience, and improving the accuracy of dynamic process decomposition. (3) The present invention can be used for batch analysis of satellite remote sensing sea surface temperature data to reveal the differences in the causes of eddy sea surface temperature anomalies in different sea areas; at the same time, it can be extended to the inversion of marine salinity, chlorophyll and other structures. Attached Figure Description

[0034] Figure 1 This is a flowchart of the dynamic process decomposition method of the present invention.

[0035] Figure 2 This is a schematic diagram of sea surface temperature anomaly field observation data according to an embodiment of the present invention.

[0036] Figure 3 This is a schematic diagram of sea surface temperature anomalies caused by advection processes in an embodiment of the present invention.

[0037] Figure 4 This is a schematic diagram of sea surface temperature anomalies caused by vertical motion in an embodiment of the present invention.

[0038] Figure 5 This is a diagram showing the reconstruction results of the vortex sea surface temperature anomaly according to an embodiment of the present invention. Detailed Implementation

[0039] The technical solution of the present invention will be further described below with reference to the accompanying drawings.

[0040] like Figure 1 As shown, the method for decomposing the dynamic process of ocean eddy sea surface temperature anomalies includes the following steps.

[0041] Step 1: Obtain input data.

[0042] Step 1.1: Obtain observation data on sea surface height and sea surface temperature in the vortex region.

[0043] Step 1.2: Calculate the sea surface height anomaly and geostrophic flow field, sea surface background temperature field, and sea surface temperature anomaly field in the vortex region.

[0044] Step 2: Simulation of advection process and vertical motion process.

[0045] 1) Considering diffusion and damping effects during advection, the eddy sea surface temperature anomaly Satisfy the following equation:

[0046] (1)

[0047] in The diffusion coefficient is... The damping coefficient is... For ocean grid points The vortex background temperature field at that location, and These represent the flow velocities in the east-west and north-south directions, respectively. For the Nabla operator, For time.

[0048] Assuming the initial eddy sea surface temperature anomaly is zero, equation (1) is integrated using a numerical method over a fixed time. Sea surface temperature anomalies caused by advection can be obtained. .

[0049] 2) Upwelling or downwelling currents at the center of the eddy cause sea surface temperature anomalies. The sea surface temperature anomalies caused by the vertical movement of seawater are modeled using the following function:

[0050] (2)

[0051] in Reaction vertical motion intensity The scope of the reaction's influence For ocean grid points The sea level is abnormal.

[0052] Step 3: Overall simulation model and parameter inversion of eddy sea surface temperature anomalies.

[0053] Step 3.1: The overall eddy sea surface temperature anomaly can be represented by horizontal current processes and vertical motion processes.

[0054] (3)

[0055] Step 3.2, parameter inversion, searching Minimize the cost function:

[0056] (4)

[0057] in For ocean grid points The observed sea surface temperature anomalies at the location were used. Numerical optimization methods, such as the least squares method, were employed for parameter inversion.

[0058] Step 4: Calculation of sea surface temperature anomalies during advection and vertical motion processes.

[0059] 1) Invert the optimal parameters , , Substituting into equation (1) yields the vortex sea surface temperature anomaly caused by the advection process. .

[0060] 2) Invert the optimal parameters , Substituting into equation (2) yields the vortex sea surface temperature anomaly caused by the vertical motion process. .

[0061] In some optional embodiments, error analysis can be performed on the reconstructed results of the eddy sea surface temperature anomaly and the observed data. The reconstructed results of the eddy sea surface temperature anomaly can be given by the following formula:

[0062] (5)

[0063] Reconstructed results and observation data Error analysis can employ commonly used statistical methods such as mean square error, root mean square error, relative error, and correlation coefficient. For example, using root mean square error:

[0064] (6)

[0065] in This represents the number of grid point data points.

[0066] In some optional embodiments, the contributions of advection and vertical motion can be quantified. The relative contributions of advection and vertical motion to the eddy sea surface temperature anomaly are given by the following equation:

[0067] (7)

[0068] (8)

[0069] The closer the contribution ratio is to 1, the more significant the impact of the dynamic process on the eddy sea surface temperature anomaly.

[0070] The method described in this invention will be verified through specific experiments below.

[0071] This experiment was conducted in the Northwest Pacific Ocean ( , An anticyclonic vortex was conducted in the region on June 18, 2015.

[0072] 1. Obtain input data

[0073] The required data includes observational and climatological data. Sea surface height observations are from AVISO, with a temporal resolution of 1 day and a spatial resolution of [missing information]. The sea surface temperature observations used the National Oceanic and Atmospheric Administration's Optimal Interpolated Sea Surface Temperature (OISST) v2.1 version, with a temporal resolution of 1 day and a spatial resolution of [missing information]. The climatological monthly mean sea surface temperature data were extracted from the World Ocean Atlas (WOA) database, and the climatological monthly mean sea surface height was calculated from the monthly mean of ocean satellite remote sensing data from January 1, 2010 to December 31, 2019 using AVISO.

[0074] Interpolate observational data and climatological data to a unified spatial resolution. The grid points were used. The method described in patent CN202510433678.7 was employed to calculate the sea surface height anomaly in the eddy region. Geostrophic flow field Sea surface background temperature field Sea surface temperature anomaly field The geostrophic flow field can be calculated from sea level height data using the geostrophic equation. Sea surface temperature anomaly field. And the identified vortex center (pentagram) and 1 radius (dashed circle) as shown Figure 2 As shown.

[0075] 2. Simulation of advection and vertical motion processes: Establishing a general model of the eddy sea surface temperature anomaly, and using the least squares method for numerical optimization to obtain the optimal parameters. , , , , .

[0076] 3. Invert the optimal parameters , , Substituting into equation (1) yields the vortex sea surface temperature anomaly caused by the advection process. The result is as follows Figure 3 As shown. The optimal parameters will be inverted. , Substituting into equation (2) yields the vortex sea surface temperature anomaly caused by the vertical motion process. The result is as follows Figure 4 As shown.

[0077] 4. The reconstruction results of the eddy sea surface temperature anomaly are calculated using equation (5), and the results are as follows: Figure 5 As shown. The reconstructed results calculated using the root mean square error are compared with the observed data. The error is, and the result is According to equations (7) and (8), the relative contributions of advection and vertical motion to the eddy sea surface temperature anomaly are measured. After calculation, the contribution of advection is 41% and the contribution of vertical motion is 59%.

[0078] The results show that vertical motion mainly contributes to the temperature anomaly in the central region, while advection mainly contributes to the asymmetric structure in the periphery. The root mean square error of the model is about 0.2°C, which verifies the effectiveness and physical explanation capability of the present invention.

[0079] The dynamic process decomposition system for ocean eddy sea surface temperature anomalies described in this invention comprises:

[0080] Dynamic process decomposition simulation unit: For vortex regions with sea surface temperature anomalies, the advection process is simulated considering diffusion and damping effects, and the vertical motion process is simulated considering the rising and sinking flows at the center of the vortex.

[0081] Overall Simulation and Parameter Inversion Unit: Based on the advection process simulation and vertical motion process simulation, an overall simulation model of the vortex sea surface temperature anomaly is established, and parameter inversion is performed to obtain the optimal inversion parameters;

[0082] Vortex sea surface temperature anomaly reconstruction unit: Based on the optimal inversion parameters, determine the vortex sea surface temperature anomalies caused by the advection process and the vertical motion process, respectively.

[0083] The computer-readable storage medium of the present invention stores a computer program, which, when executed by a processor, implements the method for decomposing the dynamic process of sea surface temperature anomalies in ocean eddies.

[0084] The computer program product of the present invention includes a computer program that, when executed by a processor, implements a method for decomposing the dynamic process of sea surface temperature anomalies in ocean eddies.

[0085] The computer-readable storage medium may include RAM, ROM, EEPROM, CD-ROM or other optical disc storage devices, magnetic disk storage devices or other magnetic storage devices, flash memory or any other medium that can be used to store program code in the form of instructions or data structures and is accessible by a computer.

[0086] The processor is used to execute a computer program stored in memory to implement the various steps in the methods described in the above embodiments.

Claims

1. A method for decomposing the dynamic processes of sea surface temperature anomalies in ocean eddies, characterized in that, Includes the following steps: For vortex regions with sea surface temperature anomalies, the advection process is simulated considering diffusion and damping effects, and the vertical motion process is simulated considering the rising and sinking flows at the center of the vortex. Based on the advection process simulation and vertical motion process simulation, an overall simulation model of vortex sea surface temperature anomaly was established, and parameter inversion was performed to obtain the optimal inversion parameters; The sea surface temperature anomalies of the eddies caused by the advection process and the vertical motion process are determined based on the optimal inversion parameters. At marine grid points At this location, the vortex of the advection process causes sea surface temperature anomalies. The simulation model is as follows: in, The diffusion coefficient is... The damping coefficient is... For ocean grid points The vortex background temperature field at that location, and These represent the flow velocities in the east-west and north-south directions, respectively. For the Nabla operator, For time; At marine grid points At this location, the vertical motion process vortex causes sea surface temperature anomalies. The simulation model is as follows: in, Reaction vertical motion intensity The scope of the reaction's influence For ocean grid points The sea level is abnormally high; Overall simulation model of vortex sea surface temperature anomalies as follows: in, The integration time is the time for the advection process model; the search is achieved through parameter inversion. , , , , The optimal parameter values.

2. The method for decomposing the dynamic process of sea surface temperature anomalies in ocean eddies according to claim 1, characterized in that, The cost function is obtained by solving the numerical optimization method. , , , , The optimal parameter values ​​and the cost function are: in, For ocean grid points Sea surface temperature anomaly observed at the location.

3. The method for decomposing the dynamic processes of sea surface temperature anomalies in ocean eddies according to claim 1, characterized in that, Based on the vortex sea surface temperature anomalies caused by advection and vertical motion processes, the relative contributions of the vortex sea surface temperature anomalies caused by advection and vertical motion processes to the vortex sea surface temperature anomalies are calculated respectively.

4. The method for decomposing the dynamic process of sea surface temperature anomalies in ocean eddies according to claim 1, characterized in that, The relative contribution of the eddy sea surface temperature anomaly caused by the advection process to the eddy sea surface temperature anomaly is: The relative contribution of the eddy sea surface temperature anomaly caused by vertical motion to the eddy sea surface temperature anomaly is: in, For ocean grid points At the location, the vortex sea surface temperature anomaly caused by the advection process is determined according to the optimal inversion parameters; For ocean grid points The vertical motion process at the location determined according to the optimal inversion parameters results in vortex sea surface temperature anomalies.

5. A dynamic process decomposition system for ocean eddies and sea surface temperature anomalies, characterized in that, include: Dynamic process decomposition simulation unit: For vortex regions with sea surface temperature anomalies, the advection process is simulated considering diffusion and damping effects, and the vertical motion process is simulated considering the rising and sinking flows at the center of the vortex. Overall Simulation and Parameter Inversion Unit: Based on the advection process simulation and vertical motion process simulation, an overall simulation model of the vortex sea surface temperature anomaly is established, and parameter inversion is performed to obtain the optimal inversion parameters; The eddy sea surface temperature anomaly reconstruction unit: determines the eddy sea surface temperature anomalies caused by advection and vertical motion processes based on the optimal inversion parameters; and reconstructs the anomalies at ocean grid points. At this location, the vortex of the advection process causes sea surface temperature anomalies. The simulation model is as follows: in, The diffusion coefficient is... The damping coefficient is... For ocean grid points The vortex background temperature field at that location, and These represent the flow velocities in the east-west and north-south directions, respectively. For the Nabla operator, For time; At marine grid points At this location, the vertical motion process vortex causes sea surface temperature anomalies. The simulation model is as follows: in, Reaction vertical motion intensity The scope of the reaction's influence For ocean grid points The sea level is abnormally high; Overall simulation model of vortex sea surface temperature anomalies as follows: in, The integration time is the time for the advection process model; the search is achieved through parameter inversion. , , , , The optimal parameter values.

6. A computer-readable storage medium storing a computer program, characterized in that, When the computer program is executed by the processor, it implements the method for decomposing the dynamic process of sea surface temperature anomalies in ocean eddies according to any one of claims 1-4.

7. A computer program product, comprising a computer program, characterized in that, When the computer program is executed by the processor, it implements the method for decomposing the dynamic process of sea surface temperature anomalies in ocean eddies according to any one of claims 1-4.

Images