the interest in studying record events in a general sense is evident as they by definition tend to be rare events of a special nature . as examples of record - breaking events
most readers might think about athletic records like a world record in 100 m dash , or long jump @xcite .
however records also play an important role in other naturally occurring and man - made phenomena .
in particular , climatic records attract much attention and are frequently mentioned in the media in the context of global warming .
the statistics of record temperatures , precipitation amounts and floods have been adressed extensively in recent work @xcite .
also in evolutionary biology record events have a certain relevance as was pointed out in @xcite , and they play an important role in the dynamics of disordered systems like spin glasses and dirty superconductors @xcite .
last but not least records can be studied in the financial sciences , e.g. when considering record - breaking values or changes in stock prices @xcite . in all these examples not
only the rate of record events is of interest but also their correlations : given that a record flood just occurred , can one relax because it will take a while before the next record flood is to be expected ?
based on our previous work in @xcite , this question will be addressed in the present paper .
most previous studies of record events and record values @xcite have been restricted to the case where the underlying random variables ( rv s ) @xmath0 are independently drawn from a common distribution @xmath1 ( i.i.d .
rv s ) . in this case
, record events can be shown to be stochastically independent , i.e. the knowledge about previous record events does not indicate whether a future record event is more or less likely .
however in many applications , the underlying random variables are not i.i.d .
for example , in recent work majumdar and ziff computed the statistics of records in a symmetric random walk @xcite ( see also @xcite ) . here , record events are obviously not independent : in a random walk the probability for a second record to occur immediately after after a previous record is simply @xmath2 , whereas the unconditional probability for a record in the @xmath3th step is @xmath4 @xcite .
in this article we consider correlations between records in sequences of rv s that are independent but not identically distributed .
an important example is the statistics of hot or cold records in weather data . in @xcite
the occurrence of record - breaking temperatures in european and american weather recordings was analyzed by modeling daily temperature measurements as uncorrelated rv s with a linear drift .
it was shown that the increase of hot records and the decrease of cold records due to global warming can be described and quantified by this model .
we will refer to this particular example of non - i.i.d .
rv s , which was originally introduced in the 1980 s in the context of athletic records @xcite , as the linear drift model ( ldm ) .
the precise definition of the ldm is as follows : for a positive constant @xmath5 , each random variable @xmath6 has an i.i.d .
part @xmath7 with density @xmath1 and an additive linear drift such that @xmath8 in other words , while the shape of the density @xmath9 of the @xmath10 rv is invariant , its overall position is shifted by a constant amount @xmath5 each time a new rv is drawn . as was already noted in @xcite , the record events under this model
are not in general independent ( see also @xcite for the related case of rv s with changing variance ) . below we shall see that the sign and magnitude of correlations depend in a highly non - trivial way on the underlying probability density @xmath1 . to quantify the correlations , consider the joint probability @xmath11 of a record occurring in the @xmath12 as well as in the @xmath13 step , @xmath14.\end{aligned}\ ] ] to normalize this joint probability , we divide by the probabilities @xmath15 and @xmath16 of records in the @xmath13 and @xmath17 steps , obtaining the key quantity under consideration throughout this paper , @xmath18     correlations between record events in the @xmath19 and @xmath20 rv for different distributions of the underlying i.i.d .
part . for the gumbel case ,
the well known stochastic independence of record events is found , while simulations for the other distributions show that generally the records are correlated .
note in particular the curve for the symmetric lvy - stable distribution with index @xmath21 , which indicates an attraction between record events . ]    the simulations shown in fig .
[ rec_dep ] illustrate two known results , the fact that record events for i.i.d .
rvs are stochastically independent @xcite ( all curves meet at unity for @xmath22 ) , and that for the ldm with gumbel distributed i.i.d . part , record events are stochastically independent for all values of @xmath5 @xcite .
however the curves for other distributions show that stochastic independence does not hold in general .
for some choices of the underlying probability density @xmath1 record events tend to repel each other ( @xmath23 ) , wheras in other cases the correlations are positive , i.e. record events attract each other ( @xmath24 ) . in particular this last , somewhat counterintuitive numerical finding triggered our interest in this problem .    in the next section
we will derive a general expression for the correlations between record events as quantified by @xmath25 .
for this we employ methods similar to those used in @xcite . in section [ examples ] , we provide some explicit examples illustrating the rich and counterintuitive behavior of the correlations . in particular , we discuss the behavior of correlations with respect to the three classes of extreme value statistics @xcite .
we find that while in the weibull class of distributions with finite support the correlations are generally negative , positive correlations mainly arise when the rv s are in the frecht class of distributions with power law tails . in the intermediate gumbel class of distributions with exponential - like tails correlations
can be both positive and negative . finally in section [ conclusions
] , we summarize our results and discuss possible applications and open questions .
consider a sequence @xmath26 of independent but not necessarily identically distributed rv s , with @xmath27 denoting the probability density of @xmath28 .
the probability that the maximum of the first @xmath3 rv s is smaller than a given value @xmath29 then factorizes into the probabilities of the rv s individually being less than @xmath29 , @xmath30=\prod_{l=1}^{n}\int^x{\mathrm{d}}x ' f_{l}(x')=\prod_{l=1}^{n}f_l(x),\ ] ] where @xmath31 denotes the cumulative distribution function for the density @xmath27 , and the omitted lower ( upper ) integral boundary is understood to be the lower ( upper ) bound of the support of @xmath32 respectively . throughout
we assume distribution functions to be continuous .    using eq .
[ factor ] the probability @xmath15 that the @xmath13 rv is a record can be written as @xmath33 and the joint probability that both @xmath34 and @xmath35 are records is @xmath36 for the ldm we have @xmath37 and @xmath38 , where @xmath1 and @xmath39 are the density and distribution function of the i.i.d .
part @xmath28 in ( [ ldm ] ) .
it then follows that @xcite @xmath40 and similarly @xmath41      the expressions ( [ joint ] ) and ( [ joint_ldm ] ) allow us to verify two known results : both for the i.i.d .
case ( @xmath22 ) and for the ldm with gumbel distributed i.i.d .
part @xmath28 , record events are independent .
we start with the i.i.d .
first note that the probability that the @xmath13 rv is a record is the probability that it is the largest among @xmath3 equivalent rv s , so by symmetry @xmath42 . for the i.i.d .
case , the subscripts in ( [ joint ] ) can be dropped . using the substitution @xmath43 and @xmath44 , one obtains @xmath45 it is possible to insert an arbitrary spacing between the two record events and still obtain this factorization to yield @xmath46 .
similarly , one can consider more than one record event .
the factorization property implies that in the i.i.d .
case , record events are stochastically independent @xcite .
next we consider the ldm with gumbel distributed i.i.d .
part , @xmath47 and @xmath48 .
using that in this case @xcite @xmath49 the substitution @xmath43 used above yields @xcite @xmath50 with @xmath51 , and @xmath52 again one can introduce an arbitrary number of spaces and intermediate records and still have the corresponding joint probability factorize .
this implies that for the ldm with gumbel distributed i.i.d .
part , as for arbitrary distributions without drift , record events are stochastically independent .      as shown above , for the i.i.d .
case @xmath53 , record events are independent . on the other hand , for very large drift , almost every rv is a record , and thus trivially @xmath54 for @xmath55 ( see fig .
[ rec_dep ] ) .
the behavior of @xmath56 in between these two limits is not as easily characterized . in particular , while for some probability densities the records seem to repel each other for intermediate @xmath5 ( @xmath57 ) , for other probability densities the records seem to _ attract _ ( @xmath24 )
.     correlations between record events as the number of rv s @xmath3 increases . for the i.i.d .
part drawn from a gaussian distribution with variance @xmath58 , the correlations are negative and for sufficiently large @xmath3 tend to a limiting form . for the symmetric lvy distributed i.i.d .
part , the correlations are positive and increase with @xmath3 . ]
figure [ no_fse ] shows examples for both types of behavior , illustrating in addition how the correlations evolve with increasing @xmath3 for the cases of lvy and gauss distributed i.i.d
. parts . clearly , the correlations are not a finite size effect , but their @xmath3-dependence is markedly different in the two cases . for the gaussian case ,
the correlations are negative and seem to approach a form independent of @xmath3 . for the lvy case , on the other hand , the correlations seem to increase with @xmath3
. we will argue later that the correlations eventually become @xmath3-independent also for heavy - tailed distributions , but this limit is approached much more slowly than in the gaussian case .
a first step towards understanding the behavior of the correlations is to determine the sign with which the curves in figs.[rec_dep ] and [ no_fse ] depart from the i.i.d .
case @xmath22 . to address this question
, we compute the record rate @xmath59 and the joint probability @xmath60 in the limit of small @xmath5 and @xmath61 , where @xmath62 is the standard deviation or some other measure of the width of the distribution . for the record rate
this problem was discussed extensively in @xcite ( see also @xcite ) . knowing that @xmath42 , a taylor expansion for small @xmath5 yields @xmath63 where we have defined the quantity used in @xcite .
the integral ( [ i(n ) ] ) also appears in the analysis of the density of near - extreme events @xcite . ]
@xmath64 in the same small @xmath5 regime we can also easily determine the joint probability of occurrence of two consecutive records @xmath65 . to leading order in @xmath5 , the expansion of eq.([joint_ldm ] ) yields @xmath66 the first term can be evaluated by partial integration to @xmath67 , and another partial integration yields @xmath68 in addition to the zeroth order term @xmath67 , we find a positive and a negative correction term that describe the effect of the drift .
depending on which one of these two terms is larger , positive or negative correlations will result .
a similar expansion can be set up for the normalized correlation @xmath25 defined in eq.([lnn ] ) .
writing @xmath69 a straightforward computation yields @xmath70\nonumber\\     & & -\frac{n^2(n-1 ) } { 2}i(n-2)-\frac{(n-1)^2(n-2)}{2}i(n-3 )   \end{aligned}\ ] ] again , both positive and negative terms occur in ( [ taylor_mid ] ) , the relative magnitude of which depends explicitly on the underlying probability density @xmath1 .
thus from ( [ taylor_mid ] ) , it is not at all obvious how the counterintuitive positive correlations observed numerically arise , how their occurrence depends on the properties of @xmath71 and how the clustering of record events behaves as a function of the number @xmath3 of rv s considered . before turning to the explicit evaluation of @xmath72 in section [ examples ] ,
we provide some heuristics on the asymptotics of the correlations in the limit of very large @xmath3 , and address the behavior of correlations between rv s that are more than one time step apart .
it is clear from eqs .
( [ pn_ldm ] ) and ( [ joint_ldm ] ) that the asymptotic behavior of the record rate @xmath73 and the joint probability @xmath74 hinges on the existence of the function @xmath75 which was rigorously examined by ballerini and resnick . in @xcite , they prove that the limit ( [ g_limit ] ) exists and is nonzero whenever the density @xmath1 has a finite first moment . for completeness
we provide a heuristic version of their argument .
taking the logarithm of eq .
( [ g_limit ] ) one has to consider the convergence of the series @xmath76 now @xmath77 and for any finite @xmath29 and @xmath78 there is a @xmath79 large enough such that @xmath80 .
therefore @xmath81 where @xmath82 is a finite sum and the logarithms in the remaining series were approximated by their taylor expansion around unity .
the condition for @xmath1 to have a finite first moment is that @xmath83 decays faster than @xmath84 , which implies convergence of the infinite series in ( [ rest ] ) and hence of @xmath85 ( in the opposite case the series diverges to @xmath86 and @xmath87 ) .
since @xmath85 is a probability , it is bounded from above by unity and it follows from ( [ pn_ldm ] ) and ( [ joint_ldm ] ) that @xmath88 and @xmath89 have finite , non - zero limits for @xmath90 .
the same statement then applies to the ratio ( [ lnn ] ) .
although the focus of this paper is on the case of a positive drift which enhances the occurrence of ( upper ) records , it is instructive to compare the asymptotic behavior for @xmath78 described above to the case @xmath91 . for negative drift
both the record rate @xmath73 and the joint probability @xmath74 vanish for @xmath90 , and they do so more rapidly than their i.i.d .
counterparts . the asymptotic behavior of the ratio ( [ lnn ] ) is then _ a priori _ undetermined .
we consider a simple example where the quantities of interest can be computed explicitly .
let the probability density of the i.i.d .
variables be given by the negative exponential distribution , @xmath92 and @xmath93 elsewhere .
then for @xmath91 the integrand in ( [ pn_ldm ] ) and ( [ joint_ldm ] ) is of the form @xmath94 , \;\;\ ; x \leq 0,\ ] ] and integration yields @xmath95 , \ ; p_{n , n-1}(c ) = \frac{1}{n(n-1)}\exp\left[\frac{c}{2}n(n-1)\right].\ ] ] for @xmath91 the record rate is suppressed superexponentially , such that the expected number of records remains finite for @xmath90 ( see also @xcite )
. however , forming the ratio ( [ lnn ] ) we see that @xmath96\ ] ] which , in contrast to the case of positive drift , grows without bound for @xmath90 .
we will see in section [ examples ] that such a behavior is not restricted to this specific example .
the non - existence of the @xmath90 limit for @xmath91 suggests that the limiting function @xmath97 may not be smooth for @xmath98 , as will be explicitly demonstrated in section [ examples ] .      before embarking on the study of specific distributions , we briefly comment on the behavior of the correlations as a function of the distance between record events . denoting by @xmath99 the joint probability for observing a record in the @xmath100th and the @xmath3th event , it was shown in @xcite that @xmath101 ^ 2 \equiv p(c)^2,\ ] ] which implies that record events become uncorrelated for @xmath102
this leads to the expectation that in a finite time series , the correlations between record events at distance @xmath103 should decay with @xmath103 , in agreement with the simulations shown in fig . [ dist_corr ] .
the correlations are seen to be maximal at @xmath104 . in the present paper we therefore focus on the case of nearest neighbor correlations , although extending the computations to @xmath105 is in principle straightforward .
correlations between record events at distance @xmath103 from time series of length @xmath106 for gaussian ( @xmath58 ) , lvy - stable ( @xmath107 ) , uniform ( on @xmath108 $ ] ) and pareto ( @xmath107 ) distributions .
* left : * full @xmath5 dependence of @xmath109 . *
right : * at the drift velocity @xmath110 where the correlations are maximal ( @xmath111 for gaussian and uniform , @xmath112 for the lvy and @xmath113 for the pareto case ) , @xmath114 is shown as function of @xmath103 to illustrate the decay of correlations.,title="fig:",scaledwidth=48.0% ]   correlations between record events at distance @xmath103 from time series of length @xmath106 for gaussian ( @xmath58 ) , lvy - stable ( @xmath107 ) , uniform ( on @xmath108 $ ] ) and pareto ( @xmath107 ) distributions . *
left : * full @xmath5 dependence of @xmath109 .
* right : * at the drift velocity @xmath110 where the correlations are maximal ( @xmath111 for gaussian and uniform , @xmath112 for the lvy and @xmath113 for the pareto case ) , @xmath114 is shown as function of @xmath103 to illustrate the decay of correlations.,title="fig:",scaledwidth=48.0% ]
whereas the statistics of record events in sequences of i.i.d .
rv s is completely universal @xcite , in the presence of drift one has to distinguish between distributions belonging to the three different universality classes of extreme value theory , the classes of weibull , gumbel and frchet @xcite . for each of these classes we will analyze the correlations in the ldm for a few exemplary distributions , starting with the weibull class , and summarize the observed behavior in a unifying scaling picture in section [ sec : unified ] .
the results presented below rely on the study of the record rate @xmath115 in the ldm presented in @xcite , as well as on the expansions for @xmath65 [ eq .
( [ pnm ] ) ] and @xmath116 [ eqs .
( [ taylor_end],[taylor_mid ] ) ] derived above in section [ sec : expansion ] .       *
left : * numerical simulations of @xmath25 for a uniform distribution of width @xmath117 . for each @xmath5 and @xmath3
we simulated @xmath118 time series .
the figure shows data for @xmath119 and @xmath120 . at @xmath121
we find a steep descent of @xmath25 . for @xmath122 and @xmath123
we show the analytic prediction ( [ prediction_uniform ] ) ( steep lines ) . with increasing @xmath5 , @xmath25 appears to become completely independent of @xmath3 . *
right : * simulation results for @xmath124 as a function of @xmath3 for @xmath125 and @xmath126 .
we simulated @xmath127 series of @xmath128 rv s in each case . for @xmath91
we find a very steep increase of @xmath124 that is in good agreement with the analytical prediction ( [ prediction_uniform ] )
. for @xmath78 , @xmath124 quickly saturates at a constant value and our analytical results are not very useful.,title="fig:",scaledwidth=50.0% ] * left : * numerical simulations of @xmath25 for a uniform distribution of width @xmath117 . for each @xmath5 and @xmath3
we simulated @xmath118 time series .
the figure shows data for @xmath119 and @xmath120 . at @xmath121
we find a steep descent of @xmath25 . for @xmath122 and @xmath123
we show the analytic prediction ( [ prediction_uniform ] ) ( steep lines ) . with increasing @xmath5 , @xmath25 appears to become completely independent of @xmath3 . *
right : * simulation results for @xmath124 as a function of @xmath3 for @xmath125 and @xmath126 .
we simulated @xmath127 series of @xmath128 rv s in each case .
for @xmath91 we find a very steep increase of @xmath124 that is in good agreement with the analytical prediction ( [ prediction_uniform ] ) .
for @xmath78 , @xmath124 quickly saturates at a constant value and our analytical results are not very useful.,title="fig:",scaledwidth=50.0% ]    the weibull class is the class of distributions with bounded support . as a very simple member of the weibull class we consider a uniform distribution on the interval @xmath129 $ ] .
we compute @xmath116 in the small @xmath5 regime by making use of ( [ taylor_end],[taylor_mid ] ) and find for @xmath130 @xmath131 the correlations are negative and their magnitude increases rapidly with @xmath3 ( but note that the leading order approximation is valid only for @xmath132 ) .
we performed numerical simulations to check this result with @xmath117 and @xmath133 .
figure [ fig : uniform ] shows how the descent of @xmath116 for @xmath78 steepens with increasing @xmath3 . in the figure
we also plot our prediction ( [ prediction_uniform ] ) for small values of @xmath3 , however the steepness of the descent at zero does not allow us to estimate the quality of our approximation reliably .
the figure shows that for large enough @xmath3 and not too small @xmath5 , @xmath116 eventually becomes independent of @xmath3 , as expected from the general arguments given in section [ sec : largen ] . in this case
the approach to the @xmath90 limit is particularly rapid , because distributions that are more than @xmath134 steps apart do not overlap .
the sharp initial drop of @xmath56 indicates that the limiting distribution @xmath135 is discontinuous at @xmath121 , which is consistent with the non - existence of a limiting function ) is equivalent to the uniform distribution , and therefore the argument of section [ sec : largen ] applies here as well . ] for @xmath91 . as a more general subset of the weibull class we next consider the class of distributions defined by @xmath136 with @xmath137 and @xmath138 . for @xmath139 this produces a uniform distribution similar to the one used above .
in @xcite we found that for @xmath140 the integral in ( [ i(n ) ] ) is given by @xmath141 using this and assuming @xmath142 we find the following expression for @xmath116 : @xmath143 and making use of the stirling approximation we finally obtain @xmath144 similar to the uniform case , the correlations between neighboring record events are always negative and their magnitude increases rapidly with increasing @xmath3 , suggesting a discontinuity of the limiting function @xmath135 at @xmath121 .
as a first example in the gumbel class we consider the exponential distribution with mean @xmath145 , @xmath146 , @xmath147 .
in @xcite it was found that in this case the effect of the linear drift on the record rate is independent of @xmath3 to leading order in @xmath5 . using this result it is straightforward to obtain @xmath148 we compared this result to numerical simulations in fig .
[ fig : exp ] and found a very good agreement .
apparently @xmath116 assumes a value independent of @xmath3 already for very small @xmath3 .
numerical simulations of @xmath25 for an exponential distributions with unit mean . for each @xmath5 and @xmath3
we simulated @xmath118 time series .
the figure shows results for @xmath119 and @xmath120 along with our @xmath3-independent analytical prediction ( [ lexp]).,scaledwidth=90.0% ]    another important member of the gumbel class is the normal distribution with mean zero and standard deviation @xmath62 .
here , the computation is usually a bit more complicated ; however , as in the previous examples , most of the work was already done in @xcite .
there we found that for @xmath130 and @xmath149 @xmath150 using ( [ taylor_mid ] ) this result allows us to estimate the effect of the drift on the correlations .
we find @xmath151 unlike the case of the exponential distribution , in the gaussian case the correlations are negative and depend strongly on @xmath3 .
the behavior is similar to that found in the weibull case , and matches the result ( [ weibull_gen ] ) ( up to logarithimic factors ) for @xmath152 .
numerical results for the gaussian distribution are shown in fig .
[ fig : gaussian ] .
similar to the uniform distribution , @xmath116 shows a steep descent at @xmath121 which becomes steeper for increasing @xmath3 .
however , the descent appears to be smoother than in the uniform case . in agreement with the considerations of section [ sec : largen ] , for @xmath78 and large
@xmath3 the correlations become independent of @xmath3 . in contrast , for @xmath91 the correlations increase monotonically with @xmath3 and do not appear to saturate . in the case of the gaussian distribution
our analytical prediction ( [ prediction_gaussian ] ) is only valid in a small regime for small @xmath78 and relatively small @xmath3 , but for @xmath91 the approximation is significantly better .     numerical simulations of @xmath25 for a gaussian distribution of width @xmath153 .
* left figure * shows @xmath25 for @xmath119 and @xmath120 . for each @xmath5 and @xmath3
we simulated @xmath127 time series .
the figure illustrates the steep descent of the correlator at @xmath121 when @xmath5 is small .
* right figure * shows @xmath25 for different , fixed values of @xmath5 together with our analytical results . here , we analyzed @xmath154 series of rv s for each drift rate .
again , we find agreement between the simulations and our analytical computations for small @xmath3 and @xmath5 .
we manually fitted curves @xmath155 .
interestingly for @xmath91 the agreement is a lot better and @xmath25 does not saturate at a constant value.,title="fig:",scaledwidth=50.0% ] numerical simulations of @xmath25 for a gaussian distribution of width @xmath153 .
* left figure * shows @xmath25 for @xmath119 and @xmath120 . for each @xmath5 and @xmath3
we simulated @xmath127 time series .
the figure illustrates the steep descent of the correlator at @xmath121 when @xmath5 is small . *
right figure * shows @xmath25 for different , fixed values of @xmath5 together with our analytical results . here , we analyzed @xmath154 series of rv s for each drift rate .
again , we find agreement between the simulations and our analytical computations for small @xmath3 and @xmath5 .
we manually fitted curves @xmath155 .
interestingly for @xmath91 the agreement is a lot better and @xmath25 does not saturate at a constant value.,title="fig:",scaledwidth=50.0% ]    to understand the marked difference between the exponential and gaussian distributions , we analyze the general class of gumbel - type distributions of the form @xmath156 with normalization @xmath157 and @xmath158 . in this case
it was found in @xcite that @xmath159 with some further computations this leads us to the following behavior of the leading order correction coefficient @xmath160 : @xmath161 where @xmath162 and @xmath163 are positive constants not depending on @xmath3 .
the exact values of @xmath162 and @xmath163 are very difficult to compute and we will not consider them in this article .
nevertheless , this expression nicely reproduces our results for the exponential and the gaussian distribution and shows how both positive and negative correlations may emerge in the gumbel class . for @xmath164 the leading first term vanishes and @xmath72 reduces to a positive constant . for all @xmath165
the first term dominates and the correlations grow with @xmath3 , with a sign determined by @xmath166 . for values @xmath167 of distributions that decay faster than the standard exponential we find negative correlations between the records , in agreement with the gaussian example , while for stretched exponential distributions ( @xmath168 ) positive , @xmath3-dependent correlations result ( see fig .  [
fig : comp ] for an example ) . in hindsight , the special role of the exponential distribution should not come as a surprise , since the gumbel distribution ( [ gumbeldist ] ) ( for which correlations are completely absent ) has an exponential tail . in the sense of extreme value theory , the exponential distribution is close to the gumbel , and the corresponding records in the ldm are therefore almost uncorrelated ( up to small residual correlations which are independent of @xmath3 ) .      as a subset of the distributions in the frchet class of extreme value statistics we consider the well known pareto distribution @xmath169 with @xmath170 and @xmath171 .
in @xcite it was found that in this case @xmath172 this allows us to compute @xmath116 for the pareto distribution .
we find @xmath173     * left : * numerical simulations of @xmath25 for pareto distributions . for each @xmath5 and @xmath174
we simulated @xmath118 time series of length @xmath123 .
the figure shows simulations for @xmath175 and @xmath176 , along with our analytical predictions for small @xmath5 . at @xmath121
we find a steep ascent of @xmath25 which gets steeper with increasing @xmath174 .
our analytical result from ( [ prediction_pareto ] ) is in good agreement with the simulations for small enough @xmath5 . *
right : * simulation results for a pareto distribution with @xmath177 and different values of @xmath5 .
we analyzed @xmath154 series of @xmath178 rv s for each @xmath5 .
our analytical work predicts a square - root behavior of @xmath179 , which is in good agreement with the simulations .
, title="fig:",scaledwidth=50.0% ] * left : * numerical simulations of @xmath25 for pareto distributions . for each @xmath5 and @xmath174
we simulated @xmath118 time series of length @xmath123 .
the figure shows simulations for @xmath175 and @xmath176 , along with our analytical predictions for small @xmath5 . at @xmath121
we find a steep ascent of @xmath25 which gets steeper with increasing @xmath174 .
our analytical result from ( [ prediction_pareto ] ) is in good agreement with the simulations for small enough @xmath5 . *
right : * simulation results for a pareto distribution with @xmath177 and different values of @xmath5 .
we analyzed @xmath154 series of @xmath178 rv s for each @xmath5 .
our analytical work predicts a square - root behavior of @xmath179 , which is in good agreement with the simulations .
, title="fig:",scaledwidth=50.0% ]    which leads to the expression @xmath180 for large @xmath3 . for positive @xmath5
, we thus expect positive correlations between neighbouring record events for all distributions of pareto form .
these correlations are increasing with @xmath3 slower than linearly .
we expect that , asymptotically , this behavior is universal for all distributions within the frchet class which show negative correlations between record events .
we have checked that this behavior is not asymptotic , and that the correlations become positive for larger values of @xmath3 . ] .
we compare eq.([prediction_pareto ] ) to numerical simulations in fig .
[ fig : pareto ] , finding good agreement for small @xmath5 , both for the @xmath5- and the @xmath3-dependence .
the agreement improves when the distribution becomes broader for smaller @xmath174 .
the strength of correlations for strongly heavy - tailed distributions of the pareto class is remarkable .
for example , for a pareto distribution with coefficient @xmath177 we found correlations @xmath181 for large @xmath3 ( see fig .
[ fig : comp ] ) . the correlations displayed in the right panel of fig.[fig : pareto ] increase ( for @xmath182 ) or decrease ( for @xmath183 ) with @xmath3 without showing any sign of saturation ,
although we know from the general considerations of section [ sec : largen ] that @xmath124 must approach an @xmath3-independent limit for @xmath182 and @xmath184 .
as we will see in the next subsection , this reflects the fact that the time scale @xmath185 at which the limit is attained is ( for a given value of @xmath5 ) particularly large for distributions of the frchet class .
we also note that because @xmath186 for @xmath183 , the normalized correlations are confined to lie between 0 and 1 in this case .
this implies that the divergence of @xmath56 for @xmath90 and @xmath183 , which was demonstrated in section [ sec : largen ] to occur for a specific example in the weibull class , can not happen here .
we therefore conjecture that , for distributions in the frchet class , the limit ( [ last ] ) exists also for @xmath183 .
the results of the preceding subsections can be summarized in a simple scaling picture .
we first recall from @xcite and @xcite that the ldm ( with @xmath78 ) contains a characteristic time scale @xmath185 at which the record rate @xmath88 crosses over from the i.i.d .
behavior @xmath187 to the limiting value @xmath188 . for the different distributions discussed above ,
this time scale diverges for @xmath98 according to @xmath189 ignoring for the moment the logarithmic factor in the gumbel case , these behaviors can be further simplified by expressing the different universality classes in terms of the generalized pareto distribution @xcite @xmath190 for @xmath191 this is of pareto type with @xmath192 , for @xmath193 it reduces to a weibull - type distribution similar to ( [ kuramaswamy ] ) with @xmath194 , and the gumbel class is represented by the exponential distribution which arises from ( [ genpareto ] ) for @xmath195 . using this representation , the different cases in ( 45 ) , ( 46 ) and ( 47 ) reduce to @xmath196
note that @xmath197 in the weibull class but @xmath198 in the frchet class , which explains the slow convergence to the @xmath90 limit in the latter case .
moreover , it was shown in @xcite that the integral ( [ i(n ) ] ) behaves for large @xmath3 as ( see also @xcite ) @xmath199 to relate this to the behavior of the correlations , we note that for large @xmath3 eq.([taylor_mid ] ) can be approximately written as @xmath200 inserting ( [ i(n)scal ] ) we thus conclude that , to leading order , @xmath201 which correctly reproduces both the sign and the order of magnitude of the correlations derived in section [ examples ] for the weibull and frchet classes : correlations are positive ( negative ) for @xmath191 ( @xmath193 ) , and they scale sublinearly ( superlinearly ) with @xmath3 in the two cases .
a similar calculation using the refined estimate ( [ igumbel ] ) of @xmath202 for the gumbel class shows that the correlations are negative ( positive ) for @xmath203 ( @xmath204 ) , in agreement with ( [ jgumbel ] ) .
numerical simulations of @xmath25 with @xmath205 for six different distributions . in each case
we averaged over @xmath118 time series of length @xmath206 .
the results for @xmath25 were binned logarithmically in order to improve the averaging for large @xmath3 . for the uniform and
the gaussian distribution @xmath25 is clearly negative for @xmath78 and it decreases with growing @xmath3 . at some @xmath207 the correlations become independent of @xmath3 .
for the special case of the exponential distribution we find a constant value of @xmath25 which is only slightly larger than unity .
both the stretched exponential ( with @xmath208 ) and the pareto distributions ( with @xmath209 ) show positive correlations .
note the slow convergence for the pareto distribution with @xmath210.,scaledwidth=90.0% ]
in this paper we analyzed the effect of a linear drift on the correlations between records drawn from series of independent rv s , as quantified by the normalized joint probability @xmath116 . in section 2
we derived general expressions ( exact and approximate ) for @xmath116 and recalled the fact that records are independent both for i.i.d .
rv s ( @xmath121 ) and in the special case where the rv s are drawn from the gumbel distribution .
our main analytic results were obtained in section 3 by way of an expansion in the small @xmath5 limit , similar to the approach developed previously in @xcite . using this approach
we were able to show that the correlations are generally negative ( ` repulsive ' ) for distributions in the weibull class , and positive ( attractive ) for distributions in the frchet class . in the gumbel class the sign of the correlations depends on the stretching exponent @xmath211 in ( [ fbeta ] ) , with the border between positive and negative correlations being given by the exponential case @xmath212 .
in contrast to all other cases , for distributions with an exponential tail the correlations are weak and independent of @xmath3 , which is consistent with the fact that this class of distributions also contains the gumbel distribution , which has no correlations at all .
simulation results illustrating the different cases are summarized in fig .
[ fig : comp ] . a special role of the exponential distribution in separating two regimes of qualitatively different behaviors has been noted previously in the related context of near - extreme events @xcite .
perhaps the most surprising and counterintuitive outcome of our work is the discovery of strong positive record correlations for distributions with a power law tail . in view of the substantial interest in detecting and explaining heavy - tailed distributions in all areas of science @xcite
, our finding suggests that drift - induced record correlations could be used as a distribution - free test for detecting power laws or streched exponentials in empirical data @xcite .
an interesting open question concerns the structure of the record correlations in the asymptotic limit @xmath90 , where the record rate approaches a nonzero constant and the record process thus becomes stationary .
based on the work of ballerini and resnick @xcite , we have argued in section [ sec : largen ] that the limit ( [ last ] ) exists for @xmath78 whenever the underlying distribution has a finite mean , whereas for @xmath183 it is possible that @xmath56 diverges for @xmath90 .
the fact that the coefficient of the leading order term in the small @xmath5 expansion generally diverges with @xmath3 indicates that the limiting function @xmath135 may be singular for @xmath98 , and we have presented numerical evidence for the occurrence of a discontinuity at @xmath121 for the weibull class . for heavy - tailed distributions in the frchet class ,
the large-@xmath3 asymptotics is difficult to ascertain numerically because of the slow convergence , but we have argued that in this case @xmath135 may exist also for @xmath91 .
rigorous work addressing these questions along the lines of @xcite would be most welcome .    this work has been supported by dfg within the _ bonn cologne graduate school of physics and astronomy _
, by friedrich ebert stiftung through a fellowship to gw , and by studienstiftung des deutschen volkes through a fellowship to jf .
spmpsci gembris , d. , taylor , j.g .
, suter d. : sports statistics : trends and random fluctuations in athletics .
nature * 417 * , 506 ( 2002 ) gembris , d. , taylor , j.g . , suter d. : evolution of athletic records : statistical effects versus real improvements . j. appl
. stat . * 34 * , 529 ( 2007 ) bassett , g.w .
: breaking recent global temperature records .
climatic change * 21 * , 303 ( 1992 ) benestad , r.e . : how often can we expect a record event ?
climate research * 25 * , 3 ( 2003 ) redner , s. , petersen , m. r. : role of global warming on the statistics of record - breaking temperatures > phys .
e * 74 * , 061114 ( 2006 ) meehl , g.a .
, tebaldi , c. , walton , g. , easterling , d. , mcdaniel , l. : relative increase of record high maximum temperatures compared to record low minimum temperatures in the u.s .
lett . * 36 * , l23701 ( 2009 ) wergen , g. , krug , j. : record - breaking temperatures reveal a warming climate .
epl * 92 * , 30008 ( 2010 ) newman , w.i . ,
malamud , b.d . ,
turcotte , d.l . : statistical properties of record - breaking temperatures .
e * 82 * , 066111 ( 2010 ) anderson , a. , kostinski , a. : reversible record breaking and variability : temperature distributions across the globe . j. appl .
49 * , 1681 ( 2010 ) de haan , l. , ferreira , a. : extreme value theory - an introduction .
springer , new york ( 2006 ) sibani , p. , brandt , m. , alstrm , p. : evolution and extinction dynamics in rugged fitness landscapes .
* 12 * , 361 ( 1998 ) krug , j. , jain , k. : breaking records in the evolutionary race .
physica a * 358 * , 1 ( 2005 ) park , s .- c .
, krug , j. : evolution in random fitness landscapes : the infinite sites model .
p04014 ( 2008 ) oliveira , l.p . ,
jensen , h.j .
, nicodemi , m. , sibani , p. : record dynamics and the observed temperature plateau in the magnetic creep - rate of type - ii superconductors .
b * 71 * , 104526 ( 2005 ) sibani , p. , rodriguez , g.f . , kenning , g.g . : intermittent quakes and record dynamics in the thermoremanent magnetization of a spin - glass .
phys . rev .
b * 74 * , 224407 ( 2006 ) wergen , g. , bogner , m. , krug , j. : record statistics for biased random walks , with an application to financial data . phys
e * 83 * , 051109 ( 2011 ) franke , j. , wergen , g. , krug , j. : records and sequences of records from random variables with a linear trend .
: theor . exp .
p10013 ( 2010 ) arnold , b.c . ,
balakrishnan , n. , nagaraja , h.n . : records .
wiley , new york ( 1998 ) nevzorov , v. b. : records : mathematical theory .
american mathematical society ( 2001 ) glick , n. : breaking records and breaking boards .
monthly * 85 * , 2 ( 1978 ) majumdar , s.n . ,
ziff , r.m . : universal record statistics of random walks and lvy flights .
* 101 * , 050601 ( 2008 ) , p. , wiese , k.j . : driven particle in a random landscape : disorder correlator , avalanche distribution , and extreme value statistics of records .
e * 79 * , 051105 ( 2009 ) sabhapandit , s. : record statistics of continuous time random walk .
epl * 94 * , 20003 ( 2011 ) ballerini , r. , resnick , s. : records from improving populations .
* 22 * , 487 - 502 ( 1985 ) ballerini , r. , resnick , s. : records in the presence of a linear trend .
* 19 * , 801 - 828 ( 1987 ) krug , j. : records in a changing world .
p07001 ( 2007 ) gumbel , e. j. : statistical theory of extreme values and some practical applications . national bureau of standards applied mathematics series * 33 * ( 1954 ) galambos , j. : the asymptotic theory of extreme order statistics .
wiley , new york ( 1987 ) sornette , d. : critical phenomena in natural sciences .
springer , berlin ( 2004 ) sabhapandit , s. , majumdar , s.n . : density of near - extreme events .
* 98 * , 140201 ( 2007 ) pickands iii , j. : statistical inference using extreme order statistics .
. stat . * 3
* , 119 ( 1975 ) christensen , k. , moloney , n.r .
: complexity and criticality . imperial college press , london ( 2005 ) franke , j. , wergen , g. , krug , j. ( in preparation )


while the inherent complexity of the solar corona and its dynamic processes has engendered great scientific interest and study over the past half century , it is this complexity that is perhaps the greatest challenge when attempting to model and understand the underlying processes of observed phenomena in a quantitative fashion @xcite . for this reason , the development and use of fully 3d magnetohydrodynamic ( mhd ) models describing the
corona and solar wind , which has undergone a veritable renaissance in the computing age , has been critical to furthering our understanding of the solar wind and eruptive events @xcite .
a particularly exciting avenue has been the development of 3d observation - based techniques , which provide the opportunity to include realistic conditions of the corona into studies of the global solar wind @xcite and to address how this governs dynamic events @xcite .
+ with its origin in the long - studied fundamental question of how the solar corona can be heated to and maintain temperatures in the million kelvin ( mk ) regime , a primary focus of this study is to address the issue of energy input into the corona in the context of a fully 3d global model . because non - mhd processes of energy transport and local conditions in the low corona are responsible for coronal heating mechanisms , it is often quite difficult and computationally expensive to include the small - scale micro physics of reconnection and turbulence responsible for coronal heating in 3d models , especially since the exact processes involved are not universally agreed upon . as a result
, heating models are often parameterized as a heating term that depends on various local magnetic and thermodynamic properties that is included in the energy equation ( e.g. @xcite ) . however , due to the large number of heating models and their ad - hoc formulations , it is crucial to use real observations to compare and constrain these models and deduce the best combination of models , both for active regions and the quiet sun , that yields the most realistic global simulations .
when using such models to study transient events in the solar corona , particularly coronal mass ejections ( cmes ) and coronal waves seen in the extreme ultraviolet regime ( euv waves , @xcite ) , it is a motivating principle that only by constraining physical theories and scenarios through as many observable manifestations available can we provide the best possible avenue to further our knowledge and precision in understanding these dynamic events ( e.g. @xcite ) . the ultimate goal being to strip away layers of empirical approximation with more fundamental physical terms and mechanisms .
however , it is equally important to validate and asses the ability of these improvements to adequately represent the global conditions in the solar corona .
+ to address this issue through a deterministic , data - driven approach , we modify the global corona model of the space weather modeling framework ( swmf ) @xcite to include the transition region between the chromosphere and the corona , where non - mhd thermodynamic terms of energy transport , such as electron heat conduction , radiative losses , and coronal heating all become important . using techniques similar to the pioneering work of @xcite , in this refined version of the model
the inner boundary is placed at the chromosphere or upper transition region rather than in the low corona and the now relevant non - mhd terms are added to the governing mhd energy equation .
+ in section [ section : model ] we discuss the new lower corona ( lc ) model and physical considerations added as part of this work . in section [ section : model_runs ] we overview our runs and results and discuss the application of this model to the study of the corona ( in particular the comparison to euv imaging observations ) . in section [ section : detailed_analysis ] we provide a detailed analysis of a particular run , and benchmark it via comparison to a previous corona model .
we conclude in section [ section : conclusion ] .
the main model and starting point for this work is the swmf solar corona ( sc ) model @xcite . as far as numerical methods ,
the code is fully parallelized and designed to be highly customizable in terms the methods , solvers , and equations used .
a complete description of the mhd equations and their implementation in the swmf can be found in @xcite as well as the addition of source terms in @xcite .
for the current version of the solar corona model , a steady state solar wind solution is obtained using a variable @xmath0 model @xcite .
the latter references use the wang - sheeley - arge ( wsa ) model @xcite as an initial condition to derive the variation in the polytropic index .
this then achieves a steady state mhd solution after nominal integration ( henceforth refered to as the `` standard sc model '' ) .
+ [ section : model_swmf ]    physically , the most important advantage of this tool and others like it lies in its ability to simulate the complete 3d environment of any event and carrington rotation ( standardized solar rotation number , abbreviated cr ) for which data for the photospheric magnetic field is available throughout the entire solar surface .
the initial magnetic configuration is extrapolated using the potential field source surface method ( pfssm ) @xcite , which uses magnetic coefficients derived from observations of the synoptic photospheric magnetic field ( typically high order maps using data from the mdi instrument aboard the soho observatory , or low order maps the wilcox solar observatory ) .
using this method , the subsequent evolution of the magnetic field towards a steady state during the simulation is no longer strictly potential .
+ while the standard version of the sc model has been very successful at reproducing the bi - modal solar wind structure and observations at 1au @xcite , which was the fundamental goal of that empirical model , the lower boundary at the solar surface is quite smooth and nearly uniform . in order to reproduce and study the fine structures of the the low corona we significantly modify the sc model to address the unique physics that take place at this boundary ( henceforth refered to as the the lower corona ( lc ) model ) .
this modification of the energy equation takes the form of : @xmath1 where , the subscript @xmath2 refers to the standard mhd terms in the swmf and the additional terms @xmath3 , @xmath4 , @xmath5 are described below . because these terms account for realistic energy transport , the polytropic index is no longer variable , and is set to a uniform value of @xmath6 .
the standard form of anisotropic electron heat conduction in the collisional limit is included as an additional energy flux term , @xmath7 , in the energy equation ( equivalent physically to adding a source term of @xmath8 .
the parallel component of spitzer conductivity @xcite , this term is especially important in the transition region and low corona and plays a critical role in determining the equilibrium density and temperatures at the base of the corona .
we use a value of @xmath9 in cgs units .
+ we also include a method for broadening the equilibrium scale length of the transition region , which is needed to resolve the entire transition region when using a chromospheric boundary ( described in detail by @xcite ) .
this reduces to both fixing @xmath10 to a constant value of @xmath11 and reducing the radiative cooling term by the inverse factor @xmath12 ) for temperatures below a transition region value , @xmath13 . in this work we use @xmath14
. + [ section : heat_conduction ]      the next term , @xmath15 , accounts for the radiative losses of the hot coronal plasma in the optically thin limit . the loss function @xmath16 is calculated using the chianti version 5 radiative loss routines @xcite and linearly interpolated as a tabulated function in the model .
the assumption of a fully ionized hydrogen gas is already implicit in this model , which gives @xmath17 .
it is important to note that the choice of abundances used in calculating the radiative losses function can change the total losses significantly in the coronal temperature regime .
for this reason we choose to use the coronal abundances file : ` sun_coronal_ext.abund ` @xcite when calculating the radiative cooling curve with chianti , and do not vary this function in this work .
the first empirical coronal heating term , @xmath19 , that we examine is equivalent in form to that used by @xcite and was originally based on work relating the surface unsigned magnetic flux of a convective envelope to the x - ray luminosity of stars @xcite .
this model is expressed as : @xmath20    [ heating_model_1 ] where @xmath21 , @xmath22 are fixed parameters , @xmath23 represents the total unsigned magnetic flux at the solar surface , @xmath24 is the local heating weighting function , and @xmath25 is a normalization constant .
this formulation calculates the total amount of energy input into the corona and distributes it via @xmath24 , which is chosen to be a function of @xmath18 .
the primary benefit of this model is its simplicity , the amount of heating at any location is given by the value of @xmath24 , and the total amount of heat input to the corona ( nearly constant ) is normalized by the global integral in the denominator , @xmath26 , which is carried out over the entire domain . in this implementation , @xmath23 is calculated by integrating the unsigned flux at the model boundary , @xmath27 , which is essentially the photospheric magnetogram used to initiate the model .
+ because both the relationship between x - ray luminosity and total power input into the corona is not well constrained , and varying magnetogram resolution does not constrain the true unsigned flux , @xmath25 represents a relatively free parameter with which to adjust the total power on a case by case basis . in our model runs we choose @xmath28 and linear weighting with the magnetic field magnitude , @xmath29 .
this gives the heating term the operational equivalent of @xmath30 , where @xmath31 is typically around @xmath32 in cgs units , and total integrated power of @xmath33 if applied uniformly .
+ one drawback of this term ( and of other @xmath18 weighted heating terms in global models ) is that the heating scale height , defined as the point along a given loop where @xmath34 , is not a free parameter and varies with the strength and spatial distribution of @xmath35 . also , note that the value of @xmath31 depends on the specific global magnetic field configuration and the radius to which the volume integral is calculated . in this work ,
we limit the range of influence of this function to the low corona by multiplying @xmath24 by an exponential envelope function with a scale height of @xmath36 .
the second heating model is a simple exponential scale height model .
using a standard exponential form we define heating model 2 as :    @xmath37}.\ ] ]    [ heating_model_2 ] where @xmath38 is the local heating rate at @xmath39 and @xmath40 is the heating scale height .
a short heating scale height is consistent with soho and trace euv observations of coronal loops @xcite , and motivated by the `` nanoflare '' class of coronal heating models , in which most of the heat is deposited in the upper chromosphere and transition region @xcite .
it is important to note that it is the high thermal conductivity of the coronal plasma that subsequently distributes this energy along the field / flow lines . for heating quiet sun regions
, we use the values @xmath41 , and @xmath42 , which give a total power of @xmath43 if applied uniformly ( surface flux @xmath44 ) .
the lack of explicit dependence on @xmath18 with an exponential scale height model allows us examine directly the effect that purely the magnetic geometry of the 3d corona has on the resultant thermodynamic equilibrium .      to maintain realistic temperatures in the solar wind and coronal holes at large distances
, we use the same form of eq . with values of @xmath45 , and @xmath46 giving additional total power @xmath47 .
this function is applied uniformly in each model run . while not entirely realistic for the acceleration profiles observed in coronal hole regions and the solar wind , we find that including this term is adequate for reproducing coronal hole emission in the absence of short scale height heating .
[ heating_model_ch ]      one interesting focus of this study is to isolate the complementary effects that non - uniform heating at the surface and the 3d magnetic topology have in determining the equilibrium structure of the corona .
for this reason , as one of the options in the model we include the ability to apply the low scale height heating model ( either 1 or 2 ) to only regions associated with closed field in the corona .
to do so , we use the pfssm extrapolation computed to initiate the model to determine at every computational element whether the field threading the location is open ( connects to the source surface and becomes radial ) or closed ( turns back down to the solar surface ) .
if the region is associated with open field , then only the long scale height `` coronal hole '' heating model is applied ( section [ heating_model_ch ] ) , while if the region is closed , both are applied .
+ we would also like to note that any of the empirical heating terms discussed above ( sections [ heating_model_1 ] - [ heating_model_oc ] ) may also be attributed to wave absorption and dissapation , probably due to the ion interaction with alfvn wave turbulence . in this case
, the heating may be specified in terms of the wave ( spectral ) energy and their damping coefficient , and the wave turbulence may contribute not only to the energy source but also to the momentum fluxes ( see e.g. @xcite and papers cited therein ) .
however , this particular study is focused primarily on the plasma conditions low corona , and we operate under the assumption that alfvn wave pressure will be much lower than the thermal pressure there .
[ heating_model_oc ]      boundary conditions at the @xmath39 surface are applied to fix constant values temperature and density at the face boundary according to either of the models described below .
each of these use a zero flow condition and fix the tangential component of the magnetic field ( similar to those used in @xcite . for the supersonic flow at the outer boundary at
@xmath48 we use a floating ( zero - gradient ) condition .
the primary boundary condition used in this model is a simple uniform `` chromospheric '' condition .
as described in @xcite , by setting chromospheric values of electron temperature , @xmath49 , and density , @xmath50 , and including the transition region broadening method described in section [ section : heat_conduction ] , it becomes possible to resolve the entire transition region everywhere in a global model .
this is a critical advantage because the resultant topological equilibrium in the coronal part of the solution ( density enhancements , varying width of the transition region , etc . )
will depend entirely on the included physics and magnetic geometry , and not on its proximity to the boundary of the model .
this fact is particularly important for the study of any sort of dynamics in the global corona , as the sharp transition from fast alfvn speed in the quiet corona ( @xmath51 ) to extremely slow speeds in the chromosphere ( @xmath52 ) , provides a realistic dispersion buffer for fast magnetosonic waves .      in order to include a realistic lower boundary condition for the corona ,
while at the same time maintaining decent computational efficiency for the global domain , we use the method of the radiative energy balance model ( reb ) outlined by @xcite ( further formalism , including coronal heating in @xcite ) . in short ,
this model fixes the base electron temperature to a high transition region value , @xmath53 , and in turn assumes an equilibrium balance throughout the rest of the transition region below to the top of the chromosphere , from which an electron density , @xmath54 , may be obtained . by rearranging the terms in the mhd energy equation and integrating over t at constant pressure ,
one obtains : @xmath55 where , the integral of the radiative loss function is carried out between @xmath56 and @xmath53 .
this allows for the base density @xmath54 to be a function of the physical conditions included in the model ( coronal heating , anisotropic heat conduction , and radiative cooling ) and thus varies spatially along the boundary .
most importantly , the numerator of shows how the interplay between magnetic field strength , @xmath57 ( via @xmath19 in the right term ) , and normal orientation at the surface , @xmath58 ( via @xmath59 on the left ) , sets the base density on the solar surface . the ability to capture these features with such a simple , easily computed boundary is a primary motivation in using this particular model .
another aspect of this work has been to address the unique geometric concerns that arise when including the transition region and below .
a typical solar corona ( sc ) simulation is carried out with an adaptive cartesian grid over a sun - centered @xmath60r@xmath61 cube with @xmath62 cell blocks .
the average cell size is smallest at the surface ( @xmath63 r@xmath64 km ) and is incrementally increased with distance to a largest size of @xmath65 r@xmath61 near the outer boundary , giving a few million cells .
while this is a relatively large surface cell size , it is adequate considering the precision of the magnetic field boundary ( an observational limitation ) .
however , if one is required to resolve the transition region , with dynamic radial spatial scales on the order of a few km , it quickly becomes infeasible to use a cartesian grid ( with every @xmath66 refinement in length the number of surface cells increases by 4 ) .
this can be partially alleviated using methods to widen the extent of the model transition region to typical scales of order @xmath67 km without affecting the coronal solution ( section [ section : heat_conduction ] ) , but the problem is still quite substantial .
this then represents a balance between height accuracy and computational efficiency ( with a cartesian grid , the number of surface cells increases by a factor of 4 with every successive @xmath66 refinement ) .
+ to address this scale and resolution issue we adapted the generalized grid capability of the swmf , which has ability to calculate mhd fluxes in arbitrary geometries , to construct a spherical ( r@xmath68 ) grid function with highly non - uniform radial scales near the transition region ( @xmath69 km ) that smoothly transitions to near equal face area at large @xmath70 ( @xmath71 km at @xmath72 ) .
because the grid function is continuous , the model maintains the block - adaptive mesh and adaptive mesh refinement ( amr ) capability for flexible local spatial resolution both initially and as the simulation progresses .
for this study , focusing on the global structure of the corona , we use a uniform spacing in the angular directions with @xmath73 at the surface , which is coarsened by a factor of two via amr beyond @xmath74 for regions within @xmath75 and @xmath76 outside ( to avoid needlessly small cell sizes at the poles ) .
the outer boundary is fixed at @xmath48 .
degenerate cells touching the polar axes are treated with azimuthal averaging to avoid discontinuity across the poles . a grid comparison near the surface between this work and an amr cartesian mesh with a similar number of total cells
is shown in figure [ fig : geometry ] .      to study the euv corona and associated dynamics in the context of real events , we include a proper treatment of integrated coronal emission in the euv range , particularly for the purpose of comparing simulated observations to existing euv observations by a current solar observatory .
this process , following the work of @xcite , involves two main steps : ( i ) characterization of detector response based on physical parameters that are calculated in the model and ( ii ) line - of - sight ( los ) integration though the 3d data set to create synthesized image
. + [ section : los - synthesis ]    the first step is to understand and calculate the instrumental response for coronal material at a given temperature and density . for simplicity and ease to compare to studied observations we chose to model the three coronal band - passes at 171 , 195 , and 284  of the eit instrument on the soho spacecraft @xcite and the almg filter of the sxt instrument @xcite . because of
the extremely high temperature and low densities , the emission line euv and soft x - ray radiation of the solar corona above  0.5 mk can be treated as optically thin to a good approximation . as such , the intensity measured by each pixel of an imaging detector can be treated as a los integral through the coronal plasma , namely : @xmath77 where @xmath78 is the pixel response measured in data numbers ( dn ) per second .
here @xmath79 follows along the line - of - sight , the square of the electron density , @xmath80 , accounts for the amount of emitting material and @xmath81 is the instrumental response of filter @xmath82 per unit emission measure as a function of both temperature and density .
this is equivalent to integrating the resolved differential emission measure ( dem ) distribution along @xmath79 . to calculate the response functions , @xmath81
, we use the chianti version 5 emission line analysis routines @xcite to generate synthetic spectra , and the eit and sxt analysis routines for detector response ( a part of the solarsoft ( ssw ) framework written in idl @xcite ) .
in this section we present the application of the lc model to cr 1913 , initialized using a synoptic mdi magnetogram centered on aug 27 , 1996 .
occuring during solar minimum , the sun at this time features a conspicuous equatorial coronal hole extending from the north to down past disk - center , as well a large active region .
this then conveniently allows us to examine the ability of a given model to describe these three basic regimes : the average quiet - sun , coronal holes , and active regions .
additionally , this rotation has been studied in detail using similar methods ( e.g. @xcite ) and thus provides a means of comparison to existing work .
+ four model runs used to demonstrate the usefulness of the lc model are summarized in table [ table : model_runs ] . as a convenient comparison of their global thermodynamic and topological properties , the los synthesis of euv and soft x - ray emission
are shown in figure [ fig : euvcomp_4runs ] .
each of the models presented use either @xmath18 weighted heating ( model 1 , section [ heating_model_1 ] ) or exponential heating ( model 2 , section [ heating_model_2 ] ) ) .
all four use coronal hole heating ( section [ heating_model_ch ] ) applied uniformly to the domain , while all but run c use the open / closed field weighting ( section [ heating_model_oc ] ) to achieve a more realistic equilibrium in these regions . in order to include strong heating in active regions ( ar ) with high magnetic field strength , run d also includes a modification to the heating function , and transitions smoothly from exponential heating to @xmath18 weighted heating ( @xmath83 ) above @xmath84 gauss .
+    cccccc a & reb & @xmath85 & yes & @xmath86 + b & chromo & @xmath85 & yes & @xmath86 + c & chromo & exp & no & @xmath87 + d & chromo & exp + ar & yes & @xmath86 +    all four model runs show similar global features in the euv and soft x - ray when compared to observations but differ on the precise details .
for example , all four runs resolve cool temperatures along the the polarity inversion lines near the surface ( seen to the east and west of the extended equatorial coronal hole near the center of the disk ) .
this is naturally due to field line orientation and loop length governing thermodynamic equilibrium .
this commonality serves to emphasize the strong role that the 3d magnetic field topology of a given carrington rotation plays in determining dynamic equilibrium in the corona . + in the comparison between run a ( reb model bc ) and run b ( chromospheric bc ) we observe more uniform emission near the surface as well as slightly higher temperature ( as seen in the longer extent of sxt emission ) for the reb case . because the boundary is near the top of the transition region in the reb model , the height of the transition region is less able to vary in response to magnitude of the downward heat conduction flux from above
also , since the reb model does not include a region with chromospheric densities ( where the radiative cooling is much higher due to the @xmath88 dependence ) a smaller fraction of the coronal heating power is lost in the transition region , leading to higher equilibrium temperatures .
however , it is clear that the reb model , which is more suited for computational efficiency due to the longer equilibrium length scales at @xmath89 , can be used to adequately describe the ambient global corona .
+ one of the most striking results is perhaps the comparison of run c to the other model runs . in this run
, the coronal heating function was applied uniformly in the computational domain ( no dependence on @xmath90 or @xmath23 ) .
while the coronal holes and significant ar emission are unsurprisingly not well preserved in the euv , still much of the basic features of the low corona are reproduced .
the ability of simple uniform heating model to reproduce obvious features on the surface , such as average quiet sun emission , and lessening near inversion lines , as well as the bi - modal structure between open and closed field regions at higher temperatures ( @xmath91 line and soft x - rays ) speaks to the importance of both the 3d magnetic topology in determining the dynamic equilibrium of the corona and the redistribution of energy via a thermodynamic energy equation including electron heat conduction . + for a more realistic empirical model in run d , we modify run c by including ar heating and the open - field cutoff to better describe heating in active regions in coronal holes .
the most interesting finding in this comparison is that by including significantly enhanced heating in the large active region , we observe a pronounced effect on the equilibrium structure of the ar associated closed field / streamer region ( eastern side of the corona in figure [ fig : euvcomp_4runs ] ) .
this produces noticeable feedback on the global structure as seen via a larger northward tilt of the heliospheric current sheet near the ar longitudes ( shown in detail figure [ fig : streamer_br_slice ] ) , and suggests a complex relationship both between the thermodynamics of the low corona and the global structure of the solar wind .
because this effect can not be fundamentally extracted from extrapolating the magnetic field alone , this provides not only a stong motivation for using thermodynamic models , but also another element with which to constrain and refine ar heating models in future studies ( the one used in in this work being a necessary but crude approximation ) .
+ another important inference about coronal heating can be gleaned from the comparison of the two heating models directly ( run b and run d ) .
although the total heating power is the same in each model , the @xmath18 weighted heating model naturally has a shorter scale height ( due to the falloff of @xmath18 with radius ) and we subsequently observe slightly enhanced pressure in the low corona and lower temperatures ( hence density scale height ) in the closed field regions .
however , the fact that the simple exponential heating model applied to the quiet sun in run d can reproduce the mean structure of the low corona with the same fidelity as a more complicated empirical model , demonstrates that perhaps complex empirical models are not as entirely necessary as some might suggest and motivates a need for more physics based models to advance further .
[ topological analysis ]    in this section we provide a more detailed analysis of run d in particular .
we examine the magnetic and thermodynamic structures realized by the lc model and provide a comparison to the standard sc model @xcite .
an important advantage of including heat conduction in a global 3d enrionment is the ability to the study complex open / closed field topologies in a self consistent manner , which allows one to study the feedback effect that the magnetic topology has on heat distribution and vice versa .
the most striking effect can be seen by examining the differences between heat conduction along open or closed field geometries . because closed field regions represent closed systems for the flow of thermal energy , they acheive higher temperatures than their open - field counterparts and shift into the soft - xray regime ( @xmath92mk ) .
this natural correspondence can be clearly seen in figure [ fig : temperature_topology ] where we show a 3d surface at fixed temperature ( @xmath93 mk ) and an los image of soft - x - ray emission .
it is also important to examine in detail the effect that the thermodynamic model has on the magnetic structure as a whole . in figure
[ fig : magnetic_topology ] we show a comparison between closed field lines at @xmath94 for the pfssm initial condition @xcite , the standard sc model @xcite and run d. immediately obvious are the drastic changes between the closed field / streamer structure of the polytropic sc model and the thermodynamic lc model . unsurprisingly , the high temperatures and thermal energy introduced by coronal heating and conduction along the closed field produces a large thermal stress in the relatively high @xmath95 regions near the streamer cusps along the current sheet .
this in turn leads to significant reorganization of the magnetic field as compared to a pfssm or polytropic mhd model .
the fact that this sort of stress is highly sensitive to both the local nature of a given coronal heating model ( figure [ fig : streamer_br_slice ] ) and the global magnetic topology ( figure [ fig : magnetic_topology ] ) is a leading motivation for using a fully 3d , thermodynamic model when studying the structure of the global corona .      in figure
[ fig : euvcomp_run10towsa ] we show los euv and soft - x - ray image synthesis comparisons between run d and the standard sc model . in figure
[ fig : eit_slice ] we show two quantitative slice comparisons cutting an average quiet sun region and the large active region .
the improved agreement of run d in multiple filter bands is indicative of the improved temperature resolution possible with the lc model .
while not an entirely fair comparison because the standard sc model does not include a full energy equation , we demonstrate the ability to quantitatively asses the relative agreement between a given model and observations in the euv regime .
this sort of comparison also allows one to study in detail various choices of heating models and , though multi - filter comparison , provide a basis to resolve the inherent degeneracies between modifying the total coronal heating power and the heating scale height within the model .
a further examination of some of the limitations imposed by euv synthesis from mhd models with finite resolution as well as future considerations can be found in appendix [ section : appendix ] .
it is also instructive to examine the global temperature and density structure as a function of height for both run d and the standard sc model . in figures [ fig : slices_temperature ] and [ fig : slices_density ]
we show electron temperature and number density at three separate heights in the low corona . important to emphasize is the fact that with a thermodynamic model , it is possible to resolve both significant changes in temperature with height and with respect to the magnetic field topology in the model .
as expected , near the top of the transition region run d resolves cooler temperatures near the inversion lines due the short loop length , while simultaneously higher in the corona , hotter temperatures ( @xmath96 ) are achieved in close field regions , a natural result of energy flowing on a closed path via heat conduction .
this greater sensitivity to field geometry ( something difficult for a polytropic model to do without extreme fine tuning of boundary conditions ) is critical when trying to correctly describe the conditions in the low corona .
in this paper we demonstrate the clear need for a self - consistent treatment of the thermodynamic energy equation and boundary conditions when studying the global properties of the low corona with an mhd model .
we observe that the interplay between coronal heating and electron heat conduction strongly governs the details of important strucures in the low corona , which is something that connot be described by extrapolation of the magnetic field on its own ( e.g. the significant stress of the streamer regions and current sheet in response to active region heating , figures [ fig : streamer_br_slice ] and [ fig : magnetic_topology ] ) .
we also showcase the ability of the lc model to examine the effectiveness and applicability of various heating functions , and find that a simple exponential heating profile does quite well in describing the average structure of the quiet sun .
+ with this foundation , we believe that the lc model can provide a more realistic vehicle with which to study a wide range of aspects of the dynamic corona and solar wind . for example , the time - dependent dynamics of any sort of transient event in the corona , such as a cme or euv wave , depend critically on the magnetic and density structure of the corona ( density scale height being highly sensitive to temperature ) and therefore it is important to use improved models when determining the global ambient medium in which they evolve .
additionally because of the general flexibility of the swmf , future coupling of the lc model to complementary solar wind models , such as the sc model based on a full description of alfvn wave turbulence currently in development ( * ? ? ?
* in press ) , or one including a multi - species treatement to resolve electron / ion temperature decoupling in the high corona , has become a feasible scenario .
ultimately , one thing is clear : with each passing month , year , and decade , our understanding of these fundamental process is constantly refined . in furthering this process
it is critical to continue to improve the existing models and validate new efforts as they emerge .
this work is primarily supported under the nasa earth and space science fellowship ( nessf ) program nasa nessf08-helio08f-0007 as well as grants : nsf atm-0639335 ( career ) and nasa nnx08aq16 g ( lws ) .
computations were performed on the computing cluster at the institute for astronomy ( ifa ) advanced technology research center ( atrc ) pukalani , hi . c. downs would like to thank the entire swmf team for their innovative research and willing collaboration .
+ soho eit images courtesy of the soho eit consortium ; soho is a project of international cooperation between esa and nasa .
yohkoh data courtesy of the nasa - supported yohkoh legacy archive at montana state university .
one critical caveat when imaging the solar corona is that remote sensing images are fundamentally ambiguous in the sense that they are 2d realizations of a fully 3d , optically thin data set . for our purposes , it is best to frame the problem with the concepts of modeling an unresolved differential emission measure ( dem ) distribution and hydrostatic weighting bias , well described in @xcite .
these issues arise due to two main facts .
( i ) the spatial resolution of current imaging instruments ( e.g. trace , soho , stereo ) all lie below the fundamental width of an isolated flow / field line . and
( ii ) any image is a line of sight integral over a range of heights , which are preferentially biased towards hotter flow lines with longer density scale heights .
thus at further distance off of the limb , the observation may be biased towards flow lines that may not be energetically important at the surface .
local dem distributions can also be directly observed through coupled filter and tomographic inversions @xcite .
+ though the los integral for image synthesis includes the range of temperatures and densities in cells along the line of sight , it can not fundamentally account for the dem distribution that is below the physical resolution of the model . in mhd
, a given cell has a unique value for its thermodynamic properties , particularly @xmath97 , and @xmath98 , while a physical picture of the low @xmath95 multi - hydrostatic corona below the resolution limit can be treated as a superposition of densities and temperatures , each confined to an isothermal flow / field line . in this sense ,
the true dem distribution of a parcel of gas that contributes to an euv image ( e.g. the eit instrument ) is instead represented as a delta function in the simulation .
one aspect of this limitation can be seen in difference between the data and synthesized limb profiles ( figure [ fig : eit_slice ] ) . because the observational weighting bias shifts to higher temperatures further with distance from the surface , the observed limb profile can not be represented in all four filters simultaneously with only the larger scale temperature equilibrium produced by the mhd model . + to further examine this issue in the context of the lc model we use filter ratios between the euv images , which are nominally sensitive to temperature differences . in figure
[ fig : demcurves ] we show filter ratio values for soho eit observations on aug 27 , 1996 and the path in the filter ratio plane spanned by all possible filter ratios of the optically thin response functions @xmath99 and @xmath100 ( section [ section : los - synthesis ] )
. unsurprisingly , the euv data spans a much larger range than that allowed by a single temperature . to simply examine the effect that an unresolved dem distribution may have on the synthesis and subsequent comparison to observational data we construct an ad hoc dem by gaussian convolution of the filter response functions with t : @xmath101 and use @xmath102 for this discussion .
while we by no means posit that this dem distribution should resemble the corona everywhere , it includes the basic temperature splitting features of those gaussian dem distributions calculated from observations ( e.g. @xcite .
we also show the effect of this modification in the right panel of figure [ fig : demcurves ] .
+ the effect that this simple induced dem distribution has on the synthesized ratios is shown in figure [ fig : dem_effect_obs ] .
it is clear comparing the limb tracks of the synthesized images to the temperature tracks in figure [ fig : demcurves ] that the los synthesis at the current model resolution does not produce a significant dem distribution on its own ( particularly at the limb ) .
however , the improvement gained in terms of the relative location with respect to observations for the modified dem synthesis suggests that this effect should not be overlooked when constructing more detailed models in the future .
high resolution studies designed to approach the fundamental scale of isothermal separation between coronal loops , particularly those involving local , time - dependent heating for example , should be able to address this fundamental issue .
41 natexlab#1#1    abbett , w.  p. 2007
, , 665 , 1469    altschuler , m.  d. , levine , r.  h. , stix , m. , & harvey , j. 1977 , , 51 , 345    arge , c. , luhmann , j. , odstrcil , d. , schrijver , c. , & li , y. 2004 , journal of atmospheric and solar - terrestrial physics , 66 , 1295    , m.  j. , alexander , d. , hurlburt , n. , newmark , j.  s. , neupert , w.  m. , klimchuk , j.  a. , & gary , g.  a. 2000 , , 531 , 1129    aschwanden , m.  j. & nitta , n. 2000 , , 535 , l59    aschwanden , m.  j. & acton , l.  w. 2001 , , 550 , 475    aschwanden , m.  j. & schrijver , c.  j. 2002 , , 142 , 269    aschwanden , m.  j. 2008 , in astronomical society of the pacific conference series , vol .
383 , subsurface and atmospheric influences on solar activity , ed .
r.  howe , r.  w. komm , k.  s. balasubramaniam , & g.  j.  d. petrie , 1    bercik , d.  j. , fisher , g.  h. , johns - krull , c.  m. , , & abbett , w.  p. 2005 , , 631 , 529    cohen , o. , sokolov , i.  v. , roussev , i.  i. , arge , c.  n. , manchester , w.  b. , gombosi , t.  i. , frazin , r.  a. , park , h. , et al .
2007 , , 654 , l163    cohen , o. , sokolov , i.  v. , roussev , i.  i. , & gombosi , t.  i. 2008 , , 113 , a03104    , o. , attrill , g. d.  r. , manchester , w.  b. , & wills - davey , m.  j. 2009 , , 705 , 587    delaboudinire , artzner , g.  e. , brunaud , j. , gabriel , a.  h. , hochedez , j.  f. , millier , f. , song , x.  y. , au , b. , et al .
1995 , , 162 , 291    frazin , r.  a. , jacob , m. , manchester , w.  b. , morgan , h. , & wakin , m.  b. 2009 , , 695 , 636    freeland , s.  l. & handy , b.  n. 1998 , , 182 , 497    groth , c. p.  t. , de  zeeuw , d.  l. , gombosi , t.  i. , & powell , k.  g. 2000 , , 105 , 25,053    , c. , roussev , i.  i. , lugaz , n. , & poedts , s. 2009 , , 695 , l171    , e. , feldman , u. , & dere , k.  p. 2002
, , 139 , 281    landi , e. , zanna , d.  g. , young , p.  r. , dere , k.  p. , mason , h.  e. , & landini , m. 2006 , , 162 , 261    lionello , r. , linker , j.  a. , , & mikic , z. 2001 , , 546 , 542    lionello , r. , linker , j.  a. , & mikic , z. 2009 , , 690 , 902    lugaz , n. , manchester , roussev , i.  i. , tth , g. , & gombosi , t.  i. 2007 , , 659 , 788    lugaz , n. , vourlidas , a. , roussev , i.  i. , jacobs , c. , manchester , & cohen , o. 2008 , , 684 , l111    lugaz , n. , vourlidas , a. , roussev , i. , & morgan , h. 2009 , , 259 , 269    , z. , linker , j.  a. , schnack , d.  d. , lionello , r. , & tarditi , a. 1999 , physics of plasmas , 6 , 2217    mok , y. , mikic , z. , lionello , r. , , & linker , j.  a. 2005 , , 621 , 1098    mok , y. , miki , z. , lionello , r. , & linker , j.  a. 2008 , , 679 , l161    moses , d. , clette , f. , delaboudinire , j.  p. , artzner , g.  e. , bougnet , m. , brunaud , j. , carabetian , c. , gabriel , a.  h. , et al .
1997 , , 175 , 571    oran , r. , sokolov , i.  v. , roussev , i.  i. , van der  holst , b. , manchester , w.  b. , & gombosi , t.  i. 2009 , in astronomical society of the pacific conference series , vol .
406 , numerical modeling of space plasma flows : astronum-2009 proceedings of the 2nd international conference , ed
. n.  v.  pogorelov , e.  audit , p.  colella , & g.  p.  zank    parker , e.  n. 1988 , , 330 , 474    powell , k. 1999 , j. comput .
phys . , 154 , 284    roussev , i.  i. , gombosi , t.  i. , sokolov , i.  v. , velli , m. , manchester , w. , dezeeuw , d.  l. , liewer , p. , tth , g. , & luhmann , j. 2003 , , 595 , l57    , i.  i. , sokolov , i.  v. , forbes , t.  g. , gombosi , t.  i. , lee , m.  a. , & sakai , j.  i. 2004 , , 605 , l73    roussev , i.  i. & sokolov , i.  v. 2006 , washington dc american geophysical union geophysical monograph series , 165 , 89    roussev , i.  i. , lugaz , n. , & sokolov , i.  v. 2007 , . , 668 , l87    schrijver , c.  j. , sandman , a.  w. , aschwanden , m.  j. , & derosa , m.  l. 2004 , , 615 , 512    sokolov , i.  v. , roussev , i.  i. , skender , m. , gombosi , t.  i. , & usmanov , a.  v. 2009 , , 696 , 261    , l. 1965 , physics of fully ionized gases , 2nd rev .
( new york : interscience publication )    tth , g. , sokolov , i.  v. , gombosi , t.  i. , chesney , d.  r. , clauer , r.  c. , de  zeeuw , d.  l. , hansen , k.  c. , kane , k.  j. , et al .
2005 , , 110 , a12226    tsuneta , s. , acton , l. , bruner , m. , lemen , j. , brown , w. , caravalho , r. , catura , r. , freeland , et al . 1991 , , 136 , 37    withbroe , g.  l. 1988 , , 325 , 442
right : the customized spherical grid with minimum @xmath103 ( part of this work ) .
all three directions do not have to be further refined at the surface to enhance the radial resolution when using spherical geometry , leading to significant savings in computational overhead .
in this example , both simulation domains have comparable numbers of total cells.,title="fig:",scaledwidth=35.0% ] .
right : the customized spherical grid with minimum @xmath103 ( part of this work ) .
all three directions do not have to be further refined at the surface to enhance the radial resolution when using spherical geometry , leading to significant savings in computational overhead . in this example
, both simulation domains have comparable numbers of total cells.,title="fig:",scaledwidth=35.0% ]        ) causes an extreme shift in the northward tilt of the heliospheric current sheet @xmath104 .
simultaneously , the region not thermally connected to the active region on the opposite side of the sun remains unchanged .
the color contours display the radial component of the magnetic field , @xmath105 .
the chosen slice is parallel to the polar axis and intersects the inversion line between the large active region seen on aug 27 , 1996 .
the 3d magnetic streamlines intersect the slice at the same location in each figure.,title="fig:",scaledwidth=42.0% ] ) causes an extreme shift in the northward tilt of the heliospheric current sheet @xmath104 .
simultaneously , the region not thermally connected to the active region on the opposite side of the sun remains unchanged .
the color contours display the radial component of the magnetic field , @xmath105 .
the chosen slice is parallel to the polar axis and intersects the inversion line between the large active region seen on aug 27 , 1996 .
the 3d magnetic streamlines intersect the slice at the same location in each figure.,title="fig:",scaledwidth=42.0% ] ) causes an extreme shift in the northward tilt of the heliospheric current sheet @xmath104 .
simultaneously , the region not thermally connected to the active region on the opposite side of the sun remains unchanged .
the color contours display the radial component of the magnetic field , @xmath105 .
the chosen slice is parallel to the polar axis and intersects the inversion line between the large active region seen on aug 27 , 1996 .
the 3d magnetic streamlines intersect the slice at the same location in each figure.,title="fig:",scaledwidth=8.0% ]     and @xmath106 with color contours of t. the black contour lines of normal radial field direction @xmath107 are shown to illustrate how the hot , closed field streamer regions overlay the inversion lines on the surface.,title="fig:",scaledwidth=80.0% ]   and @xmath106 with color contours of t. the black contour lines of normal radial field direction @xmath107 are shown to illustrate how the hot , closed field streamer regions overlay the inversion lines on the surface.,title="fig:",scaledwidth=13.0% ]     and @xmath108 in each model.,title="fig:",scaledwidth=30.0% ]   and @xmath108 in each model.,title="fig:",scaledwidth=30.0% ]   and @xmath108 in each model.,title="fig:",scaledwidth=30.0% ]   and @xmath108 in each model.,title="fig:",scaledwidth=30.0% ]   and @xmath108 in each model.,title="fig:",scaledwidth=30.0% ]   and @xmath108 in each model.,title="fig:",scaledwidth=30.0% ]   and @xmath108 in each model.,title="fig:",scaledwidth=7.0% ]                          ( top ) , @xmath109 ( middle ) , and @xmath110 ( bottom)@xmath111 for run d ( left ) and the default sc model ( right ) for cr 1913 ( black lines at @xmath107 ) .
the full thermodynamic energy equation is able to capture large changes in temperature over short changes in radius .
as expected , near the top of the transition region we resolve cooler temperatures near the inversion lines due the short loop length .
simultaneously higher in the corona , hotter temperatures ( @xmath96 ) are achieved in close field regions , a natural result of energy flowing on a closed path via heat conduction.,title="fig:",scaledwidth=44.0% ]   ( top ) , @xmath109 ( middle ) , and @xmath110 ( bottom)@xmath111 for run d ( left ) and the default sc model ( right ) for cr 1913 ( black lines at @xmath107 ) .
the full thermodynamic energy equation is able to capture large changes in temperature over short changes in radius .
as expected , near the top of the transition region we resolve cooler temperatures near the inversion lines due the short loop length .
simultaneously higher in the corona , hotter temperatures ( @xmath96 ) are achieved in close field regions , a natural result of energy flowing on a closed path via heat conduction.,title="fig:",scaledwidth=44.0% ]   ( top ) , @xmath109 ( middle ) , and @xmath110 ( bottom)@xmath111 for run d ( left ) and the default sc model ( right ) for cr 1913 ( black lines at @xmath107 ) .
the full thermodynamic energy equation is able to capture large changes in temperature over short changes in radius .
as expected , near the top of the transition region we resolve cooler temperatures near the inversion lines due the short loop length .
simultaneously higher in the corona , hotter temperatures ( @xmath96 ) are achieved in close field regions , a natural result of energy flowing on a closed path via heat conduction.,title="fig:",scaledwidth=44.0% ]   ( top ) , @xmath109 ( middle ) , and @xmath110 ( bottom)@xmath111 for run d ( left ) and the default sc model ( right ) for cr 1913 ( black lines at @xmath107 ) .
the full thermodynamic energy equation is able to capture large changes in temperature over short changes in radius .
as expected , near the top of the transition region we resolve cooler temperatures near the inversion lines due the short loop length .
simultaneously higher in the corona , hotter temperatures ( @xmath96 ) are achieved in close field regions , a natural result of energy flowing on a closed path via heat conduction.,title="fig:",scaledwidth=44.0% ]   ( top ) , @xmath109 ( middle ) , and @xmath110 ( bottom)@xmath111 for run d ( left ) and the default sc model ( right ) for cr 1913 ( black lines at @xmath107 ) .
the full thermodynamic energy equation is able to capture large changes in temperature over short changes in radius .
as expected , near the top of the transition region we resolve cooler temperatures near the inversion lines due the short loop length .
simultaneously higher in the corona , hotter temperatures ( @xmath96 ) are achieved in close field regions , a natural result of energy flowing on a closed path via heat conduction.,title="fig:",scaledwidth=44.0% ]   ( top ) , @xmath109 ( middle ) , and @xmath110 ( bottom)@xmath111 for run d ( left ) and the default sc model ( right ) for cr 1913 ( black lines at @xmath107 ) .
the full thermodynamic energy equation is able to capture large changes in temperature over short changes in radius .
as expected , near the top of the transition region we resolve cooler temperatures near the inversion lines due the short loop length .
simultaneously higher in the corona , hotter temperatures ( @xmath96 ) are achieved in close field regions , a natural result of energy flowing on a closed path via heat conduction.,title="fig:",scaledwidth=44.0% ]   ( top ) , @xmath109 ( middle ) , and @xmath110 ( bottom)@xmath111 for run d ( left ) and the default sc model ( right ) for cr 1913 ( black lines at @xmath107 ) .
the full thermodynamic energy equation is able to capture large changes in temperature over short changes in radius .
as expected , near the top of the transition region we resolve cooler temperatures near the inversion lines due the short loop length .
simultaneously higher in the corona , hotter temperatures ( @xmath96 ) are achieved in close field regions , a natural result of energy flowing on a closed path via heat conduction.,title="fig:",scaledwidth=7.0% ]     ( top ) , @xmath109 ( middle ) , and @xmath110 ( bottom)@xmath111 for run d ( left ) and the standard sc model ( right ) for cr 1913 ( black lines at @xmath107 ) .
this demonstrates the high dynamic range in density over small changes in radius resolved with the lc model.,title="fig:",scaledwidth=44.0% ]   ( top ) , @xmath109 ( middle ) , and @xmath110 ( bottom)@xmath111 for run d ( left ) and the standard sc model ( right ) for cr 1913 ( black lines at @xmath107 ) .
this demonstrates the high dynamic range in density over small changes in radius resolved with the lc model.,title="fig:",scaledwidth=44.0% ]   ( top ) , @xmath109 ( middle ) , and @xmath110 ( bottom)@xmath111 for run d ( left ) and the standard sc model ( right ) for cr 1913 ( black lines at @xmath107 ) .
this demonstrates the high dynamic range in density over small changes in radius resolved with the lc model.,title="fig:",scaledwidth=44.0% ]   ( top ) , @xmath109 ( middle ) , and @xmath110 ( bottom)@xmath111 for run d ( left ) and the standard sc model ( right ) for cr 1913 ( black lines at @xmath107 ) .
this demonstrates the high dynamic range in density over small changes in radius resolved with the lc model.,title="fig:",scaledwidth=44.0% ]   ( top ) , @xmath109 ( middle ) , and @xmath110 ( bottom)@xmath111 for run d ( left ) and the standard sc model ( right ) for cr 1913 ( black lines at @xmath107 ) .
this demonstrates the high dynamic range in density over small changes in radius resolved with the lc model.,title="fig:",scaledwidth=44.0% ]   ( top ) , @xmath109 ( middle ) , and @xmath110 ( bottom)@xmath111 for run d ( left ) and the standard sc model ( right ) for cr 1913 ( black lines at @xmath107 ) .
this demonstrates the high dynamic range in density over small changes in radius resolved with the lc model.,title="fig:",scaledwidth=44.0% ]   ( top ) , @xmath109 ( middle ) , and @xmath110 ( bottom)@xmath111 for run d ( left ) and the standard sc model ( right ) for cr 1913 ( black lines at @xmath107 ) .
this demonstrates the high dynamic range in density over small changes in radius resolved with the lc model.,title="fig:",scaledwidth=7.0% ]     ( right ) .
the blue and green points represent euv data on an off the limb ( binned to 1000 points each ) .
the curve path follows the ratio values [ @xmath99 and @xmath112 $ ] from @xmath113k , the width of the curves represents the variation of the ratio values from @xmath114 at a given temperature .
this is intended to illustrate the large variation of filter ratio values with respect to an unresolved dem distribution.,title="fig:",scaledwidth=48.5% ]   ( right ) .
the blue and green points represent euv data on an off the limb ( binned to 1000 points each ) .
the curve path follows the ratio values [ @xmath99 and @xmath112 $ ] from @xmath113k , the width of the curves represents the variation of the ratio values from @xmath114 at a given temperature .
this is intended to illustrate the large variation of filter ratio values with respect to an unresolved dem distribution.,title="fig:",scaledwidth=49.0% ]     ( right ) . while the location of the synthesized filter ratios improves with this dem modification ,
the fine details of the local dem distribution , realized by the variation of observed eit filter ratio values , are difficult to reproduce with an mhd model involving large , global scales.,title="fig:",scaledwidth=49.0% ]   ( right ) .
while the location of the synthesized filter ratios improves with this dem modification , the fine details of the local dem distribution , realized by the variation of observed eit filter ratio values , are difficult to reproduce with an mhd model involving large , global scales.,title="fig:",scaledwidth=49.0% ]