Photospheric Boundary Research Articles

Very Long-periodic Pulsations Detected Simultaneously in a White-light Flare and Sunspot Penumbra

We investigate the origin of very long-periodic pulsations in the white-light emission of an X6.4 flare on 2024 February 22 (SOL2024-02-22T22:08), which occurred at the edge of a sunspot group. The flare white-light fluxes reveal four successive and repetitive pulsations, which are simultaneously measured by the Helioseismic and Magnetic Imager and the White-light Solar Telescope. A quasi-period of 8.6−1.9+1.5 minutes, determined by the Morlet wavelet transform, is detected in the visible continuum channel. The modulation depth, which is defined as the ratio between the oscillatory amplitude and its long-term trend, is smaller than 0.1%, implying that the quasi-periodic pulsation (QPP) feature is a weak wave process. Imaging observations show that the X6.4 flare occurs near a sunspot group. Moreover, the white-light brightening is located in sunspot penumbra, and a similar quasi-period of about 8.5 −1.8+1.6 minutes is identified in one penumbral location of the nearest sunspot. The map of Fourier power distribution suggests that a similar periodicity is universally existing in most parts of the penumbra that is close to the penumbral–photospheric boundary. Our observations support the scenario that the white-light QPP is probably modulated by the slow-mode magnetoacoustic gravity wave leaking from the sunspot penumbra.

Formation of Polar Crown Filament Magnetic Fields by Supergranular Helicity Injection

To understand the magnetic fields of the polar crown filaments (PCFs) at high latitudes near polar regions of the Sun, we perform magnetofrictional numerical simulations on the long-term magnetic evolution of bipolar fields with roughly east–west polarity inversion lines (PILs) in a 3D spherical wedge domain near polar regions. The Coriolis-effect-induced vortical motions at the boundaries of several supergranular cells inject magnetic helicity from the photospheric boundary into the solar atmosphere. Supergranular-scale helicity injection, transfer, and condensation produce strongly sheared magnetic fields. Magnetic reconnections at footpoints of the sheared fields produce magnetic flux ropes (MFRs) with helicity signs consistent with the observed hemispheric helicity rule. The cross-sectional area of MFRs exhibits an uneven distribution, resembling a “foot-node-foot” periodic configuration. Experiments with different tilt directions of PILs indicate that the PCFs preferably form along PILs with the western end close to the polar region. The bending of PILs caused by supergranular flows, forming S-shape (Z-shape) PIL segments, promotes the formation of dextral (sinistral) MFRs. The realistic magnetic models we obtained can serve as starting points for the study of the plasma formation and eruption of PCFs.

Data-driven Radiative Magnetohydrodynamics Simulations with the MURaM Code

We present a method of conducting data-driven simulations of solar active regions and flux emergence with the MURaM radiative magnetohydrodynamics (MHD) code. The horizontal electric field that is derived from the full velocity and magnetic vectors is implemented at the photospheric (bottom) boundary to drive the induction equation. The energy equation accounts for thermal conduction along magnetic fields, optically thin radiative loss, and heating of coronal plasma by viscous and resistive dissipation, which allows for a realistic presentation of the thermodynamic properties of coronal plasma that are key to predicting the observational features of solar active regions and eruptions. To validate this method, the photospheric data from a comprehensive radiative MHD simulation of solar eruption (the ground truth) are used to drive a series of numerical experiments. The data-driven simulation reproduces the accumulation of free magnetic energy over the course of flux emergence in the ground truth with an error of 3%. The onset time is approximately 8 minutes delayed compared to the ground truth. However, a precursor-like signature can be identified at the correct onset time. The data-driven simulation captures key eruption-related emission features and plasma dynamics of the ground truth flare over a wide temperature span, from to . The evolution of the flare and coronal mass ejection as seen in synthetic extreme ultraviolet images is also reproduced with high fidelity. This method helps to understand the evolution of magnetic field in a more realistic coronal environment and to link the magnetic structures to observable diagnostics.

A Magnetogram-matching Method for Energizing Magnetic Flux Ropes Toward Eruption

We propose a new “helicity-pumping” method for energizing coronal equilibria that contain a magnetic flux rope (MFR) toward an eruption. We achieve this in a sequence of magnetohydrodynamics relaxations of small line-tied pulses of magnetic helicity, each of which is simulated by a suitable rescaling of the current-carrying part of the field. The whole procedure is “magnetogram-matching” because it involves no changes to the normal component of the field at the photospheric boundary. The method is illustrated by applying it to an observed force-free configuration whose MFR is modeled with our regularized Biot–Savart law method. We find that, in spite of the bipolar character of the external field, the MFR eruption is sustained by two reconnection processes. The first, which we refer to as breakthrough reconnection, is analogous to breakout reconnection in quadrupolar configurations. It occurs at a quasi-separator inside a current layer that wraps around the erupting MFR and is caused by the photospheric line-tying effect. The second process is the classical flare reconnection, which develops at the second quasi-separator inside a vertical current layer that is formed below the erupting MFR. Both reconnection processes work in tandem with the magnetic forces of the unstable MFR to propel it through the overlying ambient field, and their interplay may also be relevant for the thermal processes occurring in the plasma of solar flares. The considered example suggests that our method will be beneficial for both the modeling of observed eruptive events and theoretical studies of eruptions in idealized magnetic configurations.

Evolution of the critical torus instability height and coronal mass ejection likelihood in solar active regions

Aims. Working towards improved space weather predictions, we aim to quantify how the critical height at which the torus instability drives coronal mass ejections (CMEs) varies over time in a sample of solar active regions. Methods. We model the coronal magnetic fields of 42 active regions and quantify the critical height at their central polarity inversion lines throughout their observed lifetimes. We then compare these heights to the changing magnetic flux at the photospheric boundary and identify CMEs in these regions. Results. In our sample, the rates of CMEs per unit time are twice as high during phases when magnetic flux is increasing than when it is decreasing, and during those phases of increasing flux, the rate of CMEs is 63% higher when the critical height is rising than when it is falling. Furthermore, we support and extend the results of previous studies by demonstrating that the critical height in active regions is generally proportional to the separation of their magnetic polarities through time. When the separation of magnetic polarities in an active region increases, for example during the continuous emergence and expansion of a magnetic bipole, the critical height also tends to increase. Conversely, when the polarity separation decreases, for example due to the emergence of a new, compact bipole at the central inversion line of an existing active region or into a quiet-Sun environment, the critical height tends to decrease.

Reconstruction of Coronal Magnetic Fields Using a Poloidal–Toroidal Representation

A new method for reconstruction of coronal magnetic fields as force-free fields (FFFs) is presented. Our method employs poloidal and toroidal functions to describe divergence-free magnetic fields. This magnetic field representation naturally enables us to implement the boundary conditions at the photospheric boundary, i.e., the normal magnetic field and the normal current density there, in a straightforward manner. At the upper boundary of the corona, a source surface condition can be employed, which accommodates magnetic flux imbalance at the bottom boundary. Although our iteration algorithm is inspired by extant variational methods, it is nonvariational and requires far fewer iteration steps than most others. The computational code based on our new method is tested against the analytical FFF solutions by Titov & Démoulin. It is found to excel in reproducing a tightly wound flux rope, a bald patch, and quasi-separatrix layers with a hyperbolic flux tube.

Active region chromospheric magnetic fields

Context. A proper estimate of the chromospheric magnetic fields is thought to improve modelling of both active region and coronal mass ejection evolution. However, because the chromospheric field is not regularly obtained for sufficiently large fields of view, estimates thereof are commonly obtained through data-driven models or field extrapolations, based on photospheric boundary conditions alone and involving pre-processing that may reduce details and dynamic range in the magnetograms. Aims. We investigate the similarity between the chromospheric magnetic field that is directly inferred from observations and the field obtained from a magnetohydrostatic (MHS) extrapolation based on a high-resolution photospheric magnetogram. Methods. Based on Swedish 1-m Solar Telescope Fe I 6173 Å and Ca II 8542 Å observations of NOAA active region 12723, we employed the spatially regularised weak-field approximation (WFA) to derive the vector magnetic field in the chromosphere from Ca II, as well as non-local thermodynamic equilibrium (non-LTE) inversions of Fe I and Ca II to infer a model atmosphere for selected regions. Milne-Eddington inversions of Fe I serve as photospheric boundary conditions for the MHS model that delivers the three-dimensional field, gas pressure, and density self-consistently. Results. For the line-of-sight component, the MHS chromospheric field generally agrees with the non-LTE inversions and WFA, but tends to be weaker by 16% on average than these when larger in magnitude than 300 G. The observationally inferred transverse component is systematically stronger, up to an order of magnitude in magnetically weaker regions, but the qualitative distribution with height is similar to the MHS results. For either field component, the MHS chromospheric field lacks the fine structure derived from the inversions. Furthermore, the MHS model does not recover the magnetic imprint from a set of high fibrils connecting the main polarities. Conclusions. The MHS extrapolation and WFA provide a qualitatively similar chromospheric field, where the azimuth of the former is better aligned with Ca II 8542 Å fibrils than that of the WFA, especially outside strong-field concentrations. The amount of structure as well as the transverse field strengths are, however, underestimated by the MHS extrapolation. This underscores the importance of considering a chromospheric magnetic field constraint in data-driven modelling of active regions, particularly in the context of space weather predictions.

A Model of Homologous Confined and Ejective Eruptions Involving Kink Instability and Flux Cancellation

In this study, we model a sequence of a confined and a full eruption, employing the relaxed end state of the confined eruption of a kink-unstable flux rope as the initial condition for the ejective one. The full eruption, a model of a coronal mass ejection, develops as a result of converging motions imposed at the photospheric boundary, which drive flux cancellation. In this process, parts of the positive and negative external flux converge toward the polarity inversion line, reconnect, and cancel each other. Flux of the same amount as the canceled flux transfers to a flux rope, increasing the free magnetic energy of the coronal field. With sustained flux cancellation and the associated progressive weakening of the magnetic tension of the overlying flux, we find that a flux reduction of ≈11% initiates the torus instability of the flux rope, which leads to a full eruption. These results demonstrate that a homologous full eruption, following a confined one, can be driven by flux cancellation.

3D MHD wave propagation near a coronal null point: New wave mode decomposition approach

Context. Ubiquitous vortex flows at the solar surface excite magnetohydrodynamic (MHD) waves that propagate to higher layers of the solar atmosphere. In the solar corona, these waves frequently encounter magnetic null points. The interaction of MHD waves with a coronal magnetic null in realistic 3D setups requires an appropriate wave identification method. Aims. We present a new MHD wave decomposition method that overcomes the limitations of existing wave identification methods. Our method allows for an investigation of the energy fluxes in different MHD modes at different locations of the solar atmosphere as waves generated by vortex flows travel through the solar atmosphere and pass near the magnetic null. Methods. We used the open-source MPI-AMRVAC code to simulate wave dynamics through a coronal null configuration. We applied a rotational wave driver at our bottom photospheric boundary to mimic vortex flows at the solar surface. To identify the wave energy fluxes associated with different MHD wave modes, we employed a wave decomposition method that is able to uniquely distinguish different MHD modes. Our proposed method utilizes the geometry of an individual magnetic field-line in the 3D space to separate the velocity perturbations associated with the three fundamental MHD waves. We compared our method with an existing wave decomposition method that uses magnetic flux surfaces instead. Over the selected flux surfaces, we calculated and analyzed the temporally averaged wave energy fluxes, as well as the acoustic and magnetic energy fluxes. Our wave decomposition method allowed us to estimate the relative strengths of individual MHD wave energy fluxes. Results. Our method for wave identification is consistent with previous flux-surface-based methods and provides the expected results in terms of the wave energy fluxes at various locations of the null configuration. We show that ubiquitous vortex flows excite MHD waves that contribute significantly to the Poynting flux in the solar corona. Alfvén wave energy flux accumulates on the fan surface and fast wave energy flux accumulates near the null point. There is a strong current density buildup at the spine and fan surface. Conclusions. The proposed method has advantages over previously utilized wave decomposition methods, since it may be employed in realistic simulations or magnetic extrapolations, as well as in real solar observations whenever the 3D fieldline shape is known. The essential characteristics of MHD wave propagation near a null – such as wave energy flux accumulation and current buildup at specific locations – translate to the more realistic setup presented here. The enhancement in energy flux associated with magneto-acoustic waves near nulls may have important implications in the formation of jets and impulsive plasma flows.

Data-driven, time-dependent modeling of pre-eruptive coronal magnetic field configuration at the periphery of NOAA AR 11726

Context. Data-driven, time-dependent magnetofrictional modeling has proved to be an efficient tool for studying the pre-eruptive build-up of energy for solar eruptions, and sometimes even the ejection of coronal flux ropes during eruptions. However, previous modeling works have illustrated the sensitivity of the results on the data-driven boundary condition, as well as the difficulty in modeling the ejections with proper time scales. Aims. We aim to study the pre- and post-eruptive evolution of a weak coronal mass ejection producing eruption at the periphery of isolated NOAA active region (AR) 11726 using a data-driven, time-dependent magnetofrictional simulation, and aim to illustrate the strengths and weaknesses of our simulation approach. Methods. We used state-of-the-art data processing and electric field inversion methods to provide the data-driven boundary condition for the simulation. We analyzed the field-line evolution, magnetic connectivity, twist, as well as the energy and helicity budgets in the simulation to study the pre- and post-eruptive magnetic field evolution of the observed eruption from AR11726. Results. We find the simulation to produce a pre-eruptive flux rope system consistent with several features in the extreme ultraviolet and X-ray observations of the eruption, but the simulation largely fails to reproduce the ejection of the flux rope. We find the flux rope formation to be likely driven by the photospheric vorticity at one of the footpoints, although reconnection at a coronal null-point may also feed poloidal flux to the flux rope. The accurate determination of the non-inductive (curl-free) component of the photospheric electric field boundary condition is found to be essential for producing the flux rope in the simulation. Conclusions. Our results illustrate the applicability of the data-driven, time-dependent magnetofrictional simulations in modeling the pre-eruptive evolution and formation process of a flux rope system, but they indicate that the modeling output becomes problematic for the post-eruptive times. For the studied event, the flux rope also constituted only a small part of the related active region.

Reflection and Evolution of Torsional Alfvén Pulses in Zero-beta Flux Tubes

Abstract We model the behavior of a torsional Alfvén pulse, assumed to propagate through the chromosphere. Building on our existing model, we utilize the zero-beta approximation appropriate for plasma in an intense magnetic flux tube, e.g., a magnetic bright point. The model is adapted to investigate the connection between these features and chromospheric spicules. A pulse is introduced at the lower, photospheric boundary of the tube as a magnetic shear perturbation, and the resulting propagating Alfvén waves are reflected from an upper boundary, representing the change in density found at the transition region. The induced upward mass flux is followed by the reversal of the flux that may be identified with the rising and falling behavior of certain lower solar atmospheric jets. The ratio of the transmitted and reflected mass flux is estimated and compared with the relative total mass of spicules and the solar wind. An example is used to study the properties of the pulse. We also find that the interaction between the initial and reflected waves may create a localized flow that persists independently from the pulse itself.

Reconstructing the Coronal Magnetic Field: The Role of Cross-field Currents in Solution Uniqueness

Abstract We present a new 3D magnetohydrostatic (MHS) direct elliptic solver for extrapolating the coronal magnetic field from photospheric boundary conditions in a manner consistent with an assumed plasma distribution. We use it to study the uniqueness of the reconstructed magnetic field as a function of how significant the plasma forcing is on the force balance of the magnetic field. To this end, we consider an analytic MHS model as ground truth. The model uses two free parameters to decompose the current into two parts: a magnetic-field-aligned component and a cross-field component. We perform a comprehensive study of the 2D parameter space to understand under what conditions the ground truth can be reproduced uniquely. We find that current oriented perpendicular to the magnetic field has a smaller solution space than the same amount of current oriented parallel to the magnetic field, and so MHS regimes with larger proportions of plasma-related forcing may be a promising avenue toward finding unique magnetic field reconstructions.

Self-consistent Nonlinear Force-Free Field Reconstruction from Weighted Boundary Conditions

Vector magnetogram data are often used as photospheric boundary conditions for force-free coronal magnetic field extrapolations. In general, however, vector magnetogram data are not consistent with the force-free assumption. In this article, we demonstrate a way to deal with inconsistent boundary data, by generalizing the "self-consistency procedure" of Wheatland & Regnier (2009). In that procedure, the inconsistency is resolved by an iterative process of constructing two solutions based on the values of the force-free parameter alpha on the two polarities of the field in the boundary (the P and N polarities), and taking uncertainty-weighted averages of the boundary alpha values in the P and N solutions. When the alpha values in the P and N regions are very different, the self-consistent solution may lose high alpha values from the boundary conditions. We show how, by altering the weighting of the uncertainties in the P or N boundary conditions, we can preserve high alpha values in the self-consistent solution. The weighted self-consistent extrapolation method is demonstrated on an analytic bipole field and applied to vector magnetogram data taken by the Helioseismic and Magnetic Imager (HMI) instrument for NOAA active region AR 12017 on 2014 March 29.

A distinct magnetic property of the inner penumbral boundary

Context. Analyses of sunspot observations revealed a fundamental magnetic property of the umbral boundary: the invariance of the vertical component of the magnetic field. Aims. We analyse the magnetic properties of the umbra-penumbra boundary in simulated sunspots and thus assess their similarity to observed sunspots. We also aim to investigate the role of the plasma β and the ratio of kinetic to magnetic energy in simulated sunspots in the convective motions because these quantities cannot be reliably determined from observations. Methods. We used a set of non-gray simulation runs of sunspots with the MURaM code. The setups differed in terms of subsurface magnetic field structure and magnetic field boundary imposed at the top of the simulation domain. These data were used to synthesize the Stokes profiles, which were then degraded to the Hinode spectropolarimeter-like observations. Then, the data were treated like real Hinode observations of a sunspot, and magnetic properties at the umbral boundaries were determined. Results. Simulations with potential field extrapolation produce a realistic magnetic field configuration on the umbral boundaries of the sunspots. Two simulations with a potential field upper boundary, but different subsurface magnetic field structures, differ significantly in the extent of their penumbrae. Increasing the penumbra width by forcing more horizontal magnetic fields at the upper boundary results in magnetic properties that are not consistent with observations. This implies that the size of the penumbra is given by the subsurface structure of the magnetic field, that is, by the depth and inclination of the magnetopause, which is shaped by the expansion of the sunspot flux rope with height. None of the sunspot simulations is consistent with the observed properties of the magnetic field and the direction of the Evershed flow at the same time. Strong outward-directed Evershed flows are only found in setups with an artificially enhanced horizontal component of the magnetic field at the top boundary that are not consistent with the observed magnetic field properties at the umbra-penumbra boundary. We stress that the photospheric boundary of simulated sunspots is defined by a magnetic field strength of equipartition field value.

Torsional Alfvénic Oscillations Discovered in the Magnetic Free Energy during Solar Flares

Abstract We report the discovery of torsional Alfvénic oscillations in solar flares, which modulate the time evolution of the magnetic free energy E f (t), while the magnetic potential energy E p (t) is uncorrelated, and the nonpotential energy varies as E np (t) = E p + E f (t). The mean observed time period of the torsional oscillations is P obs = 15.1 ± 3.9 minutes, the mean field line length is L = 135 ± 35 Mm, and the mean phase speed is v phase = 315 ± 120 km s−1, which we interpret as torsional Alfvénic waves in flare loops with enhanced electron densities. Most of the torsional oscillations are found to be decay-less, but exhibit a positive or negative trend in the evolution of the free energy, indicating new emerging flux (if positive), magnetic cancellation, or flare energy dissipation (if negative). The time evolution of the free energy has been calculated in this study with the Vertical-current Approximation (Version 4) Non-linear Force-free Field code, which incorporates automatically detected coronal loops in the solution and bypasses the non-force-freeness of the photospheric boundary condition, in contrast to traditional NLFFF codes.

Force-free Field Reconstructions Enhanced by Chromospheric Magnetic Field Data

Abstract A 3D picture of the coronal magnetic field remains an outstanding problem in solar physics, particularly in active regions. Nonlinear force-free field reconstructions that employ routinely available full-disk photospheric vector magnetograms represent state-of-the-art coronal magnetic field modeling. Such reconstructions, however, suffer from an inconsistency between a force-free coronal magnetic field and a non-force-free photospheric boundary condition, from which the coronal reconstruction is performed. In this study we focus on integrating the additional chromospheric and/or coronal magnetic field data with the vector photospheric magnetograms with the goal of improving the reliability of the magnetic field reconstructions. We develop a corresponding modification of the available optimization codes described in Fleishman et al. and test their performance using a full-fledged magnetohydrodynamics model obtained from the Bifrost code by performing a “voxel-by-voxel” comparison between the reconstructed and the model magnetic fields. We demonstrate that adding even an incomplete set of chromospheric magnetic field data can measurably improve the reconstruction of the coronal magnetic field and greatly improve reconstructions of the magnetic connectivity and of the coronal electric current.

Magnetic Structures at the Boundary of the Closed Corona: Interpretation of S-Web Arcs

Abstract The topology of coronal magnetic fields near the open-closed magnetic flux boundary is important to the the process of interchange reconnection, whereby plasma is exchanged between open and closed flux domains. Maps of the magnetic squashing factor in coronal field models reveal the presence of the Separatrix-Web (S-Web), a network of separatrix surfaces and quasi-separatrix layers, along which interchange reconnection is highly likely. Under certain configurations, interchange reconnection within the S-Web could potentially release coronal material from the closed magnetic field regions to high-latitude regions far from the heliospheric current sheet, where it is observed as slow solar wind. It has also been suggested that transport along the S-Web may be a possible cause for the observed large longitudinal spreads of some impulsive, 3He-rich solar energetic particle events. Here, we demonstrate that certain features of the S-Web reveal structural aspects of the underlying magnetic field, specifically regarding the arcing bands of highly squashed magnetic flux observed at the outer boundary of global magnetic field models. In order for these S-Web arcs to terminate or intersect away from the helmet streamer apex, there must be a null spine line that maps a finite segment of the photospheric open-closed boundary up to a singular point in the open flux domain. We propose that this association between null spine lines and arc termination points may be used to identify locations in the heliosphere that are preferential for the appearance of solar energetic particles and plasma from the closed corona, with characteristics that may inform our understanding of interchange reconnection and the acceleration of the slow solar wind.

Threshold of Non-potential Magnetic Helicity Ratios at the Onset of Solar Eruptions

Abstract The relative magnetic helicity is a quantity that is often used to describe the level of entanglement of non-isolated magnetic fields, such as the magnetic field of solar active regions. The aim of this paper is to investigate how different kinds of photospheric boundary flows accumulate relative magnetic helicity in the corona and if and how well magnetic-helicity-related quantities identify the onset of an eruption. We use a series of three-dimensional, parametric magnetohydrodynamic simulations of the formation and eruption of magnetic flux ropes. All the simulations are performed on the same grid, using the same parameters, but they are characterized by different driving photospheric flows, i.e., shearing, convergence, stretching, and peripheral- and central- dispersion flows. For each of the simulations, the instant of the onset of the eruption is carefully identified by using a series of relaxation runs. We find that magnetic energy and total relative helicity are mostly injected when shearing flows are applied at the boundary, while the magnetic energy and helicity associated with the coronal electric currents increase regardless of the kind of photospheric flows. We also find that, at the onset of the eruptions, the ratio between the non-potential magnetic helicity and the total relative magnetic helicity has the same value for all the simulations, suggesting the existence of a threshold in this quantity. Such a threshold is not observed for other quantities as, for example, those related to the magnetic energy.

Comparison of methods for modelling coronal magnetic fields

Aims. Four different approximate approaches used to model the stressing of coronal magnetic fields due to an imposed photospheric motion are compared with each other and the results from a full time-dependent magnetohydrodynamic (MHD) code. The assumptions used for each of the approximate methods are tested by considering large photospheric footpoint displacements. Methods. We consider a simple model problem, comparing the full non-linear MHD, determined with the Lare2D numerical code, with four approximate approaches. Two of these, magneto-frictional relaxation and a quasi-1D Grad-Shafranov approach, assume sequences of equilibria, whilst the other two methods, a second-order linearisation of the MHD equations and Reduced MHD, are time dependent. Results. The relaxation method is very accurate compared to full MHD for force-free equilibria for all footpoint displacements, but has significant errors when the plasma β0 is of order unity. The 1D approach gives an extremely accurate description of the equilibria away from the photospheric boundary layers, and agrees well with Lare2D for all parameter values tested. The linearised MHD equations correctly predict the existence of photospheric boundary layers that are present in the full MHD results. As soon as the footpoint displacement becomes a significant fraction of the loop length, the RMHD method fails to model the sequences of equilibria correctly. The full numerical solution is interesting in its own right, and care must be taken for low β0 plasmas if the viscosity is too high.

Magnetohydrodynamic Simulations for Studying Solar Flare Trigger Mechanism

Abstract In order to understand the flare trigger mechanism, we conduct three-dimensional magnetohydrodynamic simulations using a coronal magnetic field model derived from data observed by the Hinode satellite. Several types of magnetic bipoles are imposed into the photospheric boundary of the Nonlinear Force-free Field model of Active Region (AR) NOAA 10930 on 2006 December 13, to investigate what kind of magnetic disturbance may trigger the flare. As a result, we confirm that certain small bipole fields, which emerge into the highly sheared global magnetic field of an AR, can effectively trigger a flare. These bipole fields can be classified into two groups based on their orientation relative to the polarity inversion line: the so-called opposite polarity, and reversed shear structures, as suggested by Kusano et al. We also investigate the structure of the footpoints of reconnected field lines. By comparing the distribution of reconstructed field lines and observed flare ribbons, the trigger structure of the flare can be inferred. Our simulation suggests that the data-constrained simulation, taking into account both the large-scale magnetic structure and small-scale magnetic disturbance (such as emerging fluxes), is a good way to discover a flare-producing AR, which can be applied to space weather prediction.