text: Geophysical Research Letters Generation of shock lamellae and melting in rocks by lightning-induced shock waves and electrical heating Jiangzhi Chen1 , Chiara Elmi1 , David Goldsby1 , and Reto Gieré1 1Department of Earth and Environmental Science, University of Pennsylvania, Philadelphia, Pennsylvania, USA Abstract The very rapid energy release from impact events, such as those resulting from lightning strikes or meteorites, can drive a variety of physical and chemical processes which alter rocks and result in the formation of natural glasses (i.e., fulgurites and tektites). Fulgurite is the vitrified soil, sand, or rock resulting from lightning strikes. A thunderbolt is associated with air temperatures of up to 105 K, which can heat rocks to >2000 K within tens of microseconds. The rapid fusing and subsequent quenching of the surface of the rock leaves a distinctive, thin, garbled coating composed of a glassy to fine-grained porous material. Previous studies on rock fulgurites found planar deformation features in quartz crystals within the target rock substrate, evidence of strong shock waves during fulgurite formation. In this paper, we simulated the shock pressure and temperature caused by an idealized lightning impact on rocks and compared the model results with observations on rock fulgurites from the literature. Our model results indicate that a lightning strike can cause >7 GPa pressure on the rock surface, generate a layer of fulgurite (of radius ∼9 cm), and leave a burned region (of radius ∼11 cm). The fulgurites found on rock surfaces share many features with sand fulgurites, but their spatial distribution is completely different, as sand fulgurites are hollow tube structures. Our study on rock fulgurites provides an indirect constraint on the energy of a lightning event and also demonstrates that the presence of shock features in rocks cannot be taken as unequivocal evidence for impact events. 1. Background A lightning strike is a complicated phenomenon composed of three stages. First, a "stepped leader" descends from the cloud and contacts an oppositely charged target (e.g., the ground, architectural monuments, or other structures). Second, a "dart leader" travels back to the cloud, followed by a return stroke, i.e., the main flash, which travels from the cloud and strikes the target again. The return stroke, with an energy density of 3.3 MJ/m within a channel of radius ∼1 cm [Jones et al., 1968], is associated with air temperatures of up to 105 K [Uman, 1964] and a current of up to 30 kA [Uman and Krider, 1989]. The overall energy of the lightning strike is estimated to be 109 –1010 J [Rakov and Uman, 2007]. When lightning strikes a target, such as a rock, the target can be heated to over 2000 K within tens of microseconds [Uman and Krider, 1989]. The resulting vitrified soil, sand, or rock is called a fulgurite. Pasek et al. [2012] classified fulgurites into four types, depending on whether the lightning strikes sand, clay, or rock. Sand fulgurites are usually found in deserts, featuring subsurface, branching, and hollow glass tubes [Rakov and Uman, 2007]. Occasionally, sand fulgurites can reach several meters in length (e.g., one fulgurite found in Florida was ∼5 m long [see Rakov and Uman, 2007]). Rock fulgurites, however, appear as a distinctive, thin, garbled coating composed of glassy to fine-grained porous material [Essene and Fisher, 1986; Grapes and Müller-Sigmund, 2010; Gieré et al., 2015; Elmi et al., 2017] and may include burned organic material present on the rock surface before lightning impact [Elmi et al., 2017]. In some cases, fulgurites are found on man-made structures, such as windmills [Brocklesby, 1869, see pp. 148]. Fulgurites can also be formed artificially by high-voltage electrical arcing [Williams and Johnson, 1980; Uman and Krider, 1989; Kumazaki et al., 1993; Pasek and Hurst, 2016]. Our model focuses on rock fulgurites, particularly, fulgurites from the granite surface at Les Pradals, France [Gieré et al., 2015] (Figure 1b) and Mount Mottarone, Italy [Elmi et al., 2017] (Figure 1c). Besides high temperature, the lightning strike also generates a strong shock wave in the vicinity of the strike point, evidenced by the presence of shock lamellae in the target rock, the so-called "planar deformation RESEARCH LETTER 10.1002/2017GL073843 Key Points: • A cylindrical air dynamics model was employed to estimate the shock wave pressure caused by a lightning strike on rocks • An idealized point current source model was employed to calculate the electrical heating and subsequent cooling of rock fulgurites • Model results were compared with literature data collected on natural rock fulgurites to explain the existence of shock lamellae Correspondence to: J. Chen, jiangzhi@sas.upenn.edu Citation: Chen, J., C. Elmi, D. Goldsby, and R. Gieré (2017), Generation of shock lamellae and melting in rocks by lightning-induced shock waves and electrical heating, Geophys. Res. Lett., 44, doi:10.1002/2017GL073843. Received 18 APR 2017 Accepted 7 JUL 2017 Accepted article online 14 JUL 2017 ©2017. American Geophysical Union. All Rights Reserved. CHEN ET AL. LIGHTNING SHOCK AND FULGURITE FORMATION 1 Geophysical Research Letters 10.1002/2017GL073843 Figure 1. (a) Schematic of lightning hitting granite at 90∘. The lightning causes a cylindrical shock wave in the air, leaving microscale shock features in the quartz crystals on the surface. Electrical heating causes the temperature to rise, melts the rock incongruently, and generates a thin fulgurite layer on the surface when the temperature is above the granite melting point. The organic material on the weathered granite surface within a certain distance of the strike point is burned and leaves a distinctive black mark. (b) Rock fulgurite from Les Pradals, France [Gieré et al., 2015]. (c) Rock fulgurite from Mount Mottarone, Italy. features" (PDFs) in crystalline quartz [Gieré et al., 2015]. Shock lamellae are also associated with meteorite impacts [Bohor et al., 1984; Goltrant et al., 1992; Langenhorst and Deutsch, 1994; Alvarez et al., 1995], but lightning impacts may also produce shock features in rocks [Gieré et al., 2015]. Analysis of shocked crystalline quartz suggests shock pressures>10 GPa [Grieve et al., 1996] and sometimes∼22–30 GPa in magnitude [Langenhorst and Deutsch, 1994; Carter et al., 2010a; Ende et al., 2012]. However, the temperature rise associated with the lightning strike must be relatively moderate, likely lower than the quartz annealing point (1200 ∘C [Grieve et al., 1996]), otherwise the observed PDFs would disappear. It must be emphasized that the transmission of energy from the lightning bolt to the rock during shock deformation and electrical heating, and the accompanying physical and chemical changes in the rock, absorb only a fraction of the total energy of the lightning, with the rest of the energy dissipated in the air in the form of thunder, heated air, light, and electromagnetic waves. The aim of this paper is to provide an idealized model of the shock wave, as well as the rapid electrical heating and subsequent cooling, in a granitic rock. The model estimates the pressure and temperature disturbance caused by a lightning strike and sheds light on the generation of shock lamellae and the energy transferred from the lightning to the target rock. The model is compared with field observations and experimental analyses of granite fulgurite samples from Gieré et al. [2015] and Elmi et al. [2017]. Using our model, we constrain various physical parameters associated with the lightning strike on rocks, including the energy density of the lightning bolt, and estimate the physical conditions and parameters controlling the formation of shock lamellae in fulgurites. 2. Model We employed typical parameter values for the lightning and for the granite rock, as given in Table 1. The lightning is assumed to strike with a 90∘ angle of incidence with the rock surface (Figure 1a). We treat the shock wave around the lightning as being excited from a very rapid energy release from the lightning channel, which is treated as a line source. Estimates of the energy density of the lightning have a wide range of 103 to 105 J/m obtained from radiation, electrostatics, or optics [see Rakov and Uman, 2007, Table 12.1]. Collins et al. [2012] suggested a much higher energy density range, 10–50 MJ/m, based on estimates from arcing experiments. Pasek and Hurst [2016] concluded that lightning strikes with an energy density >6 MJ/m may account for only 3–4% of all fulgurites formed annually. We take a value of 3.3 MJ/m after Jones et al. [1968]. The temperature of the lightning channel Te is set to 105 K [Uman, 1964]. CHEN ET AL. LIGHTNING SHOCK AND FULGURITE FORMATION 2 Geophysical Research Letters 10.1002/2017GL073843 Table 1. Nominal Parameter Values of the Lightning and the Granite Model Parameters Value Lightning energy density E (J/m) 3.3 × 106 channel radius a (m) 0.01 peak current intensity I 0 (A) 2 × 104 rise time tr (s) 10−6 total duration tm (s) 50 × 10−6 decay constant k (s−1) 5 × 104 temperature Te (K) 3 × 104 Air specific heat ratio 훾 1.4 ambient temperature Ta ( ∘C) 25 Rainfall density 휌w (kg/m3) 1000 specific heat Cw (J/(kg K)) 4184 precipitation rate v (cm/h) 0.76 vaporization temperature Tw ( ∘C) 100 enthalpy of vaporization ΔvH (J/kg) 2.26 × 106 Granite density 휌 (kg/m3) 2670 specific heat C (J/(kg K)) 790 melting temperature Tg (K) 973 thermal diffusivity 휅 (m2/s) 0.913 × 10−6 electrical conductivity 휎 (S/m) 10−3 Quartz annealing temperature Tq ( ∘C) 1200 Organic material combustion temperature Tc ( ∘C) 300 The air is treated as an isentropic ideal gas with constant specific heats Cp and Cv, and the specific heat ratio 훾 = Cp∕Cv = 1.4. Rock surfaces are in general somewhat weathered before lightning strikes, with typically higher thermal and electrical conductivities in the weathered part than in the unaltered rocks due to the presence of moisture. However, Elmi et al. [2017] showed that the weathered layer on granite from Mottarone appears to be thinner than 1 mm, whereas the rock itself is tens of centimeters in thickness or larger. We treated the rock as a homogeneous half-space with uniform thermal and electrical properties. The thermal diffusivity 휅 of granite is 0.913 mm2/s, the specific heat C is 790 J/(kg K), and the density 휌 is 2670 kg/m3 [Robertson, 1988]. The melting point of the granite Tg is dependent on the water content, and we used a value of 973 K for wet granite [Eppelbaum et al., 2014]. The annealing temperature Tq of quartz, beyond which the shock lamellae are destroyed, is ∼1200 ∘C [Grieve et al., 1996]. The electrical conductivity 휎 of granite is taken as 10−3 S/m [Olona et al., 2010]. The parameters are treated as constants, but we further discuss alteration of air and rock properties at high temperatures and in strong electric fields below. The combustion temperature of organic matter (e.g., lichens and other plant materials) on the rock surface is unknown. However, for wood, dehydration and charring commences at ∼100–200∘C, whereas combustion occurs at ∼200–300∘C [White and Dietenberger, 2001]. We assume a combustion temperature Tc of organic material of 300 ∘C. After the lightning strike and associated heat pulse, the rock surface cools down due to evaporation of rain drops and convection of heated air. Meteorological agencies usually use a rain water accumulation rate v of 0.76 cm/h as a boundary between moderate and heavy rain, which we adopt here as the precipitation rate, and the cooling rate per unit area due to rain water evaporation on the rock surface is v휌w[Cw(Tw − Ta)+ΔvH] where the evaporation temperature Tw = 100 ∘C, and the enthalpy of vaporization of water at Tw is ΔvH=2.26 MJ/kg [Marsh, 1987]. We ignore the advective heat transfer by rain water flowing over the rock surface, because significant water flow only exists after the rock surface is below Tw. 2.1. The Pressure Disturbance Caused by the Shock Wave Gieré et al. [2015] reported that shock lamellae were only found within the top 3 μm of the rock surface, suggesting that the most relevant pressure disturbance is from the shock wave around the lightning bolt. We ignore the bow shock at the tip of the lightning flash in our model, because the return stroke (main flash) travels at a speed of ci = 1−2 × 108 m/s [Rakov and Uman, 2007], greatly exceeding the speed of sound in air, CHEN ET AL. LIGHTNING SHOCK AND FULGURITE FORMATION 3 Geophysical Research Letters 10.1002/2017GL073843 cs≈340 m/s. After the rock is struck by the lightning, the trailing shock wave hits the surface, and body waves propagate inside the rock; however, no shocked quartz has yet been observed deep in the rock interior (i.e., >3 μm deep), indicating that the effect of these body waves can be neglected. The pressure from the magnetic field is pB = 휇0I 2∕(8휋2r 2) around the line current I, where 휇0 = 4휋 × 10−7 N/A2 is the vacuum magnetic permeability. We ignore this pressure because it is only ∼104 Pa at the radius of the lightning channel, which is much smaller than the shock pressure. 2.1.1. Cylindrical Shock Wave Around the Lightning Lightning traverses in air can extend to tens or even hundreds of kilometers in length [Rakov and Uman, 2007], and the longest-distance flash reported traveled 321 km [Lang et al., 2016]. Until the lightning hits the ground, the radial shock wave can be treated as cylindrical, and the wave only propagates in the radial direction perpendicular to the lightning channel. Assuming that the air is an inviscid ideal gas, the Navier-Stokes equation with cylindrical symmetry is Du Dt = 휕u 휕t + u 휕u 휕r = −1 휌 휕p 휕r , (1) where u(r,t) is the air flow velocity, p(r,t) is the pressure, and 휌(r,t) is the air density. The continuity equation is D휌 Dt = 휕휌 휕t + u 휕휌 휕r = −휌 (휕u 휕r + u r ) , (2) and another governing equation is the isentropic relation for air D Dt ( p 휌훾 ) = ( 휕 휕t + u 휕 휕r ) ( p 휌훾 ) = 0. (3) The undisturbed air pressure is p0, the air density is 휌0, and cs= √훾p0∕휌0 is the speed of sound in undisturbed air. The shock wave has a radius R propagating at a speed U=dR∕dt. From the Rankine-Hugoniot condition, the disturbed pressure p1, density 휌1, and velocity u1 immediately behind the shock wave satisfy 휌1 휌0 = 훾 − 1 + p1(훾 + 1)∕p0 훾 + 1 + p1(훾 − 1)∕p0 ≈ 훾 + 1 훾 − 1 (4) U2 c2 s = 훾 − 1 + p1(훾 + 1)∕p0 2훾 ≈ p1(훾 + 1)∕p0 2훾 (5) u1 U = 2(p1∕p0 − 1) 훾 − 1 + p1(훾 + 1)∕p0 ≈ 2 훾 + 1 (6) with the approximations valid for strong shocks where p1∕p0 → ∞. Following the dimensional treatment of Taylor [1950a], the shock radius can be expressed as R = S(훾)E1∕4휌 −1∕4 0 t 1∕2 (7) where S(훾) is a dimensionless function depending on the specific heat ratio 훾 only. We solve the governing equations (1) to (3) using the self-similar approach as in Lin [1954] (or see Appendix A), and the results are presented in Figure 2. Numerical evaluations show that the shock radius is R = ( 4E B휌0 )1∕4 t 1∕2 (8) where the constant B≈3.94. The pressure behind the shock front is p1 = 2E BR2(훾 + 1) = 0.212ER−2. (9) CHEN ET AL. LIGHTNING SHOCK AND FULGURITE FORMATION 4 Geophysical Research Letters 10.1002/2017GL073843 Figure 2. Normalized air pressure, normalized velocity, and normalized density as functions of radial position r/R(t) inside the shock tube. Near the shock front, as r approaches R, both the pressure and density of air increase rapidly, whereas the velocity profile remains almost linear. If we take the energy density E=3.3 MJ/m [Jones et al., 1968], and if R∼a=1 cm is the radius of the lightning channel, then the shock pressure is p1 ≈ 2E Ba2(훾 + 1) = 7.0GPa (10) which is smaller than the pressure range (≥10 GPa) required to generate shock lamellae in quartz [Grieve et al., 1996]. However, if the energy density E is higher, for example, 10 MJ/m, the lower bound suggested by Collins et al. [2012], then p1 = 0.212ER−2 is large enough to generate a shock pressure of 21 GPa at the boundary of the lightning channel, sufficient to generate shocked quartz, and similar to pressures induced by meteorite impacts. 2.1.2. Shock Wave in the Rock Due to the short duration of the lightning strike, the wave propagation in solid media resulting from the lightning strike can be treated as having been excited from a point explosion in an elastic half-space. As a simplified estimate, the pressure disturbance in the rock decays with distance R from the hit point as p ∼ 0.156 a3 R3 2E Ba2(훾 + 1) . (11) The R−3 falloff of pressure means that pressures in the interior of the rock will be far too small to form shocked quartz when R is greater than ∼4 mm. We also ignore strain heating here because the strain energy release and associated heating must be smaller than electrical heating due to the strike, as discussed below. 2.2. Thermal Evolution of the Rock During and After the Lightning 2.2.1. Electrical Heating The electrical current through the rock is assumed to follow a pulse form I(t) = I 0훽(t) [Plooster, 1971] 훽(t) = { t∕tr, for t < tr exp [ −k(t − tr) ] , for tr < t < tm (12) where the peak current intensity I 0 =20 kA, the decay constant k=5 × 104s−1, and the duration tm =50 μs. The ∼1 cm radius lightning channel is so narrow that it can be treated as a point current source. We set the origin of the cylindrical coordinate system at the lightning strike point, with the z axis pointed downward into the rock. Since the rock is assumed to be homogeneous, the current is evenly distributed in a half sphere of area 2휋(r 2 + z2), so that the electric field E only has a radial component Er = I 2휋휎(r2 + z2) . (13) CHEN ET AL. LIGHTNING SHOCK AND FULGURITE FORMATION 5 Geophysical Research Letters 10.1002/2017GL073843 The electric field has a singularity at r = z = 0 due to the point current source treatment of the lightning. One way to remove this singularity is to use a lightning channel with a finite radius of ∼1 cm. Then the electric field in the vicinity of the lightning channel is dominated by a vertical component Ez ≈ I 휋휎a2 , (14) whereas for regions far from the lightning channel, the electric field is still Er. These two expressions differ by a factor of 2 because the area changes from that of a circle, 휋a2, to that of a half sphere, 2휋a2 at the channel boundary. A more detailed explanation is provided in Appendix B. The temperature evolution during the lightning strike is given by 휕T 휕t = 휅∇2T + 휎 휌C |E| 2. (15) The boundary and initial conditions are T | | | | |r,z→0 ⇒ finite, T | | | | |r,z→∞ = Ta, 휅휌C 휕T 휕z | | | | |z=0 = Qv + Qr, T | | | | |t=0 = Ta, (16) where the heat loss rate due to rainfall vaporization is Qv = v휌w[Cw(Tw − Ta)+ΔvH], and the heat radiation rate is Qr =휎SB ( T4 − T4 a ) where 휎SB =5.67 × 10−8 W/(m2 K4) is the Stefan-Boltzmann constant. We ignore the electromagnetic effects related to a time-dependent electric field and discuss the possible effects in later sections. The equation is inhomogeneous with a nonlinear boundary condition. During the lightning strike, the distance inside the rock over which heat conduction is significant isΔr∼√휅tm ≈7 mm, and the rate of heat conduction is 휅∇2T ≈휅(Te−Ta)∕Δr 2. Similarly, sinceQr ∼ 휎SBT4 e ≫ Qv, the heat loss rate at the surface near the hit point is of order 휎SBT4 e ∕(휌CΔr). These two processes are negligible compared with the electrical heating rate 휎|E| 2∕(휌C). Therefore, the temperature after the lightning strike is dominated by electrical heating Ti ≈ ∫ tm 0 휎 휌C |E| 2dt + Ta. (17) This equation can be modified readily to remove the singularity at the origin and make the temperature continuous Ti − Ta ≈ ⎧ ⎪ ⎨ ⎪ ⎩ AI2 0tm 4휋2휌C휎R4 , for R ≥ a∕ √ 2 AI2 0tm 휋2휌C휎a4 , for R < a∕ √ 2, (18) where R2 = r 2 + z2 and the constant A = t −1 m ∫ tm 0 훽2dt ≈ 0.205. The temperature at the surface as a function of r can be compared with the size of the burned regions observed on fulgurite surfaces. The region within which the granite is expected to be partially melted is obtained by solving Ti (r) = Tg, which gives a radius of ∼9.3 cm for fulgurite presence, and the radius of burned organics is∼11.6 cm, consistent with observations on Mottarone granite [Elmi et al., 2017]. On the other hand, the quartz annealing temperature is reached inside a region of radius ∼8 cm. The total increase in the internal energy of the rock is ΔU = 1 2 휌C ⎡ ⎢ ⎢ ⎢ ⎣ ∫∫∫ R<a∕ √ 2 ( Ti − Ta ) dV + ∫∫∫ R≥a∕ √ 2 ( Ti − Ta ) dV ⎤ ⎥ ⎥ ⎥ ⎦ = 2 √ 2AI2 0tm 3휋휎a = 1.2 × 108J, (19) but a fraction of the overall energy of the lightning, ∼109 –1010 J [Rakov and Uman, 2007]. 2.2.2. Cooling of the Rock by Convective Heat Loss at the Surface After the lightning strike, the temperature in the rock evolves from Ti by heat conduction and convective heat loss at the surface. The governing equation of cooling is 휕T 휕t = 휅∇2T (20) CHEN ET AL. LIGHTNING SHOCK AND FULGURITE FORMATION 6 Geophysical Research Letters 10.1002/2017GL073843 with the boundary and initial conditions T | | | | |r,z→0 ⇒ finite, T | | | | |r,z→∞ = Ta, 휅휌C 휕T 휕z | | | | |z=0 = Qv + Qr, T | | | | |t=0 = Ta. (21) Because the surface heat loss rate 휎SBT4∕(휌CΔr) ≫ 휅∇2T ∼휅(T − Ta)∕Δr 2, the surface quenches quickly and barely any heat is conducted from the hot strike center to the cooler periphery. Therefore, the peak temperature in the rock is reached right after the lightning strike stops, after which the temperature on the surface keeps dropping until the surface solidifies. 3. Discussion We assume that the air can be treated as an ideal gas, and the ratio of specific heats 훾 remains constant, which may be unrealistic because the ionization and dissociation occurring at high temperature will change the apparent value of 훾. As Taylor [1950b] pointed out, however, when comparing the same model as ours against the data from the first atomic explosion in New Mexico, the intense radiation from the center of the explosion can cause the apparent value of 훾 to rise, while ionization and dissociation at very high temperatures can decrease 훾. The overall effect is that the whole system may behave as though the value of 훾 remains unchanged. The lightning bolt is hot ionized plasma, and the lightning channel, where strong ionization occurs, is very narrow, which suggests that radiation and ionization, facilitated by extremely high temperatures, are confined to much smaller regions than in atomic explosions. As a result, we expect that the deviation of the state of the air from ideality away from the lightning channel is insignificant. We modeled the pressure exerted on a granite surface, and the temperature evolution in the rock, due to a lightning strike. The pressure model predicts a maximum pressure of ∼7 GPa, smaller than the pressure thought to be required to generate shocked quartz, ∼10 GPa. One possible reason for this discrepancy is that the energy density of the lightning channel is higher (>10 MJ/m) than the assumed value of 3.3 MJ/m suggested by Jones et al. [1968]. There exists a large variation in estimates of the energy density of lightning. Our model gives an indirect way to quantify a lower bound for the lightning energy density. In the previous sections, our models rely on two basic parameters of the lightning, the energy density E and the lightning current I. These two parameters should be related; however, the details of that relationship remain unclear because of the difficulty in measuring them. Based on our model, if a lightning bolt strikes granite, leaving shock lamellae on its surface after cooling down, it is possible to provide some constraint on this relationship. The formation of shock lamellae requires that, on the rock surface, the pressure disturbance at a distance r from the hit point is greater than Δp ∼10 GPa [Gieré et al., 2015; Langenhorst and Deutsch, 1994; Ende et al., 2012] 2E Br2(훾 + 1) > Δp, (22) while the temperature must be lower than the quartz annealing temperature Tq for the shock lamellae to survive ΔT = Tq − Ta > AI2 0tm 4휋2휌C휎r4 . (23) The result is that E I 0 > ΔpB(훾 + 1) 4휋 √ Atm 휌C휎ΔT . (24) For the annealing temperature Tq = 1200∘C, this ratio is 1.53×104 J/(A m). Therefore, if the peak current of the lightning is in the range 2–30 kA, the energy density is 30–460 MJ/m. Although the upper limit appears to be too large, the lower limit lies in the 10–50 MJ/m range proposed by Collins et al. [2012], which would yield a shock pressure greater than 7 GPa. The large upper limit can be attributed to the high temperature caused by intense Joule heating. There are several possible mechanisms that can decrease the temperature and the size of the annealed region, thus allowing for the survival of the PDFs. CHEN ET AL. LIGHTNING SHOCK AND FULGURITE FORMATION 7 Geophysical Research Letters 10.1002/2017GL073843 Our model of temperature evolution in the rock yields high temperatures (>105K) near the strike point, assuming constant thermal and electrical properties. In reality, however, a decrease in moisture content with heating reduces the electrical conductivity, whereas after the temperature reaches the solidus of the rock, the increased electron activity in the molten rock greatly enhances electrical conductivity [Olhoeft, 1981; Roberts and Tyburczy, 1999], reducing the electrical heating effect. Also, near the lightning strike point, the modeled electric field is greater than 1011 V/m, and dielectric breakdown will occur. However, microstructural observations [Elmi et al., 2017] confirm that these processes mainly take place at the rock surface, and the ∼R−4 decay of temperature and ∼ R−2 decay of the electric field ensure that the material alteration is mostly con- fined to the close vicinity of the strike point. Our model predicts a radius of the burned region (11.6 cm) that is consistent with observed sizes of burned areas on fulgurite samples from Mottarone [Elmi et al., 2017]. Thus, neglecting these variations in thermal and electrical properties with temperature and field strength does not change the overall behavior of our model. Also as stated above, the Mottarone granite has an apparent weathered layer<1 mm="" in="" thickness="" but="" the="" weathering="" may="" actually="" extend="" much="" deeper="" into="" rock="" and="" associated="" fluid-filled="" cracks="" can="" increase="" overall="" heat="" capacity="" c="" even="" more="" significantly="" electrical="" conductivity="" olona="" et="" al="" 2010="" suggested="" that="" of="" weathered="" granite="" be="" an="" order="" magnitude="" higher="" than="" pristine="" which="" would="" reduce="" temperature="" by="" according="" to="" equation="" 18="" radii="" areas="" prone="" annealing="" quartz="" melting="" burning="" organics="" will="" shrink="" a="" factor="" 101="" 4="" we="" also="" neglect="" coupling="" between="" time-varying="" electric="" magnetic="" fields="" adopted="" pulse-form="" current="" with="" simple="" frequency="" spectrum="" real="" lightning="" events="" have="" wide="" range="" from="" 300="" khz="" data-mce-fragment="1">30 MHz [Lan et al., 2011]. Bowler [2004] solved the frequency-dependent electric field distribution in a homogeneous half-space, and the results show that the field strength has an exponentially decaying component with r and z, and consequently, electrical heating decreases faster than R−4 as in our model. These two effects— the increase in electrical conductivity in the weathered granite and coupling of the electric and magnetic fields—may lower the temperature rise during the lightning strike compared to our model predictions. This would result in a smaller radius of material around the strike point within which the annealing temperature of quartz is obtained, ensuring survival of the shock features. Rock fulgurites are different from sand fulgurites, in which air fills the pore spaces in the sand (∼40% porosity for a random loose packing scenario) and acts as an insulator for the current. In sand fulgurites, heating done by the current only exists where dielectric breakdown occurs. In comparison, the electrical conductivity of granite is 10−3 S/m, much higher than that of air (∼10−15 S/m), thus the current is able to pass through the rock and induce significant heating. Without pervasive electrical heating, thermal diffusion is the major mechanism for the formation of the melt tube [Carter et al., 2010b; Pasek et al., 2012] with a radius ∼3 cm (Appendix C), consistent with observed sand fulgurite tubes [Rogers, 1946; Gailliot, 1980; Pye, 1982; Essene and Fisher, 1986; Navarro-González et al., 2007], but much smaller than the fulgurite region on rocks from Les Pradals [Gieré et al., 2015] and Mount Mottarone [Elmi et al., 2017]. Previous samples collected from the granite surface at Les Pradals, France indicate that shock lamellae are only observed in a surficial layer of the target <3 m="U∕cs" in="" thickness="" gier="" et="" al="" 2015="" for="" an="" ideal="" flat="" rock="" surface="" the="" main="" flash="" moves="" at="" ci="" u="" and="" angle="" between="" propagation="" direction="" of="" shock="" wave="" is="" 10="" 6="" small="" ensures="" that="" a="" quartz="" crystal="" on="" vertical="" force="" exerted="" only="" horizontal="" pressure="" disturbance="" can="" affect="" thin="" layer="" before="" being="" significantly="" damped="" this="" may="" explain="" why="" shocked="" has="" been="" observed="" depths="" much="" smaller="" than="" millimeter-scale="" topographic="" fluctuations="" reality="" lightning="" not="" strike="" exactly="" 90="" to="" from="" equation="" 5="" disturbed="" p1="" p0="" 2="" m2="" 1="" where="" data-mce-fragment="1">1 is the Mach number. The Mach number will be higher for supersonic air flow passing a convex corner, which leads to higher pressure near topographic depressions. But even in that case, the pressure disturbance is still limited to a shallow depth near the rock surface. The lightning strike can introduce hot air or gaseous reaction products into the molten rock which are then preserved as bubbles. Air inclusions may contain traces of ancient air which may be helpful in understanding paleoclimate [Navarro-González et al., 2007]. High temperatures can also vaporize water and alkalies at the rock surface, and the burning of organic material generates NOx and COx gases [Elmi et al., 2017], facilitating the formation of a dense population of bubbles, preserved during quenching and responsible for the porous nature of rock fulgurites. CHEN ET AL. LIGHTNING SHOCK AND FULGURITE FORMATION 8 Geophysical Research Letters 10.1002/2017GL073843 Table 2. Calculated Values of Physical Parameters for Granite Fulgurites Calculated result Value Radius of burned region 11.6 cm Radius of fulgurite region 9.3 cm Minimum radius for PDF survival 8 cm Heating energy 1.2 × 108 J Maximum temperature ∼105K Maximum field strength ∼1011 V/m Shock pressure 7 GPa (for E = 3.3 MJ/m) to 21 GPa (for E = 10 MJ/m) Shock inclination angle ∼10−6 rad Calculated energy density of lightning 30–460 MJ/m Minimum ratio of energy density to current 1.53 × 104 J/(A m) 4. Conclusion When a lightning bolt strikes a granite target, the shock wave sweeps the rock surface, and the distributed high current heats the rock. The shock wave creates shock lamellae in the substrate near the rock surface, and the electrical heating vaporizes material at the strike point, melts the granite (which is later quenched to form fulgurite), and burns organic material on a larger area of the surface. Our temperature and pressure models demonstrate that based on the observations of the spatial distribution of shocked quartz and burned organic material, we can constrain the energy and current intensity of the lightning event. We summarize the modeled parameter values for these and other physical parameters for granite fulgurites in Table 2. Our results further demonstrate that the presence of shock features in rocks cannot be taken as prima facie evidence for meteorite impacts. Appendix A: Self-Similar Solution of the Cylindrical Shock Wave Let the normalized location 휂 = r∕R(t), and the pressure, density, and air velocity of the shock wave are assumed to be related to the undisturbed air as p p0 = U2 c2 s f(휂), 휌 휌0 = 휓(휂), u U = 휙(휂). (A1) Then equations (1) to (3) become 휙′ (휂 − 휙) = f′ 훾휓 − 휙 (A2) 휓′ 휓 = 휙′ + 휙∕휂 휂 − 휙 (A3) 2f + 휂f ′ + 훾휓′ 휓 f(휙 − 휂) = f ′ 휙 (A4) with the boundary conditions f(1) ≈ 2훾 훾 + 1 , 휙(1) ≈ 2 훾 + 1 , 휓(1) ≈ 훾 + 1 훾 − 1 . (A5) The system of ordinary differential equations (A2) to (A4) can then be solved numerically. The shock radius R is determined by E = 2휋 ∫ R 0 휌 ( CvT + u2 2 ) rdr = 2휋 ∫ R 0 ( p 훾 − 1 + 휌u2 2 ) rdr = 휋휌0U2R2 ∫ 1 0 [ 2f 훾(훾 − 1) + 휓휙2 ] 휂d휂 = B휌0U2R2, (A6) where B=3.94. CHEN ET AL. LIGHTNING SHOCK AND FULGURITE FORMATION 9 Geophysical Research Letters 10.1002/2017GL073843 It is possible to apply the method to a spherically symmetric point explosion, where the energy release E′ is related to the shock wave speed U and the shock radius R as E′ = 5.33휌0U2R3. (A7) The pressure disturbance therefore is p = 0.156E′ R−3. (A8) Appendix B: Current Injection With Finite Width in Homogeneous Medium When the medium is homogeneous, the governing equation for the electric potential in the cylindrical coordinate system is ( 휕2 휕r2 + 1 r 휕 휕r + 휕 휕z2 ) 휙 = 0 (B1) with the boundary conditions 휙 | | | | |r=0 ⇒ finite, 휙 | | | | |r,z→∞ = 0. (B2) The solving process is a routine separation of variables, and the general solution is 휙 = ∫ ∞ 0 C휆J0(휆r) exp (−휆z) d휆, (B3) where J0 is the Bessel function of first kind with order zero. The current density at the surface is − 휎 휕휙 휕z | | | | |z=0 = 휎0 ∫ ∞ 0 C휆휆J0(휆r)d휆 = I 휋a2 H(a − r). (B4) where H(x) is the Heaviside step function. Using a Hankel transform a ∫ ∞ 0 J1 (휆a)J0 (휆r) d휆 = H(a − r) (B5) the solution of the potential is 휙 = I 휋a휎 ∫ ∞ 0 exp (−휆z) 휆 J1 (휆a)J0 (휆r) d휆. (B6) As Bowler [2004] pointed out, when the channel radius a ≪ r, lima→0 J1 (휆a) 휆a = 1 2 (B7) so lima→0 휙 = I 2휋휎 ∫ ∞ 0 exp (−휆z)J0 (휆r) d휆 = I 2휋휎√ r2 + z2 (B8) which is the point source potential field. However, the approximation has singularities for the potential and electric field at r, z → 0, while direct evaluation of the electric field at the surface gives 휙(r, z = 0) = I 휋a휎 ∫ ∞ 0 1 휆 J1 (휆a) d휆 = I 휋a휎 Ez(r, z = 0) = I 휋a휎 ∫ ∞ 0 J1 (휆a) d휆 = I 휋a2휎 Er(r, z = 0) = 0. (B9) CHEN ET AL. LIGHTNING SHOCK AND FULGURITE FORMATION 10 Geophysical Research Letters 10.1002/2017GL073843 Therefore, the current density, which is the current intensity per unit surface area, has two different approximations in the vicinity of the origin, depending on whether we treat the "area" of the small region as a circle or as a half sphere. Appendix C: Thermal Diffusion in Sand Fulgurites For sand fulgurites, thermal diffusion is dominant. Carter et al. [2010b] and Pasek et al. [2012] treated sand fulgurites as generated by a line source, and the temperature evolution in the cylindrical coordinate system is 휅∇2T = 휕T 휕t . (C1) For the ideal case where the line source is instantaneous, and the temperature within the lightning channel radius a is the peak temperature Te, the solution is [Carslaw and Jaeger, 1959] T = Ta + a2(Te − Ta) 4휅t exp ( − r 2 4휅t ) . (C2) If we set the temperature to the melting temperature of quartz sand∼1650 ∘C, the maximum radius of molten sand from the lightning path is about 3 cm. References Alvarez, W., P. Claeys, and S. W. Kieffer (1995), Emplacement of Cretaceous-Tertiary boundary shocked quartz from Chicxulub crater, Science, 269(5226), 930–935. Bohor, B. F., E. E. Foord, P. J. Modreski, and D. M. Triplehorn (1984), Mineralogic evidence for an impact event at the Cretaceous-Tertiary boundary, Science, 224(4651), 867–869. Bowler, N. (2004), Analytical solution for the electric field in a half space conductor due to alternating current injected at the surface, J. Appl. Phys., 95(1), 344–348. Brocklesby, J. (1869), Elements of Meteorology: With Questions for Examination, Designed for Schools and Academies, Sheldon & Co., New York. Carslaw, H. S., and J. C. Jaeger (1959), Conduction of Heat in Solids, 2nd ed., Oxford Univ. Press, New York. Carter, E. A., M. D. Hargreaves, T. P. Kee, M. A. Pasek, and H. G. M. Edwards (2010a), A Raman spectroscopic study of a fulgurite, Philos. Trans. R. Soc. London, Ser. A, 368(1922), 3087–3097. Carter, E. A., M. A. Pasek, T. Smith, T. P. Kee, P. Hines, and H. G. M. Edwards (2010b), Rapid Raman mapping of a fulgurite, Anal. Bioanal.Chem., 397(7), 2647–2658. Collins, G. S., H. J. Melosh, and M. A. Pasek (2012), Can lightning strikes produce shocked quartz?, in Proceedings of the 43rd Lunar and Planetary Science Conference, Lunar and Planetary Inst. Technical Report, vol. 43, p. 1160, LPI Contribution No. 1659, The Woodlands, Tex. Elmi, C., J. Chen, D. L. Goldsby, and R. Gieré (2017), Mineralogical and compositional features of rock fulgurites: A record of lightning effects on granite, Am. Mineral., doi:10.2138/am-2017-5971. Ende, M., S. Schorr, G. Kloess, A. Franz, and M. Tovar (2012), Shocked quartz in Sahara fulgurite, Eur. J. Mineral., 24(3), 499–507. Eppelbaum, L, I Kutasov, and A. Pilchin (2014), Applied Geothermics, Lecture Notes in Earth Syst. Sci., Springer, Berlin. Essene, E. J., and D. C. Fisher (1986), Lightning strike fusion: Extreme reduction and metal-silicate liquid immiscibility, Science, 234(4773), 189–193. Gailliot, M. P. (1980), "Petrified lightning" a discussion of sand fulgurites, Rocks Miner., 55(1), 13–17. Gieré, R., W. Wimmenauer, H. Müller-Sigmund, R. Wirth, G. R. Lumpkin, and K. L. Smith (2015), Lightning-induced shock lamellae in quartz, Am. Mineral., 100(7), 1645–1648. Goltrant, O., H. Leroux, J.-C. Doukhan, and P. Cordier (1992), Formation mechanisms of planar deformation features in naturally shocked quartz, Phys. Earth Planet. Inter., 74(3–4), 219–240. Grapes, R. H., and H. Müller-Sigmund (2010), Lightning-strike fusion of gabbro and formation of magnetite-bearing fulgurite, Cornone di Blumone, Adamello, Western Alps, Italy, Mineral. Petrol., 99(1), 67–74. Grieve, R. A. F., F. Langenhorst, and D. Stöffler (1996), Shock metamorphism of quartz in nature and experiment: II. Significance in geoscience, Meteorit. Planet. Sci., 31(1), 6–35. Jones, D. L., G. G. Goyer, and M. N. Plooster (1968), Shock wave from a lightning discharge, J. Geophys. Res., 73(10), 3121–3127. Lan, Y., Y. Jun Zhang, W. Sheng Dong, W. tao Lu, H. Yi Liu, and D. Zheng (2011), Broadband analysis of chaotic pulse trains generated by negative cloud-to-ground lightning discharge, J. Geophys. Res., 116, D17109, doi:10.1029/2010JD015159. Lang, T. J., et al. (2016), WMO world record lightning extremes: Longest reported flash distance and longest reported flash duration, Bull. Am. Meteorol. Soc., 98(6), 1153–1168. Langenhorst, F., and A. Deutsch (1994), Shock experiments on pre-heated α- and β-quartz: I. optical and density data, Earth Planet. Sci. Lett., 125(1–4), 407–420. Lin, S.-C. (1954), Cylindrical shock waves produced by instantaneous energy release, J. Appl. Phys., 25(1), 54–57. Marsh, K. N. (1987), Recommended Reference Materials for the Realization of Physicochemical Properties, Blackwell Sci. Publ., Oxford, U. K. Kumazaki, K., K. Nakamura, K. Naito, and K. Horii (1993), Production of artificial fulgurite by utilizing rocket triggered lightning, in Proc. 8th Int. Symp. on High Voltage Engineering, pp. 269–72, Yokohama, Japan. Navarro-González, R., S. A. Mahan, A. K. Singhvi, R. Navarro-Aceves, J.-L. Rajot, C. P. McKay, P. Coll, and F. Raulin (2007), Paleoecology reconstruction from trapped gases in a fulgurite from the late Pleistocene of the Libyan desert, Geology, 35(2), 171–174. Olhoeft, G. R. (1981), Electrical properties of granite with implications for the lower crust, J. Geophys. Res., 86(B2), 931–936. Olona, J., J. A. Pulgar, G. F. Fernández-Viejo, C. L. López-Fernández, and J. M. G. González-Cortina (2010), Weathering variations in a granitic massif and related geotechnical properties through seismic and electrical resistivity methods, Near Surf. Geophys., 8, 585–599. Acknowledgments We would like to express our gratitude to H.J. Melosh and one anonymous reviewer, whose thoughtful comments and criticisms greatly improved our paper, and to the Editor Andrew Newman for highlighting our paper with a Commentary by Melosh. We especially want to thank H.J. Melosh for his engaging Commentary on our paper, which amplified our conclusions and brought our work to the attention of a broader audience. Adopted parameter values in this paper are from references cited herein, and the calculated data used in this work are the output of the theoretical model described in the paper. CHEN ET AL. LIGHTNING SHOCK AND FULGURITE FORMATION 11 Geophysical Research Letters 10.1002/2017GL073843 Pasek, M. A., and M. Hurst (2016), A fossilized energy distribution of lightning, Sci. Rep., 6, 30586. Pasek, M. A., K. Block, and V. Pasek (2012), Fulgurite morphology: A classification scheme and clues to formation, Contrib. Mineral. Petrol., 164(3), 477–492. Plooster, M. N. (1971), Numerical model of the return stroke of the lightning discharge, Phys. Fluids, 14(10), 2124–2133. Pye, K. (1982), SEM observations on some sand fulgurites from northern Australia, J. Sediment. Res., 52(3), 991–998. Rakov, V. A., and M. A. Uman (2007), Lightning: Physics and Effects, Cambridge Univ. Press, Cambridge, U. K. Roberts, J. J., and J. A. Tyburczy (1999), Partial-melt electrical conductivity: Influence of melt composition, J. Geophys. Res., 104(B4), 7055–7065. Robertson, E. C. (1988), Thermal Properties of Rocks, U.S. Dep. of the Inter., Geol. Surv., Reston, Va. Rogers, A. F. (1946), Sand fulgurites with enclosed lechatelierite from Riverside county, California, J. Geol., 54(2), 117–122. Taylor, G. (1950a), The formation of a blast wave by a very intense explosion. I. Theoretical discussion, Proc. R. Soc. London, Ser. A, 201(1065), 159–174. Taylor, G. (1950b), The formation of a blast wave by a very intense explosion. II. The atomic explosion of 1945, Proc. R. Soc. London, Ser. A, 201(1065), 175–186. Uman, M. A. (1964), The peak temperature of lightning, J. Atmos. Terr. Phys., 26(1), 123–128. Uman, M. A., and E. P. Krider (1989), Natural and artificially initiated lightning, Science, 246(4929), 457–464. White, R. H., and M. A. Dietenberger (2001), Wood products: Thermal degradation and fire, in The Encyclopedia of Materials: Science and Technology, edited by K. H. J. Buchow et al., pp. 9712–9716, Elsevier, Amsterdam. Williams, D. J., and W. Johnson (1980), A note on the formation of fulgurites, Geol. Mag., 117(3), 293–296. CHEN ET AL. LIGHTNING SHOCK AND FULGURITE FORMATION 12


text:
IT Services›
Market size outsourcing IT healthcare market worldwide 2013-2018
Projected global market size IT healthcare outsourcing from 2013 to 2018 (in billion U.S. dollars)*
by Statista Research Department, last edited Aug 31, 2013
The statistic depicts the forecasted global market size of the IT healthcare outsourcing market from 2013 to 2018. For 2017, it is projected that the global IT healthcare outsourcing market will grow to around 47 billion U.S. dollars.
Market size in billion U.S. dollars
* Figures for 2014 to 2017 were calculated using figures for 2013 and 2018 and CAGR (7.6%) provided by the source.
Emergency department attendances in Irish hospitals 2006-2017
Deaths due to pneumonia in Ireland 2010-2016
Deaths due to stroke in Ireland 2006-2017
Deaths due to cancer in Ireland 2006-2017
Statistics on "State of health in Ireland"
Health Status
Medical Consultations
Procedures conducted
Share of individuals with good or very good self reported health status in Ireland in 2018, by age groupPrevalence of individuals with good/very good health in Ireland 2018, by age
Prevalence of leading health conditions in the past 12 months in Ireland in 2017/18Leading health conditions in the last 12 months in Ireland 2017/18
Share of individuals who drank alcohol in the preceding year in Ireland in 2017/18, by ageAlcohol consumption in past 12 months in Ireland 2017/18
Share of individuals who smoke in Ireland in 2017/18, by agePrevalence of individuals who smoke in Ireland 2017/18, by age
Share of individuals who binge drink on a typical drinking occasion in Ireland in 2017/18, by gender and ageIndividuals binge drinking in Ireland 2017/18, by gender and age
Average number of consultations with a general practitioner in the past year in Ireland in 2017/18, by ageAverage number of consultations with a GP in Ireland 2017/18, by age
Share of individuals with at least one day patient hospital admission in the last 12 months in Ireland in 2017, by age groupShare of individuals who with day patient admission Ireland 2017, by age
Share of individuals with at least one dentist visit in the last 12 months in Ireland in 2018, by gender and ageIndividuals who visited dentist in last 12 moths in Ireland 2018, by gender and age
Number of emergency department attendances in Ireland from 2006 to 2017Emergency department attendances in Irish hospitals 2006-2017
Number of individuals receiving dental treatment in Ireland from 2006 to 2017Number of individuals receiving dental treatment in Ireland 2006-2017
Number of individuals receiving ophthalmic treatment in Ireland from 2006 to 2017Number of individuals receiving ophthalmic treatment in Ireland 2006-2017
Average length of hospital stay for acute care in Republic of Ireland from 2013 to 2016* (in days)Average length of stay in hospitals in Ireland 2013-2016
Number of hip replacement procedures conducted in Ireland from 2007 to 2016Hip replacement procedures conducted in Ireland 2007-2016
Number of knee replacement procedures conducted in Ireland from 2007 to 2016Knee replacement procedures conducted in Ireland 2007-2016
Number of cataract surgeries conducted in Ireland from 2006 to 2016Cataract surgeries conducted in Ireland 2006-2016
Number of hysterectomies conducted in Ireland from 2006 to 2016Hysterectomies conducted in Ireland 2006-2016
Number of caesarean sections conducted in Ireland from 2006 to 2016Caesarean sections conducted in Ireland 2006-2016
Number of mastectomies conducted in Ireland from 2006 to 2016Mastectomy procedures conducted in Ireland 2006-2016
Consumption of antibacterials for systemic use in Ireland in 2017 (in DDD per 1,000 inhabitants per day)Antibacterials for systemic use: consumption in Ireland 2017
Number of deaths due to ischaemic heart disease in Ireland from 2006 to 2017Deaths due to ischaemic heart disease in Ireland 2006-2017
Number of deaths due to stroke in Ireland from 2006 to 2017Deaths due to stroke in Ireland 2006-2017
Number of deaths due to cancer in Ireland from 2006 to 2017*Deaths due to cancer in Ireland 2006-2017
Number of deaths due to cancer of the trachea, bronchus and lung in Ireland from 2006 to 2017Deaths due to cancer of the trachea, bronchus and lung in Ireland 2006-2017
Number of deaths due to cancer of the female breast in Ireland from 2006 to 2017*Deaths due to cancer of the female breast in Ireland 2006-2017
Number of deaths due to diseases of the respiratory system in Ireland between 2008 and 2017Deaths due to diseases of the respiratory system in Ireland 2008-2017
Number of deaths due to chronic lower respiratory disease in Ireland between 2008 and 2017Deaths due to chronic lower respiratory disease in Ireland 2008-2017
Number of deaths due to pneumonia in Ireland between 2010 and 2016*Deaths due to pneumonia in Ireland 2010-2016
Number of individuals with a physical disability in Ireland from 2006 to 2017Individuals registered with a physical disability in Ireland 2006-2017
Number of individuals with multiple disabilities in Ireland from 2006 to 2017Individuals registered with multiple disabilities in Ireland 2006-2017
Number of individuals with hearing loss/deafness in Ireland from 2006 to 2017Individuals registered with hearing loss in Ireland 2006-2017
Number of individuals registered with a visual disability in Ireland from 2006 to 2017Individuals registered with a visual disability in Ireland 2006-2017
Number of individuals registered with a primary speech and language impairment in Ireland from 2006 to 2017Individuals registered with a speech and language impairment in Ireland 2006-2017
Number of individuals registered with the physical and sensory disability database in Ireland from 2006 to 2017Individuals registered with a physical or sensory disability in Ireland 2006-2017
Number of individuals attending intellectual disability services in Ireland from 2006 to 2017Individuals using intellectual disability services in Ireland 2006-2017
Number of healthcare e-commerce platforms' M&A deals China 2012-2016
Number of M&A deals in healthcare industry by type in China 2016
Out-of-pocket health spending U.S. 2014, by number of conditions and insurance
Density of neurologists by region in France 2016
Number of domestic specialized hospitals' M&A deals China 2012-2016
Value of foreign specialized hospitals' M&A deals China 2013-2015
Italy: ICT spending in the healthcare system 2014
Share of FDI in the healthcare sector in India 2011-2015
Target groups of health care marketing in the U.S. 2015
Adults receiving social care services in England 2005-2013
Who should influence government funding most for health care research 2015
Belief that the U.S. has the best health care system worldwide in U.S. adults 2019
Innovative healthcare attitudes in the United Kingdom (UK) 2015, by age
Practices electronically ordering laboratory tests in the United Kingdom (UK) 2015
Achievement of cost-containment goals in various health staffing areas in U.S. 2016
M&A: average healthcare premiums Worldwide 2010-2018, by region
Doctors in Germany
National Health Service (NHS) UK
Analytics market in India
2018 State of Cyber Resilience - Healthcare Industry
Information and communication technology (ICT) in Canada
Survey of U.S. Healthcare IT Industry Landscape 2016
Taking the Physician's Pulse - Tackling Cyber Threats in Healthcare
Accelerating Artificial Intelligence in health and care: results from a state of the natio...
Future of Patient Data
Outsourcing is Dead, Long Live Outsourcing!
IT Industry Outlook 2019
IT services and digitalization of business in Finland
Share of individuals with good or very good self reported health status in Ireland in 2018, by age group
Prevalence of leading health conditions in the past 12 months in Ireland in 2017/18
Share of individuals who drank alcohol in the preceding year in Ireland in 2017/18, by age
Share of individuals who smoke in Ireland in 2017/18, by age
Share of individuals who binge drink on a typical drinking occasion in Ireland in 2017/18, by gender and age
Average number of consultations with a general practitioner in the past year in Ireland in 2017/18, by age
Share of individuals with at least one day patient hospital admission in the last 12 months in Ireland in 2017, by age group
Share of individuals with at least one dentist visit in the last 12 months in Ireland in 2018, by gender and age
Number of emergency department attendances in Ireland from 2006 to 2017
Number of individuals receiving dental treatment in Ireland from 2006 to 2017
Number of individuals receiving ophthalmic treatment in Ireland from 2006 to 2017
Average length of hospital stay for acute care in Republic of Ireland from 2013 to 2016* (in days)
Number of hip replacement procedures conducted in Ireland from 2007 to 2016
Number of knee replacement procedures conducted in Ireland from 2007 to 2016
Number of cataract surgeries conducted in Ireland from 2006 to 2016
Number of hysterectomies conducted in Ireland from 2006 to 2016
Number of caesarean sections conducted in Ireland from 2006 to 2016
Number of mastectomies conducted in Ireland from 2006 to 2016
Consumption of antibacterials for systemic use in Ireland in 2017 (in DDD per 1,000 inhabitants per day)
Number of deaths due to ischaemic heart disease in Ireland from 2006 to 2017
Number of deaths due to stroke in Ireland from 2006 to 2017
Number of deaths due to cancer in Ireland from 2006 to 2017*
Number of deaths due to cancer of the trachea, bronchus and lung in Ireland from 2006 to 2017
Number of deaths due to cancer of the female breast in Ireland from 2006 to 2017*
Number of deaths due to diseases of the respiratory system in Ireland between 2008 and 2017
Number of deaths due to chronic lower respiratory disease in Ireland between 2008 and 2017
Number of deaths due to pneumonia in Ireland between 2010 and 2016*
Number of individuals with a physical disability in Ireland from 2006 to 2017
Number of individuals with multiple disabilities in Ireland from 2006 to 2017
Number of individuals with hearing loss/deafness in Ireland from 2006 to 2017
Number of individuals registered with a visual disability in Ireland from 2006 to 2017
Number of individuals registered with a primary speech and language impairment in Ireland from 2006 to 2017
Number of individuals registered with the physical and sensory disability database in Ireland from 2006 to 2017
Number of individuals attending intellectual disability services in Ireland from 2006 to 2017
Number of M&A deals made by healthcare e-commerce platforms in China from 2012 to 2016
Number of M&A deals in healthcare industry in China in 2016, by type
Out-of-pocket health care spending among U.S. adults in 2014, by number of chronic conditions and insurance (in U.S. dollars)
Number of neurologists per 100,000 population in France in 2016, by region
Number of M&A deals made by domestic specialized hospitals in China from 2012 to 2016
Value of M&A deals made by foreign specialized hospitals in China from 2013 to 2015 (in million yuan)
Composition of ICT spending in the Italian healthcare system in 2014* (in million euros)
Healthcare as a share of Foreign Direct Investment across India from 2011 to 2015
Most important target groups of health care marketing in the United States in 2015
Number of adults receiving local authority social care services in England from 2005/2006 to 2012/2013 (in millions)
Which one of the following do you feel should have the most influence on how government funds for research to improve health are spent?
Do you believe the United States has the best health care system in the world?
Attitudes towards innovative healthcare in the United Kingdom (UK) in 2015, by age
Share of practices routinely or occasionally using electronic ordering of laboratory tests in the UK 2015
Staffing areas in U.S. healthcare organizations with the greatest achievements of cost-containment goals as of 2016
Average healthcare merger and acquisition (M&A) premiums to four week stock price Worldwide from 2010 to 2018, by region

redpajama_set_name: RedPajamaCommonCrawl


This darling Sway mini skirt is our new go to! We just adore it's fluted frill and trim detail. This mini has a zip at the back and is true to size. Looks super sweet teamed with our NLB tops, vintage or lace. For a casual weekend look pair it with our Lola crop! Made in a linen viscose blend, which is soft but also textured and slightly more weighted.
Model wears size 6 and is 175cm tall. These run slightly bigger so we recommended sizing down if you are after a more high waisted, firm fit.


Today I finished this little vest for my niece's birthday present, although this post will be delayed so it stays a surprise until after the party. The fabric was so lovely I couldn't resist, even though it was more than I'd usually pay for fabric for a child's garment. She's only little, so it didn't take much, and she's worth it!
It has a really lovely, bumpy texture, and hangs nicely. I combined it with pink denim to ensure it had a decent weight and hopefully washes well. Let's face it, toddler clothes are going to need washing!
I pre-washed and tumbledried my fabrics, as usual. I do this with all my apparel fabrics before sewing, or actually with anything I think might need washing in the future. This way, it's less likely to shrink by a large amount when washed after wear, and the seams stay nice. It drives me crazy waiting for the fabric to be ready for sewing, but it's worth it for the end result.
I know I hate it when I buy something at the shops, and then it shrinks weirdly in the wash and never sits the same way again.
I've made 2 of these previously, and I really recommend this pattern.
Happy Birthday, munchkin. Love you.