Programming code for article «Observing electrical fields and neutral winds with EISCAT 3D»
Measurements of height dependent electric field (E) and neutral wind (u) are important governing parameters of the Earth’s upper atmosphere, which can be used to study e.g., how auroral currents close, or how energy flows between the ionized and neutral constituents. The new EISCAT 3D (E3D) incoherent scatter radar will be able to measure a three-dimensional ion velocity vector (v) at each measurement point, which will allow less stringent prior assumptions about E and u to be made when estimating them from radar measurements. This study investigates the feasibility of estimating the three-dimensional electric field and neutral wind vectors along a magnetic field-aligned altitude profile from E3D measurements, using the ion momentum equation and Maxwell’s equations. The uncertainty of ion drift measurements is estimated for a time and height resolution of 5 s and 2 km. With the most favourable ionospheric conditions, the ion wind at E region peak can be measured with an accuracy of less than 1 m/s. In the worst case, during a geomagnetically quiet night, the uncertainty increases by a factor of around ten. The uncertainty of neutral wind and electric field estimates is ound to be strongly dependent on the prior constraints imposed on them. In the lower E region, neutral wind estimates have a lower standard deviation than 10 m/s in the most favourable conditions. In such conditions, also the F region electric field can be estimated with uncertainty of about 1 mV/m. Simulated measurements of v are used to demonstrate the ability to resolve the field-aligned profile of E and u. However, they can only be determined well at the heights where they significantly influence the ion drift, that is above 125 km for E and below 115 km for u. At the other heights, the results are strongly dependent on the the prior assumptions of smoothness.
Julia, 1.6 (versions down to 1.0 should work), see also readme
Python, 3.8 (versions down to 3.0 should work), see also readme
