---
title: 'OpenGGCM: Global Geospace MHD Model'
url: https://www.emergentmind.com/topics/open-geospace-general-circulation-model-openggcm
type: topic
---

# OpenGGCM: Global Geospace MHD Model

Searching arXiv for the specified OpenGGCM-related papers to ground the article in current literature.
The Open Geospace General Circulation Model (OpenGGCM) is a global magnetohydrodynamic (MHD) model of the geospace system that, in the configurations described in recent work, represents the coupled response of Earth’s magnetosphere, ionosphere, and thermosphere to time-dependent solar wind forcing, and in one coupled framework also includes the inner magnetosphere’s kinetic ring current via the Rice Convection Model (RCM) [2403.19873]. In the literature summarized here, OpenGGCM appears in three closely related roles: as a physics-based geospace model driven by heliospheric or spacecraft solar wind inputs to predict geomagnetic activity and indices [2403.19873]; as one of several global MHD models used to compare modeled and observed horizontal ground magnetic perturbations during the May 2024 geomagnetic storm [2507.07009]; and as a reference case for examining how magnetic divergence control and finite-domain boundary terms affect Biot–Savart estimates of the magnetic field at Earth derived from MHD simulations [2508.18405]. Across these uses, a recurring theme is that OpenGGCM’s constrained-transport treatment of the magnetic field keeps $\nabla \cdot \mathbf{B}$ extremely small, while finite-domain outer-boundary effects remain non-negligible in ground magnetic field estimation [2508.18405].

## 1. Model definition and coupled geospace architecture

OpenGGCM is described as a global magnetosphere model that solves a semi-conservative form of the MHD equations on a three-dimensional, stretched Cartesian grid using a second-order predictor–corrector finite-difference scheme [2508.18405]. In the coupled EUHFORIA–OpenGGCM framework, the outer magnetosphere is solved as a single-fluid, semi-conservative MHD system on a 3D stretched Cartesian grid in GSE coordinates, and the model is electrodynamically coupled to the Coupled Thermosphere Ionosphere Model (CTIM) and, crucially, to the inner magnetosphere’s kinetic ring current via the Rice Convection Model [2403.19873].

The three subsystems identified in that configuration are the magnetosphere, the ionosphere–thermosphere system, and the inner magnetosphere ring current [2403.19873]. The ionosphere–thermosphere coupling interface is at $2.1\,R_E$, where field-aligned currents computed by the magnetosphere close through the ionosphere, while magnetospheric precipitation parameters and electric fields drive ionospheric conductances and dynamo currents [2403.19873]. CTIM returns height-integrated conductances and dynamo currents to solve the ionospheric electric potential, which feeds back to the magnetosphere [2403.19873]. The RCM is embedded to represent kinetic ring-current physics, with a dayside boundary of approximately $10\,R_E$, a nightside boundary of approximately $13\,R_E$, and flanks at $\pm 12\,R_E$; it is initialized from OpenGGCM’s MHD density and pressure and returns corrected ring-current pressure and density to the MHD solution [2403.19873].

This coupled architecture matters because the recent literature does not treat OpenGGCM merely as a magnetopause-scale MHD solver. Instead, it is used as a system model whose outputs include both magnetospheric state variables and ionospheric electrodynamic quantities, and whose ring-current treatment is sufficiently explicit to support model-derived geomagnetic indices such as Dst and auroral indices [2403.19873]. A plausible implication is that OpenGGCM’s operational and research value depends not only on the outer-MHD numerics but also on the fidelity of its electrodynamic closures and inner-magnetosphere coupling.

## 2. Governing equations, numerics, and divergence control

In the cited work, OpenGGCM solves the ideal MHD equations in semi-conservative form with source terms and numerical divergence control via constrained transport (CT) [2403.19873]. The equations reported include continuity,
$$
\frac{\partial \rho}{\partial t} + \nabla \cdot (\rho \mathbf{v}) = 0,
$$
momentum in Maxwell-stress form, induction,
$$
\frac{\partial \mathbf{B}}{\partial t}
= \nabla \times \left(\mathbf{v} \times \mathbf{B} - \eta \,\nabla \times \mathbf{B}\right),
$$
and total-energy evolution with adiabatic closure, with $\gamma$ typically $5/3$ [2403.19873]. The model uses an explicit second-order predictor–corrector scheme and a stretched Cartesian mesh [2403.19873].

A central numerical property emphasized in the literature is OpenGGCM’s treatment of the solenoidal constraint. To control magnetic divergence, OpenGGCM employs the Constrained Transport method on a staggered grid, and CT “maintains $\nabla \cdot \mathbf{B} = 0$ to round-off error” [2508.18405]. No ad hoc divergence cleaning parameters, such as diffusive cleaning coefficients or hyperbolic/parabolic cleaning speeds, are used; rather, the algorithmic choice of CT with staggered storage and edge-centered updates by the curl of the electromotive force enforces the solenoidal constraint intrinsically [2508.18405]. Historically, the main disadvantage of CT is identified as the difficulty of extending it to generalized grids, with Evans and Hawley, Gombosi et al., and Mignone et al. cited in the paper summarized by [2508.18405].

In the specific coupled configuration described in [2403.19873], the computational domain contains $481$ cells in $x$ spanning $-35\,R_E$ to $+5000\,R_E$, and $180$ cells each in $y$ and $z$ spanning $-48\,R_E$ to $+48\,R_E$, with minimum cell sizes in the inner magnetosphere of approximately $0.167\,R_E$ in $x$ and $0.25\,R_E$ in $y,z$. Output is produced as 3D MHD fields at hourly cadence, 2D magnetosphere cuts at one-minute cadence, and ionosphere electrodynamic fields at one-minute cadence [2403.19873].

The significance of these details is clearest in the magnetic-field reconstruction study. There, OpenGGCM’s CT scheme is directly linked to the finding that divergence-related corrections to ground magnetic perturbation estimates are very small [2508.18405]. This suggests that, among errors arising in finite-volume magnetic-field inference from MHD simulations, OpenGGCM suppresses non-solenoidal contamination more effectively than schemes that rely on approximate divergence transport or cleaning.

## 3. Finite-volume magnetic-field decomposition and ground perturbation estimation

A substantial recent analysis of OpenGGCM concerns the mathematical consistency of estimating the magnetic field at Earth from simulation currents using the Biot–Savart law [2508.18405]. Maxwell’s equation enforces
$$
\nabla \cdot \mathbf{B} = 0,
$$
but numerical error may violate this condition in simulations [2508.18405]. The magnetospheric contribution to the magnetic field at a point $\mathbf{r}$ on Earth from current density $\mathbf{J}(\mathbf{r}')$ is estimated as
$$
\mathbf{B}(\mathbf{r}) = \frac{\mu_0}{4\pi} \int_V \frac{\mathbf{J}(\mathbf{r}') \times (\mathbf{r}-\mathbf{r}')}{\|\mathbf{r}-\mathbf{r}'\|^3}\, d^3r',
$$
where $V$ is the magnetospheric simulation volume and $\mu_0$ is the permeability of free space [2508.18405].

The study emphasizes that Helmholtz decomposition on a finite domain requires both volume and surface terms [2508.18405]. For a point $x_0$ in the inner gap region $G$, the paper writes the Biot–Savart contribution compactly as $B_{\mathrm{BS}}(x_0):=\#1\{M\}$ and derives
$$
B_{\mathrm{BS}}(x_0) = - \#1\{M\} + \mathrm{surfInt}\{\partial M\} - \mathrm{surfInt}\{\partial G\},
$$
showing that the Biot–Savart volume integral is balanced by a volume integral involving $\nabla \cdot \mathbf{B}$ and surface integrals on the outer and inner boundaries [2508.18405]. The complementary formula for the gap region is
$$
B(x_0) = - \mathrm{surfInt}\{\partial G\},
$$
which the authors note is a second way to compute the magnetospheric contribution to $B$ on Earth when the inner region is treated as dipolar [2508.18405]. The paper further defines
$$
B_{\mathrm{div}}(x_0):= - \#1\{M\}, \qquad
B_{\mathrm{out}}(x_0):= + \mathrm{surfInt}\{\partial M\},
$$
and states the consistency condition
$$
B(x_0) - B_{\mathrm{BS}}(x_0) = B_{\mathrm{div}}(x_0) + B_{\mathrm{out}}(x_0),
$$
so that Biot–Savart alone is exact only when both the divergence volume term and outer-boundary surface term vanish [2508.18405].

Within this framework, the outer-boundary surface integral has a specific physical interpretation. $B_{\mathrm{out}}(x_0)$ is the contribution arising from the finite extent of the simulation; when the integration domain is not all-space, Helmholtz decomposition introduces a surface term over $\partial V$ that encodes how the boundary fields and their normal derivatives, through the Green’s function for the Poisson operator, transmit information from outside the magnetospheric volume into the interior solution for $B$ [2508.18405]. In the reported calculations, this term is evaluated numerically by interpolating the boundary integrands over the outer surface and integrating [2508.18405].

For OpenGGCM specifically, the current density used in these integrals is the model-provided $\mathbf{j}$ from the MHD outputs, and the analysis verifies that $\mathbf{j}=(1/\mu_0)\,\nabla \times \mathbf{B}$ in the quasi-steady MHD regime [2508.18405]. The inner boundary is a sphere of radius $3\,R_E$ centered at Earth, and the outer boundary is a rectangular prism in GSE coordinates spanning $X$ from $-350\,R_E$ to $+60\,R_E$ and $Y,Z$ from $-48\,R_E$ to $+48\,R_E$ [2508.18405]. For the $\nabla \cdot \mathbf{B}$ volume term, second-order finite-difference stencils are used to compute $\nabla \cdot \mathbf{B}$ in the magnetospheric volume [2508.18405].

The broader significance is methodological. The study shows that magnetic perturbations at Earth inferred from MHD simulations are not determined by magnetospheric current density alone unless the finite-domain correction terms are negligible. In OpenGGCM, the divergence term is nearly suppressed by construction, but the outer-boundary term is not, so finite-domain consistency remains essential [2508.18405].

## 4. Quantitative behavior of divergence and boundary terms in OpenGGCM

The quantitative findings reported for OpenGGCM are unusually specific. In the simulations examined, when $|B_{\mathrm{BS}}|$ is large, such as after a southward turning of the interplanetary magnetic field, surface-averaged ratios over Earth’s surface show $B_{\mathrm{out}}$ at about $\sim 20\%$ of $B_{\mathrm{BS}}$ and $B_{\mathrm{div}}<1\%$ of $B_{\mathrm{BS}}$ [2508.18405]. When $|B_{\mathrm{BS}}|$ is small, both ratios increase in percentage terms, although OpenGGCM’s divergence contribution remains very small because of CT [2508.18405].

The analyzed solar wind driving for the OpenGGCM run uses GSM conditions at $33\,R_E$ sunward with density $n=5\,\mathrm{cm}^{-3}$, temperature $T=10^5\,\mathrm{K}$, $V_X=-400\,\mathrm{km/s}$, and $B_Y^{\mathrm{IMF}}=B_Z^{\mathrm{IMF}}=B_X^{\mathrm{IMF}}=V_Y=V_Z=0$ except that $B_Z^{\mathrm{IMF}}$ is $+5\,\mathrm{nT}$ from 00:00–06:00 and flips to $-10\,\mathrm{nT}$ at 06:00 and remains until 20:00 [2508.18405]. Dipole tilt is fixed at $0$, the run spans 00:00–20:00 UTC, and snapshots are examined at 01:00, 04:00, 07:00, 10:00, and 16:00 UTC [2508.18405]. Before 06:00, magnetospheric contributions at the ground are small and the relative fractions of $B_{\mathrm{out}}$ and $B_{\mathrm{div}}$ are higher; after the southward turning, magnetospheric currents intensify, $B_{\mathrm{BS}}$ becomes the dominant term, and the fractional contributions fall [2508.18405].

Spatial structure is also reported. Both $B_{\mathrm{div}}$ and $B_{\mathrm{out}}$ are spatially correlated, so nearby locations on Earth show integrals of similar magnitudes [2508.18405]. Scatter plots colored by geolatitude and geolongitude reveal coherent bands, reflecting the coherence of ground magnetic perturbations [2508.18405]. In the OpenGGCM analysis, $B_{\mathrm{out}}$ versus $B_{\mathrm{BS}}$ shows the closest match to the “desired” behavior for justifying Biot–Savart, but there is observable anti-correlated behavior after the IMF flip, indicating systematic boundary contributions that partially offset or modify the interior current contributions [2508.18405].

These results are important because they separate two distinct uncertainties in Biot–Savart-only analyses. In OpenGGCM, the self-consistency uncertainty associated with nonzero $\nabla \cdot \mathbf{B}$ is negligible in the cases examined, whereas the approximation uncertainty from omitting the outer boundary term remains at roughly the $20\%$ level in surface average during strong activity and becomes relatively larger when the Biot–Savart estimate is small [2508.18405]. The paper therefore recommends including the outer boundary term or using the inner boundary formulation $-\mathrm{surfInt}\{\partial G\}$ when possible [2508.18405].

## 5. Comparisons with other MHD frameworks and implications for model interpretation

OpenGGCM is compared directly with SWMF/BATS-R-US in the boundary-integral study and with MAGE and SWMF in the May 2024 geomagnetic storm study [2508.18405; 2507.07009]. In the divergence-control comparison, OpenGGCM’s CT keeps $\nabla \cdot \mathbf{B}$ at round-off, and the measured $B_{\mathrm{div}}$ is always very small, below $1\%$ for large $|B_{\mathrm{BS}}|$ [2508.18405]. By contrast, BATS-R-US uses the eight-wave method, and the paper notes its known limitation that $\nabla \cdot \mathbf{B}$ can grow where sources are large, with measured $B_{\mathrm{div}}$ reaching $5$–$20\%$ for large $|B_{\mathrm{BS}}|$ [2508.18405]. In both models, however, the outer boundary contribution is typically larger than the divergence term; in OpenGGCM it is approximately $20\%$, whereas in BATS-R-US it is $20$–$30\%$ at times of strong activity [2508.18405].

The event-type sensitivity reported in [2508.18405] is qualitative but consistent across scenarios. In the simple IMF-flip scenario, both OpenGGCM and BATS-R-US show smaller fractional correction terms when $|B_{\mathrm{BS}}|$ is large and larger fractions when $|B_{\mathrm{BS}}|$ is small [2508.18405]. In a separate Carrington superstorm case analyzed with BATS-R-US, the same pattern holds, and the paper concludes that the trends are robust across a simple IMF step and a complex storm [2508.18405]. OpenGGCM’s specific superstorm results are not provided there.

A separate caution concerns inter-model spread in ground magnetic estimates. The paper notes that differences between MHD models can reach factors of approximately $2$ for specific values of $B$ at ground, although trends are similar and cross-model agreement improves when $|B|$ is large [2508.18405]. It also states that Biot–Savart has known limitations in strong, non-symmetric, and time-dependent situations, and that MHD solutions exhibit finite signal propagation times, so these physical limitations should be considered when interpreting fast ground magnetic changes [2508.18405].

These comparisons support an interpretation in which OpenGGCM’s main advantage in this context is not the elimination of all magnetic-field inference errors, but the strong suppression of divergence-related error. This suggests that model-to-model discrepancies in ground perturbation products may persist even when solenoidality is well controlled, because finite-volume boundary effects and broader physical-model differences remain.

## 6. Forecasting workflows, geomagnetic indices, and heliospheric coupling

OpenGGCM is also used as the geospace component of a forecasting chain in which EUHFORIA provides time-dependent solar wind and IMF forcing at Earth [2403.19873]. In that framework, OpenGGCM is driven by plasma velocity $\mathbf{v}$, magnetic field $\mathbf{B}$, proton number density $n_p$, and thermal pressure $P$ obtained from EUHFORIA simulations, with EUHFORIA providing 1D time series at Earth with 10-minute cadence [2403.19873]. A relaxation window of approximately $10$–$12$ hours is required before the main phase, and in practice ingestion starts approximately $5$–$6$ hours before the CME shock arrival so that the ejecta arrives after the relaxation period [2403.19873]. To avoid numerical instabilities from abrupt IMF changes, the simulation is initialized with a low but non-zero $B_z$ of about $1\,\mathrm{nT}$ [2403.19873].

The coupled framework validates two observed geoeffective CME events. The first is the 2012-07-12 CME, simulated in EUHFORIA v2.0 with the magnetized spheromak flux-rope CME model using insertion time 2012-07-12 19:24 UT, speed $763\,\mathrm{km/s}$, longitude $-4^\circ$, latitude $-8^\circ$, radius $16.8\,R_\odot$, density $1\times10^{-18}\,\mathrm{kg\,m^{-3}}$, temperature $0.8\times10^6\,\mathrm{K}$, helicity $+1$, tilt $135^\circ$, and toroidal magnetic flux $1\times10^{14}\,\mathrm{Wb}$ [2403.19873]. The second consists of three interacting CMEs in September 2017, also represented with spheromaks and explicitly parameterized in speed, direction, size, helicity, tilt, and toroidal flux [2403.19873]. The paper states that these parameters control the modeled IMF structure and solar wind forcing at Earth, especially the sign and duration of $B_z$ and the dynamic pressure, both of which strongly influence OpenGGCM’s geospace response [2403.19873].

Within this workflow, OpenGGCM computes a model Dst by integrating currents in the inner magnetosphere via Biot–Savart to infer ground magnetic perturbations; this quantity represents the ring-current contribution, or “Dst*” in Burton-type terminology [2403.19873]. The magnetopause current contribution is then reintroduced by inverting the correction
$$
\mathrm{Dst}^*=\frac{\mathrm{Dst}-b\sqrt{P_{\mathrm{dyn}}}+c}{a},
$$
with coefficients $a=1$, $b=7.26$, and $c=11$ from the O’Brien–McPherron AK2 model [2403.19873]. Auroral indices AU and AL are constructed from model ionospheric currents by Biot–Savart integration over the electrojet region, with $AE = AU - AL$ [2403.19873].

Model performance is assessed with dynamic time warping (DTW), using
$$
\mathrm{DTW}(i,j)=d(s_i,q_j)+\min\{\mathrm{DTW}(i-1,j),\mathrm{DTW}(i,j-1),\mathrm{DTW}(i-1,j-1)\},
$$
where $d(s_i,q_j)=|s_i-q_j|$, and a Sequence Similarity Factor
$$
\mathrm{SSF}=\frac{\mathrm{DTW}_{\mathrm{score}(O,M)}}{\mathrm{DTW}_{\mathrm{score}(O,N)}}.
$$
Lower SSF indicates better alignment than a null reference [2403.19873]. For the 2012 event, DTW scores are $15076$ for the observed-input run and $17040$ for the EUHFORIA-driven run, with $\mathrm{DTW}(O,N)=37391$, giving $\mathrm{SSF}=0.40$ and $0.45$, respectively [2403.19873]. For the 2017 event, the scores are $21117$ and $13945$ against $\mathrm{DTW}(O,N)=30823$, giving $\mathrm{SSF}=0.68$ for the observed-input run and $0.45$ for the EUHFORIA-driven run [2403.19873].

Extremes reported in the event study show that the coupled framework can reproduce storm morphology while still exhibiting amplitude biases. For the 2012 event, observed Dst is $-122\,\mathrm{nT}$, while Event1-obs gives $-151\,\mathrm{nT}$ and Event1-euh gives $-205\,\mathrm{nT}$ [2403.19873]. For the 2017 event, observed Dst is $-144\,\mathrm{nT}$, Event2-obs gives $-175\,\mathrm{nT}$, and Event2-euh gives $-156\,\mathrm{nT}$ [2403.19873]. The study’s main highlight is the use of EUHFORIA simulated time series to predict Dst and auroral indices $1$ to $2$ days in advance, compared with the $1$ to $2$ hours available from L1 observations [2403.19873].

The importance of this result is practical as well as methodological. It shows OpenGGCM functioning as the terminal geospace response model in a heliosphere-to-magnetosphere forecasting chain. At the same time, the reported overestimation tendencies and timing offsets indicate that predictive skill depends strongly on upstream $B_z$ and dynamic pressure accuracy and on careful initialization of the geospace simulation [2403.19873].

## 7. Ground magnetic perturbations, GIC context, and operational limits

During the May 10–12, 2024 geomagnetic storm, OpenGGCM was one of three global magnetosphere MHD models used to simulate horizontal ground magnetic perturbations $\Delta B_H$ for comparison with magnetometer observations across the contiguous United States [2507.07009]. In that study, OpenGGCM outputs were not used to drive geoelectric fields or compute geomagnetically induced currents (GIC); instead, the comparison focused on how well model-predicted $\Delta B_H$ tracked measured $\Delta B_H$ at mid-latitude sites during the most disturbed phases of the storm [2507.07009].

The model-predicted horizontal magnetic perturbation is defined explicitly as
$$
\Delta B_H = \sqrt{(\Delta B_x)^2 + (\Delta B_y)^2},
$$
and the validation uses model outputs at one-minute cadence, with observed one-second magnetometer data averaged to one minute to match [2507.07009]. The time interval analyzed for skill metrics is 2024-05-10T15Z to 2024-05-12T06Z, and the model–data comparison covers $17$ unique magnetometer sites where simulation results were available, consisting of all $15$ NERC magnetometer sites plus magnetometers near two TVA sites [2507.07009]. The upstream drivers for the simulation models are the interplanetary magnetic field and solar wind parameters from the NOAA DSCOVR spacecraft [2507.07009].

The skill metrics are correlation coefficient $r$ and prediction efficiency (PE), defined in the same form used for GIC,
$$
\mathrm{PE}=1-\frac{\sum_t (y_t-\hat y_t)^2}{\sum_t (y_t-\bar y)^2},
$$
or, in the GIC implementation,
$$
\mathrm{PE}=1-\frac{\sum_t (\mathrm{GIC}_p(t)-\mathrm{GIC}_m(t))^2}{\sigma_m^2}.
$$
When means and variances are equal, $\mathrm{PE}=r^2$ [2507.07009]. The abstract states that $\Delta B_H$ predicted by MAGE, SWMF, and OpenGGCM had correlations ranging from $0.4$ to $0.6$ [2507.07009]. Mean metrics are reported explicitly for MAGE and SWMF only, not for OpenGGCM: MAGE has mean $r=0.44\pm0.03$ and mean $\mathrm{PE}=-3.5\pm0.7$, whereas SWMF has mean $r=0.58\pm0.01$ and mean $\mathrm{PE}=-0.682\pm0.07$ [2507.07009]. By implication and from the abstract’s range, OpenGGCM correlations fell between $0.4$ and $0.6$, but the study does not report OpenGGCM’s mean PE, RMSE, or bias [2507.07009].

The paper’s operational conclusion is explicit: the poor prediction efficiency seen for global MHD $\Delta B_H$ implies that GIC derived from model $\Delta B_H$ would be significantly degraded relative to GIC driven by measured $\Delta B_H$ convolved with magnetotelluric transfer functions [2507.07009]. In the same study, measured and modeled GIC by the Tennessee Valley Authority at four sites achieved correlation coefficients above $0.8$ and prediction efficiencies between approximately $0.4$ and $0.7$, demonstrating the stronger performance of measured-$\Delta B$-driven methods [2507.07009].

This establishes an important boundary on OpenGGCM’s present operational use. As deployed in that event, it provides physics-based geospace context and qualitative timing and relative magnitude of $\Delta B_H$, but its quantitative variance capture for ground perturbations in the contiguous United States was insufficient for reliable operational GIC forecasting when compared with methods based on measured ground magnetic fields [2507.07009]. A plausible implication is that OpenGGCM is better suited, in its current evaluated form, to geospace state interpretation and scenario forecasting than to stand-alone induction-hazard prediction at the level required by power-system operations.

Source: https://www.emergentmind.com/topics/open-geospace-general-circulation-model-openggcm