首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 31 毫秒
1.
Mechanisms of the formation and stability of sunspots are among the longest-standing and intriguing puzzles of solar physics and astrophysics. Sunspots are controlled by subsurface dynamics, hidden from direct observations. Recently, substantial progress in our understanding of the physics of the turbulent magnetized plasma in strong-field regions has been made by using numerical simulations and local helioseismology. Both the simulations and helioseismic measurements are extremely challenging, but it is becoming clear that the key to understanding the enigma of sunspots is a synergy between models and observations. Recent observations and radiative MHD numerical models have provided a convincing explanation for the Evershed flows in sunspot penumbrae. Also, they lead to the understanding of sunspots as self-organized magnetic structures in the turbulent plasma of the upper convection zone, which are maintained by a large-scale dynamics. Local helioseismic diagnostics of sunspots still have many uncertainties, some of which are discussed in this review. However, there have been significant achievements in resolving these uncertainties, verifying the basic results by new high-resolution observations, testing the helioseismic techniques by numerical simulations, and comparing results obtained by different methods. For instance, a recent analysis of helioseismology data from the Hinode space mission has successfully resolved several uncertainties and concerns (such as the inclined-field and phase-speed filtering effects) that might affect the inferences of the subsurface wave-speed structure of sunspots and the flow pattern. It is becoming clear that for the understanding of the phenomenon of sunspots it is important to further improve the helioseismology methods and investigate the whole life cycle of active regions, from magnetic flux emergence to dissipation. The Solar Dynamics Observatory mission has started to provide data for such investigations.  相似文献   

2.
We describe a procedure for automated detection of sunspots from SoHO/MDI full-disk continuum images. The MDI Level-1.8 continuum images were first corrected for the limb darkening and stray light, and then were flat-fielded. Sunspots were extracted using a newly developed automated sunspot detection procedure, which is based on the level set, namely the selective binary and Gaussian function regularized level set (SBGFRLS) method (Zhang et al., Image Vis. Comput. 28, 668, 2010). In this method we initialize a two-dimensional level-set function and evolve it using a signed pressure force (SPF) function. For sunspot detection, the level-set function was defined twice, first for umbra and then for penumbra extraction. Using this procedure, along with the characterization of detected sunspots we have also generated tracking reports of all sunspots in a fully unsupervised manner.  相似文献   

3.
The high quality of the asteroseismic data provided by space missions such as CoRoT (Michel et al. in The CoRoT Mission, ESA Spec. Publ. vol. 1306, p. 39, 2006) or expected from new operating missions such as Kepler (Christensen-Dalsgaard et al. in Commun. Asteroseismol. 150:350, 2007) requires the capacity of stellar evolution codes to provide accurate models whose numerical precision is better than the expected observational errors (i.e. below 0.1 μHz on the frequencies in the case of CoRoT). We present a review of some thorough comparisons of stellar models produced by different evolution codes, involved in the CoRoT/ESTA activities (Monteiro in Evolution and Seismic Tools for Stellar Astrophysics, 2009). We examine the numerical aspects of the computations as well as the effects of different implementations of the same physics on the global quantities, physical structure and oscillations properties of the stellar models. We also discuss a few aspects of the input physics.  相似文献   

4.
Long-term variations of solar differential rotation and sunspot activity are investigated through re-analyzing the data on parameters of the differential-rotation law obtained by Makarov, Tlatov, and Callebaut (Solar Phys. 170, 373, 1997), Javaraiah, Bertello, and Ulrich (Astrophys. J. 626, 579, 2005a; Solar Phys. 232, 25, 2005b), and Javaraiah et al. (Solar Phys. 257, 61, 2009). Our results indicate that the solar-surface-rotation rate at the Equator (indicated by the A-parameter of the standard solar-rotation law) shows a secular decrease since Cycle 12 onwards, given by about 1?–?1.5×10?3 (deg?day?1?year?1). The B-parameter of the standard differential-rotation law seems to also show a secular decrease since Cycle 12 onwards, but of weak statistical significance. The rotation rate averaged over latitudes 0°?–?40° does not show a secular trend of statistical significance. Moreover, the average sunspot area shows a secular increase of statistical significance since Cycle 12 onwards, while a negative correlation is found between the level of sunspot activity (indicated by the average sunspot area) and the solar equatorial rotation on long-term scales.  相似文献   

5.
In local helioseismology, numerical simulations of wave propagation are useful to model the interaction of solar waves with perturbations to a background solar model. However, the solution to the linearised equations of motion include convective modes that can swamp the helioseismic waves that we are interested in. In this article, we construct background solar models that are stable against convection, by modifying the vertical pressure gradient of Model S (Christensen-Dalsgaard et al., 1996, Science 272, 1286) relinquishing hydrostatic equilibrium. However, the stabilisation affects the eigenmodes that we wish to remain as close to Model S as possible. In a bid to recover the Model S eigenmodes, we choose to make additional corrections to the sound speed of Model S before stabilisation. No stabilised model can be perfectly solar-like, so we present three stabilised models with slightly different eigenmodes. The models are appropriate to study the f and p 1 to p 4 modes with spherical harmonic degrees in the range from 400 to 900. Background model CSM has a modified pressure gradient for stabilisation and has eigenfrequencies within 2% of Model S. Model CSM_A has an additional 10% increase in sound speed in the top 1 Mm resulting in eigenfrequencies within 2% of Model S and eigenfunctions that are, in comparison with CSM, closest to those of Model S. Model CSM_B has a 3% decrease in sound speed in the top 5 Mm resulting in eigenfrequencies within 1% of Model S and eigenfunctions that are only marginally adversely affected. These models are useful to study the interaction of solar waves with embedded three-dimensional heterogeneities, such as convective flows and model sunspots. We have also calculated the response of the stabilised models to excitation by random near-surface sources, using simulations of the propagation of linear waves. We find that the simulated power spectra of wave motion are in good agreement with an observed SOHO/MDI power spectrum. Overall, our convectively stabilised background models provide a good basis for quantitative numerical local helioseismology. The models are available for download from http://www.mps.mpg.de/projects/seismo/NA4/ .  相似文献   

6.
A. Gil  M. V. Alania 《Solar physics》2013,282(2):565-578
Recently Pop (Solar Phys. 276, 351, 2012) identified a Laplace (or double exponential) distribution in the number of days with a given absolute value in the change over a day, in sunspot number, for days on which the sunspot number does change. We show this phenomenological rule has a physical origin attributable to sunspot formation, evolution, and decay, rather than being due to the changes in sunspot number caused by groups rotating onto and off the visible disc. We also demonstrate a simple method to simulate daily sunspot numbers over a solar cycle using the Pop (Solar Phys. 276, 351, 2012) result, together with a model for the cycle variation in the mean sunspot number. The procedure is applied to three recent solar cycles. We check that the simulated sunspot numbers reproduce the observed distribution of daily changes over those cycles.  相似文献   

7.
One goal of helioseismology is to determine the subsurface structure of sunspots. In order to do so, it is important to understand first the near-surface effects of sunspots on solar waves, which are dominant. Here we construct simplified, cylindrically-symmetric sunspot models that are designed to capture the magnetic and thermodynamics effects coming from about 500 km below the quiet-Sun τ 5000=1 level to the lower chromosphere. We use a combination of existing semi-empirical models of sunspot thermodynamic structure (density, temperature, pressure): the umbral model of Maltby et al. (1986, Astrophys. J. 306, 284) and the penumbral model of Ding and Fang (1989, Astron. Astrophys. 225, 204). The OPAL equation-of-state tables are used to derive the sound-speed profile. We smoothly merge the near-surface properties to the quiet-Sun values about 1 Mm below the surface. The umbral and penumbral radii are free parameters. The magnetic field is added to the thermodynamic structure, without requiring magnetostatic equilibrium. The vertical component of the magnetic field is assumed to have a Gaussian horizontal profile, with a maximum surface field strength fixed by surface observations. The full magnetic-field vector is solenoidal and determined by the on-axis vertical field, which, at the surface, is chosen such that the field inclination is 45° at the umbral – penumbral boundary. We construct a particular sunspot model based on SOHO/MDI observations of the sunspot in active region NOAA 9787. The helioseismic signature of the model sunspot is studied using numerical simulations of the propagation of f, p 1, and p 2 wave packets. These simulations are compared against cross-covariances of the observed wave field. We find that the sunspot model gives a helioseismic signature that is similar to the observations.  相似文献   

8.
The Helioseismic and Magnetic Imager (HMI) instrument onboard the Solar Dynamics Observatory (SDO) satellite is designed to produce high-resolution Doppler-velocity maps of oscillations at the solar surface with high temporal cadence. To take advantage of these high-quality oscillation data, a?time?–?distance helioseismology pipeline (Zhao et al., Solar Phys. submitted, 2010) has been implemented at the Joint Science Operations Center (JSOC) at Stanford University. The aim of this pipeline is to generate maps of acoustic travel times from oscillations on the solar surface, and to infer subsurface 3D flow velocities and sound-speed perturbations. The wave travel times are measured from cross-covariances of the observed solar oscillation signals. For implementation into the pipeline we have investigated three different travel-time definitions developed in time?–?distance helioseismology: a Gabor-wavelet fitting (Kosovichev and Duvall, SCORE’96: Solar Convection and Oscillations and Their Relationship, ASSL, Dordrecht, 241, 1997), a?minimization relative to a reference cross-covariance function (Gizon and Birch, Astrophys. J. 571, 966, 2002), and a linearized version of the minimization method (Gizon and Birch, Astrophys. J. 614, 472, 2004). Using Doppler-velocity data from the Michelson Doppler Imager (MDI) instrument onboard SOHO, we tested and compared these definitions for the mean and difference travel-time perturbations measured from reciprocal signals. Although all three procedures return similar travel times in a quiet-Sun region, the method of Gizon and Birch (Astrophys. J. 614, 472, 2004) gives travel times that are significantly different from the others in a magnetic (active) region. Thus, for the pipeline implementation we chose the procedures of Kosovichev and Duvall (SCORE’96: Solar Convection and Oscillations and Their Relationship, ASSL, Dordrecht, 241, 1997) and Gizon and Birch (Astrophys. J. 571, 966, 2002). We investigated the relationships among these three travel-time definitions, their sensitivities to fitting parameters, and estimated the random errors that they produce.  相似文献   

9.
The problem of finding nonsingular charged analogue of Schwarzschild’s interior solutions has been reduced to that of finding a monotonically decreasing function f. The models are discussed in generality by imposing reality condition on f. It is shown that the physical solutions are possible only for surface density to central density ratio greater than or equal to 2/3 i.e. $\frac{\rho_{a}}{\rho_{0}}\ge2/3$ . The unphysical nature of solutions with linear equation state has been proved. A generalization procedure has been utilized to generalize solutions by Guilfoyle (1999). Recently found solutions by Gupta and Kumar (2005a, 2005b, 2005c) are generalized by taking particular form of f and seen to have higher mass and more stable. The maximum mass is found to be 1.59482 M Θ . The models have been found to be stable once the physical requirements are established due to mass to radius less than 4/9, total charge to total mass ratio less than 1 and redshift quite low.  相似文献   

10.
We present a model for the reconstruction of spectral solar irradiance between 200 and 400?nm. This model is an extension of the total solar irradiance (TSI) model of Crouch et al. (Astrophys.?J. 677, 723, 2008) which is based on a data-driven Monte Carlo simulation of sunspot emergence, fragmentation, and erosion. The resulting time-evolving daily area distribution of magnetic structures of all sizes is used as input to a four-component irradiance model including contributions from the quiet Sun, sunspots, faculae, and network. In extending the model to spectral irradiance in the near- and mid-ultraviolet, the quiet Sun and sunspot emissivities are calculated from synthetic spectra at T eff=5750?K and 5250?K, respectively. Facular emissivities are calculated using a simple synthesis procedure proposed by Solanki and Unruh (Astron. Astrophys. 329, 747, 1998). The resulting time series of ultraviolet flux is calibrated against the data from the SOLSTICE instrument on the Upper Atmospheric Research Satellite (UARS). Using a genetic algorithm, we invert quiet Sun corrections, profile of facular temperature variations with height, and network model parameters which yield the best fit to these data. The resulting best-fit time series reproduces quite well the solar-cycle timescale variations of UARS ultraviolet observations, as well as the short-timescale fluctuations about the 81 day running mean. We synthesize full spectra between 200 and 400?nm, and validate these against the spectra obtained by the ATLAS-1 and ATLAS-3 missions, finding good agreement, to better than 3?% at most wavelengths. We also compare the UV variability predicted by our reconstructions in the descending phase of sunspot cycle 23 to SORCE/SIM data as well as to other reconstructions. Finally, we use the model to reconstruct the time series of spectral irradiance starting in 1874, and investigate temporal correlations between pairs of wavelengths in the bands of interest for stratospheric chemistry and dynamics.  相似文献   

11.
In this paper, we have investigated plane symmetric cosmological models with negative constant deceleration parameter in Barber’s (Gen. Relativ. Gravit. 14:117, 1982) second self-creation theory in presence of perfect fluid source. For this we use a special law of variation for Hubble parameter proposed by Bermann (Nuovo Cim. B 74:182, 1983) that yields a constant deceleration parameter model of the universe. Some physical properties of the models and entropy are discussed and studied.  相似文献   

12.
Subdwarf B stars (sdBs) can significantly change the ultraviolet spectra of populations at age t～1 Gyr, and have been even included in the evolutionary population synthesis (EPS) models by Han et al. (Mon. Not. R. Astron. Soc. 380:1098, 2007). In this study we present the spectral energy distributions (SEDs) of binary stellar populations (BSPs) by combining the EPS models of Han et al. (Mon. Not. R. Astron. Soc. 380:1098, 2007) and those of the Yunnan group (Zhang et al. in Astron. Astrophys. 415:117, 2004; Mon. Not. R. Astron. Soc. 357:1088, 2005), which have included various binary interactions (except sdBs) in EPS models. This set of SEDs is available upon request from the authors. Using this set of SEDs of BSPs we build the spectra of Burst, E, S0–Sd and Irr types of galaxies by using the package of Bruzual and Charlot (Mon. Not. R. Astron. Soc. 344:1000, 2003). Combined with the photometric data (filters and magnitudes), we obtain the photometric redshifts and morphologies of 1502 galaxies by using the Hyperz code of Bolzonella et al. (Astron. Astrophys. 363:476, 2000). This sample of galaxies is obtained by removing those objects, mismatched with the SDSS/DR7 and GALEX/DR4, from the catalogue of Fukugita et al. (Astron. J. 134:579, 2007). By comparison the results with the SDSS spectroscopic redshifts and the morphological index of Fukugita et al. (Astron. J. 134:579, 2007), we find that the photo-z fluctuate with the SDSS spectroscopic redshifts, while the Sa–Sc galaxies in the catalogue of Fukugita et al. (Astron. J. 134:579, 2007) are classified earlier as Burst-E galaxies.  相似文献   

13.
This paper studies the stability of infinitesimal motions about the triangular equilibrium points in the elliptic restricted three body problem assuming bigger primary as a source of radiation and the smaller one a triaxial rigid body. The perturbation technique developed by Bennet (Icarus 4:177, 1965b) has been used for determination of characteristic exponents. This technique is based on Floquet’s Theory for determination of characteristic exponents in the system with periodic coefficients. The results of the study are analytical and numerical expressions are simulated for the transition curves bounding the region of stability in the μ–e plane, accurate to O(e 2). The unstable region is found to be divided into three parts. The effect of radiation parameter is significant. For small values of e, the results are in favor with the numerical analysis of Danby (Astron. J. 69:166, 1964), Bennet (Icarus 4:177, 1965b), Alfriend and Rand (AIAA J. 6:1024, 1969). The effect of radiation pressure is significant than the oblateness and triaxiality of the primaries.  相似文献   

14.
G. J. D. Petrie 《Solar physics》2014,289(10):3663-3680
It is shown that expressions for the global Lorentz force associated with a flaring active region derived by Fisher et al. (Solar Phys. 277, 59, 2012) can be used to estimate the Lorentz-force changes for strong fields in large structures over photospheric subdomains within active regions. Gary’s (Solar Phys. 203, 71, 2001) model for the stratified solar atmosphere is used to demonstrate that in large-scale structures with typical horizontal magnetic length scale ??300 km and with strong magnetic fields (≥?1 kG at the τ=1 opacity layer at 5000 Å), the Lorentz force acting on the photosphere may be approximated by a surface integral based on photospheric boundary data alone. These conditions cover many of the sunspot fields and major neutral lines that have been studied using Fisher et al.’s (2012) expressions over the past few years. The method gives a reasonable estimate of flare-related Lorentz-force changes based on photospheric magnetogram observations provided that the Lorentz-force changes associated with the flare have a lasting effect on the observed fields, and they are not immediately erased by post-flare equilibration processes.  相似文献   

15.
We applied automatic identification of sunspot umbrae and penumbrae to daily observations from the Helioseismic Magnetic Imager (HMI) on board the Solar Dynamics Observatory (SDO) to study their magnetic flux density (B) and area (A). The results confirm an already known logarithmic relationship between the area of sunspots and their maximum flux density. In addition, we find that the relation between average magnetic flux density ( $B_{\rm avg}$ ) and sunspot area shows a bimodal distribution: for small sunspots and pores (A≤20 millionth of solar hemisphere, MSH), $B_{\rm avg} \approx 800~\mbox{G}$ (gauss), and for large sunspots (A≥100 MSH), $B_{\rm avg}$ is about 600 G. For intermediate sunspots, average flux density linearly decreases from about 800 G to 600 G. A similar bimodal distribution was found in several other integral parameters of sunspots. We show that this bimodality can be related to different stages of sunspot penumbra formation and can be explained by the difference in average inclination of magnetic fields at the periphery of small and large sunspots.  相似文献   

16.
The symmetric trace free (STF) tensor formalism, developed by Hartmann et al. (Celest Mech Dyn Astron 60:139–159. doi:10.1007/BF00693097, 1994), is a nice tool, not much used in Celestial Mechanics. It is fully equivalent to the usual spherical harmonics but permits more elegant and compact formulations. The coupling between the gravitational fields of extended bodies with this formalism has been used in Mathis and Le Poncin-Lafitte (Astron Astrophys 497:889–910. doi:10.1051/0004-6361/20079054, 2009) for binary stars or planetary systems, but not yet applied to binary asteroids. However, binary asteroids are common in the Solar System and usually their study requires a full two rigid body approach. The formulation of the two-body interaction potential in the STF formalism in the full two rigid body problem is detailed and completed in this article. An application to the binary asteroid (66391) 1999 KW4 is presented with a comparison of our results with other results of the literature for validation.  相似文献   

17.
We study the capture and crossing probabilities in the 3:1 mean motion resonance with Jupiter for a small asteroid that migrates from the inner to the middle Main Belt under the action of the Yarkovsky effect. We use an algebraic mapping of the averaged planar restricted three-body problem based on the symplectic mapping of Hadjidemetriou (Celest Mech Dyn Astron 56:563–599, 1993), adding the secular variations of the orbit of Jupiter and non-symplectic terms to simulate the migration. We found that, for fast migration rates, the captures occur at discrete windows of initial eccentricities whose specific locations depend on the initial resonant angles, indicating that the capture phenomenon is not probabilistic. For slow migration rates, these windows become narrower and start to accumulate at low eccentricities, generating a region of mutual overlap where the capture probability tends to 100 %, in agreement with the theoretical predictions for the adiabatic regime. Our simulations allow us to predict the capture probabilities in both the adiabatic and non-adiabatic cases, in good agreement with results of Gomes (Celest Mech Dyn Astron 61:97–113, 1995) and Quillen (Mon Not RAS 365:1367–1382, 2006). We apply our model to the case of the Vesta asteroid family in the same context as Roig et al. (Icarus 194:125–136, 2008), and found results indicating that the high capture probability of Vesta family members into the 3:1 mean motion resonance is basically governed by the eccentricity of Jupiter and its secular variations.  相似文献   

18.
We present a method for fast and accurate azimuth disambiguation of vector magnetogram data regardless of the location of the analyzed region on the solar disk. The direction of the transverse field is determined with the principle of minimum deviation of the field from the reference (potential) field. The new disambiguation (NDA) code is examined on the well-known models of Metcalf et al. (Solar Phys. 237, 267, 2006) and Leka et al. (Solar Phys. 260, 83, 2009), and on an artificial model based on the observed magnetic field of AR 10930 (Rudenko, Myshyakov, and Anfinogentov, Astron. Rep. 57, 622, 2013). We compare Hinode/SOT-SP vector magnetograms of AR 10930 disambiguated with three codes: the NDA code, the nonpotential magnetic-field calculation (NPFC: Georgoulis, Astrophys. J. Lett. 629, L69, 2005), and the spherical minimum-energy method (Rudenko, Myshyakov, and Anfinogentov, Astron. Rep. 57, 622, 2013). We then illustrate the performance of NDA on SDO/HMI full-disk magnetic-field observations. We show that our new algorithm is more than four times faster than the fastest algorithm that provides the disambiguation with a satisfactory accuracy (NPFC). At the same time, its accuracy is similar to that of the minimum-energy method (a very slow algorithm). In contrast to other codes, the NDA code maintains high accuracy when the region to be analyzed is very close to the limb.  相似文献   

19.
We identify 565 coronal mass ejections (CMEs) between January 2007 and December 2010 in observations from the twin STEREO/SECCHI/COR2 coronagraphs aboard the STEREO mission. Our list is in full agreement with the corresponding SOHO/LASCO CME Catalog ( http://cdaw.gsfc.nasa.gov/CME_list/ ) for events with angular widths of 45° and up. The monthly event rates behave similarly to sunspot rates showing a three- to fourfold rise between September 2009 and March 2010. We select 51 events with well-defined white-light structure and model them as three-dimensional (3D) flux ropes using a forward-modeling technique developed by Thernisien, Howard and Vourlidas (Astrophys. J. 652, 763??C?773, 2006). We derive their 3D properties and identify their source regions. We find that the majority of the CME flux ropes (82?%) lie within 30° of the solar equator. Also, 82?% of the events are displaced from their source region, to a lower latitude, by 25° or less. These findings provide strong support for the deflection of CMEs towards the solar equator reported in earlier observations, e.g. by Cremades and Bothmer (Astron. Astrophys. 422, 307??C?322, 2004).  相似文献   

20.
The Cassini spacecraft collects high resolution images of the Saturnian satellites and reveals the surface of these new worlds. Tiscareno et?al. succeeded to determine the Epimetheus rotation from the Cassini Imaging Science Subsystem data, initiating studies on the rotation of Epimetheus and its companion Janus (Tiscareno et?al., Icarus 204:254?C261, 2009; Noyelles, Icarus 207:887?C902, 2010; Robutel et?al., Icarus 211:758?C769, 2011). Especially, Epimetheus is characterized by its horseshoe shape orbit and the presence of the swap has to be introduced explicitly into rotational models. During its journey in the Saturnian system, Cassini spacecraft accumulates the observational data of the other satellites and it will be possible to determine the rotational parameters of several of them. To prepare these future observations, we built rotational models of the coorbital (also called Trojan) satellites Telesto, Calypso, Helene, and Polydeuces, in addition to Janus and Epimetheus. Indeed, Telesto and Calypso orbit around the L 4 and L 5 Lagrange points of Saturn-Tethys while Helene and Polydeuces are coorbital of Dione. The goal of this study is to understand how the departure from the Keplerian motion induced by the perturbations of the coorbital body, influences the rotation of these satellites. To this aim, we introduce explicitly the perturbation in the rotational equations by using the formalism developed by érdi (Celest Mech 15:367?C383, 1977) to represent the coorbital motions, and so we describe the rotational motion of the coorbitals, Janus and Epimetheus included, in compact form.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号