Abstract
We review the range of applications and use of multi spacecraft techniques, applicable to close formation arrays of spacecraft, focusing on spatial gradient based methods, and the curlometer in particular. The curlometer was originally applied to Cluster multi-spacecraft magnetic field data, but later was updated for different environments and measurement constraints such as the NASA MMS mission, small-scale formation of 4 spacecraft; the 3 spacecraft configurations of the NASA THEMIS mision, and derived 2-4 point measurements from the ESA Swarm mission. In general, spatial gradient based methods are adaptable to a range of multi-point and multi-scale arrays. We also review the range of other techniques based on the computation of magnetic field gradients and magnetic field topology in general, including: magnetic rotation analysis and various least squares approaches. We review Taylor expansion methodology (FOTE), in particular, which has also been applied to both Cluster and MMS constellations, as well as interpretation of simulations. Four-point estimates of magnetic gradients are limited by uncertainties in spacecraft separations and the magnetic field, as well as the presence of non-linear gradients and temporal evolution. Nevertheless, the techniques can be reliable in many magnetospheric regions where time stationarity is largely applicable, or when properties of the morphology can be assumed (for example, the expected orientation of underlying large-scale structure). Many magnetospheric regions have been investigated directly (illustrated here by the magnetopause, ring current and field-aligned currents at high and low altitudes), and options for variable numbers of spacecraft have been considered. The comparative use of plasma measurements and possible new methodology for arrays of spacecraft greater than four are also considered briefly.
1 Introduction
The four Cluster II spacecraft (Escoubet et al., 2001) allowed 3-D structure and temporal evolution to be probed through the development of multi-spacecraft techniques for the first time. These techniques specifically allow spatial gradients of key quantities to be analyzed, typically through first order approximations, or Taylor expansion around measurement points (e.g., Fu et al., 2015). Such analysis was first described in a book on collected multi-spacecraft analysis techniques in 1998 (Paschmann and Daly, 1998). The application of these multi-spacecraft methods was updated in (Paschmann and Daly, 2008), where here we focus on magnetic gradients and specifically magnetic currents (e.g., Dunlop and Eastwood, 2008; Shen and Dunlop, 2008; Vogt et al., 2008). The operation of Cluster provided a spacecraft configuration which maintained a quasi-tetrahedral formation for much of its life and Cluster is still the only space physics mission to provide fully four-point coverage over a large spatial range of scales; over a time epoch of two solar cycles (see Dunlop et al., 2021b; Escoubet et al., 2021).
The methodology, covering a wide range of analysis techniques, has been continuously developed to determine key quantities and investigate a large number of phenomena and has been applied to other missions in the 22 years since then (such as MMS, THEMIS and Swarm). The Magnetospheric Multi-Scale (MMS) mission maintained a close four spacecraft configuration on smaller separation scales (a few km) than Cluster for much of its orbit (Burch et al., 2016), while during extended operations (Angelopoulos, 2008), some of the NASA THEMIS spacecraft flew in a 3-spacecraft configuration in the magnetosphere and the ESA Swarm low orbit (LEO) polar mission provided both 2 and 3 spacecraft measurements in close formations (Friis-Christensen et al., 2008) on meso-scales (∼100 km).
The future interest in multi-spacecraft methods remains strong; particularly in their development to make best use of planned larger arrays of spacecraft capable of probing multiscale phenomena, e.g., Plasma Observatory (Retinò et al., 2022), AME (Dai et al., 2020) and Helioswarm (Klein et al., 2023).
2 The curlometer and basic concepts
2.1 Integral method
The application of the curlometer to Cluster data was reviewed by Dunlop and Eastwood (2008). More recently, its adaption to the context of the high altitude ionosphere, focusing on the determination of field-aligned currents (FAC) was covered by Dunlop et al. (2020) and Trenchi et al. (2020) [see other papers in the ISSI book on ionospheric multi-spacecraft data analysis tools (Dunlop and Lühr, 2020)]. The method has also been reviewed by (Dunlop et al., 2018; Dunlop et al., 2021a; Dunlop et al., 2021b), and was surveyed by Robert and Dunlop (2022), and its application to the Earth’s ring current region has also been recently analysed in the context of MMS data (Tan et al., 2023), who also addressed its accuracy in different regimes. The Cluster Science archive (http://www.cosmos.esa.int/web/csa/software) contains method implementations also in the technical note by Middleton and Masson (2016).
The calculation uses the integral form of Ampère’s law, i.e., neglecting the displacement current () for high electrical conductivity (Russell et al., 2016), where B and E are the magnetic and electric fields and J is the current density. The technique (Dunlop et al., 1988; Robert et al., 1998a) combines four, non-planar spatial positions to make a linear estimate of the electric current density, i.e., , where ΔBi, ΔRi are the differences in the measured magnetic field at positions (i, j) to a reference spacecraft, giving a rugged and simple formalism (see also Section 2.2). The current density normal to each face of the spacecraft tetrahedral configuration is represented by the terms on the left-hand side of the equation. One of the four normal components is redundant and can be used to check stability of the estimate (Dunlop et al., 2018; Dunlop et al., 2020), by choosing different faces in the estimate of J. For irregular spacecraft configurations, this also allows some flexibility to choose which face gives the best estimate of a component, where the relative alignment of the spacecraft configuration to the local field geometry is significant and often only one face determines a stable J component (see also the methodology in Vogt et al., 2009; Shen et al., 2012b).
A partial estimate of one component can still provide useful information if the large-scale current orientation is assumed, such as for FACs and in the case of the in situ ring current, where the azimuth component is significant (Zhang et al., 2011). For example, the three magnetospheric THEMIS spacecraft (Yang et al., 2016) can be used as shown in the left panel of Figure 1 in the ring current, but these assumptions can severely limit the stability of the estimates (Tan et al., 2023), and indeed the need to project the normal component into the ring plane means assumptions on the form of the large-scale currents are critical. The right-hand panel in Figure 1, illustrates that Swarm close configurations can also be used for partial and full current estimates with assumptions on the stationarity of the currents (over a few seconds) and that the field-aligned component is dominant or force-free (Shen et al., 2012a; Ritter et al., 2013; Vogt et al., 2013; Dunlop et al., 2015b, references in Dunlop and Lühr, 2020).
FIGURE 1
The right-hand side of Figure 1 illustrates the adaption in Dunlop et al. (2015b), where adjacent positions are used to form at least four points in space. As with the standard form of the curlometer, the convection time across the array is the relevant temporal scale for the estimates. In the case of Swarm this is typically 5–10 s for separation scales of around 100–150 km for the Swarm A,C pair of spacecraft, which fly side by side in near circular, polar orbits (∼500 km altitude). Swarm B flies at a slightly higher altitude in an orbit but is only in alignment at specific times during the mission. In the special close configuration shown in Figure 1, a series of values can be made from different combinations of the five points (A,B,C,A’,C’), providing information on any temporal changes as well as comparative estimates. For resolving the FAC component ACA’C’ provides a vertical component of J (Ritter et al., 2013). The configuration ABC provides simultaneous measurements, but suffers from the fact that the plane is not well aligned to the FACs.
Although generally robust, the relative structure scales applying affect the validity range of the estimates. For Cluster separations (>100 km), the dominant error arises from nonlinear gradients, while at MMS separations (∼5–10 km) measurement uncertainties can be important [typically these affect the estimate below a threshold |J| (Dunlop et al., 2018), i.e., for the small MMS tetrahedron scales (a few km), measurement uncertainty (∼0.1 nT in B; ∼100 m for R and millisecond timing) drives the error unless the current density is greater than several nAm−2]. Both Runov et al. (2005) and Forsyth et al. (2011), for example, have examined the effect of the characteristic scale of current structures on the use of quality estimates and indeed Tan et al. (2023) have compared MMS to Cluster results, which sample on distinct large and small spatial scales in the ring current. The linear estimator has been used extensively as an indirect quality parameter (Robert et al., 1998a; Haaland et al., 2004b), along with the constellation shape (elongation and planarity) as discussed in more detail in Section 2.2. The linear estimate of the average value of divB over the volume of the tetrahedron is an integral part of the method and is given by ,
This full combination of the four positions are needed for Q to be used which is found unreliable if the spacecraft configuration is highly irregular and not well aligned to the background magnetic field structure. The estimate of Q has also been used in qualification of the FOTE (First-Order Taylor Expansion) method (see Section 2.3).
In the magnetosphere, the effect of dipole non-linear gradients, not associated with current density [first noted in Dunlop et al. (2002), while Grimald et al. (2012) considered this in the context of the ring current], can be minimised by subtracting the dipole (or IGRF) field from the measured magnetic field to give magnetic residuals, e.g., for studies of the in situ ring current (Yang et al., 2016), where dipole gradients are significant and at low altitude orbits (LEO) where the formation of magnetic residuals is normal practice, particularly for Swarm (Ritter et al., 2013; Dunlop et al., 2020).
2.2 The influence of elongation, planarity and Q: the limiting case
At the mesoscales of Cluster separations, the accuracy of the curlometer estimate primarily depends on how the neglected non-linear gradients contribute in the context of the spacecraft constellation, i.e., its scale, shape (irregularity) and relative orientation to the measured current structure. This is also true at smaller separation scales (e.g., for MMS) but then measurement errors (in the magnetic field, position and timing) also become significant. For Cluster, therefore, the spatial sampling through the constellation shape was considered. To characterize the shape using the three eigenvalues of the volumetric tensor , , (in order of magnitude, i.e., ) the square roots a, b, and c are used to define: elongation and planarity . These parameters can be used to check the degree of irregularity of the tetrahedral shape, so that they complement the value of Q, since this because a poor indicator for irregular tretrahedral shape.
In terms of the curlometer approximation, whether the method is used on three magnetospheric spacecraft from the THEMIS mission, or the tetrahedral constellations of MMS and Cluster, current density components are calculated starting from a single plane formed by 3 spacecraft. We therefore need to calculate, to first order, the closed integral (Equation 1) of the magnetic field and divide it by the area of the triangle formed by 3 spacecraft (see also discussion of the integral forms in Section 2.3).
Under the condition of limited observation data, we calculate the integral of the magnetic field, and the area of the triangle formed by three positions, from the following Equations 2, 3 (as written above in condensed notation).
In the limiting case, the three spacecraft are collinear. Suppose that spacecraft 3 is between 1 and 2, with a distance t from 1 and a distance 1-t from 2. in Equation 4 is the non-linearity term, but it could have contributions from the nonlinearity of the magnetic field, measurement errors of magnetometer, differences between the measurement at each spacecraft and time-varying error introduced by time interpolation. Putting Equations 4–6 to Equation 2, we get Equation 7
If is not zero, which is almost certain, then in this limiting case the integral of the magnetic field is not zero, while the area of the triangle formed by the three spacecraft is exactly zero; producing an infinite current error. Thus, when the 3 spacecraft are collinear, the actual current density cannot be obtained. Even if the three are not collinear, the calculated current density will increase dramatically as they approach the collinear position. Because of the nonlinear term, the area approaches zero faster than the curve integral approaches zero. Estimating accurate error is limited by the difficulty of obtaining the exact nonlinear term, however, future work will attempt to quantify this. It should also be noted that the parameter Q cannot be well used in this limit, because becomes infinite but not . Thus, Q will remain a very low value and lose its expected function. In past analysis, therefore, the Elongation for a triangle has been used to gain an empirical reliability and exclude bad results. In other words, only those calculation results meeting a certain Elongation condition (e.g., less than 0.8) can be trusted.
For the case of the full four spacecraft tetrahedral configurations (as for MMS and Cluster), a further step is needed as three current density components are obtained from three of the four planes in a tetrahedron. The current density vector is then obtained by solving equations (as indicated above). Using subscripts 1, 2, and 3 to represent the average current density J and the normal direction N of the three planes, respectively.
By denoting , then Equation 8 is reduced to Equation 9. And Jxyz is obtained (Equation 10).
The calculation process is actually naturally stable because of the closure of the equations for a sampled volume through Ampères law, and in addition the fourth face of the tetrahedron provides a check on the estimates. Nevertheless, if matrix N is not well-conditioned (i.e., ill-conditioned or near singular), we cannot get the true current (as will occur if the constellation is near planar). As mentioned earlier, the matrix N is made up of normal vectors on three planes, and the most likely (or perhaps the only) factor for poor quality is that some normal vectors are closely to be parallel. The effect on Q in this case is not well understood, but further work is expected to clarify this through specific analysis in the future. Currently, we can conservatively choose to trust the results obtained by the tetrahedral configuration with better non-coplanar conditions. This condition has been evaluated in the past through the added use of the Planarity, i.e., typically limiting this also below 0.8 (In practice, both parameters are combined into a quality index as the square route of the sum of the squares of E and P with this value <0.8, typically).
When using the curlometer method, careful attention to the configuration parameters of the constellation must be paid, therefore, due to the existence of these factors (arising from the neglected nonlinearity of the physical quantities). Tailored analysis of specific events can reduce the risk of calculation anomalies and in the case of large amounts of data can be statistically processed. Clearly, sampling by more spacecraft (and more than four spacecraft in particular) can help stabilize the estimate through alternative choice of the planes and by selection of particular tetrahedra within the constellation, but then the tracking of the more complex constellations requires more management. This and the use of the quality indicators, have been considered, for example, for new, proposed constellation mission operations, such as Plasma Observatory (Retinò et al., 2022) and Helioswarm (Klein et al., 2023).
2.3 Advanced integral theorems: geometrical approach
The curlometer method rests on the ability to estimate linear gradients in the measured quantities (specifically the magnetic field) between spacecraft positions. Apart from the linear interpolation method (Chanteur and Harvey, 1998) and the least-squares method (Harvey, 1998), the geometrical method (integral theorems), introduced above is the third way to obtain general estimators of the linear gradients of physical quantities. Recently, Shen and Dunlop (2023) have made full use of a geometrical method (summarized below in Equations 11–26) to derive the gradient, divergence, and curl of physical quantities with the integral theorems. Furthermore, this geometrical method has the special advantage to easily derive the field gradients for observations made by a planar constellation. The errors in the estimators of the linear gradients from the geometrical method were found to enter at second-order and it was illustrated that the method based on integral theorems are equivalent to the spatial interpolation method (Chanteur, 1998; Chanteur and Harvey, 1998; Vogt et al., 2009) and the least-squares method (Harvey, 1998; De Keyser, 2008; Hamrin et al., 2008) for deriving linear gradients.
In Figure 2 the position vectors of the four spacecraft in a tetrahedral configuration are . The barycenter coordinates are chosen as . The spacecraft , , and constitute a face triangle as discussed earlier , where the vertex opposite to this is , as shown in Figure 2. The three vertices , , and of the face are defined to rotate anticlockwise around its normal . In this notation, the vector area of the face of the tetrahedron is
FIGURE 2
The volume of the tetrahedron is
For a certain arbitrary scalar field f, vector field and tensor field , the satellite of the constellation yields the scalar field , vector field and tensor field . The scalar field f, vector field and tensor field obey the following integral theorems (Bittencourt, 2004):
Starting from these integral theorems, Shen and Dunlop (2023) have obtained the estimators of the gradients of a scalar field f, vector field and tensor field , as well as the curl and divergence of the vector field , respectively, as below:
Considering the reciprocal vector as defined by Chanteur (1998)then the above estimators are identical to those from the interpolation method (Chanteur, 1998).
The integral theorems method, however, has one special advantage, it can easily derive the field gradients for measurements from a planar constellation [e.g., as for the three-spacecraft THEMIS (Friis-Christensen et al., 2006) or Swarm (Angelopoulos, 2009) configurations].
By using the following integral theorem applied to the triangle :the averaged gradient of the scalar field in the plane of the constellation is readily derived as the following formula (Shen and Dunlop, 2023)
Similarly, for a vector field , its averaged gradient of the scalar field in the plane of the constellation is
A rigid error analysis has been made for this geometric approach based on Taylor expansion (Shen and Dunlop, 2023). It is verified that the truncation error of the method is at the order of , where L is the characteristic size of the constellation tetrahedron (Robert et al., 1998b) and D is the length scale of the field structure measured. It is found that the truncation error for deriving the linear gradient with the four point measurements by Cluster and MMS is actually very small and Cluster (Escoubet et al., 1997; Escoubet et al., 2001), THEMIS (Angelopoulos, 2009) and MMS (Burch et al., 2016) are generally able to yield stable estimates of current density, charge density, curvature of magnetic field lines, and other related parameters (Shen et al., 2003; Haaland et al., 2021; Pitout and Bogdanova, 2021; Shen et al., 2021b; Robert and Dunlop, 2022).
3 Magnetic gradients and topology
The gradient and curvature terms in the dyadic of the magnetic field, B, can be linearly estimated (Chanteur, 1998; Harvey, 1998; Shen and Dunlop, 2008; Vogt et al., 2008; Shen et al., 2012a; Shen et al., 2012b), from which the current density can be obtained. Key methodology includes magnetic curvature and rotation analysis (Shen et al., 2007; Shen et al., 2012a); least squares analysis of planar reciprocal vectors (De Keyser et al., 2007; Hamrin et al., 2008; Vogt et al., 2009; Vogt et al., 2013), and a range of techniques related to Taylor expansion around the measurement points (FOTE method, see Section 3.3). Estimates based on these gradient methods basically depend on both the integrity of the spacecraft array and stationary properties (temporal dependence) of the magnetic structures, although additional constraints or assumptions can be incorporated. More recently, the polynomial reconstruction of the magnetic field topology has been explored using MMS data by Denton et al. (2020), Denton et al. (2022).
3.1 Magnetic rotation analysis applications
In addition to estimating direct gradients, Magnetic Curvature Analysis (MCA) (Shen et al., 2003) and Magnetic Rotation Analysis (MRA) (Shen et al., 2007) give the 3-D topology of the magnetic field (curvature radius, normal direction and binormal direction of the magnetic field-lines). Particular results have been obtained in the magnetotail current sheet (Shen et al., 2008a; Shen et al., 2008b; Rong et al., 2011); the Earth’s ring current (Shen et al., 2014); flux ropes and plasmoids (Zhang et al., 2013; Yang et al., 2014); reconnection regions (Lavraud et al., 2016; Zhang et al., 2016), and in the cusp and at the magnetopause (Shen et al., 2011; Xiao et al., 2018).
The curvature of field-lines (from MCA) from Cluster located within the ring current are shown on the left-hand side of Figure 3, which indicates current strength in a complimentary manner to the magnetic gradient based J calculation. The plot shows a decrease in relative curvature, which implies a growth of the implied current density which is related to increasing geomagnetic activity (e.g., SYM-H). There is also a dawn-dusk asymmetry, which is most apparent at lower activity levels. The right-hand side of Figure 3 shows that the method can resolve the form of the tail current sheet in terms of different current sheet geometry. As discussed by Rong et al. (2011) these can be classified as: normal, flattened and tilted. These distinct geometries result from the combined MCA/MRA methodology, where the key properties of curvature radius, normal direction and binormal direction are extracted to inform the sheet characteristics controlling current sheet dynamics.
FIGURE 3
3.2 Least-squares methods for multi-point gradient computation
Typically, instruments record multiple data points in the convection time needed to cover a comparable distance to the separation scales, which in principle carry information relevant for calculation of the gradient, as was already noted by Harvey (1998). This is particularly true if a certain degree of time invariance in the structures of interest can be assumed. The Gradient Analysis by Least-Squares (GALS) technique (Hamrin et al., 2008) and the Least-Squares Gradient Calculation (LSGC) technique (De Keyser et al., 2007; De Keyser, 2008) apply these ideas. The latter is based on least squares gradient calculation by approximating the measured quantity (scalar or vector) through Taylor series expansion around the measurement reference point (typically the barycenter). This expansion describes the magnetic field, for example, its spatial and temporal gradients and non-linear terms at higher-order. An assumption that the gradients are constant on certain spatial scales, allows higher-order terms to be estimated. It follows that, in practice, an iterative, weighted least-squares, procedure can be devised.
The method can provide error estimates on the results. As in the case of the curlometer, these reflect errors in both the measurement and non-linear behaviour, (to simplify, uncorrelated measurement errors as well as homogeneity parameters can be assumed for all three components). The property that divB = 0, that the parallel magnetic field gradient is zero, or static structures can be added as constraints. An important application of a higher number of spacecraft is in error control. Although other numbers of measurement points can be used, quality of the results depends on the measurement errors.
It is worth also noting here that the conventional ways to calculate the gradient involve the calculating of the inverse of volume tensor (e.g., Chanteur, 1998; Harvey, 1998; Shen et al., 2003; Shen et al., 2007). The volume tensor would become an ill-conditioned matrix, however, when Cluster tetrahedron becomes an irregular shape, e.g., plane-like or line-like, so that the direct calculation of the inverse of volume tensor would yield significant error and gradient cannot be correctly calculated in this case.
Shen et al. (2012b) avoided the problem of irregular shaped configurations of spacecraft, by introducing a procedure where transforming coordinates into the eigenvector space of volume tensor allows the gradient to be universally calculated. The gradient calculation can involve the inverse of volume tensor (e.g., Chanteur, 1998; Harvey, 1998; Shen et al., 2003; Shen et al., 2007), which is problematic for very irregular shapes. A bonus of the approach of Shen et al. (2012b) is that it can be applied to three-point magnetic field observations (e.g., as in the case of THEMIS), to give current density and the vorticity of plasma flow, and to three-point plasma measurement for the vorticity of K-H waves (as for Cluster).
3.3 FOTE methods and applications: local Taylor expansion
The FOTE method is based on the Jacobian matrix δB, which is a 3 × 3 real matrix . With four-point measurements of magnetic fields, and can be easily obtained. Theoretically, such matrix has three eigenvectors, e1, e2, e3, and correspondingly three eigenvalues, λ1, λ2, λ3. The sum of these three eigenvalues is zero (λ1+λ2+λ3≡∇⋅B = 0), because the magnetic field is “non-divergent.” This implies that either all the eigenvalues are real or one is real while the two others are conjugate complex (Fu et al., 2015; Fu et al., 2020).
The immediate application of the FOTE method is to find magnetic nulls (particularly in regions containing magnetic reconnection X-lines); complementing the analysis based on the use of the Poincare index (Xiao et al., 2006; Xiao et al., 2007) and field line reconstruction methods (He et al., 2008a; He et al., 2008b; Dunlop et al., 2009). Assuming that the magnetic field changes linearly around the spacecraft tetrahedron, the position of a magnetic null can be resolved if we perform the first-order Taylor expansion of magnetic fields around this null, , where is the spacecraft position, is the distance from spacecraft to the null, and is the magnetic field measured by each spacecraft. Applying the four-spacecraft measurements to this equation, r is easily resolved, and therefore, the null position is known (Fu et al., 2015). Notice that the null-spacecraft distance is a three-dimensional vector. Such distance, however, may involve uncertainties if the null is a quasi-2D structure (e.g., X-null). In other words, if the null is 2-D, the derived null-spacecraft distance in the “out-of-plane” direction is unreliable (Fu et al., 2019). In such situation, we only consider the 2-D null-spacecraft distance, i.e., the null-spacecraft distance in the reconnection plane.
Magnetic null types can also be identified. Since the sum of the three eigenvalues is zero (), either all the eigenvalues are real or at least one is real and the other two are conjugate and complex. The different conditions affect the type of null: it is radial when the eigenvalues are real for both A- and B-type. A combination of one positive and two negative eigenvalues gives an A-type null, while two positive and one negative eigenvalue gives a B-type. In the other case, where only one eigenvalue is real, the null is an As- and Bs-type spiral (the As null corresponds to a positive real eigenvalue and the Bs null to a negative real eigenvalue). Sometimes large instrument uncertainties or magnetic field non-linearity mean the type cannot be identified, so is labelled “unknown”. A, B, As, and Bs are all the null types in 3D regime. They are labeled by using the symbols △, ▷, ▲, ▶, respectively (Fu et al., 2015; Fu et al., 2020).
The FOTE method can also determine the dimensionality of a magnetic null. Among the three eigenvalues of the Jacobian matrix, if one eigenvalue is significantly smaller than the two others, the three-dimensional A- and B-null will degenerate into two-dimensional X-null, or in other words, the magnetic topology will have a 2-D appearance; if the real part is significantly smaller than the imaginary part, the three-dimensional As- and Bs-null will degenerate into two-dimensional O-null, which certainly has the 2-D appearance (Fu et al., 2015; Wang et al., 2020). Typically, in spacecraft measurements, we simplify A and B nulls to X null if the three eigenvalues satisfy and simplify As and Bs nulls to O null if the real and imaginary parts of the three eigenvalues satisfy . The A, B, As, and Bs nulls are 3-D structures, while the X and O nulls are 2-D structures. Such 2-D structures are characterized by negligible magnetic fields in the “out-of-plane” direction. In space plasmas, the O null (or O line) is referred to plasmoid or flux rope. Figure 4 shows an example for the application of the FOTE method to find magnetic nulls during 12 magnetic-reconnection events detected by the MMS mission at the Earth’s magnetopause, with the null-spacecraft distance (see the left-side vertical axis), null types (see the symbols), null dimensionality (see the labeling system at the bottom of the figure), and the analysis error (see the gray shade and the right-side vertical axis) exhibited (adapted from Fu et al., 2019).
FIGURE 4
During periods of magnetic reconnection, the open angle of separatrix-lines can be resolved by the method. The Jacobian matrix δB has three eigenvectors, , and correspondingly three eigenvalues, . The open angle of two separatrix-lines is determined by the two eigenvectors related to the two large eigenvalues. The angle between these two eigenvectors is the open angle of the two separatrix-lines (Chen et al., 2018). For example, if , the open angle is the angle between and . Such an open angle directly determines the reconnection rate of a reconnection process (Chen et al., 2019; Wang et al., 2020). Figure 5 is an example, showing the application of the FOTE method to resolve the open angle of separatrix-lines and deduce the reconnection rate during an unsteady reconnection at the Earth’s magnetopause (adapted from Wang et al., 2020).
FIGURE 5
The FOTE method can be used in the reconstruction of magnetic topology around magnetic nulls. In eigenvector coordinates , we trace and inverse-trace a few points around the null to obtain the magnetic field topology. The step length of trace/inverse-trace is typically set to be the local magnetic strength (Fu et al., 2016). Figure 6 illustrate the application of the FOTE method to reconstruct the topology of a radial-type magnetic null and a spiral-type magnetic null, which is also referred to the magnetic flux rope in spacecraft observations (modified from Fu et al., 2017; Wang Z. et al., 2019).
FIGURE 6
Finally, in terms of the errors during applications of the method to real data, we require the null-spacecraft distance to be less than the local ion inertial length, in order to guarantee that the null positions are accurately resolved. In addition, to guarantee that the null properties are accurately identified, we define two parameters ( and ) and require them to be smaller than 0.4. These criteria are derived from the comprehensive test of three-dimensional simulation data (Fu et al., 2015; Fu et al., 2016).
4 Application of the curlometer to currents in the magnetosphere
Due to its robust and flexible nature, the curlometer calculation is perhaps the most widely used in the magnetosphere (notably applied in: the magnetopause boundary layer (Dunlop et al., 2002; Haaland et al., 2004a, e.g., Dunlop and Balogh, 2005; Panov et al., 2006; Panov et al., 2008; Shi et al., 2019); the magnetotail (Runov et al., 2006, e.g., Nakamura et al., 2008; Narita et al., 2013); the ring current (Vallat et al., 2005; Zhang et al., 2011, e.g., Shen et al., 2014; Yang et al., 2016); field-aligned currents (e.g., Forsyth et al., 2008; Marchaudon et al., 2009; Shi et al., 2010; Shi et al., 2012; Dunlop and Lühr, 2020) and other transient signatures and in the solar wind (e.g., Eastwood et al., 2002; Xiao et al., 2004; Shen et al., 2008a; Roux et al., 2015). Some of these applications are briefly reviewed here.
4.1 Basic use of the curlometer and time stationarity
The magnetopause boundary layer (MPBL) matched well the scale size of the early Cluster mission phases (100–2,000 km spacecraft separation). Figure 7 shows examples from Haaland et al. (2004a), Dunlop and Balogh (2005). The left panel shows that J orientations are predominantly in the MP plane during a number of in/out crossings resulting from magnetopause motion (average speed ∼25 km/s; with average thicknesses ∼1,200 km), while the right panel shows a thin MPBL with high current density. Typically, the Earth’s magnetopause thickness varies from 100s of km (a few ion gyro radii) to 1,000s of km (Berchem and Russell, 1982; Paschmann et al., 2005; Panov et al., 2008), while corresponding current densities vary from 10 to 200 nA/m2. In Figure 7A, the signatures outside the main MP crossing period are magnetosheath FTEs, where the current is along the mean reconnected flux tube direction (also studied by Pu et al., 2005).
FIGURE 7
Table 1 gives a summary of typical current density values in the Earth’s environment, based on estimates of large-scale morphology and transient structure.
TABLE 1
| Feature/Region | Typical values for J |
|---|---|
| Magnetopause currents | ∼10 nA m−2 (Dunlop and Eastwood, 2008), to 100 s nA m−2 (see, e.g., Panov et al., 2008) |
| Currents in flux transfer events | ∼1 nA m−2 (Dunlop and Eastwood, 2008) up to 10 nA m−2 (Pu et al., 2005) |
| Current at the cusp boundaries | ∼20 nA m−2 (Dunlop et al., 2002) |
| Field aligned currents (FAC) | ∼2 µA m−2 at 500 km altitude and ∼20 nA m−2 at 2.5 RE altitude (Dunlop et al., 2005) |
| Magnetotail current sheet | up to ∼30 nA m−2 (Runov et al., 2006) |
| Plasma sheet boundary layer | ∼10 nA m−2, (Nakamura et al., 2004) |
| Ring current | 9–27 nA m−2 at 4–4.5 RE, (Zhang et al., 2011) |
| Solar wind current sheet | ∼10 nA m−2 (Eastwood et al., 2002) |
Typical current density values (from Dunlop et al., 2021a).
Figure 7 also illustrates that the combination of the curlometer and discontinuity analysis (which can obtain boundary orientation and motion (see Dunlop and Woodward, 1998; Dunlop et al., 2002; Haaland et al., 2004a) can confirm the thickness of the current layer and the alignment of J in the local MP plane. Broad scaling of |J| (10–50 nAm−2) depending on a range of thicknesses, ΔD, can be shown to be consistent with the effective planar current . The Cluster results tend to underestimate current for higher J and thicker boundary layers (compared to the separation). Indeed, the existence of small-scale sub-layers within the MPBL, having high intensity currents, were not often resolved by Cluster, but were seen by MMS (Dunlop et al., 2021a). The MVAJ method, referred to in the right panel of Figure 7 [see also Xiao et al. (2004), who apply the method to FTE orientations] better ties the orientation of the current sheet to J (minimum variance of J obtains the orientation of a near 1D current sheet, since div J = 0, when μ0J = curl B). The velocity of the current sheet can also be obtained (Haaland et al., 2004b), where different estimates of orientation all agree to within a few %.
A second key region demonstrating the capabilities of the curlometer became accessible after the launch of the multi-point measurements from Swarm at low-Earth orbit (LEO) altitudes. Although Swarm is a three-spacecraft mission and is not always in a close constellation, the method can be generalised by using nearby positions in time as indicated in Figure 1. This provides estimates even for only 2 or 3 closely separated spacecraft when either the dominant currents are field-aligned (FAC) or characteristic currents are locally static, [typical in this region of the high altitude ionosphere and thermosphere, e.g., Ritter and Lühr, 2006; Ritter et al., 2013; Dunlop et al., 2015a; Dunlop et al., 2015b; Dong et al., 2023, and references in Dunlop and Lühr (2020)]. Subtraction of the main, background field components prior to application of the method is essential At these low LEO altitudes (400–600 km), the main, background field must be subtracted to avoid the effect of zero current non-linear gradients [typically the IGRF or Chaos model (e.g., Olsen et al., 2014) are subtracted].
Typical convection times of ∼10–15s apply in the case of the Swarm spacecraft separations of ∼100–150 km at mid to high latitudes. Thus, the multi-spacecraft estimate depends on the FACs not being highly time dependent (e.g., ULF or Alfvén waves); but usually this is only relevant for small-scale currents, which can be easily identified by differences in field measured at the individual spacecraft. Thus, in addition to identifying smaller scale and time dependent structures, the extended methodology maps out the morphology and dynamics of larger scale current systems, e.g., region 1, region 2 and NBZ FAC systems (see review in Mcpherron et al., 1973; Shiokawa et al., 1998; Cao et al., 2010).
Figure 1B shows the Swarm configuration, as considered in Dunlop et al. (2015a) for a close conjunction of Cluster and Swarm (as shown in Figure 8). A series of 2, 3 and 4 spacecraft estimates for the FACs can be made (in the context of the time shifted positions), while 4-spacecraft positions also give very low Q estimates (a few %). Changing the choice of the selected spacecraft positions can indicate how stable the estimates are (note that the effective barycentres of each set are slightly different). These comparative estimates are shown on the right-hand side of Figure 8, which also shows (top two panels) the filtered single spacecraft FAC estimates from dB/dt (Lühr et al., 2015). The 4-point estimates of the FAC profiles match those estimated from Cluster measurements most closely (with suitable scaling). Cluster moves to higher MLAT values and then back to lower MLAT in its traversal of the region but cross the local time of the Swarm orbit within minutes of the Swarm pass. The results imply that ∼ 1–200 km Swarm currents (at LEO altitudes) can coherently map to Cluster.
FIGURE 8
The 4-point calculation can also identify any perpendicular components and estimates for this event appear to show possible, associated hall signatures consistent with wire model FACs (see Gjerloev and Hoffman, 2002; Ritter et al., 2004; Wang et al., 2006; Liang and Liu, 2007; Shore et al., 2013). Related methods dealing with FAC estimates have been cross-compared to assess key events by Trenchi et al. (2020) and a similar time-shifted Swarm configuration, tailored to the low latitude regions, has attempted to extract low latitude currents (Fillion et al., 2021).
A third key region, suitable for adaption of the curlometer is the Earth’s ring current (RC). This was studied early with Cluster (Dandouras et al., 2018; Dunlop et al., 2018), and then using the 3-spacecraft magnetospheric THEMIS spacecraft, as mentioned in Section 2.1. More recently MMS measurements have been used to estimate the ring current densities. There are contrasts between each: Cluster only accessed the mesoscale separations; the 3 THEMIS spacecraft limit the estimate to the normal current component, and MMS only has survey magnetic field data (with no thermal plasma) in the RC region. As first reported by Vallat et al. (2005), the Cluster polar orbit cuts north to south through the ring plane, allowing all local times to be scanned over the mission (Zhang et al., 2011) for a narrow range of radial distance (∼4–4.5 RE). Typically, the azimuthal component, Jϕ, can be made, where the IGRF model is subtracted to reduce the effect of non-linear dipole gradients (mentioned in Section 2.1). Full azimuthal scans can highlight local time asymmetry in the RC and a possible relation to asymmetries at the MP has been suggested by Haaland and Gjerloev (2013).
As shown on the left-hand side of Figure 1, the three magnetospheric THEMIS spacecraft I near equatorial orbits also provide a wide MLT coverage (Yang et al., 2016), but for a range of radial distances, as shown in the left panel of Figure 9, so that L-shell profile can be revealed together with the Westward-Eastward current boundary on the inner edge of the RC, as shown in the right panel of Figure 9.
FIGURE 9
MMS also samples the ring current in a near equatorial orbit, providing a comparative RC estimate covering the same region as THEMIS on small separation scales. Thus, MMS can better resolve the trends seen in both the radial and azimuthal morphology, and also can identify small scale intense currents, which are not resolved by Cluster or THEMIS. The recent studies using MMS data (Tan et al., 2023), mentioned in Section 2.1 show both large and small-scale structures can be present, extending the earlier studies with Cluster and THEMIS. The morphology of the RC (left hand side of Figure 10) is broadly consistent with previous in situ studies with strong dawn/dusk and noon midnight asymmetry, but shows more detailed structure, with Tan et al. (2023) reporting a layered structure in latitude. A partial RC (or banana current), with an inner eastward current (blue) is most clear in the noon to dusk quadrant. There is no evidence of enhancement on the dusk-side during geomagnetically active periods.
FIGURE 10
Comparison of the Swarm low-altitude, dual-satellite FAC data with mapped MMS FAC signatures measured adjacent to the RC shows that RC behaviour and R2 FACs can be investigated directly (see right hand side of Figure 10) and show consistent patterns. The statistical coverage has some overlap in the location of Swarm FACs and the mapped locations of parallel currents measured adjacent to the RC (between 60 and 70 deg). In the auroral zone (particularly inside 65 deg, shown as a dashed circle), there is a qualitative R1/R2 pattern.
4.2 Recent application to MMS configuration scales
The current density can be obtained from velocity and density moments of the 3D plasma distributions for all ion species and electrons. Cluster, for example, was limited to low cadence (since the full distribution is taken over a spacecraft spin period) and typically particle distributions are measured over a limited energy range [the use of particle moments previously can be found in Henderson et al. (2008), Petrukovich et al. (2015), and indeed (Yao et al., 2014), who also used pressure gradients to estimate the perpendicular current density]. The MMS mission can measure plasma distributions at high-time resolution (primarily in the outer magnetosphere). The MMS constellation is generally much smaller size than the prevailing ion structure, while the Cluster configuration was usually of order or greater than the MPBL, for example. In many regions of the magnetosphere there are intense small-scale currents which are missed on the separation scales of Cluster (100 s km) so typically the curlometer tends to underestimate the current density. On MMS separation scales, however, the ion-scale structures can be well resolved. A number of studies have used MMS plasma moments to estimate the currents (e.g., Lavraud et al., 2016; Phan et al., 2016) and comparing these estimates to the curlometer can reveal details of sub-structure; smaller scale current layers, and details of the current carriers. Moreover, Cluster and MMS conjunctions (for example, at different locations on the MP) provide the opportunity to compare across multiple spatial scales. Below we briefly show two aspects of this analysis [a more complete treatment d given in Dunlop et al. (2021a)].
4.2.1 Sub-structure in the MP and FTEs
The left-hand side of Figure 11 shows a typical MMS magnetopause crossing where the curlometer current densities are plotted with the current density estimated from the plasma moments [(], which can be estimated at each spacecraft position. The plasma currents closely follow the curlometer when averaged over the four spacecraft positions. The normal component (panel d) shows the most significant differences between the plasma currents and the curlometer. Panel (e) shows that for currents below ∼50 nAm−2 both the measurement errors and Q are significant. Despite the close curlometer agreement with the spatially averaged plasma current (i.e., the mean current over the tetrahedron), the estimates at each spacecraft (bottom panel) vary significantly between each position. This might suggest small scale (filamentary) structure within the magnetopause layer, where the dominant current carriers are measured by the plasma moments (this would be consistent with the intense, narrow bursts of current seen in the curlometer profiles), and appears to be typical of the magnetopause layer (Dong et al., 2017). This substructure was not well resolved by Cluster (except for the smallest separation scales) and the Cluster array tends to miss these filamentary currents.
FIGURE 11
Many recent MMS studies have also focused on the ion-scale structure of FTEs (e.g., Eastwood et al., 2016; Zhao et al., 2016; Dong et al., 2017; Teh et al., 2017; Alm et al., 2018; Hwang et al., 2018; Dong et al., 2020). As first reported by Eastwood et al. (2016), agreement between the curlometer and plasma moments can be shown for current density in ion-scale FTEs, where filamentary currents were found and the main current carriers were electrons. The right-hand side of Figure 11 shows two ion-scale FTEs and their detailed current structures (Dong et al., 2017). The currents are highly inhomogeneous and mainly located in either the centre of the flux rope or on the leading edge [this was also true in the Roux et al. (2015) study of a large-scale FTE]. Central bifurcated features in the flux ropes were also present, while as a result of force-free structures, the parallel currents dominate.
4.2.2 Carriers and sources of magnetopause current
Figure 12 shows a thin MP current layer (∼100 km) event encountered by MMS, which was studied by Dong et al. (2018). There is a comparable perpendicular current from the contributions of the ions and electrons in the boundary layer, but the parallel current appears to be dominated by the electron carriers. For the perpendicular currents, the diamagnetic current term () and the directly measured shows good agreement (red line in Figures 12B, C) while the curvature current can be neglected (green line in Figures 12B, C). When the diamagnetic current is separated into ion and electron components, the perpendicular current is dominated by the ion diamagnetic current (: 85%, : 15%, Figures 12D, E). The ions and electrons ultimately carry comparable current through the redistribution of the electric field. The electron current deviates at a narrow front layer in region 2 which suggests non-MHD behaviour (beyond Chapman-Ferraro).
FIGURE 12
Note that the matched plasma and magnetic field measurements from MMS mean that the flow vorticity can be obtained from the curlometer estimates (Dunlop et al., 2002), providing linear estimates for velocity, in addition. Electron vorticities related with electron currents have been studied in the turbulent magnetosheath (Phan et al., 2016; Chasapis et al., 2018), and coherent Alfvén vortices in the same region (Wang T. et al., 2019) can be inferred from alignment of parallel current density and ion vorticity.
4.2.3 Simultaneous MMS and cluster magnetopause crossings
A simultaneous MMS and Cluster crossing of the magnetopause (MPBL) at two different locations is shown in Figure 13, where MMS is located near the subsolar region and Cluster is located on the high latitude post-noon region (Figures 13N, O) and MMS crosses into the magnetosphere ∼20 s later than Cluster. The MMS spacecraft are ∼8 km apart, which is much smaller than the estimated current layer of ∼440 km, while the Cluster spacecraft are ∼3,500 km apart, which is large compared to the current layer at their location of ∼1,400 km. Thus, the small-scale current structure can be resolved by MMS, while the overall MP current is underestimated by Cluster. In fact, the average MMS current (shaded region, Figure 13E) is ∼220 nA/m−2; one order of magnitude larger than estimated at Cluster (Figure 13M), although the simple 1-D (Chapman-Ferraro) current sheet estimated from the change in B and the thickness only differs by a factor of 2-3.
FIGURE 13
The parallel current measured by MMS () is dominant in the magnetosheath boundary layer, where magnetic field component BL contra-rotates to a negative value (Figures 13A, F). This BL feature is also observed by Cluster (Figure 13I), suggesting a similar current structure at both MMS and Cluster locations. The main MMS current sheet has a bifurcated structure, divided by a small current region, noted by the vertical dashed line. This small region corresponds to a flat BL structure (Figure 13A), which is similar to that seen at Cluster 1, 2 and 3 (vertical dashed line, Figure 13I). This suggests the bifurcated current structure extends to the region of Cluster, but is not observed by Cluster 4 ∼15 s beforehand, implying it is also highly dynamic.
Sampling the MPBL at two locations during specific IMF conditions has shown similar overall form for the magnetic field, suggesting that the MP can have similar current structure across a wide region. In another simultaneous event, however, studied by Escoubet et al. (2020) during a high-speed jet in the magnetosheath, the magnetopause current structures at Cluster and MMS are totally different. Similar current structures across a wide region may therefore only remain during relatively stable solar wind conditions.
5 Conclusion
The multi-spacecraft estimates of current density and spatial gradients from Cluster, and the later adaptions to the circumstances of the Swarm and MMS missions, have provided key information on large and small magnetospheric current systems and related transient structures, resolving 3-D currents for a range of conditions in widely different geospace regions. The curlometer in particular, has proved to be reliable and robust. The applicability of the method is limited by certain constraints, particularly those for relatively small structures compared to the spacecraft separation distances. These constraints depend on the form of the spacecraft configuration and the presence of magnetic contributions from other (non-current) sources. Particle moments, give direct, complementary current density estimates, which can be also used to add constraints in generalised methods computing magnetic gradients beyond linear order (Shen et al., 2021a; Shen et al., 2021c).
New constellation missions can make use of further adaptations of this multi-spacecraft methodology within the framework of extended arrays of more than four satellites, particularly when distributed across distinct spatial scales, where the comparative results with MMS, Cluster and THEMIS, as illustrated here, can be re-visited. Figure 14 illustrates two early examples, where close conjunctions of Cluster, Doublestar and THEMIS were interpreted in the context of both local magnetopause structure and the larger scale surrounding phenomena, which provided evidence of the operation of multiple X-lines (see also the discussion in Section 4.2). The adaption of the various quality indicators have been considered, for example, to the new, proposed constellation missions, such as Plasma Observatory (Retinò et al., 2022) and Helioswarm (Klein et al., 2023). As discussed here, there is a balance between the use of either the Q divergence condition or the configuration shape (E and P), depending on the overall comparative spatial scale of the differential measurements. Where fast plasma measurements are available to produce high time resolution moments (n, V, T), then this information can help better resolve the spatial gradients on multiple spatial scales. Finally, we note here that the low frequency spectral response in the measurements (reflecting the time stationarity over the spacecraft arrays) will be an important consideration for applications to configurations of spacecraft which sample distinct spatial scales.
FIGURE 14
Plain text
This article is an account of analysis methods relating to spatial gradients obtained from the differential measurements taken by formations of multiple spacecraft. It focusses on the so named curlometer method, which uses four positions in space to estimate electric current density from the magnetic field. Originally, these measurements were from the set of four ESA Cluster II spacecraft. The curlometer produced many results throughout the magnetosphere, while more recently it has been applied to data from other missions of multiple spacecraft (the NASA MMS and THEMIS missions and the ESA Swarm Earth explorer mission). Although there are certain caviats on its application, the technique has proved to be rugged and stable. Other techniques, such as those estimating magnetic field geometry and gradients, have incorporated the curlometer and the application if these is also covered.
• Robust multi-spacecraft analysis method relating to spatial gradients.
• Estimates of the vector, in situ, electric current density directly.
• Wide application throughout the magnetosphere from magnetic field measurements.
Statements
Author contributions
MD: Conceptualization, Data curation, Funding acquisition, Methodology, Writing–original draft, Writing–review and editing. H-SF: Conceptualization, Investigation, Methodology, Writing–original draft. CS: Investigation, Methodology, Writing–original draft. XT: Data curation, Investigation, Methodology, Software, Writing–review and editing. X-CD: Conceptualization, Investigation, Methodology, Writing–review and editing. Y-YY: Conceptualization, Methodology, Writing–review and editing. PR: Conceptualization, Methodology, Writing–review and editing. PE: Data curation, Validation, Writing–review and editing.
Funding
The author(s) declare that financial support was received for the research, authorship, and/or publication of this article. This work is supported by the NSFC grants 42350710793, 41821003, 42174208; by NERC grants NE/P016863/1 and NE/W003309/1 and by STFC in-house research grant ST/M001083/1.
Acknowledgments
The Swarm satellite data used in this study are available from ESA at http://swarm-diss.eo.esa.int/. For MMS data visit https://lasp.colorado.edu/mms/sdc/public/. Cluster data can be found at http://www.cosmos.esa.int/web/csa/.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
1
AlmL.FarrugiaC. J.PaulsonK. W.ArgallM. R.TorbertR. B.BurchJ. L.et al (2018). Differing properties of two ion-scale magnetopause flux ropes. J. Geophys. Res. Space Phys.123, 114–131. 10.1002/2017JA024525
2
AngelopoulosV. (2008). The THEMIS mission. Space Sci. Rev.141, 5–34. 10.1007/s11214-008-9336-1
3
AngelopoulosV. (2009). “The THEMIS mission,” in The THEMIS mission. Editors BurchJ. L.AngelopoulosV. (New York, NY: Springer), 5–34. 10.1007/978-0-387-89820-9_2
4
BerchemJ.RussellC. T. (1982). The thickness of the magnetopause current layer: ISEE 1 and 2 observations. J. Geophys. Res. Space Phys.87, 2108–2114. 10.1029/JA087iA04p02108
5
BittencourtJ. A. (2004). Fundamentals of plasma physics. New York, NY: Springer. 10.1007/978-1-4757-4030-1
6
BurchJ. L.MooreT. E.TorbertR. B.GilesB. L. (2016). Magnetospheric multiscale overview and science objectives. Space Sci. Rev.199, 5–21. 10.1007/s11214-015-0164-9
7
CaoJ.-B.YanC.DunlopM.RemeH.DandourasI.ZhangT.et al (2010). Geomagnetic signatures of current wedge produced by fast flows in a plasma sheet. J. Geophys. Res. Space Phys.115. 10.1029/2009JA014891
8
CarterJ. A.DunlopM.ForsythC.OksavikK.DonovonE.KavanaghA.et al (2024). Ground-based and additional science support for SMILE. Earth Planet. Phys.8, 275–298. 10.26464/epp2023055
9
ChanteurG. (1998). “Spatial interpolation for four spacecraft: theory,” in Analysis methods for multi-spacecraft data. Editors PaschmannG.DalyP. W. (Noordwijk, Netherlands: ESA Publications Division), 349–370.
10
ChanteurG.HarveyC. C. (1998). “Spatial interpolation for four spacecraft: application to magnetic gradients,” in Analysis methods for multi-spacecraft data. Editors PaschmannG.DalyP. W. (Noordwijk, Netherlands: ESA Publications Division), 371–394.
11
ChasapisA.YangY.MatthaeusW. H.ParasharT. N.HaggertyC. C.BurchJ. L.et al (2018). Energy conversion and collisionless plasma dissipation channels in the turbulent magnetosheath observed by the magnetospheric multiscale mission. Astrophysical J.862, 32. 10.3847/1538-4357/aac775
12
ChenX. H.FuH. S.LiuC. M.CaoD.WangZ.DunlopM. W.et al (2018). Magnetic nulls in the reconnection driven by turbulence. Astrophysical J.852, 17. 10.3847/1538-4357/aa9991
13
ChenZ. Z.FuH. S.WangZ.LiuC. M.XuY. (2019). Evidence of magnetic nulls in the reconnection at bow shock. Geophys. Res. Lett.46, 10209–10218. 10.1029/2019GL084360
14
DaiL.WangC.CaiZ.GonzalezW.HesseM.EscoubetP.et al (2020). AME: a cross-scale constellation of CubeSats to explore magnetic reconnection in the solar-terrestrial relation. Front. Phys.8, 89. 10.3389/fphy.2020.00089
15
DandourasI.Rochel-GrimaldS.VallatC.DunlopM. W. (2018). “Terrestrial ring current: a review of cluster results based on the curlometer technique,” in Electric currents in geospace and beyond. American geophysical union. Editors KeilingA.MarghituO.WheatlandM. (John Wiley and Sons), 115–126. 10.1002/9781119324522.ch7
16
De KeyserJ. (2008). Least-squares multi-spacecraft gradient calculation with automatic error estimation. Ann. Geophys.26, 3295–3316. 10.5194/angeo-26-3295-2008
17
De KeyserJ.DarrouzetF.DunlopM. W.DécréauP. M. E. (2007). Least-squares gradient calculation from multi-point observations of scalar and vector fields: methodology and applications with Cluster in the plasmasphere. Ann. Geophys.25, 971–987. 10.5194/angeo-25-971-2007
18
DentonR. E.LiuY.-H.HasegawaH.TorbertR. B.LiW.FuselierS.et al (2022). Polynomial reconstruction of the magnetic field observed by multiple spacecraft with integrated velocity determination. J. Geophys. Res.127. 10.1029/2022JA030512
19
DentonR. E.TorbertR. B.HasegawaH.DorsI.GenestretiK. J.ArgallM. R.et al (2020). Polynomial reconstruction of the reconnection magnetic field observed by multiple spacecraft. J. Gephys. Res.125. 10.1029/2019JA027481
20
DongX. C.DunlopM. W.TrattnerK. J.PhanT. D.FuH. S.CaoJ. B.et al (2017). Structure and evolution of flux transfer events near dayside magnetic reconnection dissipation region: MMS observations. Geophys. Res. Lett.44, 5951–5959. 10.1002/2017gl073411
21
DongX. C.DunlopM. W.WangT. Y.CaoJ. B.TrattnerK. J.BamfordR.et al (2018). Carriers and sources of magnetopause current: MMS case study. J. Geophys. Res. Space Phys.123, 5464–5475. 10.1029/2018ja025292
22
DongX. C.DunlopM. W.WangT. Y.TrattnerK. J.RussellC. T.GilesB. (2020). MMS observation of secondary magnetic reconnection beside ion‐scale flux rope at the magnetopause. Geophys. Res. Lett.47. 10.1029/2020gl089075
23
DongX. C.DunlopM. W.XiaoC.WeiD.WangT. Y.ZhaoJ. S. (2023). Simultaneous mesoscale polar cusp field-aligned currents measured on mid- and low-altitude satellites. Geophys. Res. Lett.50, e2022GL102460. 10.1029/2022GL102460
24
DunlopM. W.BaloghA. (2005). Magnetopause current as seen by Cluster. Ann. Geophys.23, 901–907. 10.5194/angeo-23-901-2005
25
DunlopM. W.BaloghA.GlassmeierK. H.RobertP. (2002). Four‐point Cluster application of magnetic field analysis tools: the Curlometer. J. Geophys. Res. Space Phys.107. 10.1029/2001ja005088
26
DunlopM. W.DongX. C.WangT. Y.EastwoodJ. P.RobertP.HaalandS.et al (2021a). Curlometer technique and applications. J. Geophys. Res. Space Phys.126. 10.1029/2021ja029538
27
DunlopM. W.EastwoodJ. P. (2008). “The curlometer and other gradient based methods,” in Multi-spacecraft analysis methods revisited. Editors PaschmannG.DalyP. W. (Noordwijk, Netherlands: ESA Communications Keplerlaan), 17–26.
28
DunlopM. W.HaalandS.DongX. C.MiddletonH. R.EscoubetC. P.YangY. Y.et al (2018). “Multipoint analysis of electric currents in geospace using the curlometer technique,” in Electric currents in geospace and beyond. American geophysical union. Editors KeilingA.MarghituO.WheatlandM. (John Wiley and Sons), 67–80. 10.1002/9781119324522.ch4
29
DunlopM. W.LührH. (2020). Ionospheric multi-spacecraft analysis tools. Springer. 10.1007/978-3-030-26732-2
30
DunlopM. W.SouthwoodD. J.GlassmeierK. H.NeubauerF. M. (1988). Analysis of multipoint magnetometer data. Adv. Space Res.8, 273–277. 10.1016/0273-1177(88)90141-X
31
DunlopM. W.TaylorM. G. G. T.DaviesJ. A.OwenC. J.PuZ.LaaksoH.et al (2005). Coordinated cluster/double star observations of dayside reconnection signatures. Ann. Geophys.23 (8), 2867–2875.
32
DunlopM. W.WangT.DongX.HaarlandS.ShiQ.FuH.et al (2021b). “Multispacecraft measurements in the magnetosphere,” in Magnetospheres in the solar system. American geophysical union. Editors MaggioloR.AndréN.HasegawaH.WellingD. T. (John Wiley and Sons), 637–656. 10.1002/9781119815624.ch40
33
DunlopM. W.WoodwardT. I. (1998). “Multi-spacecraft discontinuity analysis: orientation and motion,” in Analysis methods for multi-spacecraft data. Editors PaschmannG.DalyP. W. (Noordwijk, Netherlands: ESA Publications Division), 271–306.
34
DunlopM. W.YangJ. Y.YangY. Y.LührH.CaoJ. B. (2020). “Multi-spacecraft current estimates at swarm,” in Ionospheric multi-spacecraft analysis tools. Editors DunlopM. W.LührH. (Springer), 83. 10.1007/978-3-030-26732-2_5
35
DunlopM. W.YangJ. Y.YangY. Y.XiongC.LührH.BogdanovaY. V.et al (2015a). Simultaneous field‐aligned currents at Swarm and Cluster satellites. Geophys. Res. Lett.42, 3683–3691. 10.1002/2015gl063738
36
DunlopM. W.YangY. Y.YangJ. Y.LührH.ShenC.OlsenN.et al (2015b). Multispacecraft current estimates at swarm. J. Geophys. Res. Space Phys.120, 8307–8316. 10.1002/2015ja021707
37
DunlopM. W.ZhangQ.-H.BogdanovaY. V.LockwoodM.PuZ.HasegawaH.et al (2011). Extended magnetic reconnection across the dayside magnetopause. Phys. Rev. Lett.107, 025004. 10.1103/PhysRevLett.107.025004
38
DunlopM. W.ZhangQ. H.XiaoC. J.HeJ. S.PuZ.FearR. C.et al (2009). Reconnection at high latitudes: antiparallel merging. Phys. Rev. Lett.102, 075005. 10.1103/PhysRevLett.102.075005
39
EastwoodJ. P.BaloghA.DunlopM. W.SmithC. W. (2002). Cluster observations of the heliospheric current sheet and an associated magnetic flux rope and comparisons with ACE. J. Geophys. Res. Space Phys.107 (SSH 9-1), 9. 10.1029/2001JA009158
40
EastwoodJ. P.PhanT. D.CassakP. A.GershmanD. J.HaggertyC.MalakitK.et al (2016). Ion-scale secondary flux ropes generated by magnetopause reconnection as resolved by MMS. Geophys. Res. Lett.43, 4716–4724. 10.1002/2016GL068747
41
EscoubetC. P.FehringerM.GoldsteinM. (2001). IntroductionThe cluster mission. Ann. Geophys.19, 1197–1200. 10.5194/angeo-19-1197-2001
42
EscoubetC. P.HwangK. J.Toledo-RedondoS.TurcL.HaalandS. E.AunaiN.et al (2020). Cluster and MMS simultaneous observations of magnetosheath high speed jets and their impact on the magnetopause. Front. Astronomy Space Sci.6, 78. 10.3389/fspas.2019.00078
43
EscoubetC. P.MassonA.LaaksoH.GoldsteinM. L.DimbylowT.BogdanovaY. V.et al (2021). Cluster after 20 years of operations: science highlights and technical challenges. J. Geophys. Res. Space Phys.126, e29474. 10.1029/2021JA029474
44
EscoubetC. P.SchmidtR.GoldsteinM. L. (1997). Cluster - science and mission overview. Space Sci. Rev.79, 11–32. 10.1023/A:1004923124586
45
FillionM.HulotG.AlkenP.ChulliatA.VigneronP. (2021). Multispacecraft current density estimates in the low- and mid-latitude F-region ionosphere using the swarm constellation. J. Geophys. Res. Space Phys.126, e28872. 10.1029/2020JA028872
46
ForsythC.LesterM.CowleyS. W. H.DandourasI.FazakerleyA. N.FearR. C.et al (2008). Observed tail current systems associated with bursty bulk flows and auroral streamers during a period of multiple substorms. Ann. Geophys.26, 167–184. 10.5194/angeo-26-167-2008
47
ForsythC.LesterM.FazakerleyA. N.OwenC. J.WalshA. P. (2011). On the effect of line current width and relative position on the multi-spacecraft curlometer technique. Planet. Space Sci.59, 598–605. 10.1016/j.pss.2009.12.007
48
Friis-ChristensenE.LührH.HulotG. (2006). Swarm: a constellation to study the Earth's magnetic field. Earth, Planets Space58, 351–358. 10.1186/BF03351933
49
Friis-ChristensenE.LührH.KnudsenD.HaagmansR. (2008). Swarm – an Earth observation mission investigating geospace. Adv. Space Res.41, 210–216. 10.1016/j.asr.2006.10.008
50
FuH. S.CaoJ. B.CaoD.WangZ.VaivadsA.KhotyaintsevY. V.et al (2019). Evidence of magnetic nulls in electron diffusion region. Geophys. Res. Lett.46, 48–54. 10.1029/2018GL080449
51
FuH. S.CaoJ. B.VaivadsA.KhotyaintsevY. V.AndreM.DunlopM.et al (2016). Identifying magnetic reconnection events using the FOTE method. J. Geophys. Res. Space Phys.121, 1263–1272. 10.1002/2015JA021701
52
FuH. S.VaivadsA.KhotyaintsevY. V.AndréM.CaoJ. B.OlshevskyV.et al (2017). Intermittent energy dissipation by turbulent reconnection. Geophys. Res. Lett.44, 37–43. 10.1002/2016GL071787
53
FuH. S.VaivadsA.KhotyaintsevY. V.OlshevskyV.AndréM.CaoJ. B.et al (2015). How to find magnetic nulls and reconstruct field topology with MMS data?J. Geophys. Res. Space Phys.120, 3758–3782. 10.1002/2015ja021082
54
FuH. S.WangZ.ZongQ.ChenX. H.HeJ. S.VaivadsA.et al (2020). “Methods for finding magnetic nulls and reconstructing field topology,” in Dayside magnetosphere interactions. American geophysical union. Editors ZongQ.EscoubetP.SibeckD.LeG.ZhangH. (John Wiley and Sons), 153–172. 10.1002/9781119509592.ch9
55
GjerloevJ. W.HoffmanR. A. (2002). Currents in auroral substorms. J. Geophys. Res. Space Phys.107 (SMP 5-1), 5–13. 10.1029/2001JA000194
56
GrimaldS.DandourasI.RobertP.LucekE. (2012). Study of the applicability of the curlometer technique with the four Cluster spacecraft in regions close to Earth. Ann. Geophys.30, 597–611. 10.5194/angeo-30-597-2012
57
HaalandS.GjerloevJ. (2013). On the relation between asymmetries in the ring current and magnetopause current. J. Geophys. Res. Space Phys.118, 7593–7604. 10.1002/2013ja019345
58
HaalandS.HasegawaH.PaschmannG.SonnerupB.DunlopM. (2021). 20 Years of cluster observations: the magnetopause. J. Geophys. Res. Space Phys.126, e2021JA029362. 10.1029/2021JA029362
59
HaalandS.SonnerupB.DunlopM.BaloghA.GeorgescuE.HasegawaH.et al (2004a). Four-spacecraft determination of magnetopause orientation, motion and thickness: comparison with results from single-spacecraft methods. Ann. Geophys.22, 1347–1365. 10.5194/angeo-22-1347-2004
60
HaalandS.SonnerupB. U. Ö.DunlopM. W.GeorgescuE.PaschmannG.KleckerB.et al (2004b). Orientation and motion of a discontinuity from Cluster curlometer capability: minimum variance of current density. Geophys. Res. Lett.31. 10.1029/2004GL020001
61
HamrinM.RönnmarkK.BörlinN.VedinJ.VaivadsA. (2008). GALS – gradient analysis by least squares. Ann. Geophys.26, 3491–3499. 10.5194/angeo-26-3491-2008
62
HarveyC. C. (1998). “Spatial gradients and the volumetric tensor,” in Analysis methods for multi-spacecraft data. Editors PaschmannG.DalyP. W. (Noordwijk, Netherlands: ESA Publications Division), 307–322.
63
HeJ. S.TuC. Y.TianH.XiaoC. J.WangX. G.PuZ. Y.et al (2008a). A magnetic null geometry reconstructed from Cluster spacecraft observations. J. Geophys. Res. Space Phys.113. 10.1029/2007JA012609
64
HeJ. S.ZongQ. G.DengX. H.TuC. Y.XiaoC. J.WangX. G.et al (2008b). Electron trapping around a magnetic null. Geophys. Res. Lett.35, L14104. 10.1029/2008GL034085
65
HendersonP. D.OwenC. J.LahiffA. D.AlexeevI. V.FazakerleyA. N.YinL.et al (2008). The relationship between j × B and ∇·Pe in the magnetotail plasma sheet: cluster observations. J. Geophys. Res. Space Phys.113. 10.1029/2007ja012697
66
HwangK. J.SibeckD. G.BurchJ. L.ChoiE.FearR. C.LavraudB.et al (2018). Small-scale flux transfer events formed in the reconnection exhaust region between two X lines. J. Geophys. Res. Space Phys.123, 8473–8488. 10.1029/2018JA025611
67
KleinK. G.SpenceH.AlexandrovaO.ArgallM.ArzamasskiyL.BookbinderJ.et al (2023). HelioSwarm: a multipoint, multiscale mission to characterize turbulence. Space Sci. Rev.219, 74. 10.1007/s11214-023-01019-0
68
LavraudB.ZhangY. C.VernisseY.GershmanD. J.DorelliJ.CassakP. A.et al (2016). Currents and associated electron scattering and bouncing near the diffusion region at Earth's magnetopause. Geophys. Res. Lett.43, 3042–3050. 10.1002/2016GL068359
69
LiangJ.LiuW. W. (2007). A MHD mechanism for the generation of the meridional current system during substorm expansion phase. J. Geophys. Res. Space Phys.112. 10.1029/2007JA012303
70
LührH.ParkJ.GjerloevJ. W.RaubergJ.MichaelisI.MerayoJ. M. G.et al (2015). Field-aligned currents' scale analysis performed with the Swarm constellation. Geophys. Res. Lett.42, 1–8. 10.1002/2014gl062453
71
MarchaudonA.CerisierJ. C.DunlopM. W.PitoutF.BosquedJ. M.FazakerleyA. N. (2009). Shape, size, velocity and field-aligned currents of dayside plasma injections: a multi-altitude study. Ann. Geophys.27, 1251–1266. 10.5194/angeo-27-1251-2009
72
McPherronR. L.RussellC. T.AubryM. P. (1973). Satellite studies of magnetospheric substorms on August 15, 1968: 9. Phenomenological model for substorms. J. Geophys. Res.78, 3131–3149. 10.1029/JA078i016p03131
73
MiddletonH. R.MassonA. (2016). The Curlometer Technique: a beginner’s guide. Available at: http://www.cosmos.esa.int/web/csa/multi-spacecraft (Accessed April 11, 2024).
74
NakamuraR.BaumjohannW.FujimotoM.AsanoY.RunovA.OwenC. J.et al (2008). Cluster observations of an ion-scale current sheet in the magnetotail under the presence of a guide field. J. Geophys. Res. Space Phys.113. 10.1029/2007JA012760
75
NaritaY.NakamuraR.BaumjohannW. (2013). Cluster as current sheet surveyor in the magnetotail. Ann. Geophys.31, 1605–1610. 10.5194/angeo-31-1605-2013
76
OlsenN.LührH.FinlayC.SabakaT.MichaelisI.RaubergJ.et al (2014). The CHAOS-4 geomagnetic field model. Geophys. J. Int.197, 815–827. 10.1093/gji/ggu033
77
PanovE. V.BüchnerJ.FränzM.KorthA.SavinS. P.FornaçonK. H.et al (2006). CLUSTER observation of collisionless transport at the magnetopause. Geophys. Res. Lett.33. 10.1029/2006GL026556
78
PanovE. V.BüchnerJ.FränzM.KorthA.SavinS. P.RèmeH.et al (2008). High‐latitude Earth's magnetopause outside the cusp: cluster observations. J. Geophys. Res. Space Phys.113. 10.1029/2006ja012123
79
PaschmannG.DalyP. W. (1998). Analysis methods for multi-spacecraft data. Noordwijk, Netherlands: ESA Publications Division.
80
PaschmannG.DalyP. W. (2008). Multi-spacecraft analysis methods revisited. Noordwijk, Netherlands: ESA Communications Keplerlaan.
81
PaschmannG.SchwartzS. J.EscoubetC. P.HaalandS. (2005). Outer magnetospheric boundaries: cluster results. Springer/ISSI. 10.1007/1-4020-4582-4
82
PetrukovichA.ArtemyevA.VaskoI.NakamuraR.ZelenyiL. (2015). Current sheets in the Earth magnetotail: plasma and magnetic field structure with cluster project observations. Space Sci. Rev.188, 311–337. 10.1007/s11214-014-0126-7
83
PhanT. D.EastwoodJ. P.CassakP. A.ØierosetM.GoslingJ. T.GershmanD. J.et al (2016). MMS observations of electron-scale filamentary currents in the reconnection exhaust and near the X line. Geophys. Res. Lett.43, 6060–6069. 10.1002/2016GL069212
84
PitoutF.BogdanovaY. V. (2021). The polar cusp seen by cluster. J. Geophys. Res. Space Phys.126, e2021JA029582. 10.1029/2021JA029582
85
PuZ. Y.ZongQ. G.FritzT. A.XiaoC. J.HuangZ. Y.FuS. Y.et al (2005). Multiple flux rope events at the high-latitude magnetopause: cluster/rapid observation on 26 january, 2001. Surv. Geophys.26, 193–214. 10.1007/s10712-005-1878-0
86
RetinòA.KhotyaintsevY.Le ContelO.MarcucciM. F.PlaschkeF.VaivadsA.et al (2022). Particle energization in space plasmas: towards a multi-point, multi-scale plasma observatory. Exp. Astron.54, 427–471. 10.1007/s10686-021-09797-7
87
RitterP.LührH. (2006). Curl-B technique applied to Swarm constellation for determining field-aligned currents. Earth, Planets Space58, 463–476. 10.1186/BF03351942
88
RitterP.LührH.RaubergJ. (2013). Determining field-aligned currents with the Swarm constellation mission. Earth, Planets Space65, 1285–1294. 10.5047/eps.2013.09.006
89
RitterP.LührH.ViljanenA.AmmO.PulkkinenA.SillanpääI. (2004). Ionospheric currents estimated simultaneously from CHAMP satelliteand IMAGE ground-based magnetic field measurements: a statisticalstudy at auroral latitudes. Ann. Geophys.22, 417–430. 10.5194/angeo-22-417-2004
90
RobertP.DunlopM. W. (2022). Use of twenty years CLUSTER/FGM data to observe the mean behavior of the magnetic field and current density of Earth’s magnetosphere. J. Geophys. Res. Space Phys.127. 10.1029/2021ja029837
91
RobertP.DunlopM. W.RouxA.ChanteurG. (1998a). “Accuracy of current density determination,” in Analysis methods for multi-spacecraft data. Editors PaschmannG.DalyP. W. (Noordwijk, Netherlands: ESA Publications Division), 395–418.
92
RobertP.RouxA.HarveyC. C.DunlopM. W.DalyP. W.GlassmeierK.-H. (1998b). “Tetrahedron geometric factors,” in Analysis methods for multi-spacecraft data. Editors PaschmannG.DalyP. W. (Noordwijk, Netherlands: ESA Publications Division), 323–348.
93
RongZ. J.WanW. X.ShenC.LiX.DunlopM. W.PetrukovichA. A.et al (2011). Statistical survey on the magnetic structure in magnetotail current sheets. J. Geophys. Res. Space Phys.116. 10.1029/2011JA016489
94
RouxA.RobertP.FontaineD.ContelO. L.CanuP.LouarnP. (2015). What is the nature of magnetosheath FTEs?J. Geophys. Res. Space Phys.120, 4576–4595. 10.1002/2015JA020983
95
RunovA.SergeevV. A.BaumjohannW.NakamuraR.ApatenkovS.AsanoY.et al (2005). Electric current and magnetic field geometry in flapping magnetotail current sheets. Ann. Geophys.23, 1391–1403. 10.5194/angeo-23-1391-2005
96
RunovA.SergeevV. A.NakamuraR.BaumjohannW.ApatenkovS.AsanoY.et al (2006). Local structure of the magnetotail current sheet: 2001 Cluster observations. Ann. Geophys.24, 247–262. 10.5194/angeo-24-247-2006
97
RussellC. T.LuhmannJ. G.StrangewayR. J. (2016). Space physics: an introduction. Cambridge: Cambridge University Press. 10.1017/9781316162590
98
ShenC.DunlopM. (2023). Field gradient analysis based on a geometrical approach. J. Geophys. Res. Space Phys.128, e2023JA031313. 10.1029/2023JA031313
99
ShenC.DunlopM.MaY. H.ChenZ. Q.YanG. Q.LiuZ. X.et al (2011). The magnetic configuration of the high-latitude cusp and dayside magnetopause under strong magnetic shears. J. Geophys. Res. Space Phys.116, n/a-n/a. 10.1029/2011ja016501
100
ShenC.DunlopM. W. (2008). “Geometrical structure analysis of the magnetic field,” in Multi-spacecraft analysis methods revisited. Editors PaschmannG.DalyP. W. (Noordwijk, Netherlands: ESA Communications Keplerlaan), 27–32.
101
ShenC.LiX.DunlopM.LiuZ. X.BaloghA.BakerD. N.et al (2003). Analyses on the geometrical structure of magnetic field in the current sheet based on cluster measurements. J. Geophys. Res. Space Phys.108. 10.1029/2002ja009612
102
ShenC.LiX.DunlopM.ShiQ. Q.LiuZ. X.LucekE.et al (2007). Magnetic field rotation analysis and the applications. J. Geophys. Res. Space Phys.112. 10.1029/2005ja011584
103
ShenC.LiuZ. X.LiX.DunlopM.LucekE.RongZ. J.et al (2008a). Flattened current sheet and its evolution in substorms. J. Geophys. Res. (Space Physics)113, A07S21. 10.1029/2007JA012812
104
ShenC.RongZ. J.DunlopM. (2012a). Determining the full magnetic field gradient from two spacecraft measurements under special constraints. J. Geophys. Res. Space Phys.117. 10.1029/2012JA018063
105
ShenC.RongZ. J.DunlopM. W.MaY. H.LiX.ZengG.et al (2012b). Spatial gradients from irregular, multiple‐point spacecraft configurations. J. Geophys. Res. Space Phys.117. 10.1029/2012ja018075
106
ShenC.RongZ. J.LiX.DunlopM.LiuZ. X.MalovaH. V.et al (2008b). Magnetic configurations of the tilted current sheets in magnetotail. Ann. Geophys.26, 3525–3543. 10.5194/angeo-26-3525-2008
107
ShenC.YangY. Y.RongZ. J.LiX.DunlopM.CarrC. M.et al (2014). Direct calculation of the ring current distribution and magnetic structure seen by Cluster during geomagnetic storms. J. Geophys. Res. Space Phys.119, 2458–2465. 10.1002/2013ja019460
108
ShenC.ZhangC.RongZ.PuZ.DunlopM. W.EscoubetC. P.et al (2021a). Nonlinear magnetic gradients and complete magnetic geometry from multispacecraft measurements. J. Geophys. Res. Space Phys.126, e2020JA028846. 10.1029/2020JA028846
109
ShenC.ZhouY.GaoL.WangX.PuZ.EscoubetC. P.et al (2021b). Measurements of the net charge density of space plasmas. J. Geophys. Res. Space Phys.126, e2021JA029511. 10.1029/2021JA029511
110
ShenC.ZhouY.MaY.WangX.PuZ.DunlopM. (2021c). A general algorithm for the linear and quadratic gradients of physical quantities based on 10 or more point measurements. J. Geophys. Res. Space Phys.126, e2021JA029121. 10.1029/2021JA029121
111
ShiJ.GuoJ.DunlopM.ZhangT.LiuZ.LucekE.et al (2012). Inter-hemispheric asymmetry of dependence of the cusp location on dipole tilt during northward IMF conditions. Ann. Geophys.30, 21–26. 10.5194/angeo-30-21-2012
112
ShiJ. K.ChengZ. W.ZhangT. L.DunlopM.LiuZ. X.TorkarK.et al (2010). South-north asymmetry of field-aligned currents in the magnetotail observed by Cluster. J. Geophys. Res. Space Phys.115. 10.1029/2009JA014446
113
ShiQ. Q.TianA. M.BaiS. C.HasegawaH.DegelingA. W.PuZ. Y.et al (2019). Dimensionality, coordinate system and reference frame for analysis of in-situ space plasma and field data. Space Sci. Rev.215, 35. 10.1007/s11214-019-0601-2
114
ShiokawaK.BaumjohannW.HaerendelG.PaschmannG.FennellJ. F.Friis-ChristensenE.et al (1998). High-speed ion flow, substorm current wedge, and multiple Pi 2 pulsations. J. Geophys. Res. Space Phys.103, 4491–4507. 10.1029/97JA01680
115
ShoreR. M.WhalerK. A.MacmillanS.BegganC.OlsenN.SpainT.et al (2013). Ionospheric midlatitude electric current density inferred from multiple magnetic satellites. J. Geophys. Res. Space Phys.118, 5813–5829. 10.1002/jgra.50491
116
TanX.DunlopM. W.DongX. C.YangY. Y.DuY. S.ShenC.et al (2023). Ring current morphology from MMS observations. J. Geophys. Res. Space Phys.128. 10.1029/2023ja031372
117
TehW. L.NakamuraT. K. M.NakamuraR.BaumjohannW.RussellC. T.PollockC.et al (2017). Evolution of a typical ion-scale magnetic flux rope caused by thermal pressure enhancement. J. Geophys. Res. Space Phys.122, 2040–2050. 10.1002/2016JA023777
118
TrenchiL.KauristieK.KäkiS.VanhamäkiH.JuusolaL.BlagauA.et al (2020). “ESA field-aligned currents—methodology inter-comparison exercise,” in Ionospheric multi-spacecraft analysis tools. Editors DunlopM. W.LührH. (Springer), 167. 10.1007/978-3-030-26732-2_8
119
VallatC.DandourasI.DunlopM.BaloghA.LucekE.ParksG. K.et al (2005). First current density measurements in the ring current region using simultaneous multi-spacecraft CLUSTER-FGM data. Ann. Geophys.23, 1849–1865. 10.5194/angeo-23-1849-2005
120
VogtJ.AlbertA.MarghituO. (2009). Analysis of three-spacecraft data using planar reciprocal vectors: methodological framework and spatial gradient estimation. Ann. Geophys.27, 3249–3273. 10.5194/angeo-27-3249-2009
121
VogtJ.PaschmannG.ChanteurG. (2008). “Reciprocal vectors,” in Multi-spacecraft analysis methods revisited. Editors PaschmannG.DalyP. W. (Noordwijk, Netherlands: ESA Communications Keplerlaan), 33–46.
122
VogtJ.SorbaloE.HeM.BlagauA. (2013). Gradient estimation using configurations of two or three spacecraft. Ann. Geophys.31, 1913–1927. 10.5194/angeo-31-1913-2013
123
WangH.MaS. Y.LührH.LiuZ. X.PuZ. Y.EscoubetC. P.et al (2006). Global manifestations of a substorm onset observed by a multi-satellite and ground station network. Ann. Geophys.24, 3491–3496. 10.5194/angeo-24-3491-2006
124
WangT.AlexandrovaO.PerroneD.DunlopM.DongX.BinghamR.et al (2019). Magnetospheric multiscale observation of kinetic signatures in the Alfvén vortex. Astrophysical J.871, L22. 10.3847/2041-8213/aafe0d
125
WangZ.FuH. S.LiuC. M.LiuY. Y.CozzaniG.GilesB. L.et al (2019). Electron distribution functions around a reconnection X-line resolved by the FOTE method. Geophys. Res. Lett.46, 1195–1204. 10.1029/2018GL081708
126
WangZ.FuH. S.VaivadsA.BurchJ. L.YuY.CaoJ. B. (2020). Monitoring the spatio-temporal evolution of a reconnection X-line in space. Astrophysical J.899, L34. 10.3847/2041-8213/abad2c
127
XiaoC.LiuW.ShenC.ZhangH.RongZ. (2018). Study on the curvature and gradient of the magnetic field in Earth's cusp region based on the magnetic curvature analysis method. J. Geophys. Res. (Space Physics)123, 3794–3805. 10.1029/2017JA025028
128
XiaoC. J.PuZ. Y.MaZ. W.FuS. Y.HuangZ. Y.ZongQ. G. (2004). Inferring of flux rope orientation with the minimum variance analysis technique. J. Geophys. Res. Space Phys.109. 10.1029/2004JA010594
129
XiaoC. J.WangX. G.PuZ. Y.MaZ. W.ZhaoH.ZhouG. P.et al (2007). Satellite observations of separator-line geometry of three-dimensional magnetic reconnection. Nat. Phys.3, 609–613. 10.1038/nphys650
130
XiaoC. J.WangX. G.PuZ. Y.ZhaoH.WangJ. X.MaZ. W.et al (2006). In situ evidence for the structure of the magnetic null in a 3D reconnection event in the Earth's magnetotail. Nat. Phys.2, 478–483. 10.1038/nphys342
131
YangY. Y.ShenC.DunlopM.RongZ. J.LiX.AngelopoulosV.et al (2016). Storm time current distribution in the inner equatorial magnetosphere: THEMIS observations. J. Geophys. Res. Space Phys.121, 5250–5259. 10.1002/2015ja022145
132
YangY. Y.ShenC.ZhangY. C.RongZ. J.LiX.DunlopM.et al (2014). The force‐free configuration of flux ropes in geomagnetotail: cluster observations. J. Geophys. Res. Space Phys.119, 6327–6341. 10.1002/2013ja019642
133
YaoZ. H.PuZ. Y.OwenC. J.FuS. Y.ChuX. N.LiuJ.et al (2014). Current reduction in a pseudo‐breakup event: THEMIS observations. J. Geophys. Res. Space Phys.119, 8178–8187. 10.1002/2014ja020186
134
ZhangQ. H.DunlopM. W.LockwoodM.HolmeR.KamideY.BaumjohannW.et al (2011). The distribution of the ring current: cluster observations. Ann. Geophys.29, 1655–1662. 10.5194/angeo-29-1655-2011
135
ZhangY. C.ShenC.LiuZ. X.RongZ. J.ZhangT. L.MarchaudonA.et al (2013). Two different types of plasmoids in the plasma sheet: cluster multisatellite analysis application. J. Geophys. Res. Space Phys.118, 5437–5444. 10.1002/jgra.50542
136
ZhangY. C.ShenC.MarchaudonA.RongZ. J.LavraudB.FazakerleyA.et al (2016). First in situ evidence of electron pitch angle scattering due to magnetic field line curvature in the Ion diffusion region. J. Geophys. Res. Space Phys.121, 4103–4110. 10.1002/2016JA022409
137
ZhaoC.RussellC. T.StrangewayR. J.PetrinecS. M.PatersonW. R.ZhouM.et al (2016). Force balance at the magnetopause determined with MMS: application to flux transfer events. Geophys. Res. Lett.43, 11941–11947. 10.1002/2016GL071568
Summary
Keywords
curlometer analysis, multi-spacecraft, analysis methods, magnetosphere, magnetic gradients and reconstruction
Citation
Dunlop MW, Fu H-S, Shen C, Tan X, Dong X-C, Yang Y-Y, Robert P and Escoubet CP (2024) Curlometer and gradient techniques: past and future applications. Front. Astron. Space Sci. 11:1422341. doi: 10.3389/fspas.2024.1422341
Received
23 April 2024
Accepted
05 December 2024
Published
24 December 2024
Volume
11 - 2024
Edited by
David Knudsen, University of Calgary, Canada
Reviewed by
Vladimir Krasnoselskikh, UMR7328 Laboratoire de Physique et Chimie de l’Environnement et de l’Espace (LPC2E), France
R. S. Pandey, Amity University, India
Updates
Copyright
© 2024 Dunlop, Fu, Shen, Tan, Dong, Yang, Robert and Escoubet.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: M. W. Dunlop, malcolm.dunlop@stfc.ac.uk
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.