# Patent

## Title
Method of forecasting a magnetic or electrical environment from an ocean volume

## Abstract

A method of ocean property forecasting. The method includes generating one of a magnetic field prediction and an electric field prediction between a first observation point and a plurality of uniform ocean cells by generating a vector sum of magnetic field contributions from a plurality of magnetic field vectors and a vector sum of electric field contributions from a plurality of electric field vectors, respectively. The method further includes generating one of a magnetic field forecast and an electric field forecast of the ocean volume at a second observation point in the ocean volume by re-using the magnetic field vectors and electric field vectors and scaling one of the plurality of magnetic field vectors and the plurality of electric field vectors by a plurality of tri-axial equivalent electrical currents with identical geometry to the first observation point.

## Background
TECHNICAL FIELD
This invention relates generally to a method for forecasting an ocean property, and, more particularly, to a method for forecasting a magnetic or electrical environment of a ocean volume having one or more currents running therethrough.
BACKGROUND ART
A major source of extremely low frequency electromagnetic variations in the ocean is caused by the motion of the highly conductive water through the earth's magnetic field. These hydrodynamic variations affect magnetic and electric field sensors in the frequency range below 1 Hz. One of the first studies of these fields was performed by Longuet-Higgins et al., whose investigations were concerned with electric fields induced by the steady motion of seawater. See, e.g., Longuet-Higgins, M. S., M. E. Stem, and H. Stommel, “The Electric Field Induced by Ocean Currents and Waves, With Applications to the Method of Towed Electrodes,” Papers in Physical Oceanography and Meteorology XIII, 1, Massachusetts Institute of Technology and Woods Hole Oceanographic Institution, 1954, incorporated herein by reference. Surface waves, internal waves, solitary waves, tides, and ocean currents all produce observable magnetic and electric fields. Larson and Sanford showed that the Florida Current volume transport could be determined from cable voltage measurements of an undersea cable. See, e.g., Larsen, J. C. and T. B. Sanford, “Florida Current Volume Transports from Voltage Measurements,” Science, 227, 302-304, 1985, incorporated herein by reference.
A long period electromagnetic response of internal waves has been observed with magnetic sensors in the deep ocean, but there is little information about the character of the electric and magnetic field generated in shallow water areas. Deep water models for internal wave-induced magnetic spectra indicate that the amplitude increases with decreasing frequency.
Internal wave-induced magnetic fields for a two-layered ocean model were treated by Beal and Weaver. See, e.g., Beal, H. T. and Weaver, J. T., “Calculations of Magnetic Variations Induced by Internal Ocean Waves,” J. Geophys. Res., 75, no. 33, 1970, incorporated herein by reference. Podney followed with a more comprehensive treatment of internal waves for an exponentially stratified ocean with a horizontally uniform Brunt-Vaisala frequency profile. See, e.g., Podney, Walter, “Electromagnetic Fields Generated by Ocean Waves”, J. Geophys. Res., 80, No. 21, 1975, incorporated herein by reference. Wasylkiwskyj used a similar approach to derive solutions for the case where the Vaisala frequency profile decreases exponentially in a manner analogous to that used by Garrett and Munk (1972). See, e.g., Wasylkiwskyj, W., “Electromagnetic Fields Induced by Ocean Currents,” IDA Paper P-1399, IDA: Arlington, Va., 1979, incorporated herein by reference, and Garrett, C. and W. Munk, “Space-Time Scales of Internal Waves,” J. Geophys. Fluid Dynamics, 2, 225-264, 1972, incorporated herein by reference. The solutions derived by Wasylkiwskyj are solved for the case of an airborne sensor moving over the ocean surface. Later on Petersen and Poehls used Podney's formulation combined with the Garrett and Munk model to generate a spectral estimate of the magnetic induction. See, e.g., Petersen, R. A. and K. A. Poehls, “Model Spectrum of Magnetic Induction Caused by Ambient Internal Waves,” J. Geophys. Res., 87, No. C1, 433-440, 1982, incorporated herein by reference. Chave derived a somewhat more general solution for internal waves that used the Garrett and Munk wave spectra. See, e.g., Chave, A. D., “On the Electromagnetic Field Induced by Ocean Internal Waves,” J. Geophys. Res., 89, No. C6, 10519-10528, 1984, incorporated herein by reference. His derivation is applicable both within the ocean and on the ocean bottom and includes the effects of self and mutual induction, which should improve the model solution at low frequencies on the seafloor. He has also computed the results for a theoretical solitary internal wave. See e.g., Chave, A. D., “The Magnetic Effects of Shallow Water Internal Solutions,” Scripps Institute of Oceanography, Reference 86-7, 1986, incorporated herein by reference.
DISCLOSURE OF THE INVENTION
An embodiment of the invention includes a method of predicting magnetic field noise and electric field noise due to ocean flow hiding a target. A first observation point is provided. A three dimensional cellular model of ocean volume is provided. A plurality of earth's magnetic field vectors around the first observation point is provided. A three dimensional cellular model of ocean flow in the ocean volume is provided. A three dimensional cellular model of ocean electrical properties in the ocean volume is provided. From the three dimensional cellular model of vector ocean volume around the first observation point, the plurality of earth's magnetic field vectors around the first observation point, the plurality of three dimensional cellular model of ocean flow in the ocean volume; and the three dimensional cellular model of ocean electrical properties in the ocean volume, a plurality of uniform ocean cells between a water surface and a water bottom within the ocean volume around the first observation point with three dimensional ocean flow and with ocean electrical properties is generated. A plurality of tri-axial equivalent electrical currents corresponding to the plurality of uniform ocean cells within the ocean volume is generated. A plurality of unit tri-axial electric bipole transmitter sources corresponding to the plurality of uniform ocean cells is generated. One of a plurality of magnetic field vectors and a plurality of electric field vectors from the plurality of unit electric bipole transmitter sources is determined, wherein the water bottom and the water surface are approximated to have infinite horizontal dimensions. One of the plurality of magnetic field vectors and the plurality of electric field vectors is scaled by the plurality of corresponding tri-axial equivalent electrical currents. One of a magnetic field prediction and an electric field prediction is generated between the first observation point and the plurality of uniform ocean cells by generating a vector sum of magnetic field contributions from the plurality of magnetic field vectors and a vector sum of electric field contributions from the plurality of electric field vectors, respectively. One of a magnetic field forecast and an electric field forecast of the ocean volume is generated at a second observation point in the ocean volume by re-using the magnetic field vectors and electric field vectors and scaling one of the plurality of magnetic field vectors and the plurality of electric field vectors by the plurality of tri-axial equivalent electrical currents with identical geometry to the first observation point.
Optionally, the plurality of uniform ocean cells comprises a plurality of layers between the water bottom and the water surface.
Optionally, the plurality of uniform ocean cells are geometrically sampled such that the plurality of uniform ocean cells comprise a higher density of cells nearer to the first observation point than further away from the first observation point.
Optionally, each cell of the plurality of cells comprises a rectangular prismatic cross-section.
Optionally, the method further includes setting a detection threshold for a target based on one of the magnetic field forecast and the electric field forecast of the ocean volume at the observation point.



BRIEF DESCRIPTION OF THE DRAWINGS
FIG. 1 is a graph of a normalized contribution to the sum of the magnetic field for bipoles spaced 1 meter apart as a function of radial range from the observation point.
FIG. 2 is a three-dimensional model of a region of an ocean volume.
FIG. 3 is a map view for the geometric sampling scheme in a horizontal layer for source current elements where elements are spaced most densely near an observation point P. In this illustrative case, the geometric factor is three and there are three expansions or levels of increasing area size shown.



BEST MODES OF CARRYING OUT THE INVENTION
An embodiment of the invention is described as follows. The embodiment of the invention includes a method of predicting magnetic field noise and electric field noise due to ocean flow hiding a target. A first observation point is provided. A three dimensional cellular model of ocean volume is provided. A plurality of earth's magnetic field vectors around the first observation point is provided. A three dimensional cellular model of ocean flow in the ocean volume is provided. A three dimensional cellular model of ocean electrical properties in the ocean volume is provided. From the three dimensional cellular model of vector ocean volume around the first observation point, the plurality of earth's magnetic field vectors around the first observation point, the plurality of three dimensional cellular model of ocean flow in the ocean volume; and the three dimensional cellular model of ocean electrical properties in the ocean volume, a plurality of uniform ocean cells between a water surface and a water bottom within the ocean volume around the first observation point with three dimensional ocean flow and with ocean electrical properties is generated. A plurality of tri-axial equivalent electrical currents corresponding to the plurality of uniform ocean cells within the ocean volume is generated. A plurality of unit tri-axial electric bipole transmitter sources corresponding to the plurality of uniform ocean cells is generated. One of a plurality of magnetic field vectors and a plurality of electric field vectors from the plurality of unit electric bipole transmitter sources is determined, wherein the water bottom and the water surface are approximated to have infinite horizontal dimensions. One of the plurality of magnetic field vectors and the plurality of electric field vectors is scaled by the plurality of corresponding tri-axial equivalent electrical currents. One of a magnetic field prediction and an electric field prediction is generated between the first observation point and the plurality of uniform ocean cells by generating a vector sum of magnetic field contributions from the plurality of magnetic field vectors and a vector sum of electric field contributions from the plurality of electric field vectors, respectively. One of a magnetic field forecast and an electric field forecast of the ocean volume is generated at a second observation point in the ocean volume by re-using the magnetic field vectors and electric field vectors and scaling one of the plurality of magnetic field vectors and the plurality of electric field vectors by the plurality of equivalent electrical currents with identical geometry to the first observation point.
Optionally, the plurality of uniform ocean cells comprises a plurality of layers between the water bottom and the water surface.
Optionally, the plurality of uniform ocean cells are geometrically sampled such that the plurality of uniform ocean cells comprise a higher density of cells nearer to the first observation point than further away from the first observation point.
Optionally, each cell of the plurality of cells comprises one of rectangular cross-section and an annular cross-section.
Optionally, the method further includes setting a detection threshold for a target based on one of the magnetic field forecast and the electric field forecast of the ocean volume at the observation point.
The Foundational EM Model
Another embodiment of the invention is described in the foundational EM model as follows. The purpose of an embodiment of this invention is to efficiently compute the magnetic and electric fields due to moving seawater from a model ocean. A method according to such an embodiment is an integral equation approach to model the electric and magnetic field from fluid motion of seawater. The method takes the fluid velocity and electrical conductivity for a fluid cell along with the earth's magnetic field at the location and compute the associated magnetic and electric fields at some observation point. Vector fluid velocity and conductivity values are obtained from a time step of the finite volume general circulation model like the NRL-MIT nonhydrostatic ocean model. A description of the NRL-MIT nonhydrostatic ocean model is found, for example, in C. Gallacher, Michael Schaferkotter, and Will Avera, “Nonhydrostatic Hindcasts of High Amplitude Internal Waves in the Mid-Atlantic Bight,” Proceedings of the OCEANS 2007 MTS/IEEE Vancouver Conference, incorporated herein by reference.
Computation of the magnetic field due to moving seawater has usually been limited to specific cases that have an analytic theoretical solution or power spectral characterization. In this embodiment of the invention, a method is provided for efficiently calculating the magnetic fields using products from a finite volume ocean flow model. Advantages of this method for computing the electric and magnetic fields from moving seawater include: utilizing Navy oceanographic circulation models, accounting for complex oceanographic conditions, accounting for complex ocean bathymetric and coastal conditions, taking into account the water bottom and surface conductivity contrast; taking into account the conductivity structure within the water, approximating infinite horizontal boundaries for the water, the bottom, and air (i.e., no electromagnetic (“EM”) reflections from the edge of the model. Finite difference EM models have a limited number of cells to approximate the computational edge of the model). This method is further computationally efficient because its algorithms do not require a computer to compute values at all of the nodes within a model space to obtain the fields for a single point whereas finite difference and finite element EM models require computing fields at all node locations.
The following discussion of this embodiment of the invention will focus on the fields produced by nonlinear internal ocean waves. However, it should be understood to those of ordinary skill in the art that the method can be applicable to all ocean flow dynamics.
The electro-magnetic (“EM”) model computations for this embodiment of the invention are dependent on the velocity and conductivity properties derived from the NRL-MIT nonhydrostatic model. The nonhydrostatic model solves for the finite volume representation of the fluid motion as a function of time and space. At each time step, there are three orthogonal velocity components associated with the fluid motion for each fluid volume element. From these velocity components, the effective electromotive force across a volume element can be computed using Eq. (1),

Ēi= Vi× B.  (1)

Where B is the vector magnetic field of the Earth, Vi is the vector velocity of the ith water volume element and Ēi is the vector electric field across the ith volume element. An overstrike line indicates a vector quantity. The current density Ji (A/m2) for the ith fluid element with conductivity σi is then,

 Ji=σi( Vi× B).  (2)

Several approximations are used to simplify the hydrodynamic model fluid element for the EM model computations. The magnetic field generated by the flow is very small compared to the Earth's main field and is assumed to not modify the primary field B. No magnetic field gradients are included such that the Earth's field is assumed to be constant across the model domain. In addition, due to the limited size of the finite volume model, assumptions about the conductivity and velocity at the horizontal boundaries are required to properly compute the EM fields near the boundaries.
For each fluid element of the model, there is an electric current across the element that is associated with the electric field generated by the fluid motion. Because the velocity is known during a time step of the hydrodynamic model, the electric current across the fluid element during the time step can be computed. In order to model the electric and magnetic fields for any field observation (measurement) point outside the fluid element, the algorithm will compute the fields based on the electric current from each fluid element near the observation point and sum the vector fields together from all the fluid elements. In this context, ‘near’ the observation point will be quantified later in the discussion. To simplify the EM model and capitalize on existing well-tested computational techniques and software, (e.g., Kong, J. A., “Electromagnetic Fields Due To Dipole Antennas Over Stratified Anisotropic Media,” Geophysics, 37, n. 6, 985-996, December 1972, incorporated herein by reference, and Xiong, Z., “Electromagnetic fields of electric dipoles embedded in a stratified anisotropic earth,” Geophysics, 54, no. 12, 1643-1646, December 1989, incorporated herein by reference), a bipole representation of the electric current in the fluid element is used to compute the EM fields. While this is not truly representative of the electric geometry for a square fluid element, the bipole representation is found to be a very good approximation provided that the current is small and scaled appropriately. In addition, problems with singularities are avoided by introducing the bipole representation at the center of the fluid elements and restricting field observation points to the edges of the fluid elements.
In EM model calculations, the fluid element can be replaced with three orthogonal bipole line sources each having an equivalent bipole moment ( M) associated with the vector current density for the fluid element,

 M=∫ Jdv  (3)

Then for a rectangular fluid element having dimensions of ΔxΔyΔz for which the current density is constant, the equivalent moment would be:

 M= JΔvΔyΔz  (4)

Each bipole line source is positioned at the center of the fluid element and having a length Δd equal to the related dimension of the fluid element. Then, for a bipole of length Δd and current I,

 Mi=ĪiΔd.  (5)

Thus, for a bipole representation of the fluid element, the current Ii in the equivalent bipole should be scaled to,

 JiΔxΔyΔz=ĪiΔd  (6)

And solving for the tri-axial (vector) equivalent electric current Ii,





  
    
      
        
          
            I
            _
          
          i
        
        =
        
          
            
              
                J
                _
              
              i
            
            ⁢
            Δ
            ⁢
            
              
            
            ⁢
            x
            ⁢
            
              
            
            ⁢
            Δ
            ⁢
            
              
            
            ⁢
            y
            ⁢
            
              
            
            ⁢
            Δ
            ⁢
            
              
            
            ⁢
            z
          
          
            Δ
            ⁢
            
              
            
            ⁢
            d
          
        
      
    
    
      
        (
        7
        )
      
    
  




The EM computation model described here uses an integral equation solution for a bipole source in a layered medium to solve for the EM fields at an observation point. The integral equation technique is described by Kong (1972) and the computational technique is described by Xiong (1989).
Using this technique reduces the computational load because the primary field and reflections of the sources due to the horizontal layers (e.g., the surface of water, the water bottom, or the water subbottom) are included in the solution. As a result, solutions need only be computed for the lateral extent of the sources within the water associated with the fluid motion and only to the extent that they have a significant contribution at the field observation point. The model configuration (e.g., air, water, bottom) described herein has a water thickness much less than the horizontal extent of the model, and no motion induced sources exist in the air or the bottom. For this configuration, it is more efficient to compute the fields at an observation point from fluid elements in layers, and the layers are summed together to compute the resulting vector field at the observation point. In this sum of sources throughout the water volume, the current in each bipole source in each fluid element is scaled by Eq. (7) for the hydrodynamic parameters of that element.
The integral equation model takes advantage of the numerical computer codes that have been developed and thoroughly tested as described above in the Kong (1972) and the Xiong (1989) references. These computer codes use a Green's function approach to solve for the electric and magnetic fields in a layered half space. For an observation point P, at a distance r from the origin, the magnetic ( H) and electric (Ē) field due to a current source J at a distance ri can be computed using the integral equation technique for a layered media. If the current source is selected to be unity (1 A/m2), then the resulting fields H1 (rp:ri) and Ē1(rp:ri) represent a scaling function for the field contributions at point P from the bipole source at ri. Combining with Eq. (7), the fields at point P from the bipole source at r, are then

 Hi(rp)= H1(rp:ri)·Īi,Ēi(rp)=Ē1(rp:ri)·Īi  (8)

The magnetic and electric field at point P is then the vector sum of all the field contributions from all the bipole sources.





  
    
      
        
          
            
              H
              _
            
            ⁡
            
              (
              
                r
                p
              
              )
            
          
          =
          
            
              ∑
              
                i
                =
                1
              
              N
            
            ⁢
            
              
            
            ⁢
            
              
                
                  H
                  _
                
                i
              
              ⁡
              
                (
                
                  
                    r
                    P
                  
                  :
                  
                    r
                    i
                  
                
                )
              
            
          
        
        ,
        
          
            
              E
              _
            
            ⁡
            
              (
              
                r
                P
              
              )
            
          
          =
          
            
              ∑
              
                i
                =
                1
              
              N
            
            ⁢
            
              
                
                  E
                  _
                
                i
              
              ⁡
              
                (
                
                  
                    r
                    P
                  
                  :
                  
                    r
                    i
                  
                
                )
              
            
          
        
      
    
    
      
        (
        9
        )
      
    
  




The contributions to the vector sum in Eq. (9) come from all the source current elements within the water volume. However, the magnetic and electric field contribution from any particular bipole decreases rapidly with distance from the source. As a result, it is adequate to only sum the contributions within some region about the observation point. If the source contributions are coherent over large scales, then the sum in Eq. (9) will increase with the distance (i.e., with increased volume) away from the observation point and will converge to a finite value. An example of this is shown in FIG. 1 for the summed magnitude relative to a uniform grid of coherent current sources plotted as a function of distance from the observation point. In these terms, coherence can be related to the scale of the fluid flow phenomena in the model. For example, tidal flow would have scales much larger than the submesoscale ocean model and would contribute in a coherent way to the fields at an observation point. So, in general, the fields from a particular ocean flow can be considered as coherent for scales less than the scale of the ocean feature. For internal wave features with scales on the order of a kilometer or less, the sum would be roughly coherent at scales less than a kilometer and less coherent for larger scales. This affects the computations by providing some guidelines as to how large an aerial coverage the sum in Eq. (9) needs to cover to get an accurate answer. Based on the observed convergence in FIG. 1 and consideration for the scales of the ocean features of interest, the sum in Eq. (9) should encompass two kilometers or more from the observation point.
As previously discussed, the focus of this embodiment of the instant invention is to calculate fields for submesoscale internal ocean flow (e.g., 10 kilometers or less) and particularly for the case where the water depth is significantly less than the horizontal scale of the model. This suggests that the aspect ratio for the model has a horizontal dimension that is on the order of ten times or more relative to the vertical dimension. Taking into consideration the previous discussion and that the water is in general less than 1 kilometer deep, then the sum in Eq. (9) should include the entire water column vertical dimension. It is computationally efficient to construct the sum in Eq. (9) as two sums: a first sum over a horizontal layer of source current elements with length scales on the order of several kilometers, and a second sum over all the layers within the water column. Breaking the computations into two sums simplifies the complexity of determining sample points for the sum. Equation (9) then becomes,




  
    
      
        
          
            
              H
              _
            
            ⁡
            
              (
              
                r
                p
              
              )
            
          
          =
          
            
              ∑
              
                j
                =
                1
              
              
                N
                ⁢
                
                  
                
                ⁢
                v
              
            
            ⁢
            
              
            
            ⁢
            
              
                ∑
                
                  i
                  =
                  1
                
                
                  N
                  ⁢
                  
                    
                  
                  ⁢
                  h
                
              
              ⁢
              
                
              
              ⁢
              
                
                  
                    H
                    _
                  
                  
                    i
                    ⁢
                    
                      
                    
                    ⁢
                    j
                  
                
                ⁡
                
                  (
                  
                    
                      r
                      P
                    
                    :
                    
                      r
                      
                        i
                        ⁢
                        
                          
                        
                        ⁢
                        j
                      
                    
                  
                  )
                
              
            
          
        
        ,
        
          
            
              E
              _
            
            ⁡
            
              (
              
                r
                P
              
              )
            
          
          =
          
            
              ∑
              
                j
                =
                1
              
              
                N
                ⁢
                
                  
                
                ⁢
                v
              
            
            ⁢
            
              
            
            ⁢
            
              
                ∑
                
                  i
                  =
                  1
                
                
                  N
                  ⁢
                  
                    
                  
                  ⁢
                  h
                
              
              ⁢
              
                
              
              ⁢
              
                
                  
                    E
                    _
                  
                  
                    i
                    ⁢
                    
                      
                    
                    ⁢
                    j
                  
                
                ⁡
                
                  (
                  
                    
                      r
                      P
                    
                    :
                    
                      r
                      
                        i
                        ⁢
                        
                          
                        
                        ⁢
                        j
                      
                    
                  
                  )
                
              
            
          
        
      
    
    
      
        (
        10
        )
      
    
  




where Nh is the number of source current elements summed within the horizontal layer and Nv is the number of horizontal layers summed within the vertical water column. Hij are the respective contributions for the individual current elements (i) within each layer and for the layers (j) representing the water volume. The summation is now discretized so that the current element contributions within each layer can be spaced more efficiently to reduce the computational load.

To reduce the number of computations while maintaining the accuracy, the source current elements can be spaced using a geometric sampling density scheme within each layer. This samples the volume close to the source at every one meter cube fluid element cell and efficiently reduces the sampling density further away from the observation point in a way to completely fill the area of the layer without overlap. The geometric scheme chosen computes the X and Y positions for the center point of each cubic cell within a layer as a function of an integer counter. The positions are computed using,




  
    
      
        
          P
          ⁢
          
            
          
          ⁢
          x
          ⁢
          
            
          
          ⁢
          
            y
            
              n
              ⁢
              
                
              
              ⁢
              j
            
          
        
        =
        
          
            
              G
              n
            
            2
          
          +
          
            
              G
              n
            
            ·
            j
          
        
      
    
    
      
        (
        11
        )
      
    
  




where Pxynj is the x or y position of the source current element relative to the observation point P, and identified in FIG. 2 as reference number 70. G is the geometric factor, n is the number of expansions or levels, and j is an integer counter with values from j=0 to j=n−1. G must be an odd number for the source current elements to always fall in the center of a cell. An example of a geometric sampling density scheme 50 is shown by way of example in FIG. 2 for the case for G=3 and n=3.

Application of the Foundational EM Model

An illustrative application of the foundational EM model discussed above is fairly straightforward. Once the fields are computed with Eq. (8), then the fields are summed together with Eq. (10) to generate the resulting field at point P. Most of the computational work is done in computing the H1 (rp:ri) and Ē1 (rp:ri) scaling functions in Eq. (8) and in prepping the oceanographic model parameters (e.g., vector velocity and conductivity) into uniform cells of similar size. An illustrative example of an ocean volume 10 of the oceanographic fluid cells of uniform size is given in FIG. 3. Fluid cell 50 is an illustrative oceanographic fluid cell in ocean volume 10. Ocean volume 10 includes a region of ocean volume subdivided into uniform fluid cells within which an arrow is shown representing one of three perpendicular bipole sources for each fluid cell. A dot 20 on the face of each fluid cells of ocean volume 10 is the projection of the bipole arrow onto the face of the visible fluid cells. An illustrative field measurement (i.e., observation) point 30 is also shown. Field measurement point 30 is a point where the magnetic and/or electric field prediction for the ocean volume 10 is desired. The lines 40 emanating from the ocean volume 10 to the field measurement point reference 30 represent the magnetic and/or electric field contribution from each of the bipoles in each of the cells for the ocean volume 10. For a given ocean volume 10, the relationship between the water velocity for a given fluid cell, the electric current field for the given fluid cell, and the magnetic field for the given fluid cell is shown illustratively by arrows 60 adjacent to the ocean volume 10. While it is not essential that the ocean volume to be divided into fluid cells of a uniform size, it greatly simplifies the computational bookkeeping. Accordingly, it should be understood that in alternative embodiments of the invention, the ocean volume is optionally divided into non-uniform fluid cells.
The functions H1(rp:ri) and Ē1(rp:ri) are computed with a layered conductivity model consisting of, for example, air, water, and ocean bottom. The water conductivity is, for example, divided into sub-layers, but the primary conductivity contrast occurs at the air and bottom interfaces. So, in an embodiment of the invention, choosing an average value for the water works best for computing these functions. In another embodiment of the invention, choosing another value for the water is also acceptable.
The function Īi in Eq. (8) is computed for each fluid element cell in the oceanographic model, described above. It is possible to evaluate the relative importance of the parameters used to calculate Īi and determine which parameters are most significant to the accuracy of the magnetic field computations. Looking at a scale analysis of the terms in Eq. (2), the conductivity σ, velocity V, and earth's primary field B are the main environmental parameters. Variations in these parameters across a region of open ocean (e.g., tens of kilometers) are on the order of or less than 1 S/m out of 5 S/m for the conductivity, on the order of up to ±1 m/s (for a vector direction) for the velocity, and on the order of a few nT to 100's of nT out of 50,000 nT for the Earth's magnetic field. Based on these order of magnitude variations, the parameter variation in Eq (2) is in general on the order of:




  
    
      
        
          
            δσ
            σ
          
          ≅
          0.2
        
        ,
        
          
            
              δ
              ⁢
              
                V
                _
              
            
            V
          
          ≅
          2
        
        ,
        
          
            
              δ
              ⁢
              
                B
                _
              
            
            B
          
          ≅
          0.0
        
      
    
    
      
        (
        12
        )
      
    
  




Therefore, B can be treated as a constant across the model in the source current element computations. However, σ and V are controlling variables for computing the equivalent bipole scaling currents.

In general, the hydrodynamic elements of the oceanographic model are sized for optimum computation of the vector fluid velocities. To generate a consistent fluid element size in the source current element computations, the σ and V parameters are, for example, gridded into uniform cells with a standard gridding routine to prep the ocean model. As mentioned above, other embodiments of the invention include non-uniform fluid cells. A grid cell size of one meter cube has proved simple to work with and is conceptually easy to check. But, of course, other grid cell sizes or shapes are consistent with the instant invention are employed in alternative embodiments of the instant invention, depending on the application desired. Once the ocean model parameters are gridded, the equivalent bipole scaling currents Iij are computed using Eq. (7) for each cell in the model, and then Eq. (10) is computed by summing up all the contributing source elements.
Illustrative alternatives to calculate the magnetic and electric fields for moving seawater are as follows. Analytic solutions discussed above used spectral estimates for the internal wave dynamics. However, instead of using spectral estimates, a simplified Biot-Savart type calculation for the fluid elements is beneficially used in an alternative embodiment of the invention. (See, e.g., Frederick J. Buesche, Introduction to Physics for Scientist and Engineers, McGraw-Hill, 1975, incorporated herein by reference). This simplified type of calculation ignores the distortions of the fields from the surface and bottom boundaries, which can amount to errors on the order of 20%. Instead of spectral estimates or simplified Biot-Savart type calculations, another embodiment of the invention employs another numerical technique using finite element electromagnetic models to calculate the fields due to moving seawater. An example of such a numerical technique is discussed in Chaillout, J. J., V. Poulbot, J. Berthier, and R. Blanpain, “Finite Element Calculation of Electromagnetic Fields Due to Moving Sea Water,” IEEE Transactions On Magnetics, V. 31, no. 3, May 1995, incorporated herein by reference In using a finite element model, an oceanographic model is divided into three elementary models and the result of such division is combined.
For completeness, one or more of the procedures described above are optionally implemented in terms of computer-implemented methods executed on a computer or a network of computers. For example, a computer-implemented method of an embodiment of the invention is stored on a computer-readable medium that is executable on such a computer or network of computers. A computer-readable medium storing such a computer-implemented method according to an embodiment of the invention is an example of a computer program product.
Obviously, many modifications and variations of the present invention are possible in light of the above teachings without departing from the true scope and spirit of the invention. It is therefore to be understood that the scope of the invention should be determined by referring to the following appended claims.

## Claims

1. A computer-implemented method of predicting magnetic field noise and electric field noise due to ocean flow hiding a target, the method comprising:
providing a first observation point;
providing a three dimensional cellular model of ocean volume;
providing a plurality of earth's magnetic field vectors around the first observation point;
providing a three dimensional cellular model of ocean flow in the ocean volume;
providing a three dimensional cellular model of ocean electrical properties in the ocean volume;
generating, from the three dimensional cellular model of ocean volume around the first observation point, the plurality of earth's magnetic field vectors around the first observation point, the three dimensional cellular model of vector ocean flow in the ocean volume; and the three dimensional cellular model of ocean electrical properties in the ocean volume, a plurality of uniform ocean cells between a water surface and a water bottom within the ocean volume around the first observation point with three dimensional ocean flow and with ocean electrical properties;
generating a plurality of tri-axial equivalent electrical currents corresponding to the plurality of uniform ocean cells within the ocean volume;
generating a plurality of unit tri-axial electric bipole transmitter sources corresponding to the plurality of uniform ocean cells;
determining one of a plurality of magnetic field vectors and a plurality of electric field vectors from the plurality of unit electric bipole transmitter sources, wherein the water bottom and the water surface are approximated to have infinite horizontal dimensions;
scaling one of the plurality of magnetic field vectors and the plurality of electric field vectors by the plurality of corresponding tri-axial equivalent electrical currents;
generating one of a magnetic field prediction and an electric field prediction between the first observation point and the plurality of uniform ocean cells by generating a vector sum of magnetic field contributions from the plurality of magnetic field vectors and a vector sum of electric field contributions from the plurality of electric field vectors, respectively;
generating one of a magnetic field forecast and an electric field forecast of the ocean volume at a second observation point in the ocean volume by re-using the magnetic field vectors and electric field vectors and scaling one of the plurality of magnetic field vectors and the plurality of electric field vectors by the plurality of tri-axial equivalent electrical currents with identical geometry to the first observation point.

2. The method according to claim 1, wherein the plurality of uniform ocean cells comprises a plurality of layers between the water bottom and the water surface.

3. The method according to claim 1, wherein the plurality of uniform ocean cells are geometrically sampled such that the plurality of uniform ocean cells comprise a higher density of cells nearer to the first observation point than further away from the first observation point.

4. The method according to claim 1, wherein each cell of the plurality of cells comprises a rectangular prismatic cross-section.

5. The method according to claim 1, further comprising:
setting a detection threshold for a target based on one of the magnetic field forecast and the electric field forecast of the ocean volume at the observation point.

6. A computer-readable medium storing instructions being executed by a computer, the instructions instructing the computer to generate a magnetic field forecast from an ocean volume, the instructions comprising:
providing a first observation point;
providing a three dimensional cellular model of ocean volume;
providing a plurality of earth's magnetic field vectors around the first observation point;
providing a three dimensional cellular model of ocean flow in the ocean volume;
providing a three dimensional cellular model of ocean electrical properties in the ocean volume;
generating, from the three dimensional cellular model of ocean volume around the first observation point, the plurality of earth's magnetic field vectors around the first observation point, the three dimensional cellular model of vector ocean flow in the ocean volume; and the three dimensional cellular model of ocean electrical properties in the ocean volume, a plurality of uniform ocean cells between a water surface and a water bottom within the ocean volume around the first observation point with three dimensional ocean flow, and with ocean electrical properties;
generating a plurality of tri-axial equivalent electrical currents corresponding to the plurality of uniform ocean cells within the ocean volume;
generating a plurality of unit tri-axial electric bipole transmitter sources corresponding to the plurality of uniform ocean cells;
determining one of a plurality of magnetic field vectors and a plurality of electric field vectors from the plurality of unit electric bipole transmitter sources, wherein the water bottom and the water surface are approximated to have infinite horizontal dimensions;
scaling one of the plurality of magnetic field vectors and the plurality of electric field vectors by the plurality of corresponding tri-axial equivalent electrical currents;
generating one of a magnetic field prediction and an electric field prediction between the first observation point and the plurality of uniform ocean cells by generating a vector sum of magnetic field contributions from the plurality of magnetic field vectors and a vector sum of electric field contributions from the plurality of electric field vectors, respectively;
generating one of a magnetic field forecast and an electric field forecast of the ocean volume at a second observation point in the ocean volume by re-using the magnetic field vectors and electric field vectors and scaling one of the plurality of magnetic field vectors and the plurality of electric field vectors by the plurality of tri-axial equivalent electrical currents with identical geometry to the first observation point.

7. The computer-readable medium according to claim 6, wherein the plurality of uniform ocean cells comprises a plurality of layers between the water bottom and the water surface.

8. The computer-readable medium according to claim 6, wherein the plurality of uniform ocean cells are geometrically sampled such that the plurality of uniform ocean cells comprise a higher density of cells nearer to the first observation point than further away from the first observation point.

9. The computer-readable medium according to claim 6, wherein each cell of the plurality of cells comprises a rectangular prismatic cross-section.

10. The computer-readable medium according to claim 6, further comprising:
setting a detection threshold for a target based on one of the magnetic field forecast and the electric field forecast of the ocean volume at the observation point.