Nebular Abundances of Fifteen Intermediate Luminosity Star-Forming Galaxies

Metal Abundances of KISS Galaxies. V. Nebular Abundances of Fifteen Intermediate Luminosity Star-Forming Galaxies


We present high S/N spectroscopy of 15 emission-line galaxies (ELGs) cataloged in the KPNO International Spectroscopic Survey (KISS), selected for their possession of high equivalent width [O iii] lines. The primary goal of this study was to attempt to derive direct-method () abundances for use in constraining the upper-metallicity branch of the  relation. The spectra cover the full optical region from [O ii]3726,3729 to [S iii]9069,9531 and include the measurement of [O iii]4363 in 13 objects. From these spectra, we determine abundance ratios of helium, nitrogen, oxygen, neon, sulfur, and argon. We find these galaxies to predominantly possess oxygen abundances in the range of 8.0 12+log(O/H)  8.3. We present a comparison of direct-method abundances with empirical SEL techniques, revealing several discrepancies. We also present a comparison of direct-method oxygen abundance calculations using electron temperatures determined from emission lines of O and S, finding a small systematic shift to lower  (1184 K) and higher metallicity (0.14 dex) for sulfur-derived  compared to oxygen-derived . Finally, we explore in some detail the different spectral activity types of targets in our sample, including regular star-forming galaxies, those with suspected AGN contamination, and a local pair of low-metallicity, high-luminosity compact objects.

Subject headings:
galaxies: abundances – galaxies: starburst – H ii regions

1. Introduction

Analysis of emission-line spectra provides the principal diagnostic tools for determining gas-phase physical characteristics of star-forming galaxies. Heavy element abundances are valuable in the analysis of star formation histories (e.g., Kennicutt, 1998; Contini et al., 2002), chemical evolution (e.g., Steidel et al., 1996; Kobulnicky & Kewley, 2004; Lamareille et al., 2006), and general galactic evolution trends such as the luminosity-metallicity relation (e.g., Rubin et al., 1984; Lee et al., 2004; van Zee et al., 2006), and the mass-metallicity relation (e.g., Lequeux et al., 1979; Tremonti et al., 2004). Determination of this gas-phase property remains a difficult task for all but the most nearby and luminous H ii regions and emission-line galaxies (ELGs), and is intrinsically limited both by our ability to accurately measure intensity ratios of nebular emission lines and by our calibration of the conversion of observational data into abundance measurements.
Metallicities determined using the temperature of the electron gas (the so-called direct- or -method), considered to be the most accurate and robust available, are calculated using measurements of often exceedingly weak temperature-sensitive auroral emission lines such as [O iii]4363 (Osterbrock & Ferland, 2006). In computing the chemical composition of photoionized nebulae via the direct-method, the limiting factor is typically the measurement precision of the electron gas temperature. Furthermore, the emission strength of the temperature-sensitive auroral lines scales with increasing , and are thus preferentially detected in the hotter nebulae of low-abundance H ii regions and ELGs. Determination of the electron temperature for high-metallicity systems is therefore a difficult observational task. Reliable observations of high quality are thus relatively rare in the literature.
For the majority of cases in which direct metallicity determinations are not possible, relations utilizing consistently detectable strong lines have been developed. These are either empirical relations that are tied to direct-method abundances, theoretical relations that are based on photoionization models, or semi-empirical relations that are a mixture of the two (e.g., Moustakas et al., 2010). The most common such strong-emission-line (SEL) method is  = ([O ii]3726,3729+[O iii]4959,5007)/H, first introduced by Pagel et al. (1979).  provides an estimate of total cooling and thus represents a proxy for global metallicity (Kewley & Dopita, 2002), but is famously double-valued at high- and low-abundance and includes an ambiguous turnaround region. Calibration of this relation relies upon comparison to accurate measurements of metallicity, derived via the direct-method or through photoionization models (e.g., Kewley et al., 2001). Many  calibrations of the latter sort are available in the literature (e.g., Edmunds & Pagel, 1984; Skillman et al., 1989; Zaritsky et al., 1994). As -method abundances are more readily determined for lower-abundance systems, the  relation is much better calibrated using direct-method data below the turnaround region than above (Pilyugin, 2000; Pilyugin & Thuan, 2005).
The primary goal of this study is to attempt to derive direct-method () abundances of star-forming galaxies for use in constraining the upper-metallicity branch of the  relation. We utilize the superior light-gathering capability of the Keck 10-m telescope to observe a specially selected set of strong-lined ELGs from the KISS sample (Salzer et al., 2000, 2001). These observations were planned with the hope of detecting and measuring the temperature-sensitive [O iii] emission line in higher-metallicity star-forming systems situated above the  relation turnaround region. With a robust calibration of , it becomes possible for the full range of ELG metallicities to be reliably studied, rather than being limited only to the brightest and most metal-poor objects. Accurate abundances for high-metallicity targets allows for the comprehensive understanding of a more representative sample of star-forming galaxies in the local Universe. Furthermore, the refinement of reliable empirical metallicity estimation methods from our observations of the total integrated light of the ELGs is particularly useful in the study of higher-redshift galaxies, where spectral observations are global by necessity and the detection of auroral lines becomes observationally difficult. Calibration of the  method using global spectra of local star-forming galaxies will provide for a more consistent approach to metallicity relations and chemical enrichment studies across cosmic time.
There are often offsets in metallicity estimates when employing direct- and SEL-methods, where in general -method abundances are low as compared to techniques utilizing ratios of strong lines (e.g., Kennicutt et al., 2003; Bresolin et al., 2005; Zurita & Bresolin, 2012). Biases in the sample selection of observational targets used to calibrate some empirical relations can be a major contributor to such abundance discrepancies. The magnitude of these offsets is typically of order a few tenths of a dex. A more comprehensively populated  relation, including -method abundances of the upper-metallicity branch, will help to reduce such discrepancies. An alternative approach to resolving this disparity involves a reexamination of the -method. In computing direct-method metallicities, the electron temperature can be estimated from ions of similar electronic structure to O, such as doubly-ionized sulfur (e.g., Garnett, 1989; Bresolin et al., 2005; Pérez-Montero et al., 2006). Models show that the S zone straddles the zones of both O and O (Garnett, 1992). The electron temperature derived from [S iii] emission lines may therefore be a better single representation of the temperature of the nebula as a whole.
In §2 we discuss the galaxy sample, the observations, and the data-reduction methods. In §3 we present the spectral data obtained using the W. M. Keck Observatory. Section 4 presents a detailed analysis of the 13 objects with abundance-quality spectra. We measure properties of the nebular emission regions such as electron temperature and electron density and calculate abundance ratios of heavy elements with respect to hydrogen. Section 5 discusses the abundance results, presenting the preliminary findings of our program to constrain the upper-metallicity branch of the  empirical strong-line relation. We do so with direct () method abundances of relatively metal-rich ELGs, using the KISS sample for target selection and the Keck telescope to deliver high-quality spectra in order to detect weak temperature-sensitive emission lines. In addition, we compare our abundance results with various strong-emission-line (SEL) abundance determination methods and calibrations available in the literature, and compare metallicity estimates computed via the direct method using electron temperatures derived from both doubly-ionized oxygen (O) and sulfur (S) to test the effect of sampling different zones of ionization within the nebula. Finally, we scrutinize our sample to better understand the physical limitations inherent to studies of this nature, compare direct-method abundances calculated using electron temperatures characteristic of different ionization zones, and explore in more detail the three categories of objects uncovered in our observations. Our results are summarized in §6.

2. Observations and Data Reduction

2.1. Sample Selection

Our spectroscopic targets for this project were potentially metal-rich galaxies identified by the KPNO International Spectroscopic Survey (KISS; Salzer et al., 2000). KISS used low-dispersion objective-prism spectra to identify ELG candidates via detection of line emission in galaxies with redshifts of less than 0.095. The “red” survey (Salzer et al., 2001; Gronwall et al., 2004b; Jangren et al., 2005) cataloged objects by means of the H line, while the “blue” survey (Salzer et al., 2002) distinguished objects through strong [O iii]5007 line emission. The KISS sample is comprised of a wide variety of ELGs, including starburst nucleus galaxies, H ii galaxies, irregular galaxies with significant star formation, and blue compact dwarfs (BCDs). In addition, the survey detected low-ionization nuclear emission-line regions (LINERs) and Seyfert galaxies as well as quasars whose redshifts place their emission lines into the objective-prism spectral bandpass.
The galaxies discussed in this paper were chosen from the KISS H-selected catalog with selection criteria based on a combination of [O iii] line strength and coarse abundance estimates (Melbourne & Salzer, 2002; Salzer et al., 2005b). We were specifically interested in objects for which [O iii]4363 might be detectable in order to obtain the nebular electron temperature , which is measured using the line ratio [O iii]4363/[O iii]4959,5007 (e.g., Izotov et al., 1994). Additionally, we are interested in galaxies with high enough oxygen abundance to be located above the turnaround region of the  relation (12+log(O/H)  8.3) as a means to constrain/calibrate the location of the upper-metallicity branch.
The selection of targets for this study was based on the emission-line data obtained from “quick-look” follow-up spectra of KISS galaxies (Wegner et al., 2003; Melbourne et al., 2004; Gronwall et al., 2004a; Jangren et al., 2005; Salzer et al., 2005a). Based on the published spectral data, we were able to generate lists of potentially interesting objects. These spectra were usually not of high enough quality for the measurement of direct abundances, but they did yield reasonably accurate strong-line ratios that gave a good estimate of the properties of each object.
Galaxies were chosen for observation based on their inclusion in one of three lists of KISS objects generated with varying selection criteria. First, our targets were selected to be among those with the highest line strengths, in order to have the best chance of detecting the weak temperature-sensitive emission lines. Using the equivalent width (EW) of the [O iii]5007 line measured in our quick-look spectra as a discriminant, list 1 requires an EW > 100 Å, list 2 requires an EW between 50 and 100 Å, and list 3 requires an EW between 5 and 50 Å. Second, our objects require higher abundances, as our goal is to probe into the upper-metallicity branch of the  relation. Using coarse abundance estimates (Melbourne & Salzer, 2002; Salzer et al., 2005b), list 1 requires a 12+log(O/H) > 8.2, list 2 requires a 12+log(O/H) > 8.3, and list 3 requires a 12+log(O/H) > 8.5. For the third list, additional selection criteria include the requirement that the Balmer decrement c < 0.7 in order to avoid the most heavily reddened objects. The number of galaxies selected to each of the three lists is as follows: List 1 includes 20 objects, list 2 includes 22 objects, and list 3 includes 28 objects. In total, four objects were observed from list 1, eight objects from list 2, and three objects from list 3, for fifteen total objects observed during our run. Images of each of the fifteen galaxies taken from the Sloan Digital Sky Survey (SDSS; Abazajian et al., 2005) are shown as Figure 1.
The galaxies selected for this study are displayed on a spectral activity diagnostic diagram (e.g., Baldwin et al., 1981) in Figure 2 using emission-line ratios taken from the original KISS survey data, overlaying data points representing all star-forming galaxies from KISS. Objects from the three target lists are indicated by symbols of different color, while those to be observed for this study with Keck are further indicated by the heavy black rings. The solid line represents a star-formation sequence derived by Dopita & Evans (1986) from theoretical models, while the dashed line represents an empirically defined demarcation between starburst galaxies and AGN from Kauffmann et al. (2003). It should be noted that all initial observations of relevant line ratios place the galaxies considered for this study firmly within the boundaries of regular star-forming galaxies (i.e., no targets were believed to exhibit spectral contamination from an AGN). In general, these criteria allowed us to identify starburst galaxy candidates of medium- to high-abundance with a moderate degree of success. In a few cases, however, we discovered examples of galaxies with spectral contamination from AGN and galaxies with anomalously low abundances and high luminosities for their redshifts. We will explore the inhomogeneity of these targets within our sample in greater detail in §5.2. A list of all observed galaxies with coordinates and basic information is presented as Table 1, including references to important information relevant to KISS from these papers.

Figure 1.— The fifteen galaxies selected from KISS observed for this study, as taken from the Sloan Digital Sky Survey (SDSS) (Abazajian et al., 2005). Targets range in redshift between = 0.00791 and = 0.09115. Image boxes are 100 arcsec per side and in some cases the objects are seen to be quite compact.
Figure 2.— Diagnostic diagram of KISS star-forming galaxies using emission-line ratios from the original survey data, with galaxies from each of the three lists of varying selection criterion overlaid as blue circles (List 1: High [O iii] EWs), green circles (List 2: Lower [O iii] EWs), and red circles (List 3: Higher Coarse Abundances). The fifteen galaxies to be observed with Keck for this study are labeled and signified with black rings. The solid line represents a star formation sequence derived from theoretical models (Dopita & Evans, 1986), while the dashed line represents an empirically defined demarcation between starburst galaxies and AGN (Kauffmann et al., 2003).
KISSR Field KISS ID RA Dec [sec] [sec] [sec]
0148 F1250 317 12:54:45.4 +28:55:31 17.06 -15.93 0.00791 1600 1600
0242 F1315 3582 13:16:03.9 +29:22:54 16.71 -19.46 0.03782 2900, 1300 2900 1300
0258 F1315 1234 13:18:40.3 +28:59:57 17.57 -20.15 0.07567 2600 2600
0451 F1405 3370 14:06:35.0 +29:27:37 17.81 -18.59 0.04158 3600 2600 1600
0475 F1410 4626 14:11:06.6 +29:22:21 17.58 -19.89 0.06641 3600 2600 1600
0512 F1415 2664 14:17:48.4 +29:03:11 17.36 -20.46 0.07719 2900, 1600 2900 1600
0590 F1455 8836 14:55:01.2 +29:33:31 19.28 -16.48 0.03072 3600 2600 1600
0653 F1510 4675 15:10:52.7 +28:46:46 17.53 -20.10 0.06990 3600 2600 1600
1056 F1635 4097 16:37:18.2 +29:51:11 19.32 -18.91 0.09115 2600 2600
1379 F1246 1011 12:48:13.4 +43:57:04 18.02 -17.67 0.02969 3600 2600 1600
1734 F1411 4039 14:10:59.2 +43:02:47 18.44 -18.97 0.06556 2600, 1300 2600 1300
1955 F1528 4743 15:27:59.4 +42:50:58 16.83 -17.88 0.01849 1900 1900
2117 F1610 7310 16:09:53.6 +44:13:22 16.58 -19.59 0.03716 3600 2600 1600
2125 F1610 5338 16:10:20.4 +43:00:35 16.58 -18.73 0.02517 11200, 1300 11200 1300
2132 F1610 3202 16:11:45.8 +43:30:44 17.78 -18.37 0.03685 2600 2600
Table 1Object Parameters and Identifiers for Observed Galaxies

2.2. Observations

Spectra of 15 KISS ELGs were obtained using the Keck I 10 m telescope on 25 May 2006 with the Low-Resolution Imaging Spectrometer (LRIS), a double spectrograph that includes a dichroic that directs light toward a red and a blue side (Oke et al., 1995). The use of two distinct dispersion elements and cameras enables simultaneous spectral coverage spanning from 3000 Å to 10,000 Å with good dispersion. The KISS catalog is comprised of galaxies characterized by a single (often central) dominant unresolved source of emission, which for our targets were mostly encompassed within our 1”.5 wide slit. For each observation, the slit was aligned with the parallactic angle to accommodate for differential atmospheric dispersion. Blue side exposures used the 600/400 grism, with a dispersion of 0.63 Å pixel, a resolution of 6 Å, and wavelength coverage between 3000 Å and 5600 Å. Red side exposures used the 900/5500 grating, with a dispersion of 0.53 Å pixel, a resolution of 3 Å, and wavelength coverage between 5700 Å and 7300 Å. Each galaxy was first observed using the blue and red sides simultaneously to determine if further observations would yield usable spectra. If the raw spectra were of substantially high quality, the red side grating was switched to observe NIR lines. This was completed for ten of the fifteen galaxies. For NIR exposures, the red side was switched to the 400/8500 grating, with a dispersion of 1.16 Å pixel, a resolution of 6 Å, and wavelength coverage between 6300 Å to 10,000 Å in order to observe the near-infrared [S iii] nebular lines at 9069 and 9531. Typical individual exposure times were 600 seconds, although some targets necessitated individual exposure times of 1200, 900, or 300 seconds. Exposure time information is listed in Table 1.
Example spectra of one galaxy from this study (KISSR 0242) seen as Figures 3 and 4 illustrates the wavelength coverage afforded by Keck LRIS, extending from the near-UV to the near-IR. Furthermore, Figure 4 reveals clear detection of, e.g., the [O iii]4363 line required to derive -method oxygen abundances. This object is of particularly high excitation and signal strength and should not be taken as typical for our study, but rather its spectral characteristics were chosen to demonstrate the capabilities of our apparatus and observing program.

Figure 3.— Optical spectrum of KISSR 0242, in all three wavelength regimes afforded by Keck LRIS, including “blue”, “red”, and “infrared”. Emission lines that are important to abundance analysis for this study are labeled.
Figure 4.— Same as Figure 3, with reduced vertical scaling to highlight emission line detail. Emission lines important for the abundance analysis for this study are labeled.

For five of the fifteen galaxies in our sample, an inspection of the first two exposures indicated that the object did not possess adequately strong emission lines to warrant further observation, and so the red side grating was not switched and thus observations of the near-infrared [S iii] lines were not acquired. Our observations included spectra of a Hg-Cd-Zn and Ar-Ne lamps to set the wavelength scale as well as observations of the spectrophotometric standard stars GD 153 (Bohlin et al., 1995), Feige 56 (Hamuy et al., 1992), Feige 66, BD+28 4211, and Feige 110 (Massey et al., 1988) for flux calibration.

2.3. Data Reduction

The data reduction was carried out with the Image Reduction and Analysis Facility (IRAF1). Processing of the two-dimensional spectral images followed standard methods. All of the reduction steps mentioned below were carried out on the red and blue spectral images independently. The mean bias level was determined and subtracted automatically from each image by the data acquisition software. A mean bias image was created by combining 11 zero-second exposures taken on the night of the observations. This image was subtracted to correct the science images for any possible two-dimensional structure in the bias. Flat-fielding was achieved using an average-combined quartz lamp image that was corrected for the wavelength-dependent response of the system. For cosmic ray rejection, we used L.A.Cosmic (van Dokkum, 2001).
One-dimensional spectra were extracted using the IRAF APALL routine. The extraction width (i.e., distance along the slit) was selected independently for each source by examination. As the LRIS spectrograph uses a dichroic beam splitter to break the spectrum into red and blue sections, in order to analyze the data it is necessary that the two sections are on a consistent flux scale with each other, accomplished by setting the extraction regions to be the same. For compact targets our observations acquired global spectra, while for more extended galaxies we analyzed the nuclear region only. In all cases, the extraction width was limited to 10 pixels or 2.15 arcseconds for the red side CCD (pixel scale of 0.215” pix). For the blue side CCD we scaled the extraction width accordingly to 15.62 pixels (from a pixel scale of 0.135” pix). Sky subtraction was also performed at this stage, with the sky spectrum being measured in regions on either side of the object extraction window. The Hg-Cd-Zn and Ar-Ne lamp spectra were used to assign a wavelength scale, and the spectra of the spectrophotometric standard stars were used to establish the flux scale. The standard star data were also used to correct the spectra of our target ELGs for telluric absorption. This is important because at the redshifts of some of our targets the [S ii] doublet falls in the strong telluric -band and would lead to a severe underestimate of the true emission-line flux.
We utilized the galaxy observations to measure any low-level flux scale differences between the red and blue sides. To do this, we measured the continuum flux of the galaxies in 50 Å bins from 5200 Å to 6000 Å. A plot of flux versus wavelength reveals any small shifts between the continuum fluxes on either side of the dichroic break. Since the spectral extraction regions covered the same angular extent on both sides, the correction factors were expected to be small. This turned out to be the case, with the derived correction factors for each observation averaging a few percent. For our instrumental setup, H and H are always on opposite sides of the dichroic. Therefore, the main effect of the correction factor is to modify the H/H ratio and the value of the reddening parameter c. This necessary procedure introduces an additional level of uncertainty into the derivation of the nebular abundances in the sense that it affects reddening corrections. Great care was taken in ensuring that the red and blue sections of the spectra were extracted consistently.

3. The Spectral Data

Emission-line flux measurements were carried out using the IRAF routine SPLOT. When multiple images of the same object were available in any of the three wavelength ranges, the spectra were reduced separately then combined into a single high S/N spectrum prior to the measurement stage. For an initial estimate of the internal reddening in each galaxy we calculate c from the H/H line ratio. The value of c is then used to correct the measured line ratios for reddening, following the standard procedure (e.g., Osterbrock & Ferland, 2006),

where ) is derived from studies of absorption in the Milky Way (using values taken from Rayo et al., 1982). The final c values used for the reddening correction were determined from a simultaneous fit to the reddening and underlying stellar absorption in the Balmer lines, using all available Balmer line ratios. Under the simplifying assumption that the same EW of underlying absorption applies to each Balmer line, an absorption line correction was applied to each spectrum that ranged in value from 0 to 5 Å. The value for the underlying absorption was varied until a self-consistent value of c for the Balmer line ratios, H/H, H/H, and H/H, was found. This process led to the determination of a characteristic value for the underlying absorption for each galaxy. The average value was typically equal to 2 Å (e.g., McCall et al., 1985; Skillman & Kennicutt, 1993); line strengths of the Balmer emission lines were adjusted accordingly. In the case of KISSR 0512, the derived value for c is slightly negative, although the formal error in c is such that the measured value is consistent with c = 0.00. We therefore adopt c = 0.00 for the computation of reddening-corrected line ratios for this object. Results of our spectroscopic observations of the Keck sample, presented as reddening-corrected line ratios relative to H, are presented in Table 2.
Determinations of c for each galaxy derived from the earlier follow-up spectroscopy of the original KISS galaxies were used as a source of comparison for our study (Wegner et al., 2003; Melbourne et al., 2004; Gronwall et al., 2004a; Jangren et al., 2005; Salzer et al., 2005a). In the majority of cases, these values agree. For some galaxies, however, disparities in value may represent cases in which the placement of the slit was different between the two sets of observations, and as such properties of different regions of the galaxy were measured (e.g., a nuclear emission region compared to a knot of H ii in a spiral arm). The spectra of our fifteen objects are presented in Figures 5, 6, and 7.

Figure 5.— Spectra of “blue” and “red” wavelength ranges of the first five galaxies studied for this project, with inset showing the H and [O iii]4363 emission lines.
Figure 6.— Spectra of “blue” and “red” wavelength ranges of the next five galaxies studied for this project, with inset showing the H and [O iii]4363 emission lines.
Figure 7.— Spectra of “blue” and “red” wavelength ranges of the last five galaxies studied for this project, with inset showing the H and [O iii]4363 emission lines.

4. Metal Abundances

4.1. Electron Density and Temperature

In determining the physical characteristics of H ii regions and ELGs, a two-zone model for the star-forming nebulae is typically assumed. The outer zone is represented by a lower-temperature region in which oxygen is singly-ionized, while the inner zone is a medium-temperature zone in which oxygen is doubly-ionized. Hydrogen is assumed to be fully ionized within the radius of the outer zone, and is assumed to be neutral outside this zone.
In each zone we would ideally like to know both the electron density and temperature. The available observational diagnostics, however, only allow for an electron density measurement in the low-temperature zone and an electron temperature measurement in the medium-temperature zone. We use the [S ii6716/6731 (Osterbrock & Ferland, 2006) ratio to determine the electron density, which in all measurable cases is roughly 100 e cm. For galaxies for which electron density was not derived from our observations, we assume a density of 100 e cm. Our abundance analysis assumes that the density does not vary from zone to zone.


The electron temperature in the medium-temperature zone is determined by the oxygen line ratio [O iii]4363/[O iii]4959,5007, which represents a comparison in doubly-ionized oxygen of electron transitions with considerably different excitation energies (Osterbrock & Ferland, 2006). As can be seen by inspection of the insets in Figures 5, 6, and 7, several of our 4363 detections are quite weak, and great care was required in their measurement. The uncertainties in the final 4363 line fluxes are reflected in the large uncertainties in the derived values of . Calculations of both the electron density and electron temperature are carried out using the program ELSA (Johnson et al., 2006), which makes use of the latest collision strengths and radiative transition probabilities.
We estimate the temperature in the low-temperature zone by using the algorithm presented in Skillman et al. (1994) based on the nebular models of Stasińska (1990),

where are temperatures measured in units of 10 K. The measured electron densities and temperatures are presented in Table 3.

Ion [Å] KISSR 0148 KISSR 0242 KISSR 0258 KISSR 0451 KISSR 0475 KISSR 0512 KISSR 0590 KISSR 0653
[O ii] 3727 3.287 0.135 2.816 0.115 3.241 0.140 3.908 0.164 3.582 0.147 1.767 0.039 3.501 0.146 2.939 0.122
H 10 3798 0.042 0.004 0.040 0.002 0.025 0.002
H 9 3835 0.065 0.004 0.064 0.002 0.064 0.008 0.052 0.003
[Ne iii] 3869 0.201 0.008 0.404 0.015 0.134 0.009 0.243 0.010 0.171 0.007 0.050 0.004 0.222 0.010 0.106 0.005
He I+H 8 3889 0.175 0.007 0.198 0.007 0.140 0.009 0.161 0.007 0.168 0.007 0.102 0.005 0.152 0.007 0.170 0.007
[Ne iii]+H 3970 0.184 0.007 0.271 0.009 0.124 0.008 0.133 0.005 0.141 0.005 0.090 0.004 0.160 0.007 0.136 0.005
He I 4026 0.013 0.001
[S ii] B 4073 0.023 0.001 0.045 0.007 0.036 0.002 0.027 0.002 0.034 0.004 0.032 0.003 0.024 0.002
H 4102 0.249 0.009 0.272 0.008 0.234 0.010 0.235 0.008 0.240 0.008 0.194 0.006 0.239 0.009 0.236 0.008
H 4340 0.456 0.013 0.491 0.013 0.426 0.014 0.455 0.012 0.452 0.012 0.387 0.010 0.456 0.013 0.455 0.012
[O iii] 4363 0.017 0.003 0.051 0.001 0.057 0.006 0.028 0.002 0.018 0.002 0.030 0.003 0.013 0.002
He I 4472 0.035 0.003 0.042 0.001 0.025 0.006 0.050 0.002 0.050 0.002 0.033 0.004 0.048 0.003 0.045 0.002
He II 4687 0.006 0.001 0.004 0.001
[Ar iv] 4711 0.006 0.001
[Ar iv] 4740
H 4861 1.000 0.022 1.000 0.021 1.000 0.024 1.000 0.022 1.000 0.021 1.000 0.023 1.000 0.023 1.000 0.021
He I 4922 0.009 0.001
[O iii] 4959 0.862 0.019 1.357 0.028 0.446 0.012 0.952 0.021 0.701 0.015 0.293 0.008 0.883 0.020 0.592 0.013
[O iii] 5007 2.589 0.057 4.091 0.086 1.294 0.030 2.873 0.062 2.091 0.045 0.867 0.020 2.644 0.059 1.774 0.038
[N i] 5199 0.006 0.003 0.017 0.001 0.034 0.002 0.028 0.003
[N ii] 5755 0.011 0.001
He I 5876 0.107 0.004 0.117 0.004 0.092 0.005 0.119 0.004 0.107 0.004 0.087 0.004 0.112 0.004 0.117 0.004
[O i] 6300 0.074 0.004 0.051 0.002 0.090 0.005 0.057 0.003 0.057 0.002 0.086 0.004 0.067 0.003 0.047 0.002
[S iii] 6312 0.020 0.002 0.012 0.001 0.010 0.001 0.010 0.001 0.013 0.002 0.009 0.001
[O i] 6364 0.024 0.002 0.016 0.001 0.031 0.003 0.023 0.001 0.015 0.001 0.024 0.003 0.014 0.002 0.014 0.001
[N ii] 6548 0.084 0.004 0.090 0.004 0.385 0.017 0.100 0.005 0.161 0.007 0.319 0.008 0.109 0.005 0.143 0.006
H 6563 2.858 0.123 2.821 0.121 2.732 0.123 2.837 0.124 2.845 0.122 2.705 0.059 2.827 0.123 2.854 0.124
[N ii] 6583 0.259 0.012 0.230 0.010 1.110 0.051 0.316 0.014 0.492 0.021 0.981 0.022 0.328 0.015 0.416 0.018
He I 6678 0.030 0.002 0.028 0.001 0.015 0.003 0.026 0.002 0.033 0.002 0.018 0.003 0.021 0.002 0.029 0.002
[S ii] 6717 0.367 0.017 0.207 0.010 0.414 0.020 0.315 0.015 0.348 0.016 0.435 0.010 0.381 0.018 0.315 0.015
[S ii] 6731 0.255 0.012 0.154 0.007 0.332 0.016 0.231 0.011 0.262 0.012 0.317 0.008 0.272 0.013 0.237 0.011
He I 7065 0.016 0.002 0.028 0.001 0.017 0.001 0.013 0.001 0.019 0.002 0.019 0.001
[Ar iii] 7136 0.082 0.005 0.071 0.004 0.056 0.003 0.018 0.001 0.049 0.003 0.079 0.005 0.060 0.003
[O ii] 7319 0.031 0.002 0.017 0.001 0.035 0.003
[O ii] B 7325 0.038 0.002 0.045 0.003 0.034 0.002 0.057 0.004 0.040 0.002
[O ii] 7330 0.008 0.001 0.017 0.001 0.022 0.002
[S iii] 9069 0.145 0.011 0.198 0.016 0.168 0.013 0.383 0.009 0.205 0.016 0.244 0.019
[S iii] 9531 0.358 0.029 0.424 0.035 0.359 0.029
Table 2Emission-Line Intensity Ratios Relative to H for Observed Galaxies (1 of 2)
Ion [Å] KISSR 1056 KISSR 1379 KISSR 1734 KISSR 1955 KISSR 2117 KISSR 2125 KISSR 2132
[O ii] 3727 3.204 0.133 3.022 0.125 1.806 0.073 3.437 0.150 3.466 0.142 3.490 0.156 4.188 0.174
H 10 3798 0.029 0.003 0.042 0.002 0.048 0.002
H 9 3835 0.058 0.004 0.048 0.003 0.065 0.003 0.069 0.003
[Ne iii] 3869 0.272 0.011 0.188 0.007 0.489 0.018 0.095 0.008 0.168 0.007 0.068 0.005 0.250 0.010
He I+H 8 3889 0.156 0.007 0.185 0.007 0.189 0.007 0.178 0.010 0.195 0.008 0.134 0.007 0.168 0.008
[Ne iii]+H 3970 0.182 0.007 0.162 0.006 0.303 0.010 0.113 0.008 0.208 0.007 0.099 0.005 0.151 0.007
He I 4026 0.016 0.001 0.015 0.001
[S ii] B 4073 0.037 0.003 0.019 0.002 0.025 0.001 0.049 0.006 0.021 0.001 0.029 0.003 0.063 0.004
H 4102 0.238 0.008 0.256 0.008 0.259 0.008 0.235 0.011 0.278 0.009 0.226 0.008 0.232 0.008
H 4340 0.452 0.013 0.478 0.013 0.473 0.012 0.439 0.015 0.492 0.013 0.461 0.013 0.450 0.013
[O iii] 4363 0.027 0.002 0.019 0.002 0.076 0.002 0.015 0.001 0.020 0.003 0.028 0.003
He I 4472 0.050 0.003 0.055 0.002 0.042 0.001 0.077 0.006 0.039 0.002 0.043 0.003 0.029 0.003
He II 4687 0.009 0.001 0.009 0.001
[Ar iv] 4711 0.004 0.001 0.009 0.001 0.009 0.005
[Ar iv] 4740 0.004 0.001
H 4861 1.000 0.022 1.000 0.021 1.000 0.021 1.000 0.025 1.000 0.022 1.000 0.021 1.000 0.022
He I 4922 0.008 0.001 0.009 0.001
[O iii] 4959 1.027 0.023 0.803 0.017 1.666 0.035 0.258 0.008 0.859 0.019 0.253 0.006 0.925 0.021
[O iii] 5007 3.051 0.067 2.389 0.051 4.969 0.105 0.795 0.020 2.593 0.056 0.744 0.016 2.776 0.061
[N i] 5199 0.023 0.001 0.022 0.005 0.021 0.001 0.044 0.002 0.026 0.003
[N ii] 5755 0.007 0.001 0.006 0.001
He I 5876 0.108 0.004 0.119 0.004 0.117 0.004 0.109 0.005 0.103 0.004 0.098 0.004 0.114 0.005
[O i] 6300 0.078 0.003 0.043 0.002 0.048 0.002 0.061 0.004 0.053 0.002 0.054 0.002 0.064 0.003
[S iii] 6312 0.011 0.001 0.015 0.001 0.011 0.003 0.009 0.001 0.005 0.001 0.008 0.002
[O i] 6364 0.028 0.002 0.011 0.001 0.015 0.001 0.016 0.003 0.016 0.001 0.019 0.001 0.015 0.002
[N ii] 6548 0.103 0.005 0.103 0.005 0.097 0.004 0.222 0.011 0.141 0.006 0.225 0.010 0.095 0.005
H 6563 2.843 0.123 2.849 0.123 2.808 0.119 2.860 0.129 2.867 0.123 2.772 0.130 2.835 0.123
[N ii] 6583 0.325 0.014 0.311 0.014 0.277 0.012 0.680 0.031 0.423 0.018 0.708 0.034 0.305 0.014
He I 6678 0.024 0.002 0.031 0.002 0.032 0.001 0.026 0.003 0.024 0.001 0.020 0.001 0.034 0.003
[S ii] 6717 0.361 0.017 0.286 0.013 0.176 0.008 0.572 0.028 0.299 0.014 0.424 0.022 0.400 0.019
[S ii] 6731 0.268 0.013 0.207 0.010 0.136 0.006 0.414 0.021 0.226 0.011 0.305 0.016 0.289 0.014
He I 7065 0.017 0.001 0.038 0.002 0.010 0.003 0.033 0.002 0.009 0.001
[Ar iii] 7136 0.073 0.004 0.053 0.003 0.070 0.005 0.081 0.004 0.032 0.002
[O ii] 7319 0.022 0.001 0.022 0.001 0.017 0.001
[O ii] B 7325 0.042 0.002 0.041 0.002 0.069 0.004 0.026 0.002
[O ii] 7330 0.020 0.001 0.018 0.001 0.013 0.001
[S iii] 9069 0.125 0.010 0.122 0.009 0.288 0.022 0.133 0.012
[S iii] 9531 0.238 0.019 0.385 0.031 0.281 0.025
Table 2Emission-Line Intensity Ratios Relative to H for Observed Galaxies (2 of 2)
KISSR [O iii [K] [S ii [cm] log(He/H)+12 log(N/H)+12 log(N/O) log(O/H)+12
0148 10120 343 100 10.91 0.05 6.94 0.05 -1.30 0.08 8.23 0.06
0242 12490 136 68.5 25.4 10.96 0.01 6.99 0.02 -0.95 0.02 7.93 0.02
0258 23730 2894 176 38.1 10.98 0.12 6.97 0.04 -0.42 0.02 7.38 0.03
0451 11390 403 48.8 24.7 10.96 0.01 6.91 0.04 -1.20 0.02 8.12 0.04
0475 10920 500 84.8 28.3 10.91 0.01 7.09 0.04 -1.00 0.02 8.10 0.06
0512 10000 49.6 43.5 10.82 0.04 7.43 0.01 -0.46 0.01 7.89 0.01
0590 12020 449 15.7 28.3 10.94 0.04 6.89 0.04 -1.13 0.02 8.02 0.05
0653 10380 335 84.1 28.7 10.95 0.01 7.08 0.03 -1.00 0.02 8.08 0.04
1056 11050 380 66.0 23.6 7.00 0.04 -1.14 0.02 8.14 0.04
1379 10730 516 39.5 28.1 10.96 0.01 6.98 0.05 -1.13 0.02 8.11 0.06
1734 13610 154 114 34.0 10.96 0.01 7.04 0.03 -0.91 0.03 7.95 0.02
1955 10000 41.4 31.4 10.92 0.04 7.18 0.01 -0.89 0.02 8.06 0.01
2117 9751 295 91.0 26.5 10.89 0.04 7.18 0.03 -1.15 0.03 8.33 0.04
2125 10000 100 10.91 0.04 6.67 0.02 -0.67 0.02 7.34 0.03
2132 11590 417 31.5 26.1 10.95 0.02 6.85 0.04 -1.27 0.02 8.12 0.04
Table 3Electron Densities, Temperatures, and Ionic Abundances (1 of 2)
KISSR 12+log(Ne/H) log(Ne/O) 12+log(S/H) log(S/O) 12+log(Ar/H) log(Ar/O)
0148 7.57 0.09 -0.67 0.04 6.81 0.07 -1.42 0.07 5.88 0.06 -2.36 0.08
0242 7.33 0.02q -0.61 0.01 6.35 0.02 -1.59 0.02 5.69 0.02 -2.25 0.01
0258 6.68 0.05 -0.70 0.04 5.75 0.05 -1.63 0.02
0451 7.45 0.05 -0.67 0.01 6.60 0.03 -1.52 0.04 5.59 0.04 -2.53 0.03
0475 7.43 0.07 -0.67 0.01 6.54 0.04 -1.55 0.04 5.14 0.06 -2.95 0.04
0512 7.11 0.04 -0.78 0.04 6.85 0.01 -1.04 0.01 5.67 0.04 -2.23 0.04
0590 7.33 0.06 -0.68 0.02 6.55 0.03 -1.47 0.04 5.69 0.06 -2.33 0.05
0653 7.28 0.05 -0.80 0.02 6.77 0.03 -1.31 0.04 5.71 0.04 -2.37 0.02
1056 7.51 0.05 -0.63 0.02 6.05 0.03 -2.09 0.02
1379 7.43 0.07 -0.68 0.02 6.38 0.04 -1.73 0.04 5.76 0.05 -2.35 0.02
1734 7.32 0.02 -0.63 0.01 6.17 0.02 -1.78 0.02 5.40 0.02 -2.55 0.02
1955 7.59 0.05 -0.48 0.04 6.70 0.03 -1.37 0.03 5.81 0.04 -2.25 0.05
2117 7.59 0.05 -0.73 0.02 6.87 0.03 -1.46 0.04 5.91 0.06 -2.42 0.06
2125 6.64 0.05 -0.70 0.03 6.06 0.02 -1.28 0.02 4.91 0.03 -2.43 0.03
2132 7.48 0.05 -0.64 0.02 6.05 0.03 -2.08 0.02
Table 3Electron Densities, Temperatures, and Ionic Abundances (2 of 2)


Determination of S temperatures requires observation of the near-infrared [S iii]9069,9531 doublet and the temperature-sensitive auroral [S iii]6312 line (Garnett, 1989). These NIR sulfur emission lines become increasingly important to nebular cooling at lower temperatures (Díaz & Pérez-Montero, 2000). While they benefit from their being bright and relatively free from reddening effects, NIR sulfur requires observations that reach to substantially longer wavelengths than are commonly used for nebular abundance work. The 6312 line is generally weaker than 4363, and can suffer blending due to proximity to the [O i] night-sky line at 6300 Å as well as nebular [O i] emission, however spectrograph CCDs are often more efficient at redder wavelengths. In addition, 4363 can be subject to significant reddening effects by the presence of dust that impedes measurement, particularly for metal-rich nebulae in which the line strength is naturally diminished. A  measurement utilizing [S iii]6312 is comparatively less influenced by dust and could remain viable for a wider range of physical conditions. Additionally, while the strength of both the 4363 and 6312 auroral lines weakens significantly at high metallicity, the studies of Bresolin et al. (2004) and Bresolin et al. (2005) have shown that the 6312 line remains a viable alternative for the determination of nebular  to somewhat higher abundances than 4363.
For ten of our fifteen galaxies, the red spectral grating was switched to observe the nebular NIR [S iii] doublet, both of which were detected in six objects. In the other four cases, the 9531 line was redshifted off of the CCD, and its flux value was estimated from the strength of the 9069 line and the theoretical ratio of 9531/9069 = 2.44 (Mendoza & Zeippen, 1982). The weak auroral 6312 line was measured in twelve of our fifteen galaxies, including nine of which possessed NIR spectra. KISSR 0512 alone possessed NIR [S iii] data without recovery of the 6312 auroral line. The [O iii]4363 line was not measured for KISSR 2125, and so a comparison of O and S -method abundances was not possible in this case. This reduces the total count of galaxies for which NIR [S iii], [S iii]6312, and [O iii]4363 were all measured to eight: KISSR 0242, KISSR 0451, KISSR 0475, KISSR 0590, KISSR 0653, KISSR 1379, KISSR 1734, and KISSR 2117.
The ELSA program was used to compute a S temperature for the objects that possessed the necessary lines. Figure 8 compares the  calculated using O and S lines. The solid line indicates the one-to-one equivalence, while the dashed line represents an empirically derived offset of 1184 K to lower  based upon the average difference between the measured O and S electron temperatures. The dotted line represents the fit from Garnett (1992), which uses photoionization models to relate electron temperatures of the different ionization zones, with functional form:

We find that in all but one case, the measured electron temperature calculated using the [S iii] lines is lower than the value determined from the [O iii] lines.
Included in Figure 8 are red points representing H ii region data from the Kennicutt et al. (2003) study of M101 and blue points representing the H ii region data from the Bresolin et al. (2009) study of NGC 300. These data cover a similar range in O  as our work. Bresolin et al. (2009) data points strongly correlate to the equality line, especially at  below 10,000 K, with an average offset of 75 K to higher electron temperatures using S. Data from the Kennicutt et al. (2003) study, however, shows fairly good agreement with our offset above 12,000 K, with an overall average offset of 564 K to lower electron temperatures using S. The work of Berg et al. (2015) revisits the study of H ii regions in M101 by Kennicutt et al. (2003), finding that an adoption of updated atomic data substantially decreases the measured  as determined with S (by roughly 10%, or 1000 K), while the O remains mostly unchanged (see their Figure 3). This result appears to agree with the results from our study, although additional investigation is required.

Figure 8.— Comparison of  calculated using O and S for galaxies where all necessary emission lines are measured. The solid line marks the one-to-one equivalence, while the dashed line represents an empirically derived offset of 1184 K to lower  based upon the average difference between the measured O and S electron temperatures. The dotted line marks the model prediction fit from Garnett (1992). The red and blue points represent H ii regions in M101 and NGC 300 from Kennicutt et al. (2003) and Bresolin et al. (2009), respectively, over a similar O  range.

4.2. Direct-Method Metal Abundances

We use the program ELSA (Johnson et al., 2006) to calculate ionic abundances relative to hydrogen. The input data include the densities and temperatures of each ionization zone, as calculated above, and the emission-line ratios previously presented as Table 3. We calculate the abundance of the following ions with respect to H: O, O, N, S, S, Ne, and Ar.
The total oxygen abundance is assumed to be given by

following the standard practice. Additional ionization states for other elements that are present in the nebula but do not emit in the optical spectrum are accounted for with ionization correction factors (ICFs). We use the prescriptions given by Izotov et al. (1994) for the ICFs for N, Ne, S, and Ar:

Izotov et al. (1994) suggest that the flux for O can be inferred by

In the few cases in which we observe He in the spectra, however, the line is very noisy, and we assume that the amount of O in the nebulae is negligible.
Using the ionic abundances and ionization correction factors given above, we calculate the He/H, N/H, N/O, O/H, Ne/H, Ne/O, S/H, S/O, Ar/H, and Ar/O ratios for each galaxy. Our results are presented in Table 3. Direct-method oxygen abundances for our galaxies are plotted against  ratios in Figure 9. Eleven of the fifteen observed galaxies appear to be regular star-forming galaxies with a range of oxygen abundances between 8.02 12+log(O/H)  8.33, inhabiting the turnaround region of . The remaining four display incongruous characteristics and will be discussed in more detail in §5.2. For two of the regular star-forming galaxies (KISSR 1955 and KISSR 2125), the [O iii]4363 line was not clearly detected. Hence we are unable to provide oxygen-derived direct-method abundances for these two targets.
For all eight galaxies with an estimate of  derived using [S iii] emission lines, a S direct-method oxygen abundance was calculated using ELSA. As expected, the metallicities determined from sulfur-derived  are higher than those determined from oxygen-derived  where the S temperature is lower than that of O (seven of eight cases). We return to the issue of metallicities from [S iii] emission lines in §5.1.2.

Figure 9.—  diagram with direct-method oxygen abundances from this study. The dotted line represents the  calibration from Edmunds & Pagel (1984), while crosses are reference points from McGaugh (1991). Filled black circles represent the nine regular star-forming galaxies with measured . Open green triangles represent possible AGN contaminated targets, while open blue squares are low-metallicity, high-luminosity compact star-forming galaxies. Both of these anomalous target types are discussed in more detail in §5.2.

4.3. Strong Emission-Line (SEL) Method Abundances

The McGaugh (1991) model grid was derived using photoionization model calculations from CLOUDY (Ferland & Truran, 1981) and employs the abundance-sensitive  line ratio alongside the excitation-sensitive line ratio, defined as [O iii]4959,5007/[O ii]3727 and incorporates an upper- and lower-metallicity branch along the axis defined by the  line ratio. Data points representing each of the fifteen galaxies from this sample are shown on Figure 10, differentiated by object type (to be discussed in §5.2.). Typical abundance uncertainties of the McGaugh model grid method are 0.05 for the low-metallicity branch and 0.10 for the high-metallicity branch, dominated by the uncertainties in the model stellar radiation fields (McGaugh, 1991). The nine regular star-forming galaxies with measured  from this study fall between the upper- and lower-metallicity branches of the grid. Abundances for these systems were calculated to be in the range of 8.07 12+log(O/H)  8.74, and are presented in Table 4. The two low-metallicity, high-luminosity compact star-forming galaxies appear to reside within the turnaround region, albeit with noticeably higher excitation. Conversely, the two galaxies for which 4363 emission was undetected exhibit comparatively low-excitation.
The coarse abundance method (Melbourne & Salzer, 2002; Salzer et al., 2005b) incorporates the [O iii]5007/H and [N ii]6583/H line ratios calibrated to a combination of  abundances and previously published -like relations. These two line ratios are observed in the optical spectrum of most local ELGs, typically have good signal-to-noise ratios, and due to the relative proximity of the individual lines, are fairly insensitive to reddening corrections. A spectral activity diagnostic diagram shows that in the low-metallicity regime, abundance increases with nitrogen line strength as a smooth, single-valued function. Conversely, for the high-metallicity regime, abundance is anti-correlated with increasing oxygen line strength. Typical uncertainties of the coarse abundance method are 0.15-0.20 (Melbourne & Salzer, 2002; Salzer et al., 2005b). While not exceptionally precise, coarse abundances can be quite useful for large-scale statistical work. For the nine regular star-forming galaxies observed for this study, the computed coarse-method metallicities range from 8.29 12+log(O/H)  8.44, and are listed in Table 4.

Figure 10.— McGaugh model abundance grid for a 60 IMF upper mass cutoff (McGaugh, 1991) including the fifteen galaxies from our sample. The filled black circles represent the nine regular star-forming galaxies for which [O iii]4363 was measured. Green open triangles have potential AGN spectral contamination. Blue open squares are low-metallicity, high-luminosity compact objects. Red open hexagons are regular star-forming galaxies that lack reliable  measurement. These will be discussed further in §5.2. While the  ratio on the x-axis is sensitive to abundance (marked by vertical lines on the grid ranging from 6.9 to 9.3), the ratio on the y-axis is sensitive to levels of excitation (marked by horizontal lines on the grid ranging from 0.0001 to 0.1). The black grid represents photoionization models calculated from CLOUDY (Ferland & Truran, 1981). Errors in the line ratio measurements are comparable in size to the plot symbols.
KISSR O  [K] S  [K] O S McGaugh Coarse Type
0148 10120 343 8.23 0.06 8.07 0.10 8.29 0.20 SF
0242 12490 136 11670 242 7.93 0.02 8.02 0.02 8.16 0.10 8.21 0.20 low-
0258 23730 2894 7.38 0.03 8.76 0.10 8.67 0.20 AGN
0451 11390 403 9485 410 8.12 0.04 8.35 0.04 8.26 0.10 8.31 0.20 SF
0475 10920 500 10150 773 8.10 0.06 8.18 0.07 8.64 0.10 8.44 0.20 SF
0512 14370 627 7.47 0.05 8.94 0.10 8.71 0.20 AGN
0590 12020 449 10500 498 8.02 0.05 8.18 0.04 8.12 0.10 8.33 0.20 SF
0653 10380 335 7952 285 8.08 0.04 8.43 0.04 8.74 0.10 8.43 0.20 SF
1056 11050 335 8.14 0.04 8.34 0.10 8.31 0.20 SF
1379 10730 516 11340 920 8.11 0.06 8.05 0.05 8.33 0.10 8.34 0.20 SF
1734 13610 154 12870 406 7.95 0.02 8.01 0.02 8.07 0.10 8.21 0.20 low-
1955 8.79 0.10 8.64 0.20 no
2117 9751 295 7851 278 8.33 0.04 8.55 0.04 8.36 0.10 8.38 0.20 SF
2125 7819 924 8.79 0.10 8.66 0.20 no
2132 11590 417 8.12 0.04 8.35 0.10 8.31 0.20 SF
Table 4Electron Temperature, Oxygen Abundance Comparison, and Spectral Activity Type

5. Discussion

This study was undertaken in effort to constrain the upper-metallicity branch of the  relation with direct-method oxygen abundances. It was understood from the start that this would be a challenging task. Galaxies for which the temperature-sensitive [O iii]4363 auroral line is strong enough to be readily measured correspond to those of highest excitation and lowest overall metallicity. The two primary selection criteria used to generate the targets for this project – strong, high-EW emission lines and high coarse abundances – are fundamentally at odds with one another: Our goal was to populate the lower portion of the upper-metallicity branch, above the turnaround region, with robustly obtained abundances and provide observational constraint for the  relation. We hoped to achieve this goal by (1) utilizing the large light-collecting area of the Keck telescope and the high sensitivity of LRIS, and (2) to pre-select optimal candidates for observation using the best targets from the deep KISS catalogs.
Unfortunately, we were unsuccessful in achieving the measurements of high enough metal abundances to populate the upper-metallicity branch of the  relation. We were, however, successful in deriving reliable direct-method -abundances for the majority of our targets. Furthermore, we have compared our results to two SEL abundance methods in order to explore the apparent discrepancy in metallicity results, and have explored a non-traditional avenue of utilizing NIR [S iii] measurements to obtain electron temperatures for our abundance analysis. Finally, our study has revealed disparity in the homology of an outwardly-appearing homogeneous set of star-forming galaxies.

5.1. Metallicity Comparison

Several studies have indicated the existence of discrepancies between oxygen abundances computed via photoionization models and the direct-method. In such cases, the auroral lines may not produce a reliable  measurement (Stasińska, 2005). There are strong exponential dependences on the strength of the weak, temperature-sensitive auroral emission lines. Small-scale electron temperature fluctuations are a likely contributor toward these discrepancies (e.g., Peimbert, 1967; Peimbert & Costero, 1969; Peimbert, 2003; Peimbert et al., 2007). In such a case, an unresolved ELG spectrum encompasses a large volume of ionized gas that includes many different individual regions of varying physical properties. The measured spectrum would therefore comprise a global average that is not representative of any specific region of the nebula. Detection of [O iii]4363 emission would therefore be weighted toward areas of highest excitation and may not necessarily coincide spatially with other regions of weaker [O iii] emission, where the temperature of the electron gas is lower. Because these spatially distinct regions of the nebula are measured as part of a single spectrum, their overall physical properties become blended together.
Additionally, Garnett (1992) explains that thermal gradients can be responsible for an underestimation of derived abundances from measured electron temperatures. Observations and photoionization models show that nebulae are not necessarily isothermal and the average  can vary for different ionic species (Garnett, 1992). In such cases the measured  is biased high, influencing the estimated metallicities through the emitting regions. Different ionic species’ ionization potentials and radiation hardness gradients with optical depth lead to inhomogeneous nebular structure. For instance, Stasińska (1980) has demonstrated an increase in electron temperature with increasing radius from the ionizing source for metal-rich nebulae.
The study of Andrews & Martini (2013) suggest that O zone temperatures, often obtained using model-derived relationships involving the measured O temperatures, tend to be overestimated, leading to underestimates in the O abundance. As the overall abundance is assumed as the sum of singly- and doubly-ionized oxygen, the undervalued O abundance results in a net deficit of the calculated direct-method metallicity. Recently, discrepancies in nebular abundances derived from electron temperature measurements have been explored with the introduction of the “-distribution” of electron energies (see Nicholls et al., 2012, Nicholls et al., 2013, and Dopita et al., 2013). An adoption of model H ii regions exhibiting high-energy (suprathermal) excess in electron energies appears to resolve some temperature and metallicity discrepancies as compared to an equilibrium described by a Maxwell-Boltzmann distribution. The study of Mendoza & Bautista (2014), however, found that the reconciliation of  measurements for higher-excitation zones (such as O) using -distributions is not significant.
Some studies (notably Andrews & Martini, 2013) utilize spectra stacked in bins of mass, averaged to achieve the needed signal-to-noise to detect 4363. While such a method does away with certain selection effect biases described above, the disparity between direct- and SEL-method abundance estimations remains. Spectral stacking, however, may introduce other uncertainties that can skew measurements. In particular, for a given galaxy mass bin, the inclusion of particularly strong [O iii] emitters can over-represent 4363 emission strength relative to the averaged 4959,5007, biasing temperatures high and thus total oxygen abundances low.
With increasing recent evidence that the standard procedure in computing oxygen abundances may be less robust than traditionally assumed, it becomes useful to explore the utility of alternative methods. A comparison of SEL method metallicities, described in §5.1.1., as well as a comparison of -method metallicities derived from O and S, described in §5.1.2, is shown as Figure 11. The results of our abundance analyses are summarized in Table 4. We find that oxygen abundances calculated using empirical relations are systematically higher for our targets as compared to direct-methods, and furthermore that metallicities derived from S  tend to be higher than as calculated using O , but are still lower than SEL-method abundances.

Figure 11.— Metallicity comparison for all 15 galaxies from this study utilizing both [O iii] and [S iii] for -method abundances (red points and black points, respectively) and the SEL empirical McGaugh model grid (blue points) and coarse abundance method (green points). For the eight galaxies with calculated S  (including six regular star-forming galaxies and two low-metallicity, high-luminosity compact galaxies), the sulfur-derived abundances are higher than oxygen-derived abundances for all but one case. The [S iii] lines used to derive the electron temperature represent a different ionization zone in the nebulae than the [O iii] lines, and when studied in conjunction with traditional -methods provides a reasonable range for the global physical characteristics. Furthermore, assuming  calculated from the zone of doubly-ionized sulfur reduces the discrepancy between direct- and SEL methods, though we stress that this technique should be considered a thought experiment in deriving a lower-limit to the overall metallicity.

Comparison Between Direct and SEL Abundances

Offsets may arise when estimating metallicities for the same object while employing direct- and SEL-methods, where generally -method abundances tend to be biased low when compared to SEL calibrations. Kennicutt et al. (2003) found that their direct-method metallicities of individual H ii regions in M101 were systematically low by 0.2–0.5 dex compared to the most widely used SEL methods. Bresolin et al. (2009) came to a similar conclusion based on a study of H ii regions in NGC 300. The study of Kewley & Ellison (2008) found an even greater discrepancy in star-forming galaxies from the Sloan Digital Sky Survey (SDSS; Abazajian et al., 2005) of up to 0.7 dex. Andrews & Martini (2013) argue that the most likely cause of the offset between theoretical- and direct-method metallicities lies in the breakdown in one or more of the assumptions in the physics of H ii regions in the stellar synthesis or photoionization models. Stellar population synthesis models typically assume that the ionizing source is a zero age main sequence starburst, which affects line fluxes as a function of age (Berg et al., 2011). Additionally, photoionization models typically invoke simplifying assumptions or treatments that reduce their applicability to real-life nebulae. These include the patently false assumption of spherical or plane-parallel H ii region geometry, poorly constrained treatment of metal depletion onto dust grains (e.g., Zurita & Bresolin, 2012 and references therein), and the assumption of smooth (rather than clumpy) gas and dust density distribution (Kewley & Ellison, 2008).
In order to quantify the amount of offset between direct- and SEL-methods and the effect this has on the calibration of the  relation, a comparison of three abundance estimates for the nine regular star-forming galaxies in this study with measured [O iii]4363 are plotted with their respective  ratios as Figure 12. Using the McGaugh model grid, the estimated oxygen abundance values appear to be systematically high by an average of 0.22 0.25 dex as compared to those calculated using the -method. Metallicities estimated using the coarse abundance method also tend to be high relative to the  abundances with an average offset of 0.21 0.11 dex. There is zero average offset between the SEL McGaugh and coarse abundance methods for the nine regular star-forming galaxies, but with a standard deviation of 0.17 dex. With a few exceptions, the SEL method abundances are greater than 12+log(O/H) = 8.2 and fall on the upper-branch of the  metallicity diagram, while direct-method abundances calculated from the same data remain in the turnaround region.
The two primary sample selection criteria for this study specify both high-EW [O iii] lines and potentially high oxygen abundances as estimated using the SEL coarse abundance method. Galaxies selected from KISS that satisfy the necessity for strong [O iii] emission will be high-excitation and have higher electron gas temperatures, indicative of inefficient cooling and thus must be relatively metal-poor. The satisfaction of the selection criterion mandating high potential metallicities for these galaxies, as estimated by the coarse method, therefore introduces a selection bias. We believe that the viable targets initially picked for this study are specifically those whose coarse-method abundance estimations occupy the high end of the scatter of the metallicity distribution. In the majority of cases, as seen in Figure 12, the different metallicity estimates agree within their respective errors, with SEL method values occupying the high side of the -method scatter. Metallicity estimates for a few of the regular star-forming galaxies, namely KISSR 2117, agree quite well between the direct- and SEL methods. This provides evidence that any offset present is one that is not entirely systematic.

Figure 12.—  diagram using oxygen -abundances, McGaugh model grid abundances, and coarse method abundances for our nine regular star-forming galaxies for which [O iii]4363 was measured. The dotted line represents the  calibration from Edmunds & Pagel (1984), while the crosses are reference points from McGaugh (1991). The SEL-method abundances (McGaugh and coarse) appear systematically high compared to the direct -method. Note that while some -method values are widely variant from SEL method values, others agree quite well, implying that this offset is not entirely a systematic effect, but can include bias in sample selection. Errors in line ratio measurements (x-axis) are comparable in size to the dots. Errors in SEL abundance values (y-axis) are the high-end estimates for either method; 0.10 for the McGaugh grid and 0.20 for the coarse method.

Comparison Between O and S  Abundances

We have derived the electron temperature for several of the ELGs from this study using [S iii] lines (§4.1.2). For all but one of these objects, the sulfur  was lower than that of oxygen, however the two temperatures agree within the errors. All eight of these galaxies are considered star-forming galaxies (e.g., none are thought to possess spectral contamination by AGN), though two fall into the classification of high-luminosity, low-metallicity compact objects (see §5.2.3).
It is possible to compare the oxygen abundance 12+log(O/H) as calculated using electron temperatures derived from both oxygen and sulfur. When the temperature-sensitive auroral [S iii]6312 line is measured in ELG spectra, the S  is usually used only to compute the sulfur abundance. The S electron temperature is representative of a slightly lower-energy ionization zone within the nebulae as compared to O, and as such the  of the S zones are not necessarily appropriate substitutes for those of the O zones. In our study, however, we have opted to explore a non-standard method of deriving oxygen abundances by employing the electron temperature estimated from the S zone. Utilizing sulfur as the nebular thermometer for abundance calculations may reduce certain biases in oxygen measurement (see §4.1.2.).
The adoption of the S  for the temperature of the high-ionization O zone gives a range in physically viable values for the calculation of metallicities. Implementing abundance estimates utilizing sulfur-derived electron temperatures should not be taken as the solution to oxygen measurement biases, but rather as a complementary method. Abundances using S  are correspondingly higher compared to the traditional method for the same galaxies. These values fall between the O -method and SEL abundance estimates, potentially offering some reconciliation of the disparity discussed in §4.3.1. We note that we chose not to include the S -substituted oxygen abundance for KISSR 2125 in our analysis because the resulting error bars on the estimated value, 12+log(O/H) = 8.06 0.41, are unreasonably large (more than double the error of the coarse method abundance).
Regular star-forming galaxies for which we obtained NIR spectral coverage have direct-method metallicities estimated to be higher by an average of 0.16 0.14 dex. Should the S ionization zone temperature be more representative of the nebula as a whole, then the metallicity derived using [O iii] emission lines will underestimate the global oxygen abundance. This small, systematic offset between O and S -method oxygen abundances reduces the average overall discrepancy between direct- and SEL-method metallicities seen in Figures 11 and 12. Adopting the oxygen abundances of galaxies estimated using sulfur temperatures reduces the average offset between direct-method and SEL-method metallicities from 0.26 0.22 to 0.11 0.22 dex (for the McGaugh model grid) and from 0.25 0.10 to 0.11 0.16 dex (for the coarse method) for our targets that possess NIR spectral coverage. Because their origin is from a lower-energy transition, S zone-derived electron temperatures could therefore be less subject to the effects of temperature fluctuations than O on the derivation of direct-method metallicities.
Esteban et al. (2014) show that within their study, giant H ii regions in dwarf star-forming galaxies exhibit the highest temperature fluctuation parameters. As the measured discrepancy between direct- and SEL-method abundances is thought to be at least in part a consequence of such fluctuations, the adoption of sulfur-derived electron temperatures in direct-method analyses could produce intermediate metallicity estimates useful in reconciling the standard techniques. While some of the discrepancy between direct- and SEL-method abundance calculations appears to be accounted for by utilizing S rather than O electron temperatures, we do not prescribe that such a surrogate is the best solution for deriving oxygen abundances from ELG spectra. Rather, this method can be used as a constraint to direct-method metallicities by providing a reasonable range of electron temperature values for use in -method abundance determinations.
For a galaxy or H ii region in which an O  cannot be derived, measurements of the [S iii] emission lines may enable an alternative direct-method metallicity. For higher-metallicity nebulae, [S iii]6312 will be the preferably-detected auroral line as the strength of [O iii]4363 diminishes (Croxall et al., 2015). The S temperature could then be used to estimate an O temperature, from which an oxygen abundance could be found. Further observations and analysis of targets with NIR spectral coverage and S  determinations will provide additional data points from which to refine this offset.

5.2. Inhomogeneity of Sample Galaxies

When we selected our objects for inclusion in this study, we tacitly assumed that all of the KISS galaxies in our lists were typical star-forming galaxies. After obtaining and analyzing our spectra, however, it became clear that not all of our sources could be categorized as simple or “regular” star-forming systems. That is, the targets of this study appear to compose an inhomogeneous set of objects. We have segregated our galaxies into three distinct categories: Regular star-forming galaxies (further distinguished into those with and without [O iii]4363 detection), those with suspected AGN spectral contamination, and low-metallicity, high-luminosity compact objects. We discuss each of these subgroups below. Figure 13 presents a line-ratio diagnostic diagram plotting the newly measured ratios and using different symbols to indicate the four classes of objects in our sample.

Figure 13.— Diagnostic diagram of galaxies we observed with Keck segregated by type, including regular star-forming galaxies for which -abundances have been derived in this study (black filled circles), those with no believable detections of the temperature-sensitive [O iii]4363 Å line and thus lower precision abundance measurements (red open hexagons), those with suspected contamination by a non-thermal photoionizing source, e.g., AGN (green open triangles), and possible nearby examples of low-metallicity, high-luminosity galaxies (blue open squares). The line ratios plotted here are those found using Keck LRIS from this study.

Regular Star-Forming Galaxies

Eleven of the 15 targets were determined to be “regular” star-forming galaxies, for which -method abundances could be calculated for nine (represented by filled circles in Figure 13). Of the latter group, the calculated metallicities all inhabit a range of 8.0 12+log(O/H)  8.3, and therefore were not successful in extending above the  turnaround regime and populating the upper-metallicity branch, as seen in Figure 12. These -method abundances are on average 0.21 0.12 dex lower than the initial coarse abundance estimates from the original KISS follow-up spectroscopy. For some galaxies, such as KISSR 2117, the direct-method metallicities agree quite well with the initial coarse method estimates, with 12+log(O/H) = 8.33 and 8.36, respectively. Others, such as KISSR 0653, are widely discrepant, with 12+log(O/H) = 8.08 and 8.46, respectively. Despite the fact that we were unsuccessful in determining direct abundances above the turnaround region, we have provided robust measurements of the oxygen abundance for several KISS galaxies. These will be added to the overall sample of KISS galaxies possessing direct abundances and will be useful, for example, in future refinements to the calibration of the SEL coarse abundance method.
The two remaining star-forming galaxies, KISSR 1955 and KISSR 2125, lacked any discernible detection of the [O iii]4363 temperature-sensitive line (represented by open hexagons in Figure 13). Therefore we were unable to calculate a -method oxygen abundance. While appearing to be regular star-forming galaxies in all other respects, the initial abundance estimate for these two targets were amongst the highest of the sample (12+log(O/H) = 8.60 and 8.65, respectively), and so it should not be surprising that they represent the technical limit of our ability to recover metallicities via the direct-method. Lacking a measurement of a direct-method oxygen abundance, these two galaxies were eliminated from further analysis in this sample.

AGN Contamination

Two of the 15 targets, KISSR 0258 and KISSR 0512, appear to possess emission-line ratios indicative of spectral contamination by a non-thermal ionizing source (represented by open triangles in Figure 13). Both of these galaxies appear near or past the Kauffmann et al. (2003) boundary for AGN activity, despite initially appearing to be regular star-forming galaxies (Figure 2) by the original KISS spectral data. A direct comparison with Figure 2 furthermore indicates that while the initial position of KISSR 0512 is nearby its location on Figure 13, there appears to be a major discrepancy in the position of KISSR 0258. While initially believed to be a regular star-forming galaxy, this target has clearly shifted beyond the demarcation line segregating AGN. A comparison with the initial KISS follow-up spectroscopy implies that the original pointings of this object placed the slit over a different region of the galaxy than was observed by Keck, e.g., an outlying knot of star formation compared to nuclear emission. Based on the location of this target on the diagnostic diagram of Figure 13, it should not have been included as a potential candidate for observation by this study. The location of KISSR 0512 on the diagnostic diagrams depicting original KISS data and our values from Keck (Figures 2 and 13, respectively) resides near the boundary separating star-forming galaxies and AGN and therefore was already suspect.
As a consequence of their placement on the diagnostic diagram, we feel that these two galaxies are likely contaminated by AGN activity, and hence cannot be counted amongst the regular star-forming galaxies. In particular, their detected [O iii]4363 lines appear to be far too strong for a purely star-forming emission-line source. These AGN-contaminated targets have major offsets due to the non-thermal nature of the photoionization sources powering the nebulae. Because of the high-energy photons emitted by matter accreting onto the galaxy’s central supermassive black hole, there is a comparatively harder ionizing radiation field that disproportionally inflates the emission of [O iii]4363 in relation to stronger lines such as [O iii]4959,5007. As a result, the temperature-diagnostic emission-line ratio produces a measured electron temperature that is skewed high.
Consequently, these two galaxies present the most disparate values of metallicity between the direct- and SEL methods. -method oxygen abundances for KISSR 0258 and KISSR 0512 were calculated to have values of 12+log(O/H) = 7.38 and 7.47, respectively, as seen in Table 4. These values are extremely low and are therefore probably spurious. Especially considering that these two galaxies are the most luminous in our sample ( = -20.15 and -20.46, respectively), their low measured metallicities are very surprising when placed on the luminosity-metallicity () relation. The disparity of the  abundances compared to their coarse method estimates (12+log(O/H) = 8.41 and 8.86, respectively) is quite drastic. For such cases, the physical assumptions involved with the star-formation models that are integral to the development and calibration of SEL method relations are no longer directly applicable.

Low-Metallicity, High-Luminosity Compact ELGs

The remaining two targets, KISSR 0242 and KISSR 1734, both exhibit very strong [O iii]4363 emission. Both were found to be highly underabundant (12+log(O/H) = 7.93 and 7.95, respectively), but also exhibit very compact morphologies and high luminosities ( = -19.46 and -18.97, respectively). This combination contrasts with the well-established relation, whereby the most luminous galaxies possess the greatest degree of enrichment. We note that the direct-method metallicities for these objects relied on  measurements utilizing very high signal-to-noise spectra, owing to the high-excitation nature of these galaxies and the exceptionally strong emission of [O iii]4363 (see Figures 5 and 7). As a result, the oxygen abundances calculated for these galaxies are very robust.
Galaxies of this nature have been previously detected in the KISS sample (Salzer et al., 2009), however those objects were all detected as strong [O iii]-lined sources redshifted such that they were detected in the KISS H filter. The redshift range of these rare [O iii]-detected KISS galaxies is = 0.29 – 0.42, corresponding to a lookback time of roughly 3–4 Gyr (assuming a cosmology of = 70 km s, = 0.27, and = 0.73) compared to a redshift range of = 0.04 – 0.07 for our two targets from the current study. To match what is typically observed in the local Universe today, these higher-redshift [O iii]-detected KISS galaxies would have necessarily been required to experience an immense degree of metal enrichment (a factor of 10) or to have undergone a massive luminosity decrease (a factor of 28). Either of these scenarios show that such luminous, metal-poor galaxies represent extreme phases of evolution. As shown in Salzer et al. (2009), these galaxies are very rare, representing less than 0.02% of all galaxies in this luminosity range ( from -19.5 to -22.0) at these lookback times.
Figure 14 illustrates the location of these higher-redshift objects in relation to the set of regular star-forming KISS galaxies on an diagram (metallicities estimated using the coarse abundance method), and includes these two possible local analogues. Solid squares represent initial abundance estimates, demonstrating that these objects are metal-poor and/or over-luminous relative to typical KISS galaxies. The degree by which they are discrepant from the rest of the KISS sample, however, is fairly moderate: Their placement on the outside edge of the main cluster of data points puts them roughly half way between the relation fit and the offset fit to higher-redshift KISS star-forming galaxies from Salzer et al. (2009). Utilizing coarse abundance estimates, then, these two objects were not expected to be extreme outliers. The open squares connected via dotted lines represent -abundances measured from this study, indicating that, modulo uncertainties in the coarse abundance method, KISSR 1734 is a 2 outlier and KISSR 0242 is nearly a 3 outlier relative to the mean relation for the full KISS sample.
Other recent studies, including Hoyos et al. (2005), Kakazu et al. (2007), and the “Green Pea” galaxies (GPs) introduced by Cardamone et al. (2009), have also reported discovery of luminous, underabundant galaxies in various redshift ranges. A study by Amorín et al. (2010) shows that GPs are found to follow a mass-metallicity () relation parallel to regular star-forming galaxies, but one that is offset by 0.3 dex to lower metallicities. These objects are therefore similar but less extreme in their divergence from typical galaxies than those from Salzer et al. (2009). Potentially, GPs and their equivalents are a subset of star-forming luminous compact galaxies (LCGs). This idea demonstrates that such objects exist at any redshift regime, but present differing observable properties inducing classification as a variety of objects (Izotov et al., 2011). Cardamone et al. (2009) provides evidence that GPs/LCGs are nearby analogues of starbursting UV-luminous galaxies (UVLGs) such as Lyman-break galaxies (LBGs) and Ly emitting galaxies (LAEs). This inference is verified by UV spectra of some GPs showing Ly emission (Jaskot & Oey, 2014). Star-formation taking place within UVLGs in a redshift range 1.9 3.4 represents a relatively large portion of the current total stellar mass (Reddy & Steidel, 2009). GPs/LCGs therefore represent relatively nearby laboratories with which to study the mechanisms at work consistent with the bulk of star-formation in the Universe.
Another related class of galaxies are the luminous blue compact galaxies (LBCGs) (Garland et al., 2004; Werk et al., 2004). These objects possess luminosities equivalent to or greater than the Milky Way, and yet are substantially more compact. Though common at intermediate and high redshifts, LBCGs are locally rare. With activity dominated by intense star formation, these blue galaxies contributed significantly to the overall star-formation rate density of the Universe. The study of Werk et al. (2004) examined 16 LBCG candidates selected from KISS, including KISSR 0242. They found that LBCGs have higher H luminosities than average, indicative of strong star-formation activity, and lower metal abundances than expected based on the relation for KISS galaxies.
The import of finding low-metallicity, high-luminosity star-forming galaxies in the very local Universe is one that speaks to the ongoing understanding of how these systems evolve. Objects that are divergent low from the relation, like those at = 0.29 – 0.42 discussed in Salzer et al. (2009), are implied to have experienced significant chemical enrichment over the subsequent 4 Gyr in order to evolve into more typical star-forming systems at = 0. These galaxies are thought to be amongst the final massive systems in the Universe that are in an early phase of gravitational collapse and initiation of star formation. But finding local representatives of these objects ( = 0.038 and 0.066 for KISSR 0242 and KISSR 1734, respectively) indicates that such an intense initial evolutionary stage may have been delayed further until even very recently.
That two galaxies of this type were found from a set of only 15 is not an indication that these are common. On the contrary, of initial the list of potential KISS targets, KISSR 0242 and KISSR 1734 possessed the highest [O iii]5007 equivalent widths. Furthermore, their estimated oxygen abundances and absolute magnitudes were also amongst the lowest of all listed galaxies. Because the selection criterion of this study mandated strong [O iii] emission, it should be no surprise that low-metallicity, high-luminosity targets would be preferentially selected from the overall sample for observation. The fact that a handful of additional local galaxies from our target list possessing similar qualities have not yet been observed leaves the exciting prospect that more systems in an equivalent evolutionary stage are available for further study. Future work promises to shed new light on this type of system.

Figure 14.— Luminosity-metallicity (-) relation diagram for 1400+ low- KISS regular star-forming galaxies (small dots) adopted from Salzer et al. (2009) with 13 [O iii]-detected star-forming galaxies (red triangles). Solid blue squares are the two possible local analogues of low-metallicity, high-luminosity compact galaxies discovered by this study. Metallicities are estimated using the coarse abundance method. Open blue squares represent -method abundances calculated in this study, connected to initial metallicity estimates with a dotted line. These abundances are more than 2 outliers to the norm. The solid line is a linear fit to the low- galaxies, while the lower dashed line has the same slope but fits the higher- galaxies with an offset of -1.1 dex.

6. Summary and Conclusions

We have presented spectral analysis of 15 ELGs observed with Keck LRIS and selected from KISS to be both strong [O iii] emitters and metal-rich in effort to find -method abundances to populate and constrain the upper-metallicity branch of . We were successful in recovering -method abundances for 13 of the 15 galaxies observed. An examination of the spectra and their location on a diagnostic diagram reveals that the spectral activity of our targets are inhomogeneous, including two objects with probable AGN contamination and two compact luminous galaxies with low abundance. Oxygen abundances calculated for nine regular star-forming galaxies of this study were found to be constrained to the turnaround region of the  diagram, corresponding to 8.0 12+log(O/H)  8.3.
Of the three lists of initial selection criteria, galaxies from the third list (modest [O iii] EW and highest estimated abundance) were bereft of measurable 4363, and therefore of a -method abundance measurement. The second list (moderate [O iii] EW and moderate estimated abundance) produced galaxies with direct-method metallicities firmly in the  turnaround region, a useful determination as a robust value of oxygen abundance but falling short of the goals for this study. The first list (highest [O iii] EW and moderate estimated abundance), however, selected galaxies that we either identified as low-metallicity, high-luminosity compact star-forming galaxies (and are therefore interesting as a separate case) or as regular star-forming galaxies, including some with fairly high measured oxygen abundances.
The highest metallicity we found overall was for a galaxy from the first list, KISSR 2117, for which we calculated an oxygen abundance of 12+log(O/H) = 8.33 0.04. While not enriched enough to populate the upper-metallicity branch of , the abundance calculated for this galaxy is approaching the necessary limit to do so. To fulfill our original project goal, we suggest that searches for targets of future observations should focus on the strength of [O iii]5007, in order to recover 4363 for determination of .
We have also reaffirmed the results of previous studies that in comparison to the -method, the McGaugh model grid and coarse abundance empirical strong-line methods result in oxygen abundance values that are systematically higher (by 0.22 0.25 dex and 0.21 0.11 dex, respectively). While typically thought to be a result of certain well-known observational biases such as temperature fluctuations and gradients, in this case such disparity is likely enhanced by a bias in the selection criteria employed for our study. Preliminary metallicity estimates based on the coarse abundance method, for objects with strong [O iii] emission, will be preferentially located in the high-end scatter of the distribution, and as such the discrepancies between our direct- and SEL-method abundances may be significantly enhanced.
A comparison of electron temperature calculations utilizing emission lines of O and S, representing overlapping ionization zones in the nebulae, yield systematic differences in direct-method oxygen abundance values, a result consistent with previous studies. For the eight cases in which all required lines were measured in the galaxy spectra, all but one showed overall lower temperatures when measured with S lines, corresponding to an increase in measured metallicity of an average of 0.14 dex. While the adoption of a  derived from sulfur for abundance calculations cannot be taken as the solution to oxygen measurement biases, a comparison with standard procedures (e.g.,  calculated from oxygen) provides a complementary value. The resulting lower measured electron temperatures correspond to higher estimated abundances, potentially offering some reconciliation of the disparities between direct- and SEL-method metallicities. Whether this is because sulfur-derived  is a better representation of the true nebular temperature or not is unknown. Furthermore, with a standardized temperature offset, it is possible to estimate an electron temperature for a zone of doubly-ionized oxygen from the measurement of emission lines of doubly-ionized sulfur. This is potentially very useful for observations of targets where reddening is a major impediment or blue sensitivity is unavailable.
Based on analysis of the spectral activity type derived from the diagnostic diagram (Figure 13), we have concluded that our targets comprise an inhomogeneous set of galaxies. The majority of objects (11 of 15) are regular star-forming galaxies, albeit for two (KISSR 1955 and KISSR 2125) we were unable to measure the temperature-sensitive [O iii]4363 auroral line and thus recover direct-method abundances. Two of our galaxies (KISSR 0258 and KISSR 0512) demonstrate elevated line-ratios indicative of non-thermal photoionizing sources, and are thus likely to be contaminated by AGN activity. Finally, two of our galaxies (KISSR 0242 and KISSR 1734) appear to be under-abundant for their luminosities, falling below the relation defined by the KISS sample. They appear to represent very local analogues of the low-metallicity, high-luminosity compact objects of Salzer et al. (2009) and the “Green Peas” found in the SDSS (Cardamone et al., 2009), with implication that they are exhibiting an extreme phase of galaxy evolution not yet heavily studied in the nearby Universe. That several additional local galaxies from our initial target list were left unobserved leaves the exciting prospect that more such systems are available for study.

We gratefully acknowledge financial support for the KISS project from an NSF Presidential Faculty Fellow to J. J. S. (NSF AST 95-53020), as well as additional support for the ongoing follow-up spectroscopy campaign (NSF AST 00-71114). The quality of this manuscript was greatly improved by the many useful suggestions made by an anonymous referee. We thank the many KISS team members who have participated in the spectroscopic follow-up observations. Finally, we wish to thank the support staff of the W. M. Keck Observatory for their excellent assistance in obtaining the spectroscopic observations that made this work possible.


  1. IRAF is the Image Reduction and Analysis Facility distributed by the National Optical Astronomy Observatory, which is operated by the Association of Universities for Research in Astronomy (AURA) under cooperative agreement with the National Science Foundation (NSF).


  1. Abazajian, K., Adelman-McCarthy, J. K., Agüeros, M. A., et al. 2005, AJ, 129, 1755
  2. Amorín, R. O., Pérez-Montero, E., & Vílchez, J. M. 2010, ApJ, 715, L128
  3. Andrews, B. H., & Martini, P. 2013, ApJ, 765, 140
  4. Baldwin, J. A., Phillips, M. M., & Terlevich, R. 1981, PASP, 93, 5
  5. Berg, D. A., Croxall, K. V., Skillman, E. D., et al. 2015, ArXiv e-prints
  6. Berg, D. A., Skillman, E. D., & Marble, A. R. 2011, ApJ, 738, 2
  7. Bohlin, R. C., Colina, L., & Finley, D. S. 1995, AJ, 110, 1316
  8. Bresolin, F., Garnett, D. R., & Kennicutt, Jr., R. C. 2004, ApJ, 615, 228
  9. Bresolin, F., Gieren, W., Kudritzki, R.-P., et al. 2009, ApJ, 700, 309
  10. Bresolin, F., Schaerer, D., González Delgado, R. M., & Stasińska, G. 2005, A&A, 441, 981
  11. Cardamone, C., Schawinski, K., Sarzi, M., et al. 2009, MNRAS, 399, 1191
  12. Contini, T., Treyer, M. A., Sullivan, M., & Ellis, R. S. 2002, MNRAS, 330, 75
  13. Croxall, K. V., Pogge, R. W., Berg, D., Skillman, E. D., & Moustakas, J. 2015, ArXiv e-prints
  14. Díaz, A. I., & Pérez-Montero, E. 2000, MNRAS, 312, 130
  15. Dopita, M. A., & Evans, I. N. 1986, ApJ, 307, 431
  16. Dopita, M. A., Sutherland, R. S., Nicholls, D. C., Kewley, L. J., & Vogt, F. P. A. 2013, ApJS, 208, 10
  17. Edmunds, M. G., & Pagel, B. E. J. 1984, MNRAS, 211, 507
  18. Esteban, C., García-Rojas, J., Carigi, L., et al. 2014, MNRAS, 443, 624
  19. Ferland, G. J., & Truran, J. W. 1981, ApJ, 244, 1022
  20. Garland, C. A., Pisano, D. J., Williams, J. P., Guzmán, R., & Castander, F. J. 2004, ApJ, 615, 689
  21. Garnett, D. R. 1989, ApJ, 345, 282
  22. —. 1992, AJ, 103, 1330
  23. Gronwall, C., Jangren, A., Salzer, J. J., Werk, J. K., & Ciardullo, R. 2004a, AJ, 128, 644
  24. Gronwall, C., Salzer, J. J., Sarajedini, V. L., et al. 2004b, AJ, 127, 1943
  25. Hamuy, M., Walker, A. R., Suntzeff, N. B., et al. 1992, PASP, 104, 533
  26. Hoyos, C., Koo, D. C., Phillips, A. C., Willmer, C. N. A., & Guhathakurta, P. 2005, ApJ, 635, L21
  27. Izotov, Y. I., Guseva, N. G., & Thuan, T. X. 2011, ApJ, 728, 161
  28. Izotov, Y. I., Thuan, T. X., & Lipovetsky, V. A. 1994, ApJ, 435, 647
  29. Jangren, A., Salzer, J. J., Sarajedini, V. L., et al. 2005, AJ, 130, 2571
  30. Jaskot, A. E., & Oey, M. S. 2014, ApJ, 791, L19
  31. Johnson, M. D., Levitt, J. S., Henry, R. B. C., & Kwitter, K. B. 2006, in IAU Symposium, Vol. 234, Planetary Nebulae in our Galaxy and Beyond, ed. M. J. Barlow & R. H. Méndez, 439–440
  32. Kakazu, Y., Cowie, L. L., & Hu, E. M. 2007, ApJ, 668, 853
  33. Kauffmann, G., Heckman, T. M., Tremonti, C., et al. 2003, MNRAS, 346, 1055
  34. Kennicutt, Jr., R. C. 1998, ApJ, 498, 541
  35. Kennicutt, Jr., R. C., Bresolin, F., & Garnett, D. R. 2003, ApJ, 591, 801
  36. Kewley, L. J., & Dopita, M. A. 2002, ApJS, 142, 35
  37. Kewley, L. J., Dopita, M. A., Sutherland, R. S., Heisler, C. A., & Trevena, J. 2001, ApJ, 556, 121
  38. Kewley, L. J., & Ellison, S. L. 2008, ApJ, 681, 1183
  39. Kobulnicky, H. A., & Kewley, L. J. 2004, ApJ, 617, 240
  40. Lamareille, F., Contini, T., Brinchmann, J., et al. 2006, A&A, 448, 907
  41. Lee, J. C., Salzer, J. J., & Melbourne, J. 2004, ApJ, 616, 752
  42. Lequeux, J., Peimbert, M., Rayo, J. F., Serrano, A., & Torres-Peimbert, S. 1979, A&A, 80, 155
  43. Massey, P., Strobel, K., Barnes, J. V., & Anderson, E. 1988, ApJ, 328, 315
  44. McCall, M. L., Rybski, P. M., & Shields, G. A. 1985, ApJS, 57, 1
  45. McGaugh, S. S. 1991, ApJ, 380, 140
  46. Melbourne, J., Phillips, A., Salzer, J. J., Gronwall, C., & Sarajedini, V. L. 2004, AJ, 127, 686
  47. Melbourne, J., & Salzer, J. J. 2002, AJ, 123, 2302
  48. Mendoza, C., & Bautista, M. A. 2014, ApJ, 785, 91
  49. Mendoza, C., & Zeippen, C. J. 1982, MNRAS, 199, 1025
  50. Moustakas, J., Kennicutt, Jr., R. C., Tremonti, C. A., et al. 2010, ApJS, 190, 233
  51. Nicholls, D. C., Dopita, M. A., & Sutherland, R. S. 2012, ApJ, 752, 148
  52. Nicholls, D. C., Dopita, M. A., Sutherland, R. S., Kewley, L. J., & Palay, E. 2013, ApJS, 207, 21
  53. Oke, J. B., Cohen, J. G., Carr, M., et al. 1995, PASP, 107, 375
  54. Osterbrock, D. E., & Ferland, G. J. 2006, Astrophysics of gaseous nebulae and active galactic nuclei
  55. Pagel, B. E. J., Edmunds, M. G., Blackwell, D. E., Chun, M. S., & Smith, G. 1979, MNRAS, 189, 95
  56. Peimbert, A. 2003, ApJ, 584, 735
  57. Peimbert, M. 1967, ApJ, 150, 825
  58. Peimbert, M., & Costero, R. 1969, Boletin de los Observatorios Tonantzintla y Tacubaya, 5, 3
  59. Peimbert, M., Peimbert, A., Esteban, C., et al. 2007, in Revista Mexicana de Astronomia y Astrofisica Conference Series, Vol. 29, Revista Mexicana de Astronomia y Astrofisica Conference Series, ed. R. Guzmán, 72–79
  60. Pérez-Montero, E., Díaz, A. I., Vílchez, J. M., & Kehrig, C. 2006, A&A, 449, 193
  61. Pilyugin, L. S. 2000, A&A, 362, 325
  62. Pilyugin, L. S., & Thuan, T. X. 2005, ApJ, 631, 231
  63. Rayo, J. F., Peimbert, M., & Torres-Peimbert, S. 1982, ApJ, 255, 1
  64. Reddy, N. A., & Steidel, C. C. 2009, ApJ, 692, 778
  65. Rubin, V. C., Ford, Jr., W. K., & Whitmore, B. C. 1984, ApJ, 281, L21
  66. Salzer, J. J., Jangren, A., Gronwall, C., et al. 2005a, AJ, 130, 2584
  67. Salzer, J. J., Lee, J. C., Melbourne, J., et al. 2005b, ApJ, 624, 661
  68. Salzer, J. J., Williams, A. L., & Gronwall, C. 2009, ApJ, 695, L67
  69. Salzer, J. J., Gronwall, C., Lipovetsky, V. A., et al. 2000, AJ, 120, 80
  70. —. 2001, AJ, 121, 66
  71. Salzer, J. J., Gronwall, C., Sarajedini, V. L., et al. 2002, AJ, 123, 1292
  72. Skillman, E. D., Kennicutt, R. C., & Hodge, P. W. 1989, ApJ, 347, 875
  73. Skillman, E. D., & Kennicutt, Jr., R. C. 1993, ApJ, 411, 655
  74. Skillman, E. D., Televich, R. J., Kennicutt, Jr., R. C., Garnett, D. R., & Terlevich, E. 1994, ApJ, 431, 172
  75. Stasińska, G. 1980, A&A, 85, 359
  76. —. 1990, A&AS, 83, 501
  77. —. 2005, A&A, 434, 507
  78. Steidel, C. C., Giavalisco, M., Pettini, M., Dickinson, M., & Adelberger, K. L. 1996, ApJ, 462, L17
  79. Tremonti, C. A., Heckman, T. M., Kauffmann, G., et al. 2004, ApJ, 613, 898
  80. van Dokkum, P. G. 2001, PASP, 113, 1420
  81. van Zee, L., Skillman, E. D., & Haynes, M. P. 2006, ApJ, 637, 269
  82. Wegner, G., Salzer, J. J., Jangren, A., Gronwall, C., & Melbourne, J. 2003, AJ, 125, 2373
  83. Werk, J. K., Jangren, A., & Salzer, J. J. 2004, ApJ, 617, 1004
  84. Zaritsky, D., Kennicutt, Jr., R. C., & Huchra, J. P. 1994, ApJ, 420, 87
  85. Zurita, A., & Bresolin, F. 2012, MNRAS, 427, 1463
Comments 0
Request comment
You are adding your first comment
How to quickly get a good reply:
  • Give credit where it’s Due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should ultimately help the author improve the paper.
  • Remember, the better we can be at sharing our knowledge with each other, the faster we can move forward.
The feedback must be of minumum 40 characters
Add comment
Loading ...
This is a comment super asjknd jkasnjk adsnkj
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test description