The Problem of the Height Dependence of Magnetic Fields in Sunspots

To understand the physics of sunspots, it is important to know the properties of their magnetic field, and especially its height stratification plays a substantial role. There are mainly two methods to assess this stratification, but they yield different magnetic gradients in the photospheric layers. Determinations based on the several spectral lines of different formation heights and the slope of their profiles result in gradients of ?2 to ?3 G?km^?1, or even steeper. This is similar for the total magnetic field strength and for the vertical component of the magnetic field. The other option is to determine the horizontal partial derivatives of the magnetic field, and with the condition \(\operatorname{div} {{\boldsymbol {B}}} = 0\) also the vertical derivative is known. With this method, gradients of ?0.5 G?km^?1 and even shallower are obtained. Obviously, these results do not agree. If chromospheric spectral lines are included, only shallow gradients around ?0.5 G?km^?1 are obtained. Shallow gradients are also found from gyro-resonance measurements in the radio wave range 300?–?2000 GHz.Some indirect methods are also considered, but they cannot clarify the total picture. An analysis of a numerical simulation of a sunspot indicates a shallow gradient over a wide height range, but with slightly steeper gradients in deep layers.Several ideas to explain the discrepancy are also discussed. With no doubts cast on Maxwell’s equations, the first one is to look at the uncertainties of the formation heights of spectral lines, but a wider range of these heights would require an extension of the solar photosphere that is incompatible with observations and the theory of stellar atmospheres. Submerging and rising magnetic flux might play a role in the outer penumbra, if the resolution is too low to separate them, but it is not likely that this effect acts also in the umbra. A quick investigation assuming a spatial small scale structure of sunspots together with twist and writhe of individual flux tubes shows a reduction of the measured magnetic field strength for spectral lines sensitive to a larger height range. However, sophisticated investigations are required to prove that the explanation for the discrepancy lies here, and the problem of the height gradient of the magnetic field in sunspots is still not solved.     

Mid-term periodicities and heliospheric modulation of coronal index and solar flare index during solar cycles 22–23

We applied fast Fourier transform techniques and Morlet wavelet transform on the time series data of coronal index, solar flare index, and galactic cosmic ray, for the period 1986–2008, in order to investigate the long- and mid-term periodicities including the Rieger (\({\sim }130\) to \({\sim }190\) days), quasi-period (\({\sim }200\) to \({\sim }374\) days), and quasi-biennial periodicities (\({\sim }1.20\) to \({\sim }3.27\) years) during the combined solar cycles 22–23. We emphasize the fact that a lesser number of periodicities are found in the range of low frequencies, while the higher frequencies show a greater number of periodicities. The rotation rates at the base of convection zone have periods for coronal index of \({\sim }1.43\) years and for solar flare index of \({\sim }1.41\) year, and galactic cosmic ray, \({\sim }1.35\) year, during combined solar cycles 22–23. In relation to these two solar parameters (coronal index and solar flare index), for the solar cycles 22–23, we found that galactic cosmic ray modulation at mid cut-off rigidity (\(\hbox {Rc} = 2.43\hbox {GV}\)) is anti-correlated with time-lag of few months.     

Effective Acceleration Model for the Arrival Time of Interplanetary Shocks driven by Coronal Mass Ejections

In a previous work (Paouris and Mavromichalaki in Solar Phys. 292, 30, 2017), we presented a total of 266 interplanetary coronal mass ejections (ICMEs) with as much information as possible. We developed a new empirical model for estimating the acceleration of these events in the interplanetary medium from this analysis. In this work, we present a new approach on the effective acceleration model (EAM) for predicting the arrival time of the shock that preceds a CME, using data of a total of 214 ICMEs. For the first time, the projection effects of the linear speed of CMEs are taken into account in this empirical model, which significantly improves the prediction of the arrival time of the shock. In particular, the mean value of the time difference between the observed time of the shock and the predicted time was equal to +3.03 hours with a mean absolute error (MAE) of 18.58 hours and a root mean squared error (RMSE) of 22.47 hours. After the improvement of this model, the mean value of the time difference is decreased to ?0.28 hours with an MAE of 17.65 hours and an RMSE of 21.55 hours. This improved version was applied to a set of three recent Earth-directed CMEs reported in May, June, and July of 2017, and we compare our results with the values predicted by other related models.     

Revisiting Field Burial by Accretion onto Neutron Stars

The surface magnetic field strength of millisecond pulsars (MSPs) is found to be about 4 orders of magnitude lower than that of garden variety radio pulsars (with a spin of \({\sim }\)0.5–5 s and \(B{\sim }10^{12}\hbox { G}\)). The exact mechanism of the apparent reduction of field strength in MSPs is still a subject of debate. One of the proposed mechanisms is burial of the surface magnetic field under matter accreted from a companion. In this article we review the recent work on magnetic confinement of accreted matter on neutron stars poles. We present the solutions of the magneto-static equations with a more accurate equation of state of the magnetically confined plasma and discuss its implications for the field burial mechanism.     

On the co-orbital motion in the planar restricted three-body problem: the quasi-satellite motion revisited

In the framework of the planar and circular restricted three-body problem, we consider an asteroid that orbits the Sun in quasi-satellite motion with a planet. A quasi-satellite trajectory is a heliocentric orbit in co-orbital resonance with the planet, characterized by a nonzero eccentricity and a resonant angle that librates around zero. Likewise, in the rotating frame with the planet, it describes the same trajectory as the one of a retrograde satellite even though the planet acts as a perturbator. In the last few years, the discoveries of asteroids in this type of motion made the term “quasi-satellite” more and more present in the literature. However, some authors rather use the term “retrograde satellite” when referring to this kind of motion in the studies of the restricted problem in the rotating frame. In this paper, we intend to clarify the terminology to use, in order to bridge the gap between the perturbative co-orbital point of view and the more general approach in the rotating frame. Through a numerical exploration of the co-orbital phase space, we describe the quasi-satellite domain and highlight that it is not reachable by low eccentricities by averaging process. We will show that the quasi-satellite domain is effectively included in the domain of the retrograde satellites and neatly defined in terms of frequencies. Eventually, we highlight a remarkable high eccentric quasi-satellite orbit corresponding to a frozen ellipse in the heliocentric frame. We extend this result to the eccentric case (planet on an eccentric motion) and show that two families of frozen ellipses originate from this remarkable orbit.     

Spectral Trends of Solar Bursts at Sub-THz Frequencies

Previous sub-THz studies were derived from single-event observations. We here analyze for the first time spectral trends for a larger collection of sub-THz bursts. The collection consists of a set of 16 moderate to small impulsive solar radio bursts observed at 0.2 and 0.4 THz by the Solar Submillimeter-wave Telescope (SST) in 2012?–?2014 at El Leoncito, in the Argentinean Andes. The peak burst spectra included data from new solar patrol radio telescopes (45 and 90 GHz), and were completed with microwave data obtained by the Radio Solar Telescope Network, when available. We critically evaluate errors and uncertainties in sub-THz flux estimates caused by calibration techniques and the corrections for atmospheric transmission, and introduce a new method to obtain a uniform flux scale criterion for all events. The sub-THz bursts were searched during reported GOES soft X-ray events of class C or larger, for periods common to SST observations. Seven out of 16 events exhibit spectral maxima in the range 5?–?40 GHz with fluxes decaying at sub-THz frequencies (three of them associated to GOES class X, and four to class M). Nine out of 16 events exhibited the sub-THz spectral component. In five of these events, the sub-THz emission fluxes increased with a separate frequency from that of the microwave spectral component (two classified as X and three as M), and four events have only been detected at sub-THz frequencies (three classified as M and one as C). The results suggest that the THz component might be present throughout, with the minimum turnover frequency increasing as a function of the energy of the emitting electrons. The peculiar nature of many sub-THz burst events requires further investigations of bursts that are examined from SST observations alone to better understand these phenomena.     

On Flare-CME Characteristics from Sun to Earth Combining Remote-Sensing Image Data with <Emphasis Type="Italic">In Situ</Emphasis> Measurements Supported by Modeling

We analyze the well-observed flare and coronal mass ejection (CME) from 1 October 2011 (SOL2011-10-01T09:18) covering the complete chain of effects – from Sun to Earth – to better understand the dynamic evolution of the CME and its embedded magnetic field. We study in detail the solar surface and atmosphere associated with the flare and CME using the Solar Dynamics Observatory (SDO) and ground-based instruments. We also track the CME signature off-limb with combined extreme ultraviolet (EUV) and white-light data from the Solar Terrestrial Relations Observatory (STEREO). By applying the graduated cylindrical shell (GCS) reconstruction method and total mass to stereoscopic STEREO-SOHO (Solar and Heliospheric Observatory) coronagraph data, we track the temporal and spatial evolution of the CME in the interplanetary space and derive its geometry and 3D mass. We combine the GCS and Lundquist model results to derive the axial flux and helicity of the magnetic cloud (MC) from in situ measurements from Wind. This is compared to nonlinear force-free (NLFF) model results, as well as to the reconnected magnetic flux derived from the flare ribbons (flare reconnection flux) and the magnetic flux encompassed by the associated dimming (dimming flux). We find that magnetic reconnection processes were already ongoing before the start of the impulsive flare phase, adding magnetic flux to the flux rope before its final eruption. The dimming flux increases by more than 25% after the end of the flare, indicating that magnetic flux is still added to the flux rope after eruption. Hence, the derived flare reconnection flux is most probably a lower limit for estimating the magnetic flux within the flux rope. We find that the magnetic helicity and axial magnetic flux are lower in the interplanetary space by ～?50% and 75%, respectively, possibly indicating an erosion process. A CME mass increase of 10% is observed over a range of \({\sim}\,4\,\mbox{--}\,20~\mathrm{R}_{\odot }\). The temporal evolution of the CME-associated core-dimming regions supports the scenario that fast outflows might supply additional mass to the rear part of the CME.     

The smoke is rising but where is the fire? Exploring effective online map design for wildfire warnings

The current study sought to offer guidance for developing effective web-based mapping tools for wildfire warnings by identifying (1) the important content for facilitating individuals’ decision-making, and (2) the optimal interface design for ensuring usability and ease of information access. A map-based warning tool was prototyped in the Australian context, followed by a usability and effectiveness evaluation through individual interviews and verbal protocol analysis to assess participants’ interaction with the mapping interface and information in response to the simulated warning scenario. The results demonstrated variations in participants’ approaches to wildfire warning response, revealing varied information needs. Specifically, most participants relied on their own assessment of the prospective threat, requiring specific wildfire-related information before eliciting a response. In contrast, the decision of a minority of the participants was motivated by response guidance from agencies, and accurate wildfire information was less important for their response. Imperative information for both types of residents therefore needs to be highlighted in a map-based warning tool to cater to a wide audience. Furthermore, a number of heuristics were identified for designing effective interactive functions to facilitate the control of, and access to, the various maps and textual information presented on the map-based warning interface.     

A low-dimensional subsurface model for saturated and unsaturated flow processes: ability to address heterogeneity

A low-dimensional model that describes both saturated and unsaturated flow processes in a single equation is presented. Subsurface flow processes in the groundwater, the vadose zone, and the capillary fringe are accounted for through the computation of aggregated hydrodynamic parameters that result from the integration of the governing flow equations from the bedrock to the land surface. The three-dimensional subsurface flow dynamics are thus described by a two-dimensional equation, allowing for a drastic reduction of model unknowns and simplification of the model parameterizations. This approach is compared with a full resolution of the Richards equation in different synthetic test cases. Because the model reduction stems from the vertical integration of the flow equations, the test cases all use different configurations of heterogeneity for vertical cross-sections of a soil-aquifer system. The low-dimensional flow model shows strong consistency with results from a complete resolution of the Richards equation for both the water table and fluxes. The proposed approach is therefore well suited to the accurate reproduction of complex subsurface flow processes.