abstract: Automatically scene understanding using machine learning algorithms has been widely applied to different industries to reduce the cost of manual labor. Nowadays, insurance companies launch express vehicle insurance claim and settlement by allowing customers uploading pictures taken by mobile devices. This kind of insurance claim is treated as small claim and can be processed either manually or automatically in a quick fashion. However, due to the increasing amount of claims every day, system or people are likely to be fooled by repeated claims for identical case leading to big lost to insurance companies.Thus, an anti-fraud checking before processing the claim is necessary. We create the first data set of car damage images collected from internet and local parking lots. In addition, we proposed an approach to generate robust deep features by locating the damages accurately and efficiently in the images. The state-of-the-art real-time object detector YOLO is modified to train and discover damage region as an important part of the pipeline. Both local and global deep features are extracted using VGG model, which are fused later for more robust system performance. Experiments show our approach is effective in preventing fraud claims as well as meet the requirement to speed up the insurance claim prepossessing.
author: Pei Li
Univeristy of Notre Dame
; BingYu Shen
University of Notre dame
; Weishan Dong
IBM Research China
bibliography: egpaper_final.bib
title: An Anti-fraud System for Car Insurance Claim Based on Visual Evidence

# Introduction

Fully or semi-automatically insurance claim processing could be very useful in the insurance industry when handling small but more frequent insurance claims that under a certain amount of rate. It increases the speed of claim investigation and loss assessment. Some of the companies employ 3G mobile techniques to accelerate the investigation process by freeing the customers from submitting complicated claim proves in order to shorten the whole claim process up to eight hours. However, fraud claims for which the payment has been claimed more than one time could potentially cause loss to insurance companies. Manually validating on large scale of claims cannot meet the speed requirement for express claim process anymore and reporting a duplicated claim candidates automatically before the settlement is needed. However, this solution remains a challenging task due to a number of factors: In pictures or video frames captured by mobile devices, the viewpoint and illumination conditions are not taken in a controlled environment; Damage cannot be localized accurate and fast enough to generate robust local features for matching; No standard data set is available for analyzing this problem. We made two major contributions to address the problem:

-   we contribute the first car damage dataset which contains samples collected from both internet and local public parking lot. In addition, manually annotation is provided.

-   We proposed an anti-fraud system prevent repeated claims happens before the settlement is issued.

The whole system consists of two parts: a real-time damage detector to provide accurate damage locations in the picture and deep feature extractor to generate combined global and local deep features for anti-fraud matching in claim history database. We fine-tuned our data on the real-time object detector YOLO which reported to have 457FPS on VOC2007 with a mAP of 63.4 with its standard model. To handle the affine transformation of the damages in images due to viewpoint changes, the possible various lighting condition and the shift caused by the detection result, we tried to employ local respond normalization layer and dropout as well as data augmentation in the training process to make the model more generalized. Our system is evaluated with both data set collected on internet from three of the biggest searching engines and local public parking lot.

# Related Works

## General Object Detection and Damage Detection

Accurate localization of the damage on a car plays an important role in our system.General object detection is a important topic in computer vision. Classical detection pipeline will first conduct a sliding-window fashion candidate region selection. Then either hand-craft feature or deep features will be extracted from each of the candidate regions for binary classification that tells if it's a detection or not. Later, Deformable parts models(DPM) is introduced to run detection with detectors trained on specific regions of the image in an exhaustive manner across all locations and scales.The exhaustive search through all possible locations and scales poses a computational challenge and big time delay on each detection. Recently, some detection frameworks replaced the sliding windows with region proposals generated either from selective search or from feature map bonding box prediction, such as . Another way of doing object detection is to treat the problem as a bounding box regression problem in general , but among all of them, YOLO is the only one that can meet real-time requirement. Damage detection is a core problem during quality control. Several damage detection approaches have been proposed applied to car body damage detection. Srimal.al propose to use 3D CAD models to handle automatic vehicle damage detection via photograph. Gontscharov.al tries to solve vehicle body damage multi sensor-data fusion. Keyence Vision proposed an industrial solution for car damage by hail by applying a high-resolution Multi-camera vision system. Cha.al.adopt image based deep learning to detect crack damages in concrete while adopted a phase-based optical flow and unscented Kalman filters.summary of the up-to-date crack detection using image processing which is similar to car damage detection.

# Dataset

We introduced a car damage dataset collected from one of three biggest searching engines: Google Bing and Baidu using keywords like 'scratch','dent'. This contains 1790 images with manual annotation. In addition, we collect vehicle images with small damage in public parking lot without including the identification information of the vehicles.These data set has 92 cars each with four or five images and are also manually annotated with bounding boxes tightly bounds the damage regions. In order to mimic customer behaviors, we use mobile devices to capture the most convincing evidence for a successful insurance claim. In each case, we capture images of a whole car body and local scratch images from the same vehicle. Each of the cases contains two images captured the same damage regions which represent the first claim and the other represents the repeated claim that treated as a fraud claim. Example images are shown in Figure and Figure with random samples from the dataset. The whole dataset is used for two purposes: one for damage detection training and evaluation, another for studying the whole system which contains both damage localization and fraud claim detection.

# The anti-fraud system for car insurance claim

The anti-fraud or fraud detection system is coupled with the claim history database. Each user has their own record in the database after they have a claim for the first time. The fraud claim could happen in two ways: the same-vehicle reclaim or cross-vehicle reclaim. Same-vehicle reclaim happens when the user tries to make a quick insurance claim by uploading a similar image which is collected during the same case. Another fraud claim would happen when the user uses a similar image taken from another car as his or her own vehicle claim. Both formats could lead to reclaim for the same cases actually has been issuing the settlement. Traditional insurance claim is handled manually by representatives which is necessary when the claim value is above a certain threshold. However, when the claim number increasing each day and some of the claims are only about small damages such as scratches dent etc. Manually handle all these claims decrease the efficiency of insurance companies' service and increase the cost of staffing. In this case, a system that can help to detect the fraud claims is needed when we want to make the system more automatically. As shown in Figure , the user is required to upload several images following the requirement in order to open a new claim. The anti-fraud system will then using algorithms to extract features to search in the enrolled history database. If the system reports a high score of suspicious of a certain claim, humans can take into the loop and manually retrieve the history to do the loss assessment. To make the system more accurate and fast, we need a robust feature which based on accurate locating the damages in the image. We adopt two state-of-the-art real-time detectors for damage detection and propose an approach to form robust features for the anti-fraud searching. We will introduce each part in the following sections.

## Damage detection

In most of the cases, users are required to upload few images of the close view of the damages on the vehicle and another picture with a farther view that contains the whole body of the car as well as the plate which shows the model and identity of the car. The first step to obtaining a robust feature descriptor for fraud detection is the damage detection. This task is challenging since it is not a traditional object detection problem and the damages would have different forms. We consider the damage mainly three types due to analyzing our dataset: scratch, dent, and crack.

Localization of a damage in an image uploaded by users is a critical part of our system as the system would obtain more stable and discriminative features when accurate detection happens. In addition, to meet our real-time requirement without sacrifice detection accuracy, we employed YOLO, which is a state-of-the-art lightweight real-time detector to be part of the system. We proposed our approach to train and evaluate on the new dataset we collect.

### YOLO Detector

YOLO is a fast, accurate object detection framework.Its base network runs at 45 frames per second with no batch processing on a Titan X GPU and a fast version runs at more than 150 fps which allows us to process streaming video in real-time with less than 25 milliseconds of latency. It is extremely useful to assist needs in auto-drive, real-time responsive services. YOLO frame object detection as a single regression problem.Image pixels are mapped to bounding box coordinates and class probabilities straightly. As damages like scratch or dent might have a different shape, YOLO overcome this by learning features through regression of four coordinates.YOLO, with its standard model, achieve 45FPS with mAP of 63.4 on VOC2007. All these performances meet the requirement of the efficient needs of an express anti-fraud system for insurance claims.

## Damage Re-Identification

After damage detection generated, robust features need to be extracted to overcome lighting conditions and different angle of views. Also, both global and local features are being used to finally determine whether a claim has ever appeared in the history database as each claim requires at least one image of the image containing the close view of the damage as well as another image shows the whole view of the body of the corrosponding vehicle . To extract local features, we employ the pre-trained VGG16 object recognition model as a feature extractor. Global feature consists of global deep features and color histogram. We fuse these feature to obtain a more discriminative feature and conduct a 1 to N match between the probe images and the images in the post-claim history database.

# Experiments and result

## Damage detection

To cope with the size of the data set, we select different pre-trained models in YOLO and conduct a transfer learning by fine-tuning with our data. This contains two sets of experiments: Training and evaluation on the web-collected dataset and evaluation on the real car damage data set from a real parking lot. For the first set of experiment, we use around 1300 images for training which come from Google and Bing and 400 images for testing which are from Baidu. We evaluated on different network structures starting with the YOLO-tiny model which is the simplest and fastest network structure among the models reported in . We later tried to add dropout layers and local respond normalization layers aiming to handle the different lighting condition and reduce overfitting on our small scale data set. From the experiment result shown in Table, we can see YOLO with LRN layer perform the best in general achieved 81.7 percent of recall and 37.96 percent of precision at threshold 0.1. The YOLO with LRN layer and dropout layer achieve better precision at 47.04 percent while lower recall. The VGG-16 has a lower performance than the YOLO original model and the Tiny YOLO model achieves the lowest performance due to the size of the net architecture.

Later, we conduct an experiment over real car damage data set which is collected in the local parking lot. We designed two protocols for the training and testing. Consider the real system is updated periodically enhancing the model, all the images recorded as history in the system could be used for training the model. As shown in . In this case, we consider this evaluation protocol subject-overlapped detection. Another evaluation protocol follows the previous experiment which the training and testing set are separated by car ID. The precision-recall curve shows in indicates that the gap between subject-overlapped and subject-disjoint detection result. With more relevant images used for training the detector, the more accurate a similar damage could be detected and as the size of the history database increase, better performance of the detector is expected .

```latex
\begin{table*}[htbp]\centering
  \caption{Precision Recall @thresholds}
    \label{pr_on_baidu}
    \begin{tabular}{|l|r|r|r|r|r|}
    \toprule
          & \multicolumn{1}{l|}{YOLO\_original} & \multicolumn{1}{l|}{YOLO\_vgg\_16} & \multicolumn{1}{l|}{YOLO\_with\_LRN} & \multicolumn{1}{l|}{YOLO\_with\_LRN+dropout} & \multicolumn{1}{l|}{YOLO\_tiny} \\
    \midrule
    precision & 12.66 & 11.1  & 14.62 & 14.27 & 9.41 \\
    \midrule
    recall & 85.6  & 82.01 & 89.72 & 78.66 & 72.49 \\
    \midrule
    precision & 25.58 & 32.75 & 37.96 & 47.04 & 17.48 \\
    \midrule
    recall & 76.61 & 57.58 & 81.75 & 63.54 & 66.84 \\
    \midrule
    precision & 40.89 & 49.07 & 38.63 & 86.25 & 26.26 \\
    \midrule
    recall & 66.32 & 27.25 & 72.49 & 17.74 & 60.15 \\
    \bottomrule
    \end{tabular}%
  \label{tab:addlabel}%
\end{table*}
```

## Fraud claim detection

In this section, we design experiment aiming to mimic the real fraud cases. Under this scenario, the history database contains at least one images that has been claimed and issued settlement. We consider the most general case which is the cross vehicle search. It means that we have to search across the whole database not only the history belong to a certain user. We random chose one damage image per vehicle as probe and the rest as gallery. To get a more robust feature for fraud detection, we exploit feature fusion combining local feature, global feature as well as histogram feature to obtain a final feature descriptor. Local feature is extracted from the damage region which is detected by the damage detector. Global feature is extracted from the corresponding full vehicle image which could potentially help to make the feature more discriminative since the damages could be very similar from different claims. We also incorporate histogram feature as a global feature expecting to enhance the final descriptor. All the features from the scratch and car body are extracted from the last FC layer of the pre-trained VGG-16 model and concatenate to form the final descriptors. Cosine distance is used to compare each probe to the gallery images.Rank-1 and Rank-10 rate are reported as shown in . We also study how different dimension of the histogram feature affect the performance. We firstly used the manually annotated damage region as input which we assume that the damage detection is perfect. After this experiment, we From the experiment result shown in we gain a 41.0 percent rank-1 rate with 8 bins in each channel of the histogram feature. When we increase the bin to 16 and 32, the performance increase to 44.6 percent and 46.73 percent. For the rank-10 performance, increasing the dimension of the color histogram feature does not always boost up the performance. When adopting the 32 bin color histogram feature, a slightly performance drop happens. However, if only employ the global feature obtained from the whole vehicle body, we only gain 19.6 percent rank-1 rate, while incorporating the feature from the scratch enhance the rank-1 rate to 56.52 percent. This indicates that the histogram feature might have negative impact on the performance when the color of the candidate vehicles are very similar.

We also explore the matching performance by substitute the detection result from the ground truth with the ones yield from the YOLO detector we trained. We discover only a tiny performance difference when directly employ the detector which indicates that the detector if trained well could yield a robust result for the damage feature extraction.

# Conclusion and discussion

In this paper, we introduce a useful data set of car damage which is the first data set for car damage detection and matching. We proposed a practical system pipeline to detect the fraud claim before issuing the settlement in order to reduce loss for the insurance company. We research on different deep architectures for better damage detection result and analyze different feature fusions and its impact on the final matching performance.Extensive experiments with different experiment settings under real scenario demonstrate the effectiveness of the proposed approach.


abstract:  Some aspects of the chemical evolution of late-type dwarf galaxies are reviewed, together with their implications on some issues of cosmological relevance. A more detailed approach to model their evolution is suggested. The importance of deriving the star formation history in these systems by studying their resolved stellar populations is emphasized.
author: Monica Tosi

Osservatorio Astronomico, Via Ranzani 1, 40127 Bologna, Italy
title: **Stellar Populations and Chemical Evolution of Late–Type Dwarf Galaxies**

# Introduction

People have been studying the evolution of late-type dwarf galaxies for two-three decades, but only in the last 5–6 years this has become a fashionable research field. This late interest is probably due to the circumstance that cosmologists have recently discovered that understanding the evolution of the gaseous and stellar constituents of these systems is a necessary step for a correct analysis of some important problems of cosmological relevance. The connection between late-type dwarfs and cosmology is at least three-fold: 1) given the high gas content and low metallicity of irregular and blue compact galaxies (hereinafter, Irrs and BCDs, respectively), they are the closest analogues to primeval galaxies; 2) they are the systems currently providing the safest empirical estimate of the primordial $^4$He abundance; 3) they have been suggested 5-6 years ago (e.g. Lilly et al. 1995, Babul & Ferguson 1996) to be the local counterparts of the blue objects in excess in deep galaxy counts at intermediate redshifts.

Here I will try to summarize what can be learnt from chemical evolution models and stellar population studies of Irrs and BCDs and what this implies for these three issues of cosmological relevance. In the following discussion, I will emphasize the open questions rather than the (important) results already achieved in the field.

# Evolution of late-type dwarfs

When a galaxy forms, after a while it begins to form stars and starts a series of cycles following in general the scheme drawn by Tinsley (1980). The stars evolve and synthesize in their interiors heavier and heavier elements, and then eject them in the surrounding medium when they lose mass and die. In this way they pollute the interstellar medium (ISM) and modify both its mass and chemical composition. In the meantime, the ISM mass and metallicity may change also for gas exchanges with adjacent regions (gas losses or accretions, or both). The next generation of stars thus form in that ISM with a different initial composition and must have therefore a slightly different evolution.

Models for galaxy chemical evolution take these phenomena into account, as well as their effects on the stellar and gaseous properties, by simplisticly parametrizing the physical mechanisms. The major parameters are the star formation (SF) law, the initial mass function (IMF), the gas flows in and out of the considered region, and the quantities involved in stellar nucleosynthesis (e.g. stellar lifetimes, yields, mass loss, opacities, treatment of convection, etc.).

Chemical evolution models for Irrs and BCDs have been computed by several groups in the last twenty years (e.g. Matteucci & Chiosi 1983, Matteucci & Tosi 1985, Pilyugin 1993, Marconi et al. 1994, Carigi et al. 1995, Larsen et al. 2000, Recchi et al. 2001). The combination of the predictions of these models with the information on the evolution of Irrs and BCDs from other kinds of studies have interesting consequences on the three issues of cosmological relevance mentioned in the Introduction.

## Irrs and BCDs as almost primeval galaxies ?

Back in 1973, Searle et al. suggested that the blue colors, low metallicities and high gas fractions of BCDs could be self-consistently explained only if these systems either are experiencing now their first episode of SF activity or have a very discontinuous SF regime, dominated by short intense bursts separated by long quiescent phases (see e.g. the top right-hand panel in Fig.).

Almost thirty years later, however, no convincing case of galaxy at its first SF experience has been found yet: all the systems for which sufficiently deep photometry exists reveal the presence of stellar populations as old as the maximum lookback time allowed by the photometric depth (see also Schulte-Ladbeck, this volume). Even IZw18, the most metal-poor galaxy ever discovered, when imaged with HST, has been recognized to contain stars at least a few hundreds Myr old. Aloisi et al. (1999) have shown that the colour-magnitude diagram (CMD) derived from optical WFPC2 photometry of IZw18 contains several faint red stars (see left panel in Fig.) which can only be objects in the asymptotic giant branch (AGB) phase of low and intermediate mass stars and must, therefore, have ages between a minimum of 200-300 Myr up to a few Gyrs. Exactly the same conclusion has been independently reached by Ostlin (2000) on the basis of near-infrared Nicmos imaging (right panel in Fig.).

Since lack of evidence is certainly not evidence of lack, this result does not exclude that really young galaxies do exist, but clearly restricts their possible number to a small fraction of the known late-type dwarfs. If Searle et al. (1973) were right, then, the vast majority of BCDs should have a bursting mode of SF. This was indeed suggested by many chemical evolution models (e.g. Matteucci & Tosi 1985, Pilyugin 1993), which confirmed that the best way to reproduce simultaneously colours, chemical abundances and gas fractions of BCDs (and possibly of Irrs) was to have only a few (7–10, at most) episodes of strong SF activity, with long quiescent intervals. This conclusion has however been questioned by more recent models (e.g. Carigi et al. 1995, Legrand 2000) suggesting that continuous SF regimes, with very low rates, could account equally well for the observed features.

The lack of unique scenarios from chemical evolution models of late-type dwarfs is due to the fact that we do not have yet sufficient observational data to properly constrain them. For what concerns the SF regime, the most direct way to infer it from observations is from the CMDs of systems close enough to be resolvable in single stars.

## The helium–metallicity relation

The primordial $^4$He abundance is nowadays one of the hottest issues in cosmology, due to the recent observational data (Tytler et al. 2000 and references therein) showing that primordial D was quite low, (D/H)$_p\simeq3\times10^{-5}$, and, hence, either $^4$He was high (in mass fraction Y$_p>$`<!-- -->`{=html}0.24), or the standard theories of Big Bang Nucleosynthesis are not self-consistent. The only method currently considered reliable to derive the primordial helium abundance is based on the extrapolation to zero-metallicity of the helium–metallicity relation derived from observations of HII regions in metal-poor galaxies. Since the heavy elements more reliably derived from HII regions are oxygen and nitrogen, what people actually use are the helium–oxygen and helium–nitrogen relations, inferred from a linear fit to the data (see e.g. Izotov et al. 1997).

One of the questions on this method is whether or not these relations can be taken as linear. He, N and O are synthesized in different stellar mass ranges: helium is produced by all stars, nitrogen almost exclusively by intermediate mass stars, and oxygen only by massive ones; and this implies that the timescales for their ISM enrichment are quite different from each other. This point was explicitely addressed by Pilyugin (1993) who showed with his chemical evolution models that the evolution of He with O (and with N) in a BCD is not linear but has typically a saw-tooth shape (see left-hand panel in Fig.).

Furthermore, all chemical evolution models of these galaxies invoke the occurrence of significant galactic winds during their life, and these winds are often required to be more effective in removing elements (like oxygen) expelled from massive stars. If these differential winds occur, this also makes the behaviour of helium with oxygen non linear during the lifetime of a dwarf, as shown both by Pilyugin (1993) and Marconi et al. (1994) and illustrated in the right-hand panel of Fig.. The questions are: Do the winds actually occur and in what percentage of Irrs and BCDs ? If they do occur, do they really remove some (which ?) elements for ever ? With what efficiency ? Several optical and X-ray observations of ionized filaments and bubbles in some Irrs and BCDs confirm that somewhere the phenomenon does take place and seem to preferentially eject alpha-elements (see e.g. Meurer et al. 1992, Sahu & Blades 1997, Hensler et al. 1998, for NGC1705). However, hydrodynamical models of the SN ejecta in dwarf galaxies do not agree yet on the timescales or on the efficiency of the gas removal (cfr. Silich & Tenorio-Tagle 1998, D'Ercole & Brighenti 1999, Mac Low & Ferrara 1999).

In spite of the above arguments, I personally believe that a linear fit to the observed helium–oxygen and helium–nitrogen distributions derived from HII regions is the correct approach to infer the primordial $^4$He. In fact, each data point in the empirical distribution corresponds to the final (i.e. present epoch) stage of the He–O and He–N behaviours in each galaxy, and to derive the mean initial $^4$He one needs to average over all the possible behaviours of the whole sample of galaxies. The observed shape of the distribution suggests indeed that this is well represented by a linear fit.

## Local counterparts of faint blue galaxies ?

Some authors (e.g. Lilly et al. 1995, Babul & Ferguson 1996) have suggested that starbursting dwarfs may be the local counterparts of the faint blue galaxies usually found in excess in deep galaxy counts at redshift of 0.7–1.0. To verify whether this is a likely possibility, one needs to check if most of the dwarfs have been forming stars at the epochs corresponding to that redshift (hence about 6–7 Gyr ago) and if, in case, that SF activity was strong enough (i.e. $\sim$`<!-- -->`{=html}1 M$_{\odot}$yr$^{-1}$) to guarantee the observed luminosity levels.

We have already mentioned above that all the galaxies examined so far appear to contain stars as old as the lookback time corresponding to the faint mag limit of the available photometry. This makes it extremely probable that most of them have already been active several billion years ago. The tough question is on the intensity of the SF. Van Zee (2001) has very recently surveyed the SF histories of a large sample of distant, isolated dwarf Irrs and concludes that their SF activities are too modest and it is very unlikely that these systems could significantly contribute to the faint blue galaxy excess.

This conclusion is inevitably based on integrated properties of these distant objects. To obtain a more quantitative information, it is necessary to examine nearby systems where the SF history can be derived directly from the resolved stellar populations and a more detailed evaluation of the SF rates at different epochs can be performed.

# SF histories from stellar populations

As discussed in the previous sections, for a correct approach to study the evolution of Irrs and BCDs and for several issues of cosmological relevance, it is of crucial importance to derive the SF history of a large number of galaxies of these types. The most direct estimate of the epochs and of the intensities of the SF activities is derivable from the CMDs, since their morpholoy is the consequence of the system evolutionary conditions. Fig. shows a few examples of how the CMD morphology is affected by the SF history and IMF. We thus proposed several years ago (Tosi et al. 1989, 1991) the method of synthetic CMDs to derive from the observational CMDs as much information as possible, taking into account all the incertainties related to both observations and theories. Similar procedures have been developed also by other groups (e.g. Aparicio et al. 1996, Tolstoy & Saha 1996) and are now widely applied (see also the contributions by Grebel, Schulte-Ladbeck, Tolstoy, this volume).

The numerical procedure creates synthetic CMDs, via MonteCarlo extractions on homogeneous sets of stellar evolution tracks (e.g. the Geneva or the Padova sets) of various metallicities, taking into account all the involved theoretical parameters and uncertainties (age, metallicity, IMF, SF law, stochastic effects due to small number statistics, etc.). The synthetic CMDs have the observed number of stars of the examined region and its observational uncertainties (photometric errors, incompleteness factors, blending). The latter is a very important point which requires a special accuracy during the photometric analysis. A model can be considered satisfactory only if it reproduces all the features of the observational CMD and luminosity function (LF).

By comparing all the simulated cases with the empirical CMD and LF we select those which are more consistent with them and thus derive the number, starting epoch, duration and intensity of the SF episodes, their IMF, the durations of significant quiescent intervals (if any), and hints on the most likely metallicity of the various stellar generations. In most, if not all, the cases, the results are not unique, but we can sensibly reduce the range of acceptable values of the various parameters, thus obtaining a good indication on the most likely evolution of the examined region.

We have applied the method to Local Group Irrs with ground-based observations and to more distant late-type dwarfs observed with HST (see Table 1 for a summary of the results). Fig. shows two cases of synthetic CMDs for IZw18: it is apparent that a single recent SF burst cannot account for the observed red stars (both faint and bright) and predicts too many too bright objects, despite the rather steep IMF. The case with two episodes is instead consistent with the observed features. The most recent application is being performed on the BCD NGC1705 and its preliminary results are described by Annibali et al. (this volume). Of the 8 dwarfs examined so far by our group, only NGC1569 shows evidence for a very intense and short burst (Greggio et al. 1998), in addition to other longer, less striking episodes. All the others appear to have had a gasping SF regime, with episodes of much more moderate activity. None of them shows any evidence for quiescent phases longer than 10$^8$ yr.

If we combine the results on the SF history of the late-type dwarfs studied by our group with those obtained by other groups (see e.g. Grebel 1998 for a comprehensive review of Local Group galaxies and Schulte-Ladbeck, this volume, for BCDs outside the Local Group) it is apparent that these galaxies, with very few exceptions, have had qualitatively similar SF histories. Fig. shows the SF histories derived in 5 galaxies: NGC1569 is the only one which appears to follow the bursting scheme. The data available so far have not allowed yet to accurately derive its earlier activity, but the new HST photometry in the near infrared by Aloisi et al. (2001) should allow for a lookback time of, at least, a few Gyrs. In summary:

-   Both Irrs and BCDs seem to have a roughly continuous or gasping SF regime, with various episodes of moderate activity and no clear sign of long periods of total inactivity.

-   Only the few cases of the type of NGC1569 have bursts with SF rates high enough to allow for the luminosity attributed to faint blue galaxies at redshift 0.7–1.

-   No evidence has been found yet for dwarfs at their first SF activity.

```latex
\begin{table*}
\begin{center}
\caption{Dwarf irregulars and BCDs in our program. Listed for each
galaxy are: name, coordinates, derived distance
modulus, number of resolved stars with small photometric error, and SF mode
derived from the synthetic CMD method.}
\begin{tabular}{lccccc}
\hline\hline
Name & R.A. & DEC&  (m$-$M)$_o$ & sel obj & SF mode\\
\hline
DDO ~~70 (Sex B)   &09 57 23&+05 34 07& 25.6 &1300 & gasps \\
DDO 209 (NGC 6822) &19 42 07&$-$14 55 01& 23.5 &1772& gasps \\
DDO 210 (Aquarius) &20 44 08&$-$13 02 00& 28 ?& 633& gasps \\
DDO 221 (WLM)      &23 59 23&$-$15 44 06& 25.0 & 2000& gasps \\
DDO 236 (NGC 3109) &10 00 48&$-$25 55 00& 25.7 &2605& gasps \\
\hline
NGC 1569    & 04 26 04 & +64 44 29   & 26.7 & 801& gasps+burst \\
IZw18       & 09 34 02 & +55 14 19   & 30.0 & 150& gasps \\
NGC 1705    & 04 54 13 & $-$53 21 40 & 28.6 & $17000$& gasps \\
\hline
\end{tabular}
\end{center}
\end{table*}
```

# Future prospects

From the above presentation, it is apparent that further studies on several issues are still required to obtain a reliable scenario for the evolution of late-type dwarfs. In particular, the SF histories and the major characteristics of the gas flows triggered by SN explosions should be examined in more details. For this reason I believe that, at present, the correct approach to model the chemical evolution of these systems would be to concentrate on single representative cases, rather than to model the overall features of a large sample of galaxies. For the chosen cases I would:

-   Try to get as much information as possible on the evolutionary parameters from observations (e.g. IMF and SF from populations synthesis and synthetic CMD; gas distribution and flows from multiwavelength observations; etc.);

-   Compute the hydrodynamics of the ISM and SN ejecta in the conditions of the examined galaxy derived from the first point;

-   Only at this point, model the chemical evolution of the examined galaxy including all the results of the previous points as input data.

**Acknowledgements**

I warmly thank the organizers of this enjoyable and interesting meeting for the invitation and the financial support. I sincerely congratulate them for the success of their Graduate School. This work has been partly supported by the Italian MURST through Cofin2000.


author: B.T. Gänsicke$^1$, D. Koester$^2$, J. Girven$^1$, T.R. Marsh$^1$, D. Steeghs$^1$
$^{1}$ Department of Physics, University of Warwick, Coventry CV4 7AL, UK
$^{2}$ Institut für Theoretische Physik und Astrophysik, University of Kiel, 24098 Kiel, Germany
title: Two white dwarfs with oxygen-rich atmospheres

> **Stars with masses in the range $7-10\,M_\odot$ end their lives either as massive white dwarfs or weak type II supernovae, and there are only limited observational constraints of either channel. Here we report the detection of two white dwarfs with large photospheric oxygen abundances, implying that they are bare oxygen-neon cores and that they may have descended from the most massive progenitors that avoid core-collapse.**

White dwarfs represent the endpoints of stellar evolution for the overwhelming majority of all stars. Most white dwarfs in the Galaxy have carbon-oxygen (CO) core compositions, being descendants of low to intermediate mass stars that underwent hydrogen and helium core burning. Stars with initial masses $7M_\odot\le M \le 10M_\odot$ will reach sufficiently high core temperatures to proceed to carbon burning, and produce either oxygen-neon (ONe) core white dwarfs, or undergo a core-collapse supernova (SN II) via electron capture on the products of carbon burning . The exact outcome of stellar evolution in this mass range depends critically on the detailed understanding of mass loss , the relevant nuclear reaction rates and the efficiency of convective mixing in the stellar cores . Some observational constraints on stellar models come from analysis of SN II progenitors , which suggest a lower limit on the progenitor masses of $8^{+1.0}_{-1.5}M_\odot$. A linear extrapolation of the observed relation between the masses of white dwarfs and their progenitors up to the maximum mass of white dwarfs leads to a broadly consistent result.

Information about the core compositions of white dwarfs, in particular of those descending from the highest possible masses, would help to test and improve the theory of stellar evolution. Unfortunately, almost all white dwarfs have hydrogen and/or helium envelopes that, while low in mass, are sufficiently thick to shield the core from direct view. Asteroseismology has the potential to unveil their core composition , but observational studies attempting to exploit this potential remain ambiguous \[see e.g. the case of GD 358, \]. Core material can be directly detected in the photospheres of a small number of stars that underwent a very late thermal pulse during their asymptotic giant branch (AGB) evolution, ejecting a large fraction of the envelope and leaving a white dwarf with only a thin layer of helium. Examples of hydrogen-deficient post-AGB stars are the hot PG1159 objects, white dwarfs with helium-rich (DB) atmospheres and cool carbon-rich (DQ) white dwarfs. The most extreme case known to date is 1H1504+65, which is hydrogen and helium deficient, representing a hot, naked CO surface .

Recently, the Sloan Digital Sky Survey (SDSS), through its comprehensive spectroscopic snapshot of the Galactic stellar population, has revealed a small class of hydrogen-deficient white dwarfs . Their spectra are consistent with nearly pure carbon atmospheres . It has been suggested that these "hot" DQ white dwarfs represent the evolutionary link between objects such as PG1159 and 1H1504+65, the DB white dwarfs, and the cool DQ white dwarfs .

Until now, all known white dwarfs for which photospheric oxygen and carbon abundances have been determined have abundance ratios $\mathrm{O/C}\le1$. Some stellar models predict that the most massive stars avoiding core-collapse will result in ONe white dwarfs with very low carbon abundances . Should such a core lose its hydrogen envelope, an extremely oxygen-rich spectrum would be expected. Here we present the results of a search for white dwarfs with large photospheric oxygen abundances within the SDSS spectroscopic Data Release (DR) 7 .

We selected all spectroscopic objects in DR7 that fall within the $ugriz$ color space of the white dwarfs from . We then subjected these spectra to an automatic measurement of the equivalent widths \[see for details of the method\] of the O I 616 nm and 778 nm multiplets, which are the strongest O I lines expected in a cool, oxygen-rich atmosphere. Another strong O I multiplet is located at 845 nm, however, this region can be affected by residuals from the night sky line subtraction, and we therefore did not include it in our procedure. As an additional constraint, we required the spectra to have a signal-to-noise ratio $\ge10$ in the regions around the two O I lines, resulting in the analysis of 25639 SDSS spectra .

We found 1000 spectra with a formal $>4\sigma$ detection of both O I lines, and visual inspection of these spectra showed that 998 were erroneously flagged by our automated routine because of poor data quality. Only two objects with genuine O I 616 nm and 778 nm absorption, SDSS J092208.19+292810.9 and SDSS J110239.69+205439.4, survived our scrutiny. Both have high proper motions, 0.275"/yr and 0.163"/yr, respectively , consistent with them being nearby low-luminosity blue stars. Aside from the O I 616 nm and 778 nm absorption lines, the spectrum of SDSS 0922+2928 (Fig. 1) is typical of a cool DQ white dwarf, with strong C$_2$ Swan bands in the blue, superimposed by a number of atomic carbon lines. In contrast SDSS 1102+2054 is a unique white dwarf with a photospheric spectrum totally dominated by absorption lines of O I (Fig. 1). The only other noticeable features are some weak lines of atomic C I as well as H$\alpha$ and H$\beta$.

We determined the effective temperatures of the two white dwarfs by fitting their SDSS $ugriz$ photometry with two grids of theoretical atmosphere models: one for a pure He composition \[Supporting Online Material (SOM), Sect. 1\], and one including H, C, and O as well (SOM, Sect. 2). The white dwarf temperatures and chemical abundances determined from these fits depend only mildly on the assumed surface gravity $\log g$ and, because their masses are not constrained by the available observations, we fixed $M=0.6M_\odot$, corresponding to the canonical value of $\log g=8$ \[, SOM Sect. 2\].

Independent of the details of the atmospheric composition, the photometry suggests that both white dwarfs have temperatures in the range 8000–9000 K (Table 2). The implied distances depend strongly on the assumed surface gravity via the radius, here we used the mass-radius relation of . We used the effective temperatures (and assumed surface gravity) to calculate synthetic spectra for a final determination of atmospheric abundances from the absorption lines in the SDSS spectroscopy (SOM, Sect. 2).

For SDSS 1102+2054, we adopted $T_\mathrm{eff}=8150$ K and $\log
g=8$ and iterated the abundances of H, C, and O, until we obtained a reasonable fit to the line strengths. Because this was not possible for all oxygen lines we concentrated on the strong lines at 778 and 845 nm, which originate from the lowest lying levels of all optical lines. These lines depend only weakly on temperature and are best suited for a determination of abundances. The best compromise was achieved for the abundances in Table 3.

However, for this model the carbon multiplet near 712 nm and the 834 nm line, H$\alpha$, and many O I lines from higher levels are much too broad, indicating too high a pressure in the atmosphere. In addition, O I lines from the higher levels are much too weak. The most likely solution of this discrepancy is to assume a higher temperature, which increases the strength of the high-excitation lines compared to the strong multiplets (SOM, Sect. 2). A reasonable compromise was achieved with $T_\mathrm{eff}=10500$ K, $\log g = 8$ (Table 3). The remaining differences in the absorption lines between the model and observation are likely within the uncertainties of the atomic data and in particular the van der Waals broadening constants (SOM, Sect. 2).

This model, however, is in significant disagreement with the temperature derived from the SDSS photometry. One plausible cause for the discrepancy between the temperatures determined from our fits to the photometry and to the spectroscopy are uncertainties in the ultraviolet opacities of an atmosphere of such peculiar composition, which strongly affect the slope of the optical continuum \[, SOM Sect. 1\].

The atmospheric parameters of SDSS 1102+2054 thus remain somewhat uncertain. Nevertheless, when varying the parameters within the range demanded by photometry, line strengths and widths, the C and O abundances did not change more than factors of $\sim3$. Thus the O/C abundance ratio in SDSS 1102+2054 is clearly much larger than one, which makes this star unique among the many thousands of known white dwarfs.

Because, according to the photometry, the effective temperature for the two white dwarfs cannot be very different, the carbon abundance in SDSS 0922+2928 must be higher in order to produce the strong C$_2$ Swan bands. On the other hand, the hydrogen abundance must be lower, because no Balmer lines are visible. We fixed the H/He ratio at $10^{-5}$, the upper limit allowed by H$\beta$. Further reduction did not substantially influence the models because the major electron donors are C and O. Using $T_\mathrm{eff}=8270$ K from the fit to the photometry, and assuming $\log g=8.0$, we obtained a good fit to the spectrum with the oxygen abundance decreased by $\sim0.2$ dex and a carbon abundance increased by $\sim0.6$ dex with respect to SDSS 1102+2054 (Table 3). Therefore, also in SDSS 0922+2928, $\mathrm{O/C}>1$. With the exception of the strong O I multiplets 778 and 845 nm most of the higher excitation O I lines and the C I lines show the same problem as in SDSS 1102+2054: they are too weak and too broad. This is again much improved at a hotter temperature, e.g. 10000 K, which is, however, in conflict with the photometry (SOM, Sect. 1).

The low H abundance suggests that SDSS 0922+2928 and SDSS 1102+2054, similar as the other classes of hydrogen-deficient white dwarfs, underwent a late shell-flash, leaving a He-dominated atmosphere. At the low temperatures of these stars, the He convection zone extends sufficiently deep to dredge core material up into the atmospheres . Given the large age of the two stars ($>500$ Myr), gravitational diffusion will unavoidably lead to a larger concentration of carbon in the envelope. The only plausible explanation for the observed O/C abundances is that these two white dwarfs have overall very low carbon mass fractions, and hence represent naked ONe cores. As such, they are distinct from 1H1504+65 and the "hot" DQ white dwarfs.

Most stellar models that produce ONe cores predict a layer of CO surrounding the core that should be sufficiently thick to avoid upward diffusion of large amounts of oxygen. However, a sequence of stellar evolution calculations approaching the mass boundary of stars forming ONe white dwarfs and those undergoing electron-capture SN II show that the mass of the CO layer decreases with increasing initial stellar mass . Thus SDSS 0922+2928 and SDSS 1102+2054 may have descended from the most massive stars avoiding core-collapse, in which case they would be expected to be very massive themselves. Our current data is insufficient to provide any unambiguous measure of the masses of SDSS 0922+2928 and SDSS 1102+2054, however their C/He abundance is at the top of the range observed in the carbon-rich sequence of "cool" DQ white dwarfs, some of which have parallax measurements implying that they are of high mass.

We explored spectral models for SDSS 1102+2054 with the abundances kept fixed as in Table 3, but adopting surface gravities of $\log
g=8.5$ and 9.0, which correspond to masses of $0.9\,M_\odot$ and $1.2M_\odot$ (SOM, Sect. 2). Broadly similar fits to the absorption lines can be achieved for higher surface gravities if the temperature is increased by 1000 to 2000 K as well. For $\log g=9.0$, the strongest O I lines become somewhat too broad compared to the observations, and we conclude that the currently available data is consistent with masses of up to $\sim1M_\odot$.

Initial models of the evolution of intermediate mass stars predict that ONe cores should also contain significant amounts of magnesium. Updated nuclear reaction rates have lead to a substantial downward revision of the Mg abundances . For SDSS 1102+2054, which has the better quality spectrum of the two white dwarfs presented here, we can place an upper limit on the magnesium abundance of $\log[\mathrm{Mg/He}]<-6.1$ from the absence of the Mg II 448 nm line.

<embed src="1180228Fig1.eps" style="width:10cm" />

**Fig. 1. Sloan Digital Sky Survey spectroscopy.** The observed spectra of SDSS 0922+2928 (bottom) and SDSS 1102+2054 (top) are shown as thick gray lines. Superimposed as thin black lines are models for $T_\mathrm{eff}=10000$ K and $\log g=8$ (bottom) and $T_\mathrm{eff}=10500$ K and $\log g=8$ (top). The wavelengths of the strongest OI and CI absorption lines, as well as of H$\alpha$ and H$\beta$ are indicated by tick marks.

|           |  SDSS 0922+2928  |  SDSS 1102+2054  |
|:----------|:----------------:|:----------------:|
| RA(2000)  |   09 22 08.19    |   11 02 39.69    |
| Dec(2000) |   +29 28 10.9    |   +20 54 39.4    |
| $u$       |  $18.77\pm0.02$  |  $17.58\pm0.02$  |
| $g$       |  $18.52\pm0.02$  |  $17.24\pm0.01$  |
| $r$       |  $18.47\pm0.02$  |  $17.26\pm0.01$  |
| $i$       |  $18.58\pm0.02$  |  $17.34\pm0.01$  |
| $z$       |  $18.71\pm0.03$  |  $17.46\pm0.02$  |
| p.m.      | 0.275" yr$^{-1}$ | 0.163" yr$^{-1}$ |

|  |  |  |  |  |
|:---|:--:|:--:|:--:|:--:|
| System | $T_\mathrm{eff}$ \[K\] |  | distance \[pc\] |  |
|  | pure He | with metals | pure He | with metals |
| SDSS 0922+2928 | $8720\pm260$ | $8270\pm320$ | $141\pm6$ | $122\pm7$ |
| SDSS 1102+2054 | $8820\pm110$ | $8150\pm150$ | $81\pm2$ | $67\pm2$ |

| System ($T_\mathrm{eff}/\log g$) | $\log$(H/He) | $\log$(C/He) | $\log$(O/He) | $\log$(O/C) |
|:---|:--:|:--:|:--:|:--:|
| SDSS 1102+2054 (8150/8.0) | -3.2 | -3.2 | -1.8 | 1.4 |
| SDSS 1102+2054 (10500/8.0) | -4.0 | -3.6 | -1.8 | 1.8 |
| SDSS 0922+2928 (8270/8.0) | -5.0 | -2.6 | -2.0 | 0.6 |

*Supporting Online Material for*
Two white dwarfs with oxygen-rich atmospheres

B.T. Gänsicke$^1$, D. Koester$^2$, J. Girven$^1$, T.R. Marsh$^1$, D. Steeghs$^1$

$^{1}$ Department of Physics, University of Warwick, Coventry CV4 7AL, UK

$^{2}$ Institut für Theoretische Physik und Astrophysik, University of Kiel, 24098 Kiel, Germany

We describe here in more detail the white dwarf model atmospheres that we have computed for the analysis of the SDSS data of SDSS 0922+2928 and SDSS 1102+2054, and outline the approach of our spectral fits. Two grids of atmosphere models were calculated. The first grid consisted of pure helium models, used to provide a first estimate of the effective temperatures of both white dwarfs from fitting their SDSS photometry. The second grid, including H, C, and O, was used to refine the temperature estimate and to determine the atmospheric abundances from fitting the observed absorption lines in the SDSS spectra.

# He-model fits to the SDSS photometry

In order to determine a first estimate of the effective temperatures of the two white dwarfs, we fitted the He-models to the observed $ugriz$ magnitudes, which were corrected to the AB scale by using the corrections suggested by the SDSS project (-0.04, 0.00, 0.00, 0.01, 0.02)[^1]. Because the mass of the two white dwarfs is currently unconstrained by observations, we fixed the surface gravity to the canonical value of $\log g=8.0$. Good fits were achieved for both objects, implying temperatures in the range 8000–9000 K (Table 2, Fig. S1). We note that for both objects the SDSS data exhibits a discrepancy between the $u$-band flux implied by the photometry, and the extrapolation of the flux-calibrated spectroscopy. As part of our search for white dwarfs with oxygen-rich white dwarfs, we have inspected the SDSS data of a large number of objects, and find that such a discrepancy is rare, but not truly exceptional, and based on the available data, it is not possible to give preference to either the photometry or spectroscopy. Taking the two spectral energy disributions at face value, a fit to the photometry is bound to result in a lower temperature compared to a fit to the spectroscopy. We will come back to this point below.

# Fitting the SDSS spectroscopy with variable H,C, and O abundances

Since He is practically neutral at the temperatures estimated above, the three other elements (H, C, O) provide a significant electron density, which correspondingly changes the opacity and thermal structure of the atmospheres. The equation of state included 6 diatomic molecules (H$_2$, CH, OH, C$_2$, CO, O$_2$). In particular the very stable molecule CO is important, because it reduces the atomic carbon abundance considerably. Also included in the second grid was the continuous absorption of C and O using cross sections obtained from the TOPBASE database of the Opacity Project . The necessary O I and C I line data were obtained from the VALD database and from the Kurucz list (CDROM23). Many of the O I lines from Kurucz list that were prominent in the model were not visible or much weaker in the observed spectra, indicating that some of this line data is wrong. We therefore gave preference to the VALD data, where the problematic lines from the Kurucz list are either absent, or have $gf$ values several orders of magnitude smaller. Recent calculations confirm the data for atomic oxygen in the databases we used.

There are a number of O I lines where the core electrons are not in the ground state ${}^4S{}^\mathrm{o}$ configuration, but in the ${}^2D{}^\mathrm{o}$ or ${}^2P{}^\mathrm{o}$ state. In these cases the core electrons (without the excited electron which makes the transition) are not in the configuration of the O II ground state and the ionization therefore would occur to an excited level of the O II ion. This has to be taken into account in our approximate calculation of the van der Waals broadening constants, which depend on the distance of the levels to the ionization limit. We use the ionization limit of that series, that is 16.942 or 18.636 eV instead of the normal 13.618 eV. From that limit and the term levels we estimate the broadening constant using standard approximations . Similar approximations are also used for the Stark broadening constants.

We re-fitted the SDSS $ugriz$ photometry with H, C, and O abundances close to the final values, which resulted in slightly lower temperatures for the two stars compared to the pure He-models (Table 2). In a final step, we varied the abundances of H, C, and O while keeping the effective temperatures fixed (8270 K for SDSS 0928+2928 and 8150 K for SDSS 1102+2054) and assuming $\log
g=8.0$. In this process, we focussed on the strong O I lines at 778 and 845 nm, which originate from the lowest lying levels of all optical lines and which only mildly depend on the prevailing effective temperature. It turned out that no abundance pattern would provide a fully satisfying fit to all observed lines, and we consider the abundances listed in Table 3 to represent the best compromise.

A specific problem with this model is that the O I lines from the higher levels are much too weak. One way to reduce the pressure in the atmosphere models is to reduce the amount of He. We have calculated a series of models between 7000 and 11000 K with the O/He ratio increased up to values $>10$. In all cases the O I lines are much too strong, and we rule out this possibility. The most likely solution to this problem is then to assume a higher temperature, which increases the strength of the high-excitation lines compared to the strong multiplets. Assuming, as before, a canonical $\log g = 8$, a reasonably good fit to the absorption lines was achieved with $T_\mathrm{eff}=10500$ K (Fig. 1). However, the slope of this model is substantially bluer than the SDSS $ugriz$ photometry (Fig. S1). A fundamental uncertainty in modelling white dwarf spectra with such unusual abundances is our limited knowledge on the sources of ultraviolet opacity in their atmospheres. Any flux absorbed at ultraviolet wavelengths will be re-distributed to longer wavelengths, and hence result in a steepening of the slope in the optical bandpass. An example of the marked differences in the amount of ultraviolet opacity between pure He and C atmospheres has been illustrated by . Our models predict strong ultraviolet absorption edges of O I and C I (which we include with the values obtained from TOPBASE), and a strong L$\alpha$ line broadened mostly by neutral helium (which we describe with the theory developed by our group). The currently available data does not provide quantative constraints on the ultraviolet characteristics of SDSS 0922+2928 and SDSS 1102+2054, and hence a systematic uncertainty on their effective temperatures remains. We conclude that, within the framework of the spectroscopic analysis, it is reasonable to prefer the spectroscopic solution of the higher effective temperatures, as it is impossible to satisfactorily model the observed ratio of O I lines with high and intermediate excitation levels at the low (photometric) temperature.

As discussed in the main manuscript, the large O/C abundance ratios found for these two white dwarfs strongly suggests that they represent bare oxygen-neon cores that may have descended from massive progenitor stars. Under this assumption, their masses based on any of the current initial mass-final mass relations are expected to be $\sim1\,M_\odot$, corresponding to a likely range for the surface gravity of $8.5\le\log g\le9.0$. In order to explore the effect of a higher gravity on the spectral fit, we have calculated a grid of models spanning the range $\log g=8.0$ to 9.0. Since the spectrum is largely dominated by the ionization balance of He and the other elements, a change of $\log g$ can to some extent be compensated by a change in temperature, and to first order the abundances do not change significantly. In fact, we do find that broadly similar fits can be achieved when increasing $T_\mathrm{eff}$ along with $\log g$, but keeping the abundances fixed at values given in Table 3. Figure S2 illustrates the fits with our canonical model ($T_\mathrm{eff}=10\,500$ K, $\log g=8.0$), and models for higher surface gravities, ($T_\mathrm{eff}=11\,500$ K, $\log g=8.5$) and ($T_\mathrm{eff}=12\,500$ K, $\log g=9.0$). For $\log g=9.0$, the strongest O I lines become somewhat too broad compared to the observations. We conclude that the currently available data are consistent with a broad range of white dwarf masses, including relatively massive ($\sim1\,M_\odot$) stars, and that the main conclusion of a large O/C abundance ratio is independent of the mass assumed for the spectral fitting.

<embed src="1180228FigS1.eps" style="width:12cm" />

**Supporting Figure 1.** SDSS spectroscopy (gray lines) and photometry (green dots) of SDSS 1102 +2054 (bottom panel) and SDSS 0922+2928 (top panel). The uncertainties in the fluxes from the SDSS photometry are comparable to the size of the dots. In both cases, the $u$-band flux from the photometry lies below the extrapolation of the spectroscopy. The best-fit models to the photometry for an assumed $\log g=8.0$ are shown as blue lines ($T_\mathrm{eff}=8270$ K for SDSS 0922+2928 and $T_\mathrm{eff}=8150$ K for SDSS 1102+2054; see Tables 2 and 3). Shown as red lines are the best-fit models to the absorption lines for an assumed $\log g=8.0$ ($T_\mathrm{eff}=10\,000$ K for SDSS 0922+2928 and $T_\mathrm{eff}=10\,500$ K for SDSS 1102+2054). As part of the fitting the absorption lines, the models were normalised to the observed SDSS spectra. Here, we show the best-fit models scaled to the $r$-band magnitude without that normalisation to illustrate the systematic uncertainty in the optical slope/effective temperature.

<embed src="1180228FigS2.eps" style="width:12cm" />

**Supporting Figure 2.** Model fits to the observed absorption lines of SDSS 1102+2054, illustrating the effect of varying the surface gravity. The central panel shows the full wavelength range of the SDSS spectroscopy, with a model for ($T_\mathrm{eff}=10\,500$ K, $\log
g=8.0$) and the corresponding abundances from Table 3 superimposed on the SDSS data. Note that the continuum slope is normalised to the observed data, and that only the absorption lines are fitted. Offset by two and four flux units are spectra with ($T_\mathrm{eff}=11\,500$ K, $\log g=8.5$) and ($T_\mathrm{eff}=12\,500$ K, $\log g=9.0$), respectively, keeping the abundances fixed to those given in Table 3. The top panels show close-ups of the normalised SDSS spectrum around several O I and C I lines, as well as of H$\alpha$. Line identifications are indicated by tickmarks. A model for ($T_\mathrm{eff}=10\,500$ K, $\log g=8.0$) is superimposed on the observed spectrum, models with ($T_\mathrm{eff}=11\,500$ K, $\log
g=8.5$) and ($T_\mathrm{eff}=12\,500$ K, $\log g=9.0$) are offset by $-0.4$ and $-0.6$ flux units, respectively.

[^1]: http://www.sdss.org/dr7/algorithms/fluxcal.html#sdss2ab


abstract: Traceroute is the main tools to explore Internet path. It provides limited information about each node along the path. However, Traceroute cannot go further in statistics analysis, or *Man-Machine Interface (MMI)*.
  .
  Indeed, there are no graphical tool that is able to draw all paths used by IP routes. We present a new tool that can handle more than 1,000 Traceroute results, map them, identify graphically MPLS links, get information of usage of all routes (in percent) to improve the knowledge between countries' links. rTraceroute want to go deeper in usage of atomic traces. In this paper, we will discuss the concept of rTraceroute and present some example of usage.
author:

bibliography: metrology.bib
title: rTraceroute: Réunion Traceroute Visualization

# Introduction

## Context

Since the beginning of the Internet, the network has evolved from a simple graph with few vertices to a wide and a complex graph with uncountable number of nodes.

To improve our understanding of the Internet routing and its topology, active metrology is the only way possible with the injection of packets in the network. Traceroute is one of the most widely used network measurement tools. It reports an IP address for each network-layer device along the path from a source to a destination host in an IP network. Network operators and researchers rely on Traceroute to diagnose network problems and to infer properties of IP networks, such as the topology of the Internet.

Reunion Island is a small territory located in the Indian Ocean, near South Africa, off the eastern coast of Madagascar. However it is also one of the overseas departments of France. For Reunion Island, this issue is more important than other places due to its particular connectivity to the Internet based on two submarine cables, called *South Africa Far-East (SAFE)* and *Lower Indian Ocean Network (LION) II*, respectively in green and in red in figure . If in one way, we know the physical path of the Internet connection, in the other way, no one knows the logical path of our connection. This is why we want to graph the links between Reunion Island and the world-wide destination. This new representation will improve the global knowledge of the local topology.

## Project

The council of Region Reunion have financed a metrology project, with *European Regional Development Funds (ERDF)* to identify the difficulties of the Reunion Island connectivity. The requirements for this project are:

1.  To identify the first hop after any connection leave Reunion Island and the last hop before it reaches us. This will help on the characterization of the logical exits and entries for Reunion Island.

2.  To detect the *MultiProtocol Label Switch (MPLS)* link between Reunion Island and the different hops. It will help to understand how are the exits and entries for Reunion Island.

3.  To focus on the regional peering problem. This item will allow us to see if the physical cable and the routing policies could be superposed.

4.  To identify the most encountered node after (resp. before) leaving (resp. reaching) Reunion Island. This could improve our knowledge of the routing rules of all *Internet Service Providers (ISP)* present on the territory.

5.  To measure the minimal delay associated with each links between Reunion Island and the next (resp. previous) hop. It could detect prioritized link based on the final (resp. original) destination (resp. source).

However identifying these points are not easily obtained with small data-set. We need to analyze quite big data-set, more than one probe for each *ISP* present on the Island. Usually, each network tool generates its own measurements before creating an output in text and/or image format. It is difficult to combine data from several sources (resp. destination) from (resp. to) a same country without a platform measurement.

Moreover, before running an analysis we need to proceed methodically as follows:

1.  Cleaning data-sets.

2.  Matching data-sets from different tools, with different standards.

3.  Adding Geo-localization.

4.  Showing result on a map.

To help solve this problem, we propose *rTraceroute*, a tool which is "easy to use" but can provide a map together with statistics by analyzing data from different Traceroute sources. This tool take only one parameter: a folder which contains data from Traceroute tools. It gives us then two output: the first one is a statistic file about each node meet and the second is a map with the different link between countries.

Thus, rTraceroute allows us to save time and to focus on the analyzing of outputs. The rest of this paper is organized as follows. At first, we will describe the some existing tools in section . The section  describes rTraceroute. An example of how it works is presented in  and it will be followed by a conclusion in section .

# Existing tools

As things stand today, no tool is available to fill our needs. However, we can find some tools that could answer some of them. Indeed, there are tools for the identification of the path from a probe to a destination. This is the reason why we have created *rTraceroute* to fill the need from this project.

In this section, we will present several tools that could help us. We have categorized them as follows:

1.  Tools that give raw results, without any statistics or map.

2.  Platforms measurement that can take advantage of the previous tools and use them together.

3.  Mapping tools that can draw their result on a map.

## Raw Results

These tools would provide us data from one measurement, but they have not been built to generate a large raw data without a script.

### Traceroute

Traceroute  is the most popular tool to measure path length and to determine the routes that packets follow from a source node to a destination node. With the growth of the Internet, many paths can be found with Traceroute. But the main flaw of this tool is that it could only generate raw data (without any statistic or mapping).

### Paris-Traceroute

Paris-Traceroute  is a Traceroute-like tool, but it is less sensitive to the load-balancing phenomena. The output of this tool announces some errors in the paths, like *Unreachable Network (!N)* or *Unreachable Host (!H)*. The using of *explicit MultiProtocol Label Switch (MPLS)*  link is also provided by the tool. Using Paris-Traceroute in our work could only identify the MPLS link but it has the same flaws as Traceroute.

### TraIXceroute

TraIXceroute  is also a Traceroute-like tool with the discovering of *Internet eXchange Point (IXP)* in the paths. This version of Traceroute combines two databases, *PeeringDB*  and *Packet Clearing House* , for the identification of these particular nodes. This information could help for the reconstruction of the country crossed. Despite an accuracy of 92-93%, this tool is not adapted for our study because it can only detect the regional routing problem with the IXP databases .

### Reverse-Traceroute

Asymetric path is frequent in the Internet, and it's very difficult to identify the reverse path. To avoid this,  has presented *Reverse-Traceroute*. This tool aims to identify the return path of a classic Traceroute, with an accuracy of 87% for hops identification. It's mainly based on two IP options, *IP Record-Route* and *IP Timestamp option*. Like the original Traceroute, reverse-Traceroute can only generate Traceroute raw data but for the return way. Another method to obtain the reverse path is using a measurement platform.

## Measurement Platform

According to , an Internet measurement platform is an infrastructure of dedicated probes that periodically run network measurement tests on the Internet.

### Atlas RIPE NCC

Atlas RIPE NCC  is an active measurement platform, which allowed to use several active tools, like *ping, DNS, HTTP* or *NTP* for example. There are two different tools for path measurement. The default one is Paris-Traceroute. The second one is the original Traceroute. Atlas platform return the results in *JavaScript Object Notation (JSON)* format. It is possible also to map the path with *OpenIPMap* . Despite the fact Atlas could use Paris-Traceroute, made mapping with geo-localization and compare large data produced by its measurement, this platform is not adapted to our project because it lacks the statistics part.

### Planet-Lab

In 2003, some American researchers deployed a test-bed platform called Planet-Lab . In 2008, a European portion of Planet-Lab  have been created. On Planet-Lab nodes, researchers can install the tool they need for their works. With this politic, this platform could combine several tools. Nevertheless, it is unable to identify the last (resp. first) hop after (resp. before) reaching (resp. leaving) Reunion Island. Moreover, the management of large data is not allowed on PlanetLab.

### Archipelago

*Center for Applied Internet Data Analysis (CAIDA)* has deployed its own platform in 2007, called Archipelago . The probe connected to this platform allows us to make five main measurements, including *ping* and *Traceroute*. But Archipelago, as Traceroute, could only generate raw data without any analysis.

## Mapped results

### GTrace

The first tool created to graph Traceroute data was GTrace . This tool generates its own Traceroute data and draw it on a map. It works with city name abbreviations or airport codes, lookup client and two IP databases to validate each IP addresses location. This tool was not maintained for a long time and doesn't support the x86 architecture. Due to this, we were unable to test it.

### Topology Visualization Tool

In , a tool centralized on analysis topology in Africa was presented. This tool can generate its own Traceroute measurement and graph it after requesting information about nodes from several databases. The design is largely based on existing visualization tools, such as *OpenIPMap project*, and focuses only on the African continent. The source code of the tool was not indicated in the article, thus no test has been made.
The table  resumes the possibilities of each tool presented previously. The cells in gray represent the available options for each tool, when the white cells mean that the option is not present.

```latex
\begin{table*}[ht!]\centering
\caption{Tools option availability}
\label{tab:tools}
\resizebox{\textwidth}{!}{ %
\begin{tabular}{|c|c|c|c|c|c|c|c|c|}
\hline
\textbf{Need} &  &  &  &  &  &  &  &  \\ \cline{1-1}
\textit{Tools} & \multirow{-2}{*}{\textbf{Last Hop identification}} & \multirow{-2}{*}{\textbf{MPLS Mark}} & \multirow{-2}{*}{\textbf{Statistics}} & \multirow{-2}{*}{\textbf{Geo-localization}} & \multirow{-2}{*}{\textbf{Mapping}} & \multirow{-2}{*}{\textbf{Database}} & \multirow{-2}{*}{\textbf{Regional Routing problem}} & \multirow{-2}{*}{\textbf{Large data}} \\ \hline
\textit{Traceroute} &  &  &  &  &  &  &  &  \\
\hline
\textit{Paris-Traceroute} &  & \cellcolor[HTML]{9B9B9B}YES &  &  &  &  &  &  \\
\hline
\textit{TraIXroute} &  &  &  &  &  & \cellcolor[HTML]{9B9B9B}YES & \cellcolor[HTML]{9B9B9B}YES &  \\
\hline
\textit{Reverse-Traceroute} &  &  &  &  &  &  &  &  \\
\hline
\textit{Atlas RIPE NCC} & \cellcolor[HTML]{9B9B9B}YES & \cellcolor[HTML]{9B9B9B}YES &  & \cellcolor[HTML]{9B9B9B}YES & \cellcolor[HTML]{9B9B9B}YES & \cellcolor[HTML]{9B9B9B}YES &  & \cellcolor[HTML]{9B9B9B}YES \\
\hline
\textit{PlanetLab} & \cellcolor[HTML]{9B9B9B}YES & \cellcolor[HTML]{9B9B9B}YES &  &  &  & \cellcolor[HTML]{9B9B9B}YES & \cellcolor[HTML]{9B9B9B}YES &  \\
\hline
\textit{Archipelago} &  &  &  &  &  &  &  &  \\
\hline
\textit{Gtrace} & \cellcolor[HTML]{9B9B9B}YES &  &  & \cellcolor[HTML]{9B9B9B}YES & \cellcolor[HTML]{9B9B9B}YES & \cellcolor[HTML]{9B9B9B}YES &  &  \\
\hline
\textit{African Visual Route} & \cellcolor[HTML]{9B9B9B}YES &  &  & \cellcolor[HTML]{9B9B9B}YES & \cellcolor[HTML]{9B9B9B}YES & \cellcolor[HTML]{9B9B9B}YES &  &  \\
\hline
\textit{rTraceroute} & \cellcolor[HTML]{9B9B9B}YES & \cellcolor[HTML]{9B9B9B}YES & \cellcolor[HTML]{9B9B9B}YES & \cellcolor[HTML]{9B9B9B}YES & \cellcolor[HTML]{9B9B9B}YES & \cellcolor[HTML]{9B9B9B}YES & \cellcolor[HTML]{9B9B9B}YES & \cellcolor[HTML]{9B9B9B}YES \\
\hline
\end{tabular}%
}
\end{table*}
```

The header line of the table show our needs:

-   The identification of the last hop corresponds to the country where the node is hosted before leaving (resp. joining) a specific country. Only a graph tool can identify this point.

-   MPLS mark detect the *explicit MPLS* link between two nodes. Only Paris-Traceroute can detect these, when our rTraceroute with the graph part can detect *invisible MPLS* route.

-   Statistics part needs a large amount of data for the analyze. In general, the raw results generator does not provide any analysis. Even graphic tools can't do it because it redraws on the map on each new test. rTraceroute can analyze a large dataset and make some statistics one each meeting node.

-   The geo-localization of each IP address is based on the declaration of the RIR plus the ISP. It's well known that the Geo-localization is not accurate science. rTraceroute propose a database which could be improved by the user analysis.

-   In the general case, mapping tools need to generate their own data before draw them on a map. The first problem concerns the choice of the map and the possibility to download the final result. With rTraceroute, you can choose your own map as input and download the results for any publication or presentation. Mapping data allows us to identify the different countries pass through and detect the potential errors. With rTraceroute, these phenomena could be easily detected by the mapping and the statistics part.

-   The database is mainly used for Geo-localization of IP address or detection of the IXP. Our database is, for now, focus on Reunion Island with the help of delay analysis.

-   The regional routing problem means the identification of peering mistakes, like boomerang-routing phenomenon .

-   Each tool provide to generate its own raw data one by one. rTraceroute propose to analyze a large dataset for the statistics and mapping part.

Due to the imperfections of these tools (relatively to our need), we have made the choice to develop a new tool: *rTraceroute*.

# rTraceroute

## The concept

rTraceroute have been created to simplify the Traceroute data file analysis.

On the one hand, we must process a huge number (more than 1 million) of Traceroute files. On the other hand, we are looking for performance. The choice to handle this is to write the program with standard C, keeping in mind that it should be compatible (*Linux, OS X*, ...), without any platform-specific code. We also want the program to be easy to compile (this is why the code is a monolithic bloc, for less than three thousand lines), with a very simple Makefile. Today, the executable run on a dedicated server with *Linux Ubuntu 14.0.4 LTS*.
The figure  schematize the concept of rTraceroute. The inputs used by rTraceroute are :

-   Several files of Traceroute data, whether it comes from Paris-Traceroute  or Atlas platform , in text format.

-   The URI of the database used for geo-localization.

-   A Map.

Then, rTraceroute will proceed as follows:

1.  **Parsing and Cleaning**: As the input data from Paris-Traceroute  or Atlas platform  are text files, the program must parse JSON files as well as plain text files. The parse procedure for JSON files is made with NXJson . To be more efficient, the program has a "cleaning procedure" (as it has been explained in  ) that eliminates useless (or corrupted) files. The pre-parsing subroutine uses the following exclusion pattern:

    -   JSON corrupted.

    -   JSON file contains 3 last RTT as `0` (zero) \[unreachable\].

    -   Traceroute corrupted \[file contains more than one Traceroute\].

    -   Traceroute contains 3 stars on the last hop \[host unreachable\].

    -   Traceroute contains `!H` or `!N` or `WARN` string \[unreachable\].

2.  **Statistics**: During parsing files, rTraceroute create an internal linked-chain memorizing all information. Each IP address and hop are marked with an on the fly calculated RTT; later completed with the localization (country) from the MySQL internal database.

3.  **Geo-localization**: From this step, we need to geo-localize IP address. This is why the program needs a MySQL database (connection information is "hard-coded" in main.h). . Two tables are used: the first one to geo-localize IP address and the second one to place the country on a map.
    Then, to be able to draw each trace, *gdlib* has been used.
    To store all information obtained during the process, three kinds of tuples are used:

    -   `{hop, ip address, rtt, occurrences}` is added when reading files.

    -   `{mapx1, mapy1, mapx2, mapy2, ip_address1, ip_address2, link}` is created when drawing is done, to handle all MPLS links.

    -   `{xdeb, ydeb, xfin, yfin, occurences, rtt, mpls}` is also handled when the program draws.

4.  **MPLS Detection**: *Paris-Traceroute* can detect *explicit MPLS* links and marked them in its raw data. When rTraceroute find the marks, it memorized the link to plot them in a different coloring at the end of the mapping.

5.  **Mapping**: Some useful options have been implemented as the program can handle html color's name and create statistics:

    -   It can draw only the last links (the last hop before coming in one country).

    -   It can make the thickness of the link proportional of its use (with the `-redraw` option).

    -   It colorize MPLS links.

The program produces three result files:

-   A simple one that contains only 5 columns: `hop_position, IP address, occurrence, average_delay, country`, used for a later statistical purpose.

-   A very verbose one called `trace.txt` that explicitly details all internal operations. From this text file, it is possible to extract and produce more results (as statistics in percent of use's link). for examples:

    -   to extract all MPLS links: ` $ grep 'lien.*MPLS' trace.txt ` and we obtain this kind of result: ` // lien <-> MPLS: 22.214.171.124 188.8.131.52 `

    -   to find information on each drawing line (x1,y1 to x2,y2), occurrences, occurrence in percent between all lines, minimal time, geo-localization, in text, for all segments: ` grep '__ trajet' trace.txt ` and we obtain this kind of result: ` __ trajet: 3626 1638 -> 2581 1582: 105 (2.81%) [AU - RE et temps min 505.39 (2.81%)] `

-   Two graphical files (maps): all segments of the input files on the first one and only the segments of the last hop for a destination point (x,y) on the last map.

rTraceroute is the most adapted tool for our project. Indeed, its capacity of mapping help us to answer the first need of the project, which is the identification of the first (resp. lat) hop after (resp. before) leave (resp. reach) Reunion Island. The map function also help us on the third item of the project: the regional peering problem. Paris-Traceroute allow to detect *explicit MPLS* link. With rTraceroute, even the *invisible MPLS*  link is shown with the maps option. The statistics part of rTraceroute permits to answer the last two points, which are respectively the identification of the most encountered node after (resp. before) leaving (resp. reaching) Reunion Island and the minimal delay associated with each links between Reunion Island and the next (resp. previous) hop.

# Application of rTraceroute on Reunion Island

This paper has begun with a quick presentation of different active tools or platforms. This part will show how rTraceroute used the data provided by these tools for analysis path from a country to several destinations. The section  explain how rTraceroute can be used with direct Paris-Traceroute data, when  used Paris-Traceroute JSON formatted from Atlas RIPE NCC platform as data. The results have been previously analyzed in  and are available on the website of the authors.

## Paris-Traceroute dataset

In this example, we will test rTraceroute on the data-set publicly available and coming from . This large data was obtained during one month, between third of July and third of August 2016. With 27 Raspberry Pi used as measurement probes, each day we have tried to reach $10,000$ IPv4 spread around the world as destination. At the end of the measurement campaign we have obtained a total of $6,458,025$ raw traces.
Combined with rTraceroute, the raw data have been reduced to $1,015,180$ available for path analysis. The using of rTraceroute on Paris-Traceroute data have these advantages:

-   The first one concern the raw data. The tools can filter only the data which fill the selection criteria presented in the section

-   Mapping the path makes the path analysis between countries easier. Along with mapping, it is possible to print only the first (resp. last) node before (resp. after) leaving (resp. reaching) a country with the path and the delay associated.

-   Lastly, rTraceroute can change all *explicit MPLS* link's color found by Paris-Traceroute. The mapping would permit the detection of *invisible MPLS* links.

The map  represent the results for the first hop when data leave Reunion Island.

In this figure, we can see the first hop after Reunion Island. Despite the presence of two physical connections with Asia and Africa , we can notice that most of the first hops are not connected to these continents but to Europe.

Another interesting point is the presence of a direct link between Reunion and USA, sign of *invisible MPLS* link. We can also see that more of 96% of our data going through France before joining the destination, even if the destination is close to Reunion Island. A detailed analysis of the results are available in .
The analysis of the same data with only a *BASH* script has taken more than one day. With rTraceroute, around 10 minutes is sufficient.

## Atlas RIPE NCC dataset

In this case study, we would analysis $300,000$ measurement towards Reunion Island from all world. Our trace includes measurement performed from the third of July to the third of August 2016. $1,000$ atlas probes would reach ten Raspberry Pi  deployed over Reunion Island. The atlas probe was selected with the same distribution as the IPv4 addresses for the Paris-Traceroute case.
Without rTraceroute, we would have been constraint by using OpenIPMap . The tools present some limits: it requires to add manually the other path measurements to compare the results and the maps can't be downloaded. The "cleaning step" of rTraceroute have reduced our data to $38,714$ traces. The advantage of combine rTraceroute with Atlas is the same as with the other Traceroute tools. But in this context, we can highlight three benefits:

-   First we can easily combine data from several Atlas math measurements for the analysis. A script could download the JSON file and add it in a folder which could also be analyzed by rTraceroute.

-   The maps are hosted on the system which run rTraceroute. Maps can be downloaded and used after.

-   Eventually, the statistics file created by rTraceroute is very important. Each line of the file contains the position of an IP address meet during the measurement, the IP address, the occurrence of this couple, the meaning delay and the country associated with the IP address.

The map  is an example of all path existing to join Reunion Island using Atlas and rTraceroute.

On the figure , we can see a direct link between Paraguay and Reunion Island without a submarine cable existing. This proves us that an *invisible MPLS* link exists. Without rTraceroute, this analysis will be more difficult to be found. We can also notice the use of the section of the submarine cable from Asia and Africa. The percentage of using the exit point for submarine cable in the two continent is less than 0.5%, when more than 99% going through France, with a minimal delay of $183.30$ms.
This campaign of measurement from the world towards Reunion Island has provided $300,000$ JSON raw data for a time analysis of 3 minutes with rTraceroute. A *BASH* script developed for the analysis of these raw data has taken more than half a day.

# Conclusion

Studying the Internet connectivity of Reunion Island is important for the economic development of the island. An *ERDF* project has been funded to identify the difficulties of the Internet connectivity of the Island. Several tools are available, but none of them could cover our needs. Nearly no tools can handle a large number of Traceroute-files and draw all the path on a map. This why we have implemented a new tool, called *rTraceroute*. It also offers us a new approach to observe the Internet topology with a different point of view. As for today, there is no IP geo-localization database reliable enough to plot every IP found in the traces. The next step will be to improve IP geo-localization of each point, based on delay between to point and a system of triangulation. IPv6 deployment is in its infancy: the routes could not be optimized. Adding IPv6 in rTraceroute could also help the deployment of this new protocol in the area, but some adaptations need to be done for IPv6 because the actual work is based on IPv4. Exportation of maps in *Encapsulated PostScript (EPS)* or *Scalable Vector Graphic (SVG)* format is also planned in the future (actually, the maps result are only in PNG format).

A version of rTraceroute is in public access at <http://lim.univ-reunion.fr/rTraceroute>