Observations of narrow bipolar events in East China

Narrow bipolar events (NBEs) are a distinct class of intra-cloud lightning discharge. In this paper we present observations of 10 negative and 67 positive such events in East China. Positive NBEs occurred at 7–12 km altitude above mean sea level (MSL) with a mean altitude of 9.5 km, and negative NBEs occurred at 14–16 km altitude. Electrical/channel characteristics of these events were derived from NBE pulse waveforms based on the transmission-line model. On average, the peak current moment and the charge moment change of a NBE event is 15 kA km, and 0.12 C km, respectively. The mean time for the propagation of current front along the channel is 2.2 μs. The upper limit on channel length for NBEs in this study is 510–1060 m, the lower limit on discharge current amplitude is 12.5–43.2 kA, and the minimum charge transfer is 0.1–0.3 C.     

Modeling the responses of the middle latitude ionosphere to solar flares

In this paper, we investigate the solar flare effects of the ionosphere at middle latitude with a one-dimensional ionosphere theoretical model. The measurements of solar irradiance from the SOHO/Solar EUV Monitor (SEM) and GOES satellites have been used to construct a simple time-dependent solar flare spectrum model, which serves as the irradiance spectrum during solar flares. The model calculations show that the ionospheric responses to solar flares are largely related to the solar zenith angle. During the daytime most of the relative increases in electron density occur at an altitude lower than 300 km, with a peak at about 115 km, whereas around sunrise and sunset the strongest ionospheric responses occur at much higher altitudes (e.g. 210 km for a summer flare). The ionospheric responses to flares in equinox and winter show an obvious asymmetry to local midday with a relative increase in total electron content (TEC) in the morning larger than that in the afternoon. The flare-induced TEC enhancement increases slowly around sunrise and reaches a peak at about 60 min after the flare onset.     

Variability of equatorial ionospheric electron density at fixed heights below the F2 peak

A study on variability of the equatorial ionosphere was carried out at fixed heights below the F2 peak for two different levels of solar activity. The study covered height range of 100 km up to the peak of F2 layer using a real height step increase of 10 km. The variability index used is the percentage ratio of standard deviation over the average value for the month. Daytime minimum variability of between 3% and 10% was observed at height range of about 150–210 km during low solar activity and between 2% and 7% at height range of 160–220 km during high solar activity. The nighttime maximum of between 70% and 187% was observed at height range of about 210–250 km during low solar activity and between 42% and 127% at height range of 210–250 km during high solar activity. The height range at which daytime minimum was observed falls within the F1 height of the ionosphere. The result obtained is consistent with previous works carried out in the low latitude locations for American sector.     

Atmospheric temperature profiling by joint Raman,Rayleigh and Fe Boltzmann lidar measurements

Simultaneous and complete temperature profiles from near ground to about 100 km are essential for studying the dynamical coupling between different atmospheric layers. They are acquired by combining three different lidar techniques at Wuhan, China (30.5°N, 114.4°E). The atmospheric temperatures from about 3 to 25 km are calculated from the nitrogen molecule density profiles obtained from the N₂ vibrational Raman backscatter, while the atmospheric temperatures between 30 and ∼75 km are calculated by the standard Rayleigh scattering method. The temperatures in the 80–100 km altitude region are derived from the Fe Boltzmann technique. The temperature profiles measured by our lidar systems exhibit good agreement when compared with the radiosonde and satellite data, as well as the model. A Lomb–Scargle spectral analysis of the normalized temperature perturbations in the altitude range from 4 to 60 km shows that the spectral slopes of the vertical wave number spectra tended to −3 for large vertical wave numbers. This is consistent with the model predictions of saturated gravity wave spectra.     

The polar wind: Recent observations

The polar wind is an ambipolar outflow of thermal plasma from the high-latitude ionosphere to the magnetosphere, and it primarily consists of H⁺, He⁺ and O⁺ ions and electrons. Statistical and episodic studies based primarily on ion composition observations on the ISIS-2, DE-1, Akebono and Polar satellites over the past four decades have confirmed the existence of the polar wind. These observations spanned the altitude range from 1000 to ∼50,500 km, and revealed several important features in the polar wind that are unexpected from “classical” polar wind theories. These include the day–night asymmetry in polar wind velocity, which is 1.5–2.0 times larger on the dayside; appreciable O⁺ flow at high altitudes, where the velocity at 5000–10,000 km is of 1–4 km/s; and significant electron temperature anisotropy in the sunlit polar wind, in which the upward-to-downward electron temperature ratio is 1.5–2. These features are attributable to a number of “non-classical” polar wind ion acceleration mechanisms resulting from strong ionospheric convection, enhanced electron and ion temperatures, and escaping atmospheric photoelectrons. The observed polar wind has an averaged ion temperature of ∼0.2–0.3 eV, and a rate of ion velocity increase with altitude that correlates strongly with electron temperature and is greatest at low altitudes (<4000 km for H⁺). The rate of velocity increase below 4000 km is larger at solar minimum than at solar maximum. Above 4000 km, the reverse is the case. This suggests that the dominant polar wind ion acceleration process may be different at low and high altitudes, respectively. At a given altitude, the polar wind velocity is highly variable, and is on average largest for H⁺ and smallest for O⁺. Near solar maximum, H⁺, He⁺, and O⁺ ions typically reach a velocity of 1 km/s near 2000, 3000, and 6000 km, respectively, and velocities of 12, 7, and 4 km/s, respectively, at 10,000 km altitude. Near solar minimum, the velocity of all three species is smaller at high altitudes. Observationally it is not always possible to unambiguously separate an energized “non-polar-wind” ion such as a low-energy “cleft ion fountain” ion that has convected into a polar wind flux tube from an energized “polar-wind” ion that is accelerated locally by “non-classical” polar-wind ion acceleration mechanisms. Significant questions remain on the relative contribution between the cleft ion fountain, auroral bulk upflow, and the topside polar-cap ionosphere to the O⁺ polar wind population at high altitudes, the effect of positive spacecraft charging on the lowest-energy component of the H⁺ polar wind population, and the relative importance of the various classical and non-classical ion acceleration mechanisms. These questions pose several challenges in future polar wind observations: These include measurement of the lowest-energy component in the presence of positive spacecraft potential, definitive determination and if possible active control of the spacecraft potential, definitive discrimination between polar wind and other inter-mixed thermal ion populations, measurement of the three-dimensional ion drift velocity vector and the parallel and perpendicular ion temperatures or the detailed three-dimensional velocity distribution function, and resolution of He⁺ and other minor ion species in the polar wind population.     

Global distribution and climatological features of the 5–6-day planetary waves seen in the SABER/TIMED temperatures (2002–2007)

The paper is focused on the global spatial structure, seasonal and interannual variability of the ～5-day Rossby (W1) and ～6-day Kelvin (E1) waves derived from the SABER/TIMED temperature measurements for 6 full years (January 2002–December 2007). The latitude structure of the ～5-day W1 wave is related to the gravest symmetric wave number 1 Rossby wave. The vertical structure of the ～5-day Rossby wave amplitude consists of double-peaked maxima centred at ～80–90 km and ～105–110 km. This wave has a vertically propagating phase structure from the stratosphere up to 120 km altitude with a mean vertical wavelength of ～50–60 km. The ～6-day E1 wave is an equatorially trapped wave symmetric about the equator and located between 20°N and 20°S. Its seasonal behaviour indicates some equinoctial and June solstice amplifications, while the vertical phase structure indicates that this is a vertically propagating wave between 20–100 km altitudes with a mean vertical wavelength of ～25 km.     

Statistical wind analysis for near-space applications

Statistical wind models were developed based on the existing observational wind data for near-space altitudes between 60 000 and 100 000 ft (18–30 km) above ground level (AGL) at two locations, Akon, OH, USA, and White Sands, NM, USA. These two sites are envisioned as playing a crucial role in the first flights of high-altitude airships. The analysis shown in this paper has not been previously applied to this region of the stratosphere for such an application. Standard statistics were compiled for these data such as mean, median, maximum wind speed, and standard deviation, and the data were modeled with Weibull distributions. These statistics indicated, on a yearly average, there is a lull or a “knee” in the wind between 65 000 and 72 000 ft AGL (20–22 km). From the standard statistics, trends at both locations indicated substantial seasonal variation in the mean wind speed at these heights. The yearly and monthly statistical modeling indicated that Weibull distributions were a reasonable model for the data. Forecasts and hindcasts were done by using a Weibull model based on 2004 data and comparing the model with the 2003 and 2005 data. The 2004 distribution was also a reasonable model for these years. Lastly, the Weibull distribution and cumulative function were used to predict the 50%, 95%, and 99% winds, which are directly related to the expected power requirements of a near-space station-keeping airship. These values indicated that using only the standard deviation of the mean may underestimate the operational conditions.     

A methodology for using optimal MSIS parameters retrieved from SSULI data to compute satellite drag on LEO objects

A key application to be derived from Space Weather research will be to forecast atmospheric drag on low Earth orbit (LEO) satellites with significantly better accuracy than is attainable today. The recently launched STP P91-1 ARGOS mission will serve as a testbed for the use of future operational ultraviolet remote-sensing data to achieve such an improvement. This paper describes the associated methodology, which uses discrete inverse theory in conjunction with the data to derive correction factors in near real time for the MSISE-90 empirical thermospheric model. To simulate the application of this technique to orbit prediction, we use the Jacchia-71 operational model to generate an evolving “ground truth” upper atmospheric state over a 48 h time period. This permits a state-of-the-art Satellite Tool Kit orbit propagator to synthesize a corresponding “ground truth” orbit on a standard LEO test object at 350 km altitude. Our tests show that, for orbit prediction, the “data-enhanced” MSIS density specification can provide significant improvement over the uncorrected MSIS specification. However, for orbit prediction, the results are sensitive to the strategy selected for applying the correction factors. We contrast our results for orbit prediction with those of Marcos et al. (1998. Astrodynamics. Vol. 97(1). AAS, San Diego, pp. 501–513) for precision orbit determination. An important result in the context of Space Weather is that the Jacchia and MSIS models can show significant point-to-point disagreement, which has major implications for operational specification of thermospheric drag.     

Crustal structure of the Reykjanes Ridge near 62°N,on the basis of seismic refraction and gravity data

Explosion deep seismic sounding data sections of high quality had been obtained with RV Meteor in the Reykjanes Iceland Seismic Project (RRISP77 [Angenheister, G., Gebrande, H., Miller, H., Goldflam, P., Weigel, W., Jacoby, W.R., Pálmason, G., Björnsson, S., Einarsson, P., Pavlenkova, N.I., Zverev, S., Litvinenko, I.V., Loncarecic, B., Solomon, S., 1980. Reykjanes Ridge Iceland Seismic Experiment (RRISP 77). J. Geophys. 47, 228–238]) which close an information gap near 62°N. Preliminary results were presented by Weigel [Weigel, W., 1980. Aufbau des Reykjanes Rückens nach refraktionsseismischen Messungen. In: Weigel, W. (Ed.), Reykjanes Rücken, Island, Norwegischer Kontinentalrand. Abschlusskolloquium, Hamburg zur Meteor-Expedition, vol. 45. DFG, Bonn, pp. 53–61], and here we report on the data and results of interpretation. Clear refracted phases to 90 km distance permit crustal and uppermost mantle structure to be modelled by ray tracing. The apparent P-wave velocities are around 4.5, 6–6.5, 7–7.6 and 8.2–8.7 km/s, but no wide-angle reflections have been clearly seen. Accompanying sparker reflection data reveal thin sediment ponds in the axial zone and up to 400 m thick sediments at 10 Ma crustal age. Ray tracing reveals the following model below the sediments: (1) a distinct, 1–2 km thick upper crust (layer 2A) with Vp increasing with age (to 10 Ma) from <3.4 to 4.9 km/s and with a vertical gradient of 0.1–0.2 km/s/km, (2) a lower crust or layer 3 beginning at depths of 2 (axis) to 4 km (10 Ma age) below sea level with 6.1–6.8 km/s and similar vertical gradients as above, (3) the lower crust bottoms at 5.2–9.5 km depth below sea level (0–10 Ma) with a marked discontinuity, underneath which (4) Vp rises from about 7.5–7.8 km/s (0–10 Ma) with a positive vertical gradient of, again, 0.1–0.2 km/s/km such that 8 km/s would be reached at 12 km and deeper near the axis. Our preferred interpretation is that the mantle begins at the distinct discontinuity (“Moho”), but a deeper “Moho” of Vp ≈ 8 km/s cannot be excluded. From Iceland southward to 60°N several experiments show a decrease of crustal thickness from 14 to 8 km. Velocity trends with age across the ridge reflect cooling and filling of cracks, and thickness trends probably suggest volcanic productivity variations as previously suggested.Gravity inversion concentrates on a profile across the ridge with the above seismic a priori information; with 0.2–0.5 km depth uncertainty it leads to a good fit (±2.5 mGal where seismic data exist). Best fitting densities are (in kg/m³) for sediments, 2180; upper crust, 2450–2570; lower crust, 2850–2940; mantle lithosphere, 3215–3240 with a deficit for an asthenospheric wedge of no more than −100 kg/m³. The morphological ridges and troughs superimposed on the SE ridge flank are partly correlated, partly anti-correlated with the Bouguer anomaly and suggest that variable crustal density variations accompany the morphology variations.     

Validation of temperature measurements from MIPAS-ENVISAT with balloon observations obtained by MIPAS-B

The Michelson Interferometer for Passive Atmospheric Sounding (MIPAS), a limb sounding Fourier transform interferometer in the mid-infrared band, on board the polar-orbiting ENVIronmental SATellite (ENVISAT) was launched in March 2002 by the European Space Agency (ESA). For the MIPAS data validation, three balloon flights with MIPAS-B, a balloon version of MIPAS, were carried out on the 24th and 25th of September 2002 from Aire sur l’Adour (France 44°N, 0°E), on the 20th and 21th of March, and on the 2nd and 3rd of July 2003 from Esrange, Kiruna (Sweden 68°N, 21°E). The MIPAS operational data version 4.61 for the temperature vertical profiles were compared with the correlative MIPAS-B measurements via the coincident comparison and the trajectory comparison approaches, respectively. The precision of the MIPAS temperature was estimated to be 1.04–2.48 K in the region 123–8 hPa (15–33 km) and beyond the expectation. The systematic difference shows that the absolute accuracy of the MIPAS temperature is within the total combined errors between 356 and 3 hPa (8–39 km). The agreements in the middle stratosphere are better than in the lower stratosphere and upper troposphere. The maximum difference of 1.9 K (0.9%) was found at 228 hPa (11 km). The biases of MIPAS measurements are ?0.39 K (?0.19%) and +0.25 K (+0.12%) (“+” positive bias, “?” negative bias) corresponding to the altitude regions 143–31 hPa (14–24 km) and 31–5 hPa (24–37 km), respectively.     

Imaging Pn velocities beneath the Pannonian basin

The Pannonian depression is an extensional back-arc basin in central Europe and is an integral part of the Alpine–Carpathian orogenic mountain belts. It can be characterized by thinned lower crust, shallow Moho discontinuity, high surface heat flow and Moho temperature, implying recent active tectonic processes. Imaging the velocity structure of the upper mantle may help us to better understand the structure and formation of the Pannonian region.In this paper, P_n traveltimes from regional earthquakes are used to tomographically image the lateral velocity variations in the uppermost mantle beneath the Pannonian basin. The set of linear tomographic equations, built up of the time term equation for each source–receiver pair, is solved by a truncated singular value decomposition algorithm. The explicit computation of the generalized inverse of the tomographic equations makes it possible to deduce both the resolution matrix and the model covariance matrix, allowing us to estimate the resolution and reliability of the solution.The mean compressional wave velocity in the uppermost mantle beneath the Pannonian basin is 7.9 km/s, substantially lower than the average continental P_n velocity of 8.1 km/s. It is mostly due to the high Moho temperature having values on average 400–500 °C more than those in the surrounding areas. The velocity anomalies range from −0.3 to 0.3 km/s relative to the mean velocity of 7.9 km/s. Due to high Moho temperature, below the North Hungarian range low (7.6–7.7 km/s) velocities can be found. High-velocity anomalies of around 8.1 km/s can be detected along the W-SW boundaries of Hungary and at the junction of the Pannonian basin and the Southern Carpathians. The Great Hungarian Plain shows average (7.9 km/s) P_n velocities.     

Solar eclipse of the 29 March 2006: Results of the optical measurements by MORTI over Almaty (43.03°N, 76.58°E)

The observations of the upper mesosphere region (∼95 km altitude) in the period of 27–30 March 2006 using mesopause oxygen rotational temperature imager (MORTI) at Almaty (43.03°N, 76.58°E) are presented in this report to illustrate the mesosphere response to the solar eclipse (SE) event, which occurred on 29 March 2006. The nighttime volume emission rates and rotational temperatures, obtained from MORTI measurements, show appreciable differences in the pattern of wave-like oscillations observed during the period of interest. These oscillations are possibly due to the SE. Using a periodogram method the spectra of the observed wave-like oscillations, observed in the mesosphere, are examined. A physical mechanism is proposed to interpret the effects observed in terms of the mesosphere response to the total SE.     

On the sensitivity of total electron content (TEC) to upper atmospheric/ionospheric parameters

The total electron content (TEC) is a key ionospheric parameter for various space weather applications. Over the last decade an extensive database of TEC measurements has become available from both space- and ground-based observations, and these measurements have established the general morphology of the global TEC distributions. In particular, the TOPEX TEC measurements have shown strong longitudinal variations of TEC in addition to the observed day-to-day variabilities. To better understand the observed TEC variations and to better guide its modeling, we have studied the sensitivity of quiet-time TEC to the following key atmospheric and ionospheric parameters: neutral density, neutral wind, plasma temperatures, plasmaspheric flux, and the O⁺–O collision frequency. These parameters are often only roughly known and can cause large uncertainties in model results. For this study, we have developed a numerical mid-latitude ionospheric model, which solves the momentum and continuity equations for the O⁺ density and a simplified set of equations for the H⁺ density. To obtain TEC, the calculated ion densities have been integrated from the bottom altitude (100 km) to the altitude of the TOPEX satellite (1336 km). Our study shows that during the day the neutral wind and the neutral composition have the most important effect on TEC. In particular, the zonal component of the neutral wind can have a large effect on TEC in the southern hemisphere where the magnetic declination angle is large. During the night, most of the above-mentioned parameters can play a significant role in the TEC morphology, except for the plasma temperature, which has only a small effect on TEC. Finally, the TEC varies roughly linearly with respect to all of the parameters except for the neutral wind.     

Hybrid model for nowcasting and forecasting the K index

A new empirical model nowcasting and predicting a proxy to the geomagnetic K index is developed, which is based on the combined use of solar wind parameters and ground-based magnetic data. The present approach implements the previously developed solar wind-based MAK model, calibrating its values with magnetogram-derived K index. The new model is named as Hybrid Dourbes K (HDK) model. The HDK nowcast model provides the quantity Kdf, obtained by solar wind-based Ksw and corrected with a combination of differences between several past values of Kd and Ksw. The model error of the nowcast Kdf is found to be 0.38 KU, or nearly twice less than that of the MAK model. Kdf has a good predictability. Prediction made by weighted extrapolation 6 h ahead carries an error of 1.0 KU, while for the first 1 h the error is 0.58 KU only.     

Effect of photoelectric emission on blunt probe conductivity measurements in the stratosphere

Two identical planar blunt probes of stainless steel material, biased with a bipolar ramp voltage, are used to measure the stratospheric polar conductivities to altitudes of 34 km. One probe (DP) is mounted closer to the gondola, looking downwards and shielded from sunlight, while the other (SP) is mounted looking sideways, away from the gondola. The daytime observations of positive ions in the 29–34 km altitude range with SP, and of negative ions at 34 km with DP, show photoelectric contaminations induced by solar UV radiations in the 190–230 nm band. These contaminations are found to be due to photoemissions from the SP probe steel surface and from the carbon paint that coats the surface of the gondola, respectively. It is found that, a segment of the photocurrent contaminated I–V curve, recorded with SP at higher negative probe potentials, is linear, and it can give the ambient positive polar conductivity.     

The statistical and numerical study of the global distribution of coronal plasma and magnetic field near 2.5 Rs over a 10-year period

The basic characteristics of the global distribution for the corona plasma and magnetic field near 2.5 Rs are analyzed with the statistical and numerical methods for 136 Carrington Rotations (CRs) covering four different phases of solar activity. By using the observational data and the velocity distribution model in the corona, the statistical average distribution of the magnetic field, density and the coronal mass outputs are analyzed for the four different phases. Then, a numerical study of the global distribution near 2.5 Rs has been made by solving a self-consistent MHD system. Finally, the solar wind speed at 1 AU is given by mapping the speed at 2.5 Rs to that near 1 AU, and the comparison of the numerical results with the observational measurements and the simulation result of the Wang–Sheeley–Arge (WSA) model are made during more than 5 years. The numerical results indicate that the global distributions on the source surface of 2.5 Rs at different phases of solar activity could be used to predict the change of the solar wind in interplanetary space.     

Crustal structure variation along 30°N in the eastern Tibetan Plateau and its tectonic implications

We determined crustal structure along the latitude 30°N through the eastern Tibetan Plateau using a teleseismic receiver function analysis. The data came mostly from seismic stations deployed in eastern Tibet and western Sichuan region from 2004 to 2006. Crustal thickness and Vp/Vs ratio at each station were estimated by the H–k stacking method. On the profile, the mean crustal thickness and Vp/Vs ratio were found to be 62.3 km and 1.74 in the Lhasa block, 71.2 km and 1.79 near the Bangong–Nujiang suture, 66.3 km and 1.80 in the Qiangtang block, 59.8 km and 1.81 in the Songpan–Garze block, and 42.9 km and 1.76 in the Yangtze block, respectively. The estimated crustal thicknesses are consistent with predictions based on the topography and the Airy isostasy, except near the Bangong–Nujiang suture and in the Qiangtang block where the crust is 5–10 km thicker than predicted, indicating that the crust may be denser, possibly due to mafic underplating. We also inverted receiver functions for crustal velocity structure along the profile, which reveals a low S-wave velocity zone in the lower crust beneath the eastern Tibetan Plateau, although the extent of the low-velocity zone varies considerably. The low-velocity zone, together with previous results, suggests limited partial melting and localized crustal flow in the lower crust of the eastern Tibetan Plateau.     

Nighttime E region plasma irregularities over an equatorial station Trivandrum

In situ measurements of electron density were made over Trivandrum (8.5°N, 76.9°E) during nighttime to study E-region plasma density irregularities. Irregularities, with vertical scale sizes from a few km to 15 cm, were detected during rocket ascent and descent. Electron density profiles during ascent and descent of an earlier nighttime rocket flight from Trivandrum are also presented. Some of the important results are as follows: (i) horizontal gradients in electron density exist in 110–120 km region with horizontal scale size of at least 40 km, (ii) based on the presence/absence of electron density structures during ascent and descent of both flights, the horizontal distance over which the gradient drift instability operates is found to be at least 80 km and 90 km, for both the flights, (iii) observed irregularities in regions of negative density gradient are suggested to be produced through the gradient drift instability (GDI) driven by vertical polarization electric field as well as by electric field produced through wind shears and those in positive gradient regions by wind driven GDI, (iv) largest irregularity amplitude (≈30%) was associated with steepest gradients and so was the presence of smallest vertical scale sizes (12 m to 15 cm), which were absent at other altitudes, (v) the spectral index of irregularities was in the range of ?2.2±0.2 for large scales (few kilometers>λ>50 m), ?3.25±0.25 for medium scales (50 m>λ>10 m) and ?2.6±0.1 for smaller scales (10 m>λ>1 m) and (vi) irregularities in large and medium scales are expected to be produced directly through GDI and the small and sum-meter scales through non-linear GDI.     

Deep pooling of low degree melts and volatile fluxes at the 85°E segment of the Gakkel Ridge: Evidence from olivine-hosted melt inclusions and glasses

We present new analyses of volatile, major, and trace elements for a suite of glasses and melt inclusions from the 85°E segment of the ultra-slow spreading Gakkel Ridge. Samples from this segment include limu o pele and glass shards, proposed to result from CO₂-driven explosive activity. The major element and volatile compositions of the melt inclusions are more variable and consistently more primitive than the glass data. CO₂ contents in the melt inclusions extend to higher values (167–1596 ppm) than in the co-existing glasses (187–227 ppm), indicating that the melt inclusions were trapped at greater depths. These melt inclusions record the highest CO₂ melt concentrations observed for a ridge environment. Based on a vapor saturation model, we estimate that the melt inclusions were trapped between seafloor depths (～ 4 km) and ～ 9 km below the seafloor. However, the glasses are all in equilibrium with their eruption depths, which is inconsistent with the rapid magma ascent rates expected for explosive activity. Melting conditions inferred from thermobarometry suggest relatively deep (25–40 km) and cold (1240°–1325 °C) melting conditions, consistent with a thermal structure calculated for the Gakkel Ridge. The water contents and trace element compositions of the melt inclusions and glasses are remarkably homogeneous; this is an unexpected result for ultra-slow spreading ridges, where magma mixing is generally thought to be less efficient based on the assumption that steady-state crustal magma chambers are absent in these environments. All melts can be described by a single liquid line of descent originating from a pooled melt composition that is consistent with the aggregate melt calculated from a geodynamic model for the Gakkel Ridge. These data suggest a model in which deep, low degree melts are efficiently pooled in the upper mantle (9–20 km depth), after which crystallization commences and continues during ascent and eruption. Based on our melting model and the assumption that CO₂ is perfectly incompatible, we show that the highest CO₂ concentrations of the melt inclusions (～ 1600 ppm) are consistent with the calculated CO₂ concentrations of primary undegassed melts. The highest measured CO₂/Nb ratio (443) of Gakkel Ridge melt inclusions predicts a mantle CO₂ content of 134 ppm and would result in a global ridge flux of 2.0 × 10¹² mol CO₂/yr.     

Ground-based Pc5 ULF wave power: Solar wind speed and MLT dependence

Using over 20 years of ground-based magnetometer data from the CANOPUS/CARISMA magnetometer array, we present a statistical characterisation of Pc5 ultra-low frequency (ULF) power in the 2–10 mHz band as a function of magnetic local time (MLT), L-shell, and solar wind speed. We examine the power across L-shells between 4.2 and 7.9, using data from the PINA, ISLL, GILL and FCHU stations, and demonstrate that there is a significant MLT dependence in both the H- and D-component median 2–10 mHz power during both fast (>500 km/s) and slow (<500 km/s) solar wind speeds. The H-component power consistently dominates over D-component power at all MLTs and during both fast and slow solar wind. At the higher-L stations (L>5.4), there are strong MLT power peaks in the morning and midnight local time sectors; the morning sector dominating midnight during fast solar wind events. At lower L-shells, there is no evidence of the midnight peak and the 2–10 mHz power is more symmetric with respect to MLT except during the fastest solar wind speeds. There is little evidence in the ground-based power of a localised MLT peak in ULF power at dusk, except at the lowest L-shell station, predominantly in the H-component. The median 2–10 mHz power increases with an approximate power law dependence on solar wind speed, at all local times across the L-shell domain studied in both components. The H-component power peaks at the latitude of the GILL station, with significantly lower power at both higher and lower L-shells. Conversely, the D-component power increases monotonically. We believe that this is evidence for 2–10 mHz power accumulating at auroral latitudes in field line resonances. Finally, we discuss how such ULF wave power characterisation might be used to derive empirical radiation belt radial diffusion coefficients based on, and driven by, the solar wind speed dependence of ULF wave power.