06.07.2013 Views

Scientific Theme: Advanced Modeling and Observing Systems

Scientific Theme: Advanced Modeling and Observing Systems

Scientific Theme: Advanced Modeling and Observing Systems

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

<strong>Theme</strong> report: <strong>Advanced</strong> <strong>Modeling</strong> <strong>and</strong> <strong>Observing</strong> <strong>Systems</strong><br />

The vertical extension of the model domain has required the incorporation of additional physical processes <strong>and</strong><br />

substantial changes in the dynamical core <strong>and</strong> in its interface with the model physics. As is well known in upperatmospheric<br />

modeling, temporal <strong>and</strong> spatial variations of specific heat along air parcel trajectories have to be<br />

explicitly accounted for in the energy equation to represent correctly the effects of a highly variable thermospheric<br />

composition. These variations are commonly not accounted for in lower-atmospheric GCMs, <strong>and</strong> instead rely on the<br />

approximation of ―virtual temperature‖ <strong>and</strong> assume that the composition deviations from the uniform ―dry air,‖<br />

primarily caused by the presence of water vapor, are small. A new spectral model formulation has been developed<br />

<strong>and</strong> implemented in the GFS, using specific enthalpy cpT, where cp is specific heat at constant pressure <strong>and</strong> T is<br />

temperature, as a prognostic variable instead of the traditional virtual temperature.<br />

Molecular dissipative processes such as viscosity, heat conduction, <strong>and</strong> diffusion become so fast in the upper<br />

thermosphere that they can no longer be treated within the st<strong>and</strong>ard vertical-column interface. Horizontal molecular<br />

transport of momentum, heat, <strong>and</strong> constituents along pressure surfaces have been incorporated into the model as<br />

well. This physical dissipation eliminates the need for excessive numerical damping commonly used to stabilize<br />

numerical models in vertically extended domains. Additional physical processes incorporated in the extended model<br />

domain include UV <strong>and</strong> EUV radiative heating, infrared radiative cooling with the breakdown of local<br />

thermodynamic equilibrium, <strong>and</strong> non-orographic gravity waves. In the normal fully-coupled configuration ion drag,<br />

Joule <strong>and</strong> particle heating are calculated self-consistently within GIP. In a st<strong>and</strong>-alone WAM, these processes can be<br />

parameterized using an empirical ionosphere model.<br />

Global Ionosphere Plasmasphere (GIP) model <strong>and</strong> electrodynamics. The plasma processes required in IDEA are<br />

provided by the Global Ionosphere<br />

Plasmasphere (GIP) module, which includes<br />

a self-consistent global electrodynamic<br />

solver. GIP was developed by extracting the<br />

plasma processes embedded in the Coupled<br />

Thermosphere Ionosphere Plasmasphere<br />

electrodynamics (CTIPe) model. In earlier<br />

versions of CTIPe, the Earth‘s magnetic<br />

field was assumed to be a tilted dipole. In<br />

reality, however, the low-latitude<br />

geomagnetic field departs from the dipole<br />

equator by more than ten degrees over the<br />

Atlantic Ocean <strong>and</strong> Africa. Furthermore,<br />

ionospheric dynamo electric fields depend<br />

on the strength of the magnetic field, since<br />

the wind-driven currents depend on the<br />

field-line integration of the ionospheric<br />

conductivity <strong>and</strong> neutral wind.<br />

Figure 2. Wavelet amplitude spectrum of zonal wavenumber 2<br />

oscillation in the WAM zonal wind at 20 S <strong>and</strong> 110 km altitude as a<br />

function of period <strong>and</strong> day number in March <strong>and</strong> April.<br />

GIP now explicitly includes the International Geomagnetic Reference Field (IGRF), an accurate representation of<br />

the Earth‘s magnetic field, in order to introduce true longitude dependence. The global ionospheric dynamics <strong>and</strong><br />

energetics are solved in the magnetic APEX coordinate system. GIP receives the neutral atmosphere fields from<br />

WAM, computes the Joule heating, ion-drag, <strong>and</strong> auroral particle heating <strong>and</strong> passes these parameters back to<br />

WAM. The neutral <strong>and</strong> plasma domains are completely self-consistent <strong>and</strong> interactive.<br />

Self-consistent coupling between the WAM <strong>and</strong> GIP is completed by solving the global electrodynamics <strong>and</strong> by<br />

feeding the electric fields back into GIP. The electric fields are computed using a global potential solver in the same<br />

APEX coordinate system. The electrodynamics module receives the dynamics from WAM <strong>and</strong> the field-aligned<br />

conductivities from GIP, computes the electric fields, <strong>and</strong> provides the information back to both modules. The<br />

resultant ionospheric dynamo electric fields reflect the forcing from the lower atmosphere <strong>and</strong> realistic ionospheric<br />

structure through the use of the IGRF.<br />

Early results from WAM provide the first confirmation that the model is capable of realistically reproducing the<br />

dynamics <strong>and</strong> variability in the upper atmosphere, <strong>and</strong> that terrestrial weather features can affect the thermosphere-<br />

39

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!