21 March 2014 Porcine skin damage thresholds for 0.6 to 9.5 cm beam diameters from 1070-nm continuous-wave infrared laser radiation
Author Affiliations +
There is an increasing use of high-power fiber lasers in manufacturing and telecommunications industries operating in the infrared spectrum between 1000 and 2000 nm, which are advertised to provide as much as 10 kW continuous output power at 1070 nm. Safety standards have traditionally been based on experimental and modeling investigations with scant data available for these wavelengths. A series of studies using 1070-nm infrared lasers to determine the minimum visible lesion damage thresholds in skin using the Yucatan miniature pig (Sus scrofa domestica) for a range of beam diameters (0.6, 1.1, 1.9, 2.4, 4.7, and 9.5 cm) and a range of exposure durations (10 ms to 10 s) is presented. Experimental peak temperatures associated with each damage threshold were measured using thermal imaging. Peak temperatures at damage threshold for the 10-s exposures were ∼10°C lower than those at shorter exposures. The lowest and highest experimental minimum visible lesion damage thresholds were found to have peak radiant exposures of 19 and 432  J/cm 2 for the beam diameter-exposure duration pairs of 2.4 cm, 25 ms and 0.6 cm, 10 s, respectively. Thresholds for beam diameters <2.5  cm had a weak to no effect on threshold radiant exposure levels for exposure times ≤0.25  s , but may have a larger effect on thresholds for exposures ≥10  s .



The number of lasers in the operational range of 1000 to 2000 nm continues to increase in medical, industrial, communications, and military applications, while relatively few safety studies for skin damage from continuous-wave (CW) exposures at these wavelengths can be found in the literature. High-power fiber lasers utilize advanced technologies to combine active optical fibers with semiconductor diodes. These active fibers allow for an extremely high-intensity light out of a very small core, making it possible to produce multikilowatt lasers. For these lasers, powers >10kW are routinely used in manufacturing processes and the telecommunications industry. CW outputs can be obtained with fiber lengths >100m, with small beam divergence and high power conversion efficiencies. Experimental data used in setting the ANSI standards have not kept pace with the rapid advancement of these lasers. Furthermore, the maximum permissible exposure (MPE) recommended through the ANSI Z136.1-20071 is for a limiting, or measurement, aperture diameter of 0.35 cm in this wavelength region for exposures to skin with little empirical evidence that the standard recommendations are adequate for larger beam diameters, especially for highly penetrating wavelengths such as 1070 nm. Irradiances can easily exceed 10kW/cm2 and little or no minimum damage threshold data are available for these high levels. Therefore, this study was initiated to answer critical questions regarding the appropriateness of the standards for MPE in terms of relative penetration depth (wavelength) and for CW lasers (10 ms to 10 s) with beam diameters greater than the ANSI limiting aperture of 0.35 cm in the near-infrared (IR) region at 1070 nm.

To determine if an MPE is sufficiently safe, a series of experiments are conducted to evaluate the damage threshold for a minimum visible lesion (MVL). An MVL is considered to be any subtle, slight, morphological, or coloration, change observed at the exposure site post laser radiation at various time points. The MPE is typically several times less than the MVL, to ensure that exposures at or below the MPE would have no damaging effects.

The Yucatan miniature pig has been established as a model for human skin damage studies due to the morphological and physiological similarities to human skin.2 Skin threshold studies were conducted with Yucatan miniature pigs at wavelengths of 1070 and 1940 nm with varying spot sizes and exposure times. Results for the 1940-nm measurements have been reported3 while only portions of the 1070-nm data presented here have appeared previously.4 Also, MVL skin thresholds have been reported with these miniature pigs for various spot sizes and exposure times for wavelengths at 1318 and 1540 nm.56.7

A summary of threshold data obtained from the literature for skin exposure at 1060 nm are listed in Table 1. No published data were found for skin exposure at 1070 nm, which was used in this study, but optical properties of the skin are remarkably similar between 1060 and 1070 nm, justifying direct comparison.8,9 The time of inspection used to assess damage endpoints is assumed to be 24 h post exposure; however, the exact time of evaluation for many of the studies was left ambiguous in the text and could range from immediate response to 48 h post exposure.

Table 1

Literature summary of skin damage threshold data at the 1060-nm laser wavelength.

ReferenceSubjectExposure time (s)1/e2 spot size (cm)Damage threshold (J/cm2)
10White haired pig1.00.559.4
10Human—light complexion Chinese1.00.565.5
10Human—yellow complexion Chinese1.00.561.0
10Human—dark complexion Chinese1.00.552.3
14Guinea pig—dark pigmented10×10−90.251.0


Damage threshold calculated as the average for the radiant exposure of the lowest positive and highest negative responses 24 h post laser exposure.

There are several chromophores in the skin related to the absorption, μa, of laser radiation at near-IR wavelengths. These chromophores include melanin, water, fat, hair, and blood. The values for μa have been reported by numerous sources and are summarized in Fig. 1.

Fig. 1

Comparison of the absorption coefficient (μa) for various chromophores present in skin.8,9,1617.18.19.,20


The experiments reported in this study were designed to supplement and extend the body of knowledge summarized in Table 1 in the following ways: (1) To determine the MVL damage thresholds with beam spot size and exposure duration dependency, (2) to examine the thermal response on the surface of the skin over a variety of exposure conditions, and (3) to better understand the impact of the relative absorption coefficient (μa) of the skin on damage threshold for the deeply penetrating 1070-nm wavelength.




Animal Use

The animals involved in this study were procured, maintained, and used in accordance with the Federal Animal Welfare Act, “Guide for the Care and Use of Laboratory Animals,” prepared by the Institute of Laboratory Animal Resources National Research Council, and DoD Regulation 40-33 Secnavinst 3900.38C AFMAN 40-401(1) DARPAINST 18 USUHSINST 3203 “The Care and Use of Laboratory Animals in DOD Programs.” The Air Force Research Laboratory at Joint-Base Sam Houston, Texas, has been fully accredited by the Association for Assessment and Accreditation of Laboratory Animal Care, International since 1967.

A total of 22 female Yucatan miniature pigs were used for this study. Subjects had food withheld for 12 h prior to the procedure. Miniature pigs (30 to 50 kg) were sedated with an intramuscular (IM) injection of Tiletamine/Zolazepam (4 to 6mg/kg) IM and Xylazine (2.2mg/kg) IM. Subjects had an intravenous ear catheter (23 to 18 gauge) placed for the administration of lactated ringers. Animals were then orotracheally intubated and anesthesia was maintained with isoflurane (1 to 3%) throughout the procedure. Heart rate, oxygen saturation, core body temperature, and respiration were continuously monitored. After animals reached a surgical plane of anesthesia, electric clippers were used to remove hair from the flank. The flank was cleaned with soap and warm water, rinsed, and dried with a lint-free disposable cloth. The core body temperature of subjects was maintained using a T/Pump™ circulating water blanket (Gaymar Inc., Orchard Park, NY).

A grid pattern was drawn on the flank of the pigs using a black permanent marker with grid spacing of 1.5, 2, 3, 4, 8, or 12 cm for the nominal laser beam diameters of 0.6, 1, 2, 2.5, 5, and 10 cm, respectively. Skin biopsy punches (2 to 8 mm) were collected from select exposure sites at 1 and/or 24 h post exposure and stitched closed with nonabsorbable suture. Subjects were sedated again for 24 h post laser assessment and biopsies. Buprenorphine (0.05 to 0.1mg/kg) was administered IM as a postprocedural analgesic.

For laser exposures, anesthetized pigs were placed in a customized pig sling (BH Inc., Wheatland, WY) or on a hydraulic table in ventral recumbence. In this manner, each area of the grid, which had been drawn on the flank of the pig, was sequentially addressed with a laser exposure, according to a predetermined randomized selection of energy levels distributed over a range encompassing the anticipated threshold. Energy readings for each laser exposure and any remarkable observations were recorded on data sheets following each exposure.


Experimental Setup and Data Collection

Two different laser systems from IPG Photonics (Oxford, MA) were employed as laser sources for exposure. For the nominal beam diameters of 0.6, 1, and 2 cm, an IPG Photonics YLR-3000 laser with a maximum CW power of 3 kW was used. The experimental configuration for this system has been previously reported.3 Output of this laser from the fiber optic was approximately Gaussian in spatial energy distribution with a 1/e2 diameter of 1.3cm. The size of the beam at target was controlled with the use of a beam expanding/reducing telescope. Power at target was adjusted via electronic control within the laser head via preset PC programming. For each laser exposure, reference energy from a low-power sample of the beam was recorded with a Molectron PM10 power meter (Coherent Inc., Santa Clara, CA). Measurements acquired with the PM10 in the reference beam were compared with measurements obtained with a PM150 placed near the target plane before testing, at the scheduled midway time point during testing and after each sequence of exposures to assure integrity of the high-energy portion of the beam train and to confirm accurate quantification of energy delivered to target. The laser was co-aligned with a helium neon laser (Meredith Instruments, Glendale, AZ) to facilitate beam placement at the target. Beam size at target was measured with a Pyrocam III (Spiricon, Logan, UT) at the time of setup for the 2-cm beam configuration, while burn paper was used to measure the smaller 0.6- and 1-cm-diameter beams. Output from the telescopes was collimated to maintain beam diameter over the expected range of target distances determined from the anatomical confirmation of the animals. Consistency of the beam size was verified during testing with burn paper (ZAP-IT™, Kentek Corporation, Pittsfield, NH). A Sony (San Diego, CA) CCD-Iris color video camera was also utilized in this arrangement for target position verification. This video camera was mounted on a custom translation stage and was co-focused at a common region of interest of the thermal camera.

For the nominal beam diameters of 2.5, 5, and 10 cm, an IPG Photonics YLR-50000 laser was used. This laser system was limited to a maximum CW power of 30 kW for this study. A telescope was used to adjust the beam’s diameter at the target located 305cm from the fiber head. An optical wedge placed 122cm in front of the fiber head was used to reflect a portion of the beam to a Spectralon screen 183cm from the wedge. A camera (Alpha Near Infrared InGaS ICI NIR) was aimed and focused at the Spectralon screen to capture the beam’s profile using Spiricon software. This allowed for visualization of the beam’s profile equidistant to the target without damaging equipment. A Coherent PM150x read by a LabMax power meter was placed near the Spectralon screen to detect a small portion of the beam to capture the reference arm power with each exposure. To calibrate the energy delivered to the target, a series of energy measurements were taken at the start and end of each day using a Black Ball 3 DELE-S-GP-3 and Coherent PM5K with an EPM2000 power meter. A FLIR SC660 longwave IR camera was used in this setup to monitor the placement of the subject in addition to providing initial skin surface temperatures.

Pulse widths were verified once at the start of the setup for all of the experiments for both the utilized laser systems. This was done with a silicone photodiode and a TDS3054B or a TDS220 digital oscilloscope (Tektronix Inc., Beaverton, OR). Once it was determined, the laser consistently gave the programmed pulse width; it was not measured again. The oscilloscope showed rise and fall times for the IPG lasers were <100μs, providing better than ±0.5ms accuracy on the exposure times. Repeated shots were measured with the Coherent PM series thermopile detectors (PM5K, PM10, or PM150 depending on beam size and power configurations) and one of the following Coherent meters operating in energy mode: EPM1000, EPM2000, 3Sigma, or a LabMax.

All energy meters and detectors were verified to have been calibrated by the manufacturer within the year prior to the experiments. Uncertainty on the energy measurements was determined from the specifications from the manufacturer (Coherent) for each of the power meters and power detectors used in the calibration of the beam configuration. This was determined to be ±5.5% for beam diameters <9cm. Since the largest beam diameter exceeded the detector sizes available, the Black Ball was used to calibrate the energy at the target plane, yielding an uncertainty in delivered energy of ±7.2% for the largest beam diameter condition in the study. The uncertainty associated with the beam diameter was determined by taking the average of the measured standard deviations from multiple beam profiles analyzed from each beam configuration finding a ±2.8% (min and max were 1.1 and 5.2% found for the 9.5 and 0.6 cm beam diameter configurations, respectively).

To capture the temperature rise during and after laser exposure, both laser setups used an SC4000 (IR) high frame rate thermal camera (FLIR Systems, Boston, MA), which is sensitive to radiation in the 3- to 5-μm thermal emission band. This camera sensitivity was specifically selected to be out-of-band for the 1070-nm wavelength under test so that surface temperature changes during laser exposure could be obtained without backscatter from the target distorting the measurements. The thermal camera was triggered for data collection to begin at least 100 ms prior to laser exposure of the pig skin. IR camera frame rates were 400 fps for experiments where beam diameters were <2cm and 200 fps for experiments where the beam diameters were >2cm. The camera collected data for a period of time between 2 and 30 s beyond the completion of each laser exposure to capture the trajectory of thermal decay of the laser-heated tissue. Course focus of the thermal camera was verified with visible feedback from a monitor connected to a video camera.


Statistical Analysis

The Probit procedure21 was used to estimate the ED50 dose (i.e., damage threshold for a lesion in 50% of the observations) for creating an MVL for all laser configurations used in the skin studies at both 1- and 24-h time points. Observers were instructed that a lesion was constituted as any change in the appearance of the surface of the skin to include any erythema, no matter how slight, blistering, or change in pigmentation as viewed by eye. An exposed area was deemed a lesion if at least two out of three trained observers, blinded from each other, agreed there was a lesion. A minimum of three animals were used for each laser configuration investigated. In addition, 95% confidence intervals for dose response were calculated from the data. An adequate number of data points were collected to ensure that the magnitude of upper and lower fiducial limits at the ED50 level varied no more than 50%. In order to reduce the number of subjects needed, this requirement for the fiducial limits being within 50% of the ED50 was disregarded for the beam diameters >2cm. Furthermore, a Probit-curve slope greater than two, with aggregated 24-h observation data, was set as a minimal experimental endpoint for the MVL.



A total of 1306 skin exposures were completed on both flanks of 22 Yucatan miniature pigs with the 1070-nm laser. Sixteen different configurations of laser parameters were used as listed in Tables 2 and 3 along with correlating ED50’s, Probit slopes, average irradiance, and peak radiant exposure data. It should be noted that the beam diameter is listed as 4σ and not the usual 1/e or 2σ because the 1/e2 diameter must be used to evaluate the hazards of laser beams on an average irradiance basis. The radiant exposure is calculated using peak power, 1/e diameter, not the average which uses 1/e2 diameter.

Table 2

Probit data results from observations for lesions 1 h after exposure from the 1070-nm laser.

1-h Probit
Number of exposures (Yes, No)1/e2 diameter (cm)Time (s)ED50 (95% Conf. Inter.) (J)Probit slopeIrradiance (W/cm2)Peak radiant exposure (J/cm2)
21, 1010.60.016.6 (5.7 to 8.3)5.9249850
27, 270.60.0712.7 (11.4 to 14.1)18.168796
32, 1120.61058.4 (53.6 to 67.1)7.622442
45, 1461.10.0117.3 (15.7 to 19.9)6.0199840
74, 671.10.122.7 (21.2 to 24.4)10.726252
73, 711.11077.4 (71.9 to 82.7)10.39179
37, 171.90.0558.2 (51.2 to 63.8)12.140240
47, 831.90.2577.9 (74.9 to 81.9)15.510854
46, 471.910162.0 (152.0 to 171.0)18.06112
19, 222.40.0153.6 (29.6 to 92.3)3.4116523
23, 132.40.02540.5 (29.8 to 54.3)6.535218
27, 122.40.2554.8 (40.5 to 68.5)9.54824
11, 184.70.025315.0 (241.0 to 660.0)8.472336
19, 114.70.1283.0 (233.0 to 318.0)18.016232
33, 74.70.25325.0 (229.0 to 372.5)12.07537
12, 69.50.251142.5 (227.0 to 1367.5)11.46533

Table 3

Probit data results from observations for lesions 24 h after exposure from the 1070-nm laser. Note the ratio of the ED50 to maximum permissible exposure (MPE) is taken as Ef/MPE, where Ef is the effective irradiance computed from Eq. (1).

24-h Probit
Number of exposures (Yes, No)1/e2 diameter (cm)Time (s)ED50 (95% Conf. Inter.) (J)Probit slopeIrradiance (W/cm2)Peak radiant exposure (J/cm2)Ratio of ED50 to MPE
52, 710.60.014.0 (3.8 to 4.3)13.315183013
31, 220.60.0711.1 (9.7 to 12.0)11.76008422
38, 1050.61057.1 (52.0 to 65.9)6.62243232
46, 1431.10.0116.1 (15.1 to 17.4)9.518593721
71, 691.10.122.6 (20.5 to 24.6)6.72615217
67, 781.11082.0 (76.0 to 87.9)9.3918919
43, 91.90.0548.6 (42.5 to 53.9)16.33363413
56, 731.90.2574.2 (71.8 to 77.1)18.91035113
44, 471.910154.0 (142.0 to 166.0)10.8510611
24, 172.40.0146.6 (38.4 to 53.7)15.710132012
22, 142.40.02543.0 (32.3 to 58.0)7.0374199
27, 122.40.2555.3 (42.3 to 69.3)11.548246
21, 94.70.025213.3 (146.3 to 260.0)7.64902411
10, 304.70.1267.0 (214.0 to 305.0)13.41533110
35, 54.70.25295.0 (185.0 to 345.0)11.468349
13, 59.50.251125.0 (500.0 to 1310.0)24.664328

Lesions created from exposure durations <100ms were typically absent at 1 h and present at 24 h post exposure. Exposure durations 100ms were typically observed at the 1-h read and persisted at 24 h. This general observation is supported by comparing the number of yes/no lesions (column 1 in Tables 2 and 3) at 1 and 24 h post exposure and the ratio of the 1-h ED50 to the 24-h ED50.

To calculate the ratio of ED50MPE reported in Table 3, the ANSI Z.136.1-2007 standard for the Safe Use of Lasers recommends that for a skin exposure, the limiting aperture should be referenced to a standard 0.35 cm. Thus to compare the experimentally determined ED50’s to the MPE, the irradiance values were scaled to an effective irradiance, Ef, by


where Φ0 is the power in Watts, DL is the 1/e laser beam diameter in centimeter, and Df is the limiting aperture of 0.35 cm. It is worth noting that for beam diameters >1cm (1/e2 diameter), the exponential portion of Eq. (1) begins to approach zero, showing that Ef and the experimentally determined ED50 are nearly the same for these larger beams. Graphical comparison of Ef versus exposure time is contrasted to the MPE in Fig. 2. Excel (Microsoft Office 2007) was used to determine the equation to best fit the data for Ef versus t, for each beam diameter and the power fit equation Ef=atm finding a=143.1, 102.7, 65.9, 25.6, and 41.8 and m=0.63, 0.76, 0.79, 0.94, and 0.86 for the 0.6, 1.1, 1.9, 2.4, and 4.7-cm beam diameters, respectively. Note that for this wavelength, the current MPE=5.5t.75 between 0.01 and 10 s then is a constant 1W/cm2 for exposure durations 10s.1

Fig. 2

Damage threshold data converted to the effective irradiance contrasted to the maximum permissible exposure (solid line).


For exposures very near or above the ED50, a raised appearance known as a wheal (also spelled weal) or flare was commonly observed for lesions at all durations and spot sizes. The raised area on the skin appeared more subtle for exposures near the ED50 and more pronounced for exposures significantly higher than the ED50. The wheal and flare appeared to be proportional to the size of the beam, which was easily visualized by illuminating the lesion at a glancing angle to the skin. The wheal was not always easily visualized when illuminated under general overhead lighting conditions. An example of skin-lesion appearance over observational endpoints is shown in Fig. 3. Black marker grid lines drawn on the skin in Fig. 3 were 4×4cm apart. Though there is no pictorial documentation available, the exposed area appeared red immediately after the laser exposure, but had resolved completely by 1 h. No annotations were made of when the wheal and flare began to appear in any exposures, just that it was observed at 24 h post exposure. Note that the wheal was not used as an indicator for an MVL, but was a general observation reserved for discussion.

Fig. 3

Example of skin lesions observed in this study before and after laser exposure (10 ms exposure, 2.4 cm beam diameter, 86J±5.5%, 1.8 times higher than the 24-h ED50). Both the 1- and 24-h post laser exposure reads were positive for a lesion. The white arrows indicate a very slight brown discoloration on the epidermis; no erythema was observed at 1 h. The white bracket indicates the extent of the observed erythema at 24 h. The dashed, white circle approximates the footprint of the observed raised area on the skin 24 h post exposure.


IR camera data proved exceptionally difficult to process for this wavelength due to hot spots from the absorption of laser energy in the hair and steam rising from the skin due to phase changes in the water inside the skin. IR camera data for the 0.6- and 1.1-cm beam diameters with exposure durations <10s were processed, but had enormous variability during laser heating due to these superficial hot spots and were not suitable for rate-of-heating analysis.

IR camera data were processed by selecting a region of interest (ROI) (typically 9 to 20 pixel diameters) placed in the central area of exposure to examine the averaged temperature in the ROI at each time t during the laser exposure (0<texposure duration) and after the laser was turned off (t>exposure duration). Since hair burns very quickly at this wavelength, care was taken to avoid hot spots believed to be follicles, or in some cases, steam, sometimes requiring the ROI to be parsed into two smaller 3 by 3 pixel ROIs. The initial temperature, T0, of the exposed skin surface was determined from IR camera data acquired before the laser turned on (t0). The change in temperature, ΔT, was then determined by the difference between T0 and the peak temperature, Tpeak, at the end of the exposure. IR camera data were processed to give the linear relationships, with a forced zero intercept, for the peak temperature rise, ΔT, as a function of power P. The ΔT versus P data for the 1.9-cm beam diameter are shown in Fig. 4. The slopes, mT, were recorded from each fit to processed thermographic data for ΔT=mTP and are reported in Table 4.

Fig. 4

The peak temperature rise from numerous exposure powers for the 1.9 cm beam diameter as tested at 0.05, 0.25, and 10 s exposure durations. Units on the slope fit parameter, mT, are °C/W.


Table 4

The linear slopes, mT, (°C/W) from the fitted thermal camera data P(ΔT) and calculated ΔTMVL to achieve the minimum visible lesion at 24-h ED50’s.

Exposure duration, t (s)1/e2 diameter, d (cm)mT (°C/W)Goodness of linear fit to data, R2[(mTπd2)/4t]°C/(J/cm2)T0 (°C)ΔTMVL=mT(°C/W)×ED5024-h (W)±0.530.3

The average and standard deviation for T0 at each exposure condition shown in Fig. 4 are also reported in Table 4. The slopes in Table 4 are also corrected for radiant exposure by mT×(πd2/4t) [°C/(J/cm2)] for later discussion.

Time-temperature history for each of the three beam diameters in the 10-s exposures near ED50 are shown in Fig. 5(a). The periodic ripples in the 10-s exposures stem from the subject’s breathing; care was taken to keep the ROI in the central area of the beam to minimize the variations induced by the breath (i.e., hair follicles/hot spots moving periodically through the ROI).

Fig. 5

Time-temperature history for cases nearest ED50 resulting in a lesion at 24 h for (a) 10-s exposures, various beam diameters, and (b) largest beam diameter (9.5 cm) contrasted to the shortest exposure (10 ms) condition tested.


Time-temperature history for tested cases nearest ED50 for the largest beam diameter (250 ms, 9.5 cm) and shortest exposure (10 ms, 2.4 cm) are shown in Fig. 5(b). Notice that the experimentally determined mT(°C/W) in Table 4 for these two exposure duration-diameter pairs [depicted in Fig. 5(b)] are the same in °C/W, but very different in °C/(J/cm2). All time-temperature curves presented in Fig. 5 were positive for a lesion at 24 h. It should be noted that the uncertainty of the energy measurements was ±5.5% for all beam diameters <9cm and ±7.2% for the 9.5-cm beam reported in Figs. 5 and 6.

Fig. 6

The thermal image captured at the end of laser exposure of cases nearest 24-h ED50 for (a) 2.4 cm diameter, 10 ms, 46.5±2.6J (Video 1, MP4, 43 kB) [URL:  http://dx.doi.org/10.1117/1.JBO.19.3.035007.1] and (b) 9.5 cm diameter, 250 ms, 1067±77J (Video 2, MP4, 137 kB) [URL:  http://dx.doi.org/10.1117/1.JBO.19.3.035007.2]. The laser comes on at t=0. Temperatures are absolute T (°C). The 1-cm scale bar applies to both (a) and (b).


The thermal video files corresponding to the time-temperature histories depicted in Fig. 5(b) for the 2.4- and 9.5-cm cases are shown in Video 1 and Video 2. Figure 6 shows still frames of the supplemental video files for the thermal images captured nearest to the laser’s cutoff time.

Analyzed thermal data for the 2.4-cm beam at the 10-ms pulse did not necessarily capture the peak temperature (200 fps, 1.5 ms integration time); thus these data were analyzed differently than all other thermal data by using the two frames where the laser was on, selecting an ROI in frame 1 to find the temperature rise ΔT1 at t1 and then the same ROI in fame 2 to find ΔT2 at t2. A linear fit was then applied to the thermal data to identify ΔT as function of time and solving the linear fit for an assumed ΔTmax at t=10ms (having assumed ΔT=0 at t=0). Most of these 10-ms data only had two frames showing increasing temperature, but one 2.4-cm, 10-ms thermal video captured three frames showing increasing temperature, giving some experimental evidence of the peak temperatures achieved nearest the end of this laser pulse. This uniquely captured event is presented as supplemental data Video 1.mp4 (43 kB); the time stamp in the video was created with the assumption the peak temperature profile occurred at t=10ms. Discussion about when the peak temperature may occur from a 10-ms pulse is reserved for later.

Histology sections of Yucatan miniature pig skin punch biopsy samples were prepared with haematoxylin and eosin staining. Punch biopsies were collected from the central grid area to coincide with the center of the laser exposure site. Biopsies taken from unexposed areas of skin on the flank were used to measure skin thicknesses. Thicknesses of the stratum corneum (SC), epidermis (not including SC), and dermis ranged from 14 to 47, 52 to 82, and 2350 to 3860 μm with averages and standard deviations of 29±10, 59±17, and 2953±153μm, respectively. Histological evaluations of the damage lesions are reserved for a future paper.



MVLs are considered to be any subtle, slight, morphological or coloration change observed at the exposure site post laser radiation at 1- and 24-h time points. These changes include the faint brown discoloration or reddening (depicted in Fig. 3), or more obvious changes such as blistering or charring. Blistering and/or charring of the skin were not observed at the ED50 level for any exposure condition in this study. The MVL does not necessarily describe persistent, irreparable damage to the tissue, but rather provides evidence for radiant energy levels that will begin to cause noticeable changes to the tissue, with some changes being transient in nature.

Laser-induced erythema and wheal and flare appear to have a transient, and perhaps related, nature that is not well understood. It was common to observe immediate erythema, which disappeared at 1 h, but reappeared at 24 h post exposure for the 10-ms cases (Fig. 3). This is in contrast to 250-ms cases of the same beam diameter where disappearing erythema at 1 h usually remained resolved at 24 h. Foster et al.22 reported the curious ability of erythema caused by 1320/1440-nm laser radiation to resolve after 24 h, then reappear upon subtle thermal stimulation (i.e., hot shower) 48 h after exposure, and suggested this recall ability may be related to the wheal and flare response. Wheal and flare was observed at exposure levels at or above damage threshold for all laser beam diameters and exposure times in this study. The observation of the wheal and flare from exposed areas at or above damage threshold is not uncommon for laser exposures to skin in the visible or IR.2324.25.26 The wheal and flare is believed to be a neurogenic response caused by the stimulation of a neuropeptide protein called Substance P (SP).23,24 The exact mechanism of how and why the laser stimulates SP is not known, though one hypothesis proposed by Algermissen et al.26 suggests this involves the thermal destruction of the nerve endings in the dermis and/or epidermis from the laser exposure. The wheal and flare has not been used as an indicator of the MVL, though it should be included in future studies as part of the observational evidence lending insight into the tissue response at damage threshold exposure levels.

Comparing the MPE to the effective irradiances, Ef, used in this study found that the margin for safety between the Ef and MPE ranges from a ratio of 6.2 to >32.3 (last column in Table 3). The data trends in Fig. 2 indicate that the damage threshold may have a weak dependence on beam diameter much larger than the 0.35-cm limiting aperture used in the ANSI.1 For the 0.25-s exposures, the 1.9-cm beam had a statistically significant higher threshold (51.26J/cm2) than the larger beam diameters, indicating that the threshold is dependent on beam diameter, at least for d<2cm (Fig. 2). The 2.4-, 4.7-, and 9.5-cm beam diameters converged toward a peak radiant exposure threshold level of 30J/cm2 as indicated by their significantly overlapping 95% confidence intervals. This suggests that for 0.25-s exposures at 1070 nm, the peak radiant exposure threshold converges to 30J/cm2, settling to a safety margin (ratio of EfMPE) of 8, for d>2cm. The 0.1-s did not have overlapping 95% confidence intervals for peak radiant exposure thresholds of 52.2 and 30.65J/cm2 for the 1.1- and 4.7-cm beam diameters, respectively. If the evidence for convergence to 30J/cm2 in the 0.25-s data is considered, it is likely that the threshold for the 4.7-cm case represents a value for convergence between threshold and increasing beam diameters for 0.1-s exposures, giving a safety margin of 10 between threshold Ef and MPE for very large beam diameters. The 10-s exposure data, however, were 432.24, 189.4, and 106.38J/cm2 for the 0.6-, 1.1-, and 1.9-cm beam diameters, respectively, and exhibit no overlap in 95% confidence intervals, suggesting that convergence of the threshold with respect to beam diameter was not achieved. The larger spread in damage thresholds for the 10-s data also suggest that damage threshold is more dependent on beam diameter for very long exposures compared to exposures <10s. The 0.01-s exposure thresholds for the 0.6- and 1.1-cm diameters have overlapping confidence intervals, but the larger 2.4-cm beam diameter had a peak-radiant exposure damage threshold 33% less than the smaller beams, suggesting some dependency on beam diameter even for this very short exposure. These trends in beam diameter dependency for MVL threshold data are suspected to be caused by different damage mechanisms for very short exposures (<0.05s) compared to the longer exposure durations.

The nature of skin lesions generated with 1070-nm laser radiation is, in general, markedly different from those generated with longer IR wavelengths, such as 1940 nm. This difference is readily observed at shorter-duration exposures. Lesions generated from the 1070-nm wavelength for short (10-ms) exposure durations and smaller diameters (0.6 and 1.1 cm) seemed to be multifocal in nature and involved a desiccation process as evidenced by a white, flaky appearance to the lesion, whereas MVL studies of lesions from the 1940-nm wavelength do not report this characteristic.3,27 The 1070- and 1940-nm laser radiation wavelengths have significantly different damage thresholds in skin at comparable conditions. For example, Oliver et al.3 reported a 1940-nm wavelength with a 1.8-cm, 10-s exposure had a 24-h ED50 of 13.3J/cm2 on Yucatan mini-pig skin, as compared to the 1070-nm wavelength using a 1.9-cm, 10-s exposure, which had an ED50 of 106.3J/cm2, eight times that of a comparable setting at the longer 1940-nm wavelength. These differences in lesion appearance and damage thresholds at comparable conditions between the 1940- and 1070-nm laser radiations can be explained by differences in optical absorption and scattering of melanin and water in the skin.

Based upon the literature cited in Table 1, damage thresholds at 1070 nm should be highly dependent upon pigmentation. The absorption of pigment should be considered up to wavelengths of 1400nm, where absorption in bulk tissue should become substantially less significant than water absorption due to the low fractional content of pigment as compared to the surrounding water content. The consistency of pigmentation between Yucatan miniature pigs used in this aspect of the study likely helped minimize variance of data. Rockwell and Goldman15 and Baozhang et al.10 have reported thresholds for darker skin to be as much as 20 to 25% lower than lighter-skin subjects when using a 1064-nm wavelength. This is in stark contrast to studies at the 1940-nm wavelength by Chen et al.,28 which do not report significant differences in thresholds between light- and dark-pigmented skin. If the 24-h ED50 MVL is 25% lower for darker skin exposed to 1070-nm laser radiation, then the safety margin between Ef and MPE would reduce significantly to a range of 1.5 up to 8.1. Adequacy of the animal model should be considered when applying these findings to all pigmentation levels expected across the human race and sufficient safety margins should take this variation into account.

The 1070-nm laser radiation wavelength also penetrates deeper into tissue, heating a larger volume, as compared to the 1940-nm wavelength. When comparing the thermal burn created by a 1064-nm laser radiation wavelength (1 cm beam diameter, 40 W delivered for 8 s) to a copper-brass iron, Zhang et al.29 reported tissue damage through the skin into internal organs of a Wistar rat model where lesions had a coagulative burn appearance on the skin surface. Though the dorsal skin of the Wistar rat was significantly thinner and had far less subcutaneous tissue than the Yucatan mini-pig model, the observations reported in Zhang et al.29 demonstrated the deeply penetrating nature of this laser radiation wavelength region. It is this very nature of these deeply penetrating wavelengths that is believed to give rise to the MVL threshold dependency on beam diameter. This wavelength heats a larger volume of tissue, compared to less penetrating wavelengths, enabling radial conduction to dominate over axial thermal diffusion.

Vapor or smoke was observed for some of the higher-energy exposures near or above threshold, particularly for the 10-ms exposures, providing evidence for water phase transitions occurring in deeper layers of the skin. On occasion, dry skin debris would explode off the surface (popcorn effect) when exposed to short exposures (<100ms). However, these seemingly violent interactions did not always result in grossly observable damage to the skin surface. Due to better penetration at 1070 nm, the heat source generated by this laser within the tissue has a significantly larger volume than that generated at longer wavelengths where water absorption dominates. It is possible that the peak temperature of the skin surface occurs a few milliseconds after a 10-ms, 1070-nm laser exposure has concluded. This is feasible because the tissue beneath the surface heats rapidly, and potentially, to a higher temperature than the surface due to multiple scattering events and elevated subsurface fluence levels. Under such conditions, the thermal diffusion time required for heat to propagate to the surface would cause the peak temperature observed at the surface of the skin to occur after a sufficiently short exposure has ceased. The thermal confinement time at this wavelength can be quite long, especially for large-diameter beams. As a result, the build-up of pressure from steam generated in the volume of tissue can cause superficial dry skin to explode off the surface. This observation may reveal that there is a transition in surface damage processes with lasers operating near 1070 nm for exposure durations much longer (milliseconds versus microseconds) than typically associated with photomechanical damage.

Analysis of thermal video data revealed the thermal slopes, mT (Table 4), which described the linear relationship between the temperature rise, ΔT, per Watt delivered to the exposure cite (Fig. 4) for a given exposure duration and beam diameter. Variance in the thermal data can be attributed to the angle of incidence from laser exposure (i.e., curved surface of flank), and water and melanin content in the skin. This linear relationship assists in understanding the tissue response. For instance, the mT for the 10-ms, 2.4-cm group and 250-ms, 9.5-cm group were both 0.0065°C/W (Table 4) and had similar peak temperatures at damage threshold (ΔTMVL30°C), but when the thermal slope was corrected for radiant exposure, the two exposure conditions were 2.99 and 1.82°C/(J/cm2), respectively. Furthermore, the thermal slopes for the 10-s exposures revealed that the ΔTMVL was 20°C, significantly lower than the 30°C temperatures found for most other exposures <10s. This is consistent with thermal data analysis from a 1319-nm laser radiation skin damage study, which also found ΔTMVL for 10-s exposures to be 20°C.7 These thermography data may provide evidence for the rate of reaction as a function of peak temperature when combined with mathematical modeling.

Computational models of laser-tissue interaction traditionally use the Arrhenius equation to simulate the reaction rate of the processes with the model forced to estimate damage at time points at or very near the end of the laser exposure.30 However, the experimental damage threshold is assessed at 1 and 24 h, much beyond the exposure duration of the laser. The appearance of lesions at the 24-h end point is complicated by numerous poorly understood pathways resulting in such sequelae as neurogenic inflammation and disappearing/reappearing erythema. When combined with computational models of laser-tissue interaction, these thermographic data (Table 4) will help to better understand the processes involved in lesion formation. Computational models of the data presented here are in progress.

Thermography data showed that regardless of exposure duration and beam diameter, hair follicles always reached higher temperatures than surrounding skin, but observation of lesions did not center around, nor seem directly related to, the follicle cites or the observation of vapor or smoke. The Yucatan mini-pigs used in this study had black colored hair, which always appeared darker in contrast to subjects’ skin. Hair that had not been effectively removed by clipping appeared to burn upon exposure with high irradiance levels. For example, consider the time-temperature histories presented in Fig. 5(b) where thermal camera images (Video 1 and Video 2) showed follicular cites reaching absolute temperatures in excess of 130+°C (camera saturation) while the skin in between these hot follicles peaked near 70°C. For the 9.5-cm, 250-ms case, the lowest energy tested was about half the ED50 (thermal data not shown) and still caused the hair follicles to reach 120+°C, but no lesion was observed. For a 2.4-cm, 10-ms case using a power at 980 W (a factor of 4.7 times less than the ED50 of 4660 W), slight follicular hot spots were observed to reach 41°C, only 4°C higher than the skin peak temperature. However, a 2.4-cm, 10-ms case at 1236 W (a factor of 3.7 times less than the ED50) showed follicular hot spots beginning to appear with a follicle peaking at 54°C while the skin peaked at 40.0°C. When the 10-ms, 2.4-cm exposure was about half of the ED50, some follicles peaked above 90°C while the skin peaked at 49°C. Thermal data from these lower-power 10-ms, 2.4-cm cases indicate the hair residing in the follicle had a thermal slope, mT, between three and five times greater than the mT for skin in this study.



For 1070-nm skin exposures, beam diameters >2.5cm seem to have a weak to no effect on threshold radiant exposure levels for exposure times between 0.01 and 0.25 s, but may have a larger effect on thresholds for exposures 10s. The MVL thresholds beam diameter dependency for 0.01 s is not well understood, but is suspected to involve photomechanical damage mechanisms confounded by the large volume, especially for larger beam diameters, of tissue heated very quickly by this wavelength. At 1070 nm, damage threshold values determined in this study indicate that a buffer of at least a factor of six in the MPE as defined by the ANSI Standard 2007 (Ref. 1) exists, with higher levels of protection afforded by the standard at smaller beam diameters. These margins of safety are consistent with those reported in the literature for the pigmented pig skin at other near-IR wavelength bands between 1000 and 2000 nm.


The authors would like to thank the staff of the Air Force Research Laboratory Veterinary Sciences Branch (711 HPW/RHDV) for their expert animal care and handling service and Robert Ulibarri and Shane Johnson from the Directed Energy Effects Branch (AFRL/RDLE) for training support and the use of their laser and diagnostics assets. The authors would also like to thank J. Michael Rickman for his help in making video files, Justin Zhoner, Doug Goddard, Michael Foltz, and Victor Villavicencio for helping with logistics for the experimental setup and data collection, and Michael Wyche, Isaac Noojin, and Andrew Wharmby for thermal video processing. Research performed by TASC Inc. was conducted under USAF Contract Numbers F41624-02-D-7003 and FA8650-08-D-6930.


1. American National Standards Institute “American National Standard for safe use of lasers,” in ANSI Z136.1, pp. 62–66, 76–78, Laser Institute of America, Orlando, FL (2007). Google Scholar

2. T. A. Egglestonet al., “Comparison of two porcine (Sus scrofa domestica) skin models for in vivo near-infrared laser exposure,” Comp. Med. 50(4), 391–397 (2000).COMEFT1532-0820 Google Scholar

3. J. W. Oliveret al., “Infrared skin damage thresholds from 1940-nm continuous-wave laser exposures,” J. Biomed. Opt. 15(6), 065008 (2010).JBOPFO1083-3668 http://dx.doi.org/10.1117/1.3523622 Google Scholar

4. J. W. Oliveret al., “Skin damage thresholds with continuous-wave laser exposures at near infrared wavelengths,” in Int. Laser Safety Conf., Laser Institute of America, Orlando, FL (2011). Google Scholar

5. C. P. Cainet al., “Porcine skin visible lesion thresholds for near-infrared lasers including modeling at two pulse durations and spot sizes,” J. Biomed. Opt. 11(4), 041109 (2006).JBOPFO1083-3668 http://dx.doi.org/10.1117/1.2338815 Google Scholar

6. C. P. Cainet al., “Visible lesion thresholds with pulse duration, spot size dependency, and model predictions for 1.54-microm, near-infrared laser pulses penetrating porcine skin,” J. Biomed. Opt. 11(2), 024001 (2006).JBOPFO1083-3668 http://dx.doi.org/10.1117/1.2187987 Google Scholar

7. J. W. Oliveret al., “Infrared skin damage thresholds from 1319-nm continuous wave laser exposures,” J. Biomed. Opt. 18(12), 125002 (2013).JBOPFO1083-3668 http://dx.doi.org/10.1117/1.JBO.18.12.125002 Google Scholar

8. T. L. TroyS. N. Thennadil, “Optical properties of human skin in the near infrared wavelength range of 1000 to 2200 nm,” J. Biomed. Opt. 6(2), 167–176 (2001).JBOPFO1083-3668 http://dx.doi.org/10.1117/1.1344191 Google Scholar

9. E. Salomatinaet al., “Optical properties of normal and cancerous human skin in the visible and near-infrared spectral range,” J. Biomed. Opt. 11(6), 064026 (2006).JBOPFO1083-3668 http://dx.doi.org/10.1117/1.2398928 Google Scholar

10. M. W. Baozhanget al., “Study of injury thresholds of CW Nd:YAG laser light for human skin,” Chin. J. Lasers 12(10), 582–585 (1985).ZHJIDO0258-7025 Google Scholar

11. W. Y. Tingbiet al., “Research on injury threshold of Chinese yellow skin by CW Nd:YAG laser beam irradiation,” Chin. J. Lasers 12(10), 626–628 (1985).ZHJIDO0258-7025 Google Scholar

12. W. Y. Tingbiet al., “Research on injury threshold of Chinese yellow skin irradiated with 300 us pulsed neodymium glass laser light,” Chin. J. Lasers 12(10), 628–630 (1985).ZHJIDO0258-7025 Google Scholar

13. L. Zhaozhang, “ Laser safety, protection standards reviewed,” Chin. J. Lasers 6(3), 141–144 (1986).ZHJIDO0258-7025 Google Scholar

14. R. R. Andersonet al., “Selective photothermolysis of cutaneous pigmentation by Q-switched Nd: YAG laser pulses at 1064, 532, and 355 nm,” J. Investig. Dermatol. 93(1), 28–32 (1989).JIDEAE0022-202X http://dx.doi.org/10.1111/jid.1989.93.issue-1 Google Scholar

15. R. J. Rockwell Jr.L. Goldman, Research on Human Skin Laser Damage Thresholds, pp. 1–148, USAF School of Aerospace Medicine, Cincinnati, OH (1974). Google Scholar

16. R. L. P. van Veenet al., “Determination of VIS- NIR absorption coefficients of mammalian fat, with time- and spatially resolved diffuse reflectance and transmission spectroscopy,” in Biomedical Topical Meeting, 14 April 2004, Miami Beach, Florida, Paper SF4, Optical Society of America, Washington, DC (2004). Google Scholar

17. D. J. Segelstein, “The complex refractive index of water,” MS Thesis, University of Missouri, Kansas City, Missouri (1981). Google Scholar

18. H. E. ArimotoM. EgawaY. Yamada, “Depth profile of diffuse reflectance near-infrared spectroscopy for measurement of water content in skin,” Skin Res. Technol. 11(1), 27–35 (2005).0909-752X http://dx.doi.org/10.1111/srt.2005.11.issue-1 Google Scholar

19. A. Rogganet al., “Optical properties of circulating human blood in the wavelength range 400–2500 nm,” J. Biomed. Opt. 4(1), 36–46 (1999).JBOPFO1083-3668 http://dx.doi.org/10.1117/1.429919 Google Scholar

20. S. L. JacquesR. D. GlickmanJ. A. Schwartz, “Internal absorption coefficient and threshold for pulsed laser disruption of melanosomes isolated from retinal pigment epithelium,” Proc. SPIE 2681, 468–477 (1996). http://dx.doi.org/10.1117/12.239608 Google Scholar

21. D. J. Finney, Probit Analysis, Cambridge University Press, Cambridge, UK (1971). Google Scholar

22. K. W. FosterE. F. FincherR. L. Moy, “Heat-induced ‘recall’ of treatment zone erythema following fractional resurfacing with a combination laser (1320nm/1440nm),” Arch. Dermatol. 144(10), 1398–1399 (2008).ARDEAC0003-987X http://dx.doi.org/10.1001/archderm.144.10.1398 Google Scholar

23. D. Roostermanet al., “Neuronal control of skin function: the skin as a neuroimmunoendocrine organ,” Physiol. Rev. 86(4), 1309–1379 (2006).PHREA70031-9333 http://dx.doi.org/10.1152/physrev.00026.2005 Google Scholar

24. M. Steinhoffet al., “Modern aspects of cutaneous neurogenic inflammation,” Arch. Dermatol. 139(11), 1479–1488 (2003).ARDEAC0003-987X http://dx.doi.org/10.1001/archderm.139.11.1479 Google Scholar

25. J. M. Handley, “Adverse events associated with nonablative cutaneous visible and infrared laser treatment,” J. Am. Acad. Dermatol. 55(3), 482–489 (2006).JAADDB0190-9622 http://dx.doi.org/10.1016/j.jaad.2006.03.029 Google Scholar

26. B. Algermissenet al., “Laser-induced weal and flare reactions: clinical aspects and pharmacological modulation,” Br. J. Dermatol. 146(5), 863–868 (2002).BJDEAZ1365-2133 http://dx.doi.org/10.1046/j.1365-2133.2002.04723.x Google Scholar

27. B. Chenet al., “Histological and modeling study of skin thermal injury to 2.0 micron laser irradiation,” Lasers Surg. Med. 40(5), 358–370 (2008).LSMEDI0196-8092 http://dx.doi.org/10.1002/(ISSN)1096-9101 Google Scholar

28. B. Chenet al., “Effect of pigmentation density upon 2.0 micron laser irradiation thermal response,” Health Phys. 93(4), 273–278 (2007).HLTPAO0017-9078 http://dx.doi.org/10.1097/01.HP.0000264936.67385.04 Google Scholar

29. Y. Zhanget al., “Comparative study of 1,064-nm laser-induced skin burn and thermal skin burn,” Cell Biochem. Biophys. 67(3), 1005–1014 (2013).BJDEAZ1365-2133 http://dx.doi.org/10.1007/s12013-013-9596-6 Google Scholar

30. S. L. Jacques, “Ratio of entropy to enthalpy in thermal transitions in biological tissues,” J. Biomed. Opt. 11(4), 041108 (2006).JBOPFO1083-3668 http://dx.doi.org/10.1117/1.2343437 Google Scholar

Biographies of the authors are not available.

© The Authors. Published by SPIE under a Creative Commons Attribution 3.0 Unported License. Distribution or reproduction of this work in whole or in part requires full attribution of the original publication, including its DOI.
Rebecca L. Vincelette, Rebecca L. Vincelette, Gary D. Noojin, Gary D. Noojin, Corey A. Harbert, Corey A. Harbert, Kurt J. Schuster, Kurt J. Schuster, Aurora D. Shingledecker, Aurora D. Shingledecker, David J. Stolarski, David J. Stolarski, Semih S. Kumru, Semih S. Kumru, Jeffrey W. Oliver, Jeffrey W. Oliver, } "Porcine skin damage thresholds for 0.6 to 9.5 cm beam diameters from 1070-nm continuous-wave infrared laser radiation," Journal of Biomedical Optics 19(3), 035007 (21 March 2014). https://doi.org/10.1117/1.JBO.19.3.035007 . Submission:

Back to Top