Publications

The equatorial jets dominating the dynamics of the Jovian planets exhibit two distinct types of zonal flows: strongly eastward in the gas giants (superrotation) and strongly westward in the ice giants (subrotation). Existing theories propose different mechanisms for these patterns, but no single mechanism has successfully explained both. However, the planetary parameters of the four Solar System giant planets suggest that a fundamentally different mechanism is unlikely. In this study, we show that convection-driven columnar structures can account for both eastward and westward equatorial jets, framing the phenomenon as a bifurcation. Consequently, both superrotation and subrotation emerge as stable branches of the same mechanistic solution. Our analysis of these solutions uncovers similarities in the properties of equatorial waves and the leading-order momentum balance. This study suggests that the fundamental dynamics governing equatorial jet formation may be more broadly applicable across the Jovian planets than previously believed, offering a unified explanation for their two distinct zonal wind patterns. A bifurcation in equatorial jet formation can explain the differences between the gas and ice giants.

Saturn is characterized by strong zonal winds at the cloud level, with a wide equatorial west-to-east flow, two retrograde jets on its flanks, and a series of mid-to-high latitude jets in both hemispheres. Based on the Cassini mission gravity measurement up to zonal harmonic J₁₀, it was shown that these winds penetrate to approximately 9000 km. However, these analyses relied fully on the low, even harmonics dominated by an interior structure that is not known to high accuracy. Here, we present a Cassini-based gravity analysis, resolving the gravity harmonics up to J₂₀, by constraining the surface gravity in the polar regions. Using the extended harmonics and surface gravity, we show evidence that the wind structure poleward of latitude 45° must be shallower than 3000 km, while equatorward, the wind must extend along cylinders and decayed radially at a depth of approximately 11,000 km. We also find that nearly doubling the wind strength, while keeping its latitudinal structure, allows for an excellent fit to all 20 measured gravity harmonics, suggesting that the zonal jets reach their maximum strength below the observed cloud level.

Context. NASA's Juno mission provided exquisite measurements of Jupiter's gravity field that together with the Galileo entry probe atmospheric measurements constrains the interior structure of the giant planet. Inferring its interior structure range remains a challenging inverse problem requiring a computationally intensive search of combinations of various planetary properties, such as the cloud-level temperature, composition, and core features, requiring the computation of 109 interior models. Aims. We propose an efficient deep neural network (DNN) model to generate high-precision wide-ranged interior models based on the very accurate but computationally demanding concentric MacLaurin spheroid (CMS) method. Methods. We trained a sharing-based DNN with a large set of CMS results for a four-layer interior model of Jupiter, including a dilute core, to accurately predict the gravity moments and mass, given a combination of interior features. We evaluated the performance of the trained DNN (NeuralCMS) to inspect its predictive limitations. Results. NeuralCMS shows very good performance in predicting the gravity moments, with errors comparable with the uncertainty due to differential rotation, and a very accurate mass prediction. This allowed us to perform a broad parameter space search by computing only 104 actual CMS interior models, resulting in a large sample of plausible interior structures, and reducing the computation time by a factor of 105. Moreover, we used a DNN explainability algorithm to analyze the impact of the parameters setting the interior model on the predicted observables, providing information on their nonlinear relation.

Jupiter's equatorial eastward zonal flows reach wind velocities of ∼100 m s⁻¹, while on Saturn they are three times as strong and extend about twice as wide in latitude, despite the two planets being overall dynamically similar. Recent gravity measurements obtained by the Juno and Cassini spacecraft uncovered that the depth of zonal flows on Saturn is about three times greater than on Jupiter. Here we show, using 3D deep convection simulations, that the atmospheric depth is the determining factor controlling both the strength and latitudinal extent of the equatorial zonal flows, consistent with the measurements for both planets. We show that the atmospheric depth is proportional to the convectively driven eddy momentum flux, which controls the strength of the zonal flows. These results provide a mechanistic explanation for the observed differences in the equatorial regions of Jupiter and Saturn, and offer new understandings about the dynamics of gas giants beyond the Solar System.

The atmospheric dynamics of Jupiter are dominated by strong zonal winds engulfing the planet. Since the first gravity measurements taken by Juno at Jupiter, the low-degree gravity harmonics (J ₃J ₁₀) have been used to determine the depth and structure of the zonal winds observed at the cloud level, limiting inferences on the deep flows to the wide latitudinal structure of these harmonics. Here, using constraints on the dynamical contribution to gravity at high latitude, we present the gravity harmonics up to J _40. We find an excellent correlation between these measurements and the gravity harmonics resulting from the observed cloud-level winds extending inwards cylindrically to depths of ~10⁵ bar (3,000 km). These measurements provide direct evidence that the flows penetrate inwards along the direction of the spin axis, confirming the cylindrical nature of the flow, which has been postulated theoretically since the 1970s. Furthermore, this detailed new gravity spectrum allows us to quantify the contribution of the various jets to the gravity signal, showing the dominance of the strong zonal flows around 20° latitude in both hemispheres.

Context. The Juno mission has provided measurements of Jupiter's gravity field with an outstanding level of accuracy, leading to better constraints on the interior of the planet. Improving our knowledge of the internal structure of Jupiter is key to understanding its formation and evolution but is also important in the framework of exoplanet exploration. Aims. In this study, we investigated the differences between the state-of-the-art equations of state and their impact on the properties of interior models. Accounting for uncertainty on the hydrogen and helium equation of state, we assessed the span of the interior features of Jupiter. Methods. We carried out an extensive exploration of the parameter space and studied a wide range of interior models using Markov chain Monte Carlo simulations. To consider the uncertainty on the equation of state, we allowed for modifications of the equation of state in our calculations. Results. Our models harbour a dilute core and indicate that Jupiter's internal entropy is higher than what is usually assumed from the Galileo probe measurements. We obtain solutions with extended dilute cores, but contrary to other recent interior models of Jupiter, we also obtain models with small dilute cores. The dilute cores in such solutions extend to ∼20% of Jupiter's mass, leading to better agreement with formation-evolution models. Conclusions. We conclude that the equations of state used in Jupiter models have a crucial effect on the inferred structure and composition. Further explorations of the behaviour of hydrogen-helium mixtures at the pressure and temperature conditions in Jupiter will help to constrain the interior of the planet, and therefore its origin.

The Juno spacecraft has been collecting data to shed light on the planets origin and characterize its interior structure. The onboard gravity science experiment based on X-band and Ka-band dual-frequency Doppler tracking precisely measured Jupiters zonal gravitational field. Here, we analyze 22 Junos gravity passes to investigate the gravity field. Our analysis provides evidence of new gravity field features, which perturb its otherwise axially symmetric structure with a time-variable component. We show that normal modes of the planet could explain the anomalous signatures present in the Doppler data better than other alternative explanations, such as localized density anomalies and non-axisymmetric components of the static gravity field. We explain Juno data by p-modes having an amplitude spectrum with a peak radial velocity of 1050 cm/s at 9001200 μHz (compatible with ground-based observations) and provide upper bounds on lower frequency f-modes (radial velocity smaller than 1 cm/s). The new Juno results could open the possibility of exploring the interior structure of the gas giants through measurements of the time-variable gravity or with onboard instrumentation devoted to the observation of normal modes, which could drive spacecraft operations of future missions.

Context. While Jupiter's massive gas envelope consists mainly of hydrogen and helium, the key to understanding Jupiter's formation and evolution lies in the distribution of the remaining (heavy) elements. Before the Juno mission, the lack of high-precision gravity harmonics precluded the use of statistical analyses in a robust determination of the heavy-element distribution in Jupiter's envelope. Aims. In this paper, we assemble the most comprehensive and diverse collection of Jupiter interior models to date and use it to study the distribution of heavy elements in the planet's envelope. Methods. We apply a Bayesian statistical approach to our interior model calculations, reproducing the Juno gravitational and atmospheric measurements and constraints from the deep zonal flows. Results. Our results show that the gravity constraints lead to a deep entropy of Jupiter corresponding to a 1 bar temperature that is 515 K higher than traditionally assumed. We also find that uncertainties in the equation of state are crucial when determining the amount of heavy elements in Jupiter's interior. Our models put an upper limit to the inner compact core of Jupiter of 7 MEarth, independently of the structure model (with or without a dilute core) and the equation of state considered. Furthermore, we robustly demonstrate that Jupiter's envelope is inhomogeneous, with a heavy-element enrichment in the interior relative to the outer envelope. This implies that heavy-element enrichment continued through the gas accretion phase, with important implications for the formation of giant planets in our Solar System and beyond.

Jupiter's atmosphere is dominated by multiple jet streams which are strongly tied to its 3D atmospheric circulation. Lacking a rigid bottom boundary, several models exist for how the meridional circulation extends into the planetary interior. Here, we show, collecting evidence from multiple instruments of the Juno mission, the existence of midlatitudinal meridional circulation cells which are driven by turbulence, similar to the Ferrel cells on Earth. Different than Earth, which contains only one such cell in each hemisphere, the larger, faster rotating Jupiter can incorporate multiple cells. The cells form regions of upwelling and downwelling, which we show are clearly evident in Juno's microwave data between latitudes (Formula presented.) and (Formula presented.). The existence of these cells is confirmed by reproducing the ammonia observations using a simplistic model. This study solves a long-standing puzzle regarding the nature of Jupiter's subcloud dynamics and provides evidence for eight cells in each Jovian hemisphere.

Jupiters atmosphere has a system of zones and belts punctuated by small and large vortices, the largest being the Great Red Spot. How these features change with depth is unknown, with theories of their structure ranging from shallow meteorological features to surface expressions of deep-seated convection. We present observations of atmospheric vortices using the Juno spacecrafts Microwave Radiometer. We found vortex roots that extend deeper than the altitude at which water is expected to condense, and we identified density inversion layers. Our results constrain the three-dimensional structure of Jupiters vortices and their extension below the clouds.

Juno microwave radiometer (MWR) observations of Jupiter's midlatitudes reveal a strong correlation between brightness temperature contrasts and zonal winds, confirming that the banded structure extends throughout the troposphere. However, the microwave brightness gradient is observed to change sign with depth: the belts are microwave-bright in the p<5 bar range and microwave-dark in the p>10 bar range. The transition level (which we call the "jovicline") is evident in the MWR 11.5 cm channel, which samples the 5-14 bar range when using the limb-darkening at all emission angles. The transition is located between 4 and 10 bars, and implies that belts change with depth from being NH3-depleted to NH3-enriched, or from physically warm to physically cool, or more likely a combination of both. The change in character occurs near the statically stable layer associated with water condensation. The implications of the transition are discussed in terms of ammonia redistribution via meridional circulation cells with opposing flows above and below the water condensation layer, and in terms of the "mushball" precipitation model, which predicts steeper vertical ammonia gradients in the belts versus the zones. We show via the moist thermal wind equation that both the temperature and ammonia interpretations can lead to vertical shear on the zonal winds, but the shear is similar to 50x weaker if only NH3 gradients are considered. Conversely, if MWR observations are associated with kinetic temperature gradients then it would produce zonal winds that increase in strength down to the "jovicline", consistent with Galileo probe measurements; then decay slowly at higher pressures.<br />

During the past few years, both the Cassini mission at Saturn and the Juno mission at Jupiter provided measurements with unprecedented accuracy of the gravity and magnetic fields of the two gas giants. Using the gravity measurements, it was found that the strong zonal flows observed at the cloud level of the gas giants are likely to extend thousands of kilometres deep into the planetary interior. However, the gravity measurements alone, which are by definition an integrative measure of mass, cannot constrain with high certainty the exact vertical structure of the flow. Taking into account the recent Cassini magnetic field measurements of Saturn, and past secular variations of Jupiter's magnetic field, we obtain an additional physical constraint on the vertical decay profile of the observed zonal flows on these planets. Our combined gravity-magnetic analysis reveals that the cloud-level winds on Saturn (Jupiter) extend with very little decay, i.e. barotropically, down to a depth of around 7000 km (2000 km) and then decay rapidly in the semiconducting region, so that within the next 1000 km (600 km) their value reduces to about 1 per cent of that at the cloud level. These results indicate that there is no significant mechanism acting to decay the flow in the outer neutral region, and that the interaction with the magnetic field in the semiconducting region might play a central role in the decay of the flows.

The asymmetric gravity field measured by the Juno spacecraft has allowed the estimation of the depth of Jupiter's zonal jets, showing that the winds extend approximately 3,000 km beneath the cloud level. This estimate was based on an analysis using a combination of all measured odd gravity harmonics, J 3, J 5, J 7, and J 9, but the wind profile's dependence on each of them separately has yet to be investigated. Furthermore, these calculations assumed the meridional profile of the cloud-level wind extends to depth. However, it is possible that the interior jet profile varies somewhat from that of the cloud level. Here we analyze in detail the possible meridional and vertical structure of Jupiter's deep jet streams that can match the gravity measurements. We find that each odd gravity harmonic constrains the flow at a different depth, with J 3 the most dominant at depths below 3,000 km, J 5 the most restrictive overall, whereas J 9 does not add any constraint on the flow if the other odd harmonics are considered. Interior flow profiles constructed from perturbations to the cloud-level winds allow a more extensive range of vertical wind profiles, yet when the meridional profiles differ substantially from the cloud level, the ability to match the gravity data significantly diminishes. Overall, we find that while interior wind profiles that do not resemble the cloud level are possible, they are statistically unlikely. Finally, inspired by the Juno microwave radiometer measurements, assuming the brightness temperature is dominated by the ammonia abundance, we find that depth-dependent flow profiles are still compatible with the gravity measurements.

The Hadley circulation is a key element of the climate system. It is traditionally defined as the zonally averaged meridional circulation in the tropics, therefore treated as a zonally symmetric phenomenon. However, differences in temperature between land and sea cause zonal asymmetries on Earth, dramatically affecting the circulation. This longitudinal dependence of the meridional circulation evokes questions about where and when the actual large-scale tropical circulation occurs. Here, we look into the connection between the longitudinally dependent meridional circulation, and the actual large-scale transport of air in the tropics using a coupled Eulerian and Lagrangian approach. Decomposing the velocity field into rotational and divergent components, we identify how each component affects the actual circulation. We propose an alternative definition for the circulation, which describes the actual path of air parcels in the tropics, as a tropical atmospheric conveyor belt.

Jupiter's internal flow structure is still not fully known, but can be now better constrained due to Juno's highprecision measurements. The recently published gravity and magnetic field measurements have led to new information regarding the planet and its internal flows, and future magnetic measurements will help to solve this puzzle. In this study, we propose a new method to better constrain Jupiter's internal flow field using the Juno gravity measurements combined with the expected measurements of magnetic secular variation. Based on a combination of hydrodynamical and magnetic field considerations we show that an optimized vertical profile of the zonal flows that fits both measurements can be obtained. Incorporating the magnetic field effects on the flow better constrains the flow decay profile. This will get us closer to answering the persistent question regarding the depth and nature of the flows on Jupiter.

One of Jupiter's most prominent atmospheric features, the Great Red Spot (GRS), has been observed for more than two centuries, yet little is known about its structure and dynamics below its observed cloud level. While its anticyclonic vortex appearance suggests it might be a shallow weather-layer feature, the very long time span for which it was observed implies it is likely deeply rooted, otherwise it would have been sheared apart by Jupiter's turbulent atmosphere. Determining the GRS depth will shed light not only on the processes governing the GRS, but on the dynamics of Jupiter's atmosphere as a whole. The Juno mission single flyby over the GRS (PJ7) discovered using microwave radiometer measurements that the GRS is at least a couple hundred kilometers deep. The next flybys over the GRS (PJ18 and PJ21), will allow high-precision gravity measurements that can be used to estimate how deep the GRS winds penetrate below the cloud level. Here we propose a novel method to determine the depth of the GRS based on the new gravity measurements and a Slepian function approach that enables an effective representation of the wind-induced spatially confined gravity signal, and an efficient determination of the GRS depth given the limited measurements. We show that with this method the gravity signal of the GRS should be detectable for wind depths deeper than 300 km, with reasonable uncertainties that depend on depth (e.g., +/- 100 km for a GRS depth of 1000 km).

Jupiter's atmosphere is rotating differentially, with zones and belts rotating at speeds that differ by up to 100 metres per second. Whether this is also true of the gas giant's interior has been unknown(1,2), limiting our ability to probe the structure and composition of the planet(3,4). The discovery by the Juno spacecraft that Jupiter's gravity field is north-south asymmetric(5) and the determination of its non-zero odd gravitational harmonics J(3), J(5), J(7) and J(9) demonstrates that the observed zonal cloud flow must persist to a depth of about 3,000 kilometres from the cloud tops(6). Here we report an analysis of Jupiter's even gravitational harmonics J(4), J(6), J(8) and J(10) as observed by Juno(5) and compared to the predictions of interior models. We find that the deep interior of the planet rotates nearly as a rigid body, with differential rotation decreasing by at least an order of magnitude compared to the atmosphere. Moreover, we find that the atmospheric zonal flow extends to more than 2,000 kilometres and to less than 3,500 kilometres, making it fully consistent with the constraints obtained independently from the odd gravitational harmonics. This depth corresponds to the point at which the electric conductivity becomes large and magnetic drag should suppress differential rotation(7). Given that electric conductivity is dependent on planetary mass, we expect the outer, differentially rotating region to be at least three times deeper in Saturn and to be shallower in massive giant planets and brown dwarfs.

The depth to which Jupiter's observed east-west jet streams extend has been a long-standing question(1,2). Resolving this puzzle has been a primary goal for the Juno spacecraft(3,4), which has been in orbit around the gas giant since July 2016. Juno's gravitational measurements have revealed that Jupiter's gravitational field is north-south asymmetric(5), which is a signature of the planet's atmospheric and interior flows(6). Here we report that the measured odd gravitational harmonics J(3), J(5), J(7) and J(9) indicate that the observed jet streams, as they appear at the cloud level, extend down to depths of thousands of kilometres beneath the cloud level, probably to the region of magnetic dissipation at a depth of about 3,000 kilometres(7,8). By inverting the measured gravity values into a wind field(9), we calculate the most likely vertical profile of the deep atmospheric and interior flow, and the latitudinal dependence of its depth. Furthermore, the even gravity harmonics J(8) and J(10) resulting from this flow profile also match the measurements, when taking into account the contribution of the interior structure(10). These results indicate that the mass of the dynamical atmosphere is about one per cent of Jupiter's total mass.

The Cassini measurements of Saturn's gravity field during its Grand Finale might shed light on a long-standing question regarding the flow on Saturn. While the cloud-level winds are well known, little is known about whether these winds are confined to the outer layers of the planet or penetrate deep into the interior. An additional complexity is added by the uncertainty in the exact rotation period of Saturn, a key factor in determining the cloudlevel winds, with an effect on the north-south symmetric part of the winds. Using Saturn's cloud-level winds we relate the flow to the gravity harmonics. We give a prediction for the odd harmonicsJ(3), J(5), J(7), and J(9) as a function of the flow depth, identifying three ranges of depths. Since the odd harmonics depend solely on the flow, and are not influenced by Saturn's shape and static density distribution, any measured value of the odd harmonics by Cassini can be used to uniquely determine the depth of the flow. We also discuss the flow-induced even harmonics Delta J(2), Delta J(4),...., Delta J(12) that are affected by Saturn's rotation period. While the high-degree even harmonics might also be used to determine the flow depth, the lower-degree even harmonics serve as uncertainties for analysis of the planet's interior structure and composition. Thus, the gravity harmonics measured during the Cassini Grand Finale may be used to get a first-order estimate of the flow structure and to better constrain the planet's density structure and composition.

The close-by orbits of the ongoing Juno mission allow measuring with unprecedented accuracy Jupiter's low-degree even gravity moments J₂, J₄, J₆, and J₈. These can be used to better determine Jupiter's internal density profile and constrain its core mass. Yet the largest unknown on these gravity moments comes from the effect of differential rotation, which gives a degree of freedom unaccounted for by internal structure models. Here considering a wide range of possible internal flow structures and dynamical considerations, we provide upper bounds to the effect of dynamics (differential rotation) on the low-degree gravity moments. In light of the recent Juno gravity measurements and their small uncertainties, this allows differentiating between the various models suggested for Jupiter's internal structure.

Observations of the flow on Jupiter exists essentially only for the cloud-level, which is dominated by strong east-west jet-streams. These have been suggested to result from dynamics in a superficial thin weather-layer, or alternatively be a manifestation of deep interior cylindrical flows. However, it is possible that the observed wind is indeed superficial, yet there exists a completely decoupled deep flow. To date, all models linking the wind, via the induced density anomalies, to the gravity field, to be measured by Juno, consider only flow that is a projection of the observed cloud-level wind. Here we explore the possibility of complex wind dynamics that include both the shallow weather-layer wind, and a deep flow that is decoupled from the flow above it. The upper flow is based on the observed cloud-level flow and is set to decay with depth. The deep flow is constructed to produce cylindrical structures with variable width and magnitude, thus allowing for a wide range of possible scenarios for the unknown deep flow. The combined flow is then related to the density anomalies and gravitational moments via a dynamical model. An adjoint inverse model is used for optimizing the parameters controlling the setup of the deep and surface-bound flows, so that these flows can be reconstructed given a gravity field. We show that the model can be used for examination of various scenarios, including cases in which the deep flow is dominating over the surface wind, and discuss the uncertainties associated with the model solution. The flexibility of the adjoint method allows for a wide range of dynamical setups, so that when new observations and physical understanding will arise, these constraints could be easily implemented and used to better decipher Jupiter flow dynamics. (C) 2017 Elsevier Inc. All rights reserved.

The upcoming Juno and Cassini gravity measurements of Jupiter and Saturn, respectively, will allow probing the internal dynamics of these planets through accurate analysis of their gravity spectra. To date, two general approaches have been suggested for relating the flow velocities and gravity fields. In the first, barotropic potential surface models, which naturally take into account the oblateness of the planet, are used to calculate the gravity field. However, barotropicity restricts the flows to be constant along cylinders parallel to the rotation axis. The second approach, calculated in the reference frame of the rotating planet, assumes that due to the large scale and rapid rotation of these planets, the winds are to leading order in geostrophic balance. Therefore, thermal wind balance relates the wind shear to the density gradients. While this approach can take into account any internal flow structure, it is limited to only calculating the dynamical gravity contributions, and has traditionally assumed spherical symmetry. This study comes to relate the two approaches both from a theoretical perspective, showing that they are analytically identical in the barotropic limit, and numerically, through systematically comparing the different model solutions for the gravity harmonics. For the barotropic potential surface models we employ two independent solution methods - the potential-theory and Maclaurin spheroid methods. We find that despite the sphericity assumption, in the barotropic limit the thermal wind solutions match well the barotropic oblate potential-surface solutions.

Jupiter's Great Red Spot (GRS) is the most dominant and long-lived feature in Jupiter's atmosphere. However, whether this is a shallow atmospheric feature or a deeply rooted vortex has remained an open question. In this study, we assess the possibility of inferring the depth of the GRS by the upcoming Juno gravity experiment. This is achieved by an exploration of the possible gravitational signature of the vortex by systematically extending the surface winds into the interior and analyzing the resulting gravity signal. The gravity anomaly is then compared to the expected accuracy in the retrieval of the surface gravity at the GRS location obtained with numerical simulations of the Doppler data inversion based on the expected trajectory of the spacecraft. Starting from observations of the atmospheric velocity at the cloud level, we project the wind using a decay scale height along coaxial cylinders parallel to the spin axis and explore a wide range of decay scale heights in the radial direction. Assuming the large scale vortex dynamics are geostrophic, and therefore thermal wind balance holds, the density anomaly distribution due to Jupiter's winds can be derived from the velocity maps. The novelty of this approach is in the integration of thermal wind relations over a three-dimensional grid, and in the inclusion of the observed meridional velocity as measured during the Cassini flyby of Jupiter. The perturbations in the mean zonal flow give rise to non zero tesseral spherical harmonics in Jupiter's gravitational potential. We provide an estimate of this asymmetric gravity coefficients for different values of the wind decay scale height. We conclude that the mass anomaly associated with the GRS is detectable by the Juno gravity experiment if the vortex is deep, characterized by a vertical height larger than 2,000 km below the cloud level of Jupiter, and that the large mass involved with deep winds does not render much the ability to measure the feature.

There are many sources of very low frequency (VLF3-30kHz) and extremely low frequency (ELF3-3000Hz) radiation in the Earth-ionosphere waveguide (e.g., lightning and ELF/VLF communication transmitters). At distances of thousands of kilometers from these sources, the vertical component of the ELF/VLF AC magnetic fields is expected to be very weak and several orders of magnitude lower than the horizontal magnetic components. However, measurements in Israel show a relatively strong vertical magnetic component in both the ELF and VLF bands, at the same order of magnitude as the horizontal components. Our measurements suggest that the real Earth-ionosphere waveguide might often be very different from the theoretical waveguide used in model calculations. In addition, our results imply that using only the horizontal components for direction finding or the absolute magnetic field strength may result in errors, since often a significant fraction of the magnetic field energy hides in the vertical component.

Estimates of global thunderstorm activity have been made predominately by direct measurements of lightning discharges around the globe, either by optical measurements from satellites, or using ground-based radio antennas. In this paper we propose a new methodology in which thunderstorm clusters are constructed based on the lightning strokes detected by the World Wide Lightning Location Network (WWLLN) in the very low frequency range. We find that even with low lightning detection efficiency on a global scale, the spatial and temporal distribution of global thunderstorm cells is well reproduced. This is validated by comparing the global diurnal variations of the thunderstorm cells, and the currents produced by these storms, with the well-known Carnegie Curve, which represents the mean diurnal variability of the global atmospheric electric circuit, driven by thunderstorm activity. While the Carnegie Curve agrees well with our diurnal thunderstorm cluster variations, there is little agreement between the Carnegie Curve and the diurnal variation in the number of lightning strokes detected by the WWLLN. When multiplying the number of clusters we detect by the mean thunderstorm conduction current for land and ocean thunderstorms (Mach et al 2011 J. Geophys. Res. 116 D05201) we get a total average current of about 760 A. Our results show that thunderstorms alone explain more than 90% in the variability of the global electric circuit. However, while it has been previously shown that 90% of the global lightning occurs over continental landmasses, we show that around 50% of the thunderstorms are over the oceans, and from 00-09UTC there are more thunderstorm cells globally over the oceans than over the continents. Since the detection efficiency of the WWLLN system has increased over time, we estimate that the lower bound of the mean number of global thunderstorm cells in 2012 was around 1050 per hour, varying from around 840 at 03UTC to 1150 storms at 19UTC.

The FLASH project was implemented from 2006 to 2010 under the EU FP6 framework. The project focused on using lightning observations to better understand and predict convective storms that result in flash floods. As part of the project 23 case studies of flash floods in the Mediterranean region were examined. For the analysis of these storms lightning data from the ZEUS network were used together with satellite derived rainfall estimates in order to understand the storm development and electrification. In addition, these case studies were simulated using mesoscale meteorological models to better understand the meteorological and synoptic conditions leading up to these intense storms. As part of this project tools for short term predictions (nowcasts) of intense convection across the Mediterranean and Europe, and long term forecasts (a few days) of the likelihood of intense convection were developed. The project also focused on educational outreach through our website http://flashproject.org supplying real time lightning observations, real time experimental nowcasts, forecasts and educational materials. While flash floods and intense thunderstorms cannot be prevented as the climate changes, long-range regional lightning networks can supply valuable data, in real time, for warning end-users and stakeholders of imminent intense rainfall and possible flash floods. (C) 2011 Elsevier Ltd. All rights reserved.

Thunderstorms are often the cause of severe and disastrous flash floods. Lightning activity within these storms can be detected and monitored continuously from thousands of kilometers away and can therefore be very useful in improving forecasts and now-casts of severe thunderstorms. An improvement in the now-casting of such storms can assist in reducing damages and saving lives.Using the ZEUS ground-based VLF lightning detection network and the Warning Decision Support System-Integrated Information (WDSS-II) software, we performed now-casting simulations using 1 year of lightning data over the Mediterranean area. Thousands of thunderstorms were observed and now-cast 30, 60, 90 and 120 min ahead. Statistical analysis was then done by calculating the hit, miss and false-alarm rates, as well as the POD, FAR and CSI scores in order to determine the success of the now-casting.The results show that the algorithm is overall successful in now-casting the location of the lightning clusters, especially when applied to strong and consistent lightning events (it is these events which also have the strongest connection to flash floods). The probability of detection (POD) values range between 0.46 for 30 minute now-casts and 0.25 for 120 minute now-casts. The critical success index (CSI) values are quite similar, but slightly lower. The now-casting has a low false-alarm rate, 0.03 for 30 minute now-casts, which is also beneficial for operational use.This method has been implemented for the use in real-time now-casting and is used in the EU FLASH project to seek and track areas of thunderstorm risk according to lightning intensity. The experimental now-casts appear on the project website (www.flashproject.org/). (C) 2010 Elsevier B.V. All rights reserved.

An experimental ENSO prediction system is presented, based on an ocean general circulation model (GCM) coupled to a statistical atmosphere and the adjoint method of 4D variational data assimilation. The adjoint method is used to initialize the coupled model, and predictions are performed for the period 1980-99. The coupled model is also initialized using two simpler assimilation techniques: forcing the ocean model with observed sea surface temperature and surface fluxes, and a 3D variational data assimilation (3DVAR) method, similar to that used by the National Centers for Environmental Prediction (NCEP) for operational ENSO prediction. The prediction skill of the coupled model initialized by the three assimilation methods is then analyzed and compared. The effect of the assimilation period used in the adjoint method is studied by using 3-,6-, and 9-month assimilation periods. Finally, the possibility of assimilating only the anomalies with respect to observed climatology in order to circumvent systematic model biases is examined. It is found that the adjoint method does seem to have the potential for improving over simpler assimilation schemes. The improved skill is mainly at prediction intervals of more than 6 months, where the coupled model dynamics start to influence the model solution. At shorter prediction time intervals, the initialization using the forced ocean model or the 3DVAR may result in a better prediction skill. The assimilation of anomalies did not have a substantial effect on the prediction skill of the coupled model. This seems to indicate that in this model the climatology bias, which is compensated for by the anomaly assimilation, is less significant for the predictive skill than the bias in the model variability, which cannot be eliminated using the anomaly assimilation. Changing the optimization period from 6 to 3 to 9 months showed that the period of 6 months seems to be a near-optimal choice for this model.

One of the major factors determining the strength and extent of ENSO events is the instability state of the equatorial Pacific coupled ocean-atmosphere system and its seasonal variations. This study analyzes the coupled instability in a hybrid coupled model of the Indo-Pacific region, using the adjoint method for sensitivity studies. It is found that the seasonal changes in the ocean-atmosphere instability strength in the model used here are related to the outcropping of the thermocline in the east equatorial Pacific. From July to December, when the thermocline outcrops over a wide area in the east Pacific, there is a strong surface-thermocline connection and anomalies that arrive as Kelvin waves from the west along the thermocline can reach the surface and affect the SST and thus the coupled system. Conversely, from February to June, when the thermocline outcropping is minimal, the surface decouples from the thermocline and temperature anomalies in the thermocline depth range do not affect the surface and dissipate within the thermocline. The role of vertical mixing rather than upwelling in linking vertical thermocline movements to SST change is emphasized. It is therefore suggested that the seasonal ocean-atmosphere instability strength in the equatorial Pacific is strongly influenced by the thermocline outcropping and its seasonal modulation, a physical mechanism that is often neglected in intermediate coupled models and that can be represented properly only in models that employ the full dynamics of the mixed layer.