A Quantitative Analysis of Bulgarian Dialect Pronunciation∗

Petya Osenova,† Wilbert Heeringa and John Nerbonne Humanities Computing, University of Groningen [email protected], [email protected], [email protected]

January 4, 2007

Abstract We apply a computational measure of pronunciation dierence to a database of 36 word pronunciations from 490 sites throughout Stoykov's Bulgarian Dialect Atlases. The result is a comprehensive view of the aggregate pronunciation dierences among the 490 sites. This study aims to contribute therefore to Bulgarian dialectology, as well as to the development and testing of the computational technique, now implemented in a software package (L04).

1 Introduction

In recent years computational techniques enable the incorporation of large amounts of dialectal material into studies of language variation. Phonetic measurements of Germanic and Romance dialects (Dutch, Norwegian, Sardinian, German, American) have been conducted successfully using Levenshtein distance, also known as `(string) edit distance', as a basis. This line of work began in 1995 when Kessler introduced the use of the Levenshtein distance as a tool for measuring linguistic distances among the pronunciations of language varieties. Leven- shtein distance is a string edit distance measure, and Kessler applied this algorithm to the comparison of Irish dialects. Later the same technique was successfully applied to Dutch (Nerbonne at. al. 1996; Heeringa 2004, pp. 213278), Sardinian (Bolognesi and Heeringa, 2002), Norwegian (Gooskens and Heeringa, 2004), German (Nerbonne and Siedle, 2005), and American English (Nerbonne, to appear). Other analyses are underway. These studies have applied Levenshtein distance to large numbers of words as they are pronounced in many dierent data collection sites. They have, by and large, vindicated traditional dialectological divisions, even while systematizing

∗This work is funded by NWO, Project Number 048.021.2003.009, P. I. J. Nerbonne, Groningen, and also a grant from the Volkswagenstiftung Measuring Linguistic Unity and Diversity in Europe, P.I. E. Hinrichs, T¤ubingen. †Also at the Bulgarian Language Division, Faculty of Slavonic languages, St. Kl. Ohridski University, Soa, Bulgaria

1 the technique of comparison, enormously enlarging the base of data on which analyses are systematically based, and providing novel means of understanding the dialect landscapes to which they have been applied. It is a challenge to see how well the methods developed primarily for other language families perform for Bulgarian, a Slavic language. The task was chal- lenging: rst, because there was no digitized data for Bulgarian dialects at our disposal; second, because Bulgarian dialects had not been processed earlier with computational tools, and third, due to the fact that some language specic features had to be taken into account for the rst time. The structure of the paper is as follows: the next section describes the traditional divisions of Bulgarian dialects, providing background. Section 3 focuses on the data source and the preparation of the data. In Section 4 the Levenshtein distance measurement is presented. Section 5 discusses the analytical procedures applied to the distances computed in Section 4, namely clustering and multidimensional scaling. Section 6 illustrates the dialectal divisions by providing maps showing the distribution of three word pronunciations. Section 7 focuses on some of the dialect groups in more detail. Section 8 presents the relation of the dialects to standard Bulgarian pronunciation. Section 9 comments on the results obtained on an extended set of data, and Section 10 proposes conclusions.

2 Bulgarian Dialect Scholarship

Scholarship has oered several ideas about the geographical distribution of Bul- garian dialects, mostly relying on phonetic criteria. One prominent traditional division follows the pronunciation of the old Bulgarian vowel `yat'. It divides the Bulgarian dialects into western, where `yat' has only the reection `e' (for example bel `white'- beli `white-pl') and eastern, where `yat' has both reections, `e' and `ya' (for example bjal `white'-beli `white-pl'). This one characteristic is not by itself enough for consistent generalization, but the distinction remains one of the most important features for the comparison of the dialects. Another phonetics-based classication of Bulgarian dialects reects the realizations of the old Bulgarian `big nosovka' (ãîëÿìà íîñîâêà), a nasal vowel. It divides Bulgarian dialects into ve groups: @-dialects (northeastern and northwestern Bulgaria and the eastern part of southeastern Bulgaria); a-dialects (western Bulgaria and the eastern dialect of Pirdop); 6-dialects (the Rodopi mountain); æ-dialects (the Teteven region and two villages in eastern Bul- garia, Kazichino and Golitsa); and u-dialects (western Bulgarian areas near the Bulgarian-Serbian border). This classication is admirably simple but en- counters numerous exceptions, which mean that it cannot divide the Bulgarian dialects in a satisfactory way. There have been other attempts at the phonetic classication of Bulgarian dialects as well, but we are not going to present all of them here. The interested reader is referred to Stoykov (2002, pp. 8890). According to most morphological and lexical research Bulgaria is divided into a central part (northeastern and central Bulgaria) and a peripheral part (northwestern, southwestern and southeastern Bulgaria).

2 Because of the instability and conicting nature of various linguistic criteria Stoykov (2002) suggests a classication of Bulgarian dialects which respects geographical continuity, as well. In his standard work he distinguishes six, rather than ve areas, concluding that:

1. Bulgarian dialects are not separated categorically, but they rather form a continuum.

2. Within Bulgarian dialects there is a central (typical) area and peripheral (transition) areas. Similar situations have been observed for other languages as well, e.g. Dutch (Hoppenbrouwers and Hoppenbrouwers 2001, pp. 6667).

3. The most striking distinction of Bulgarian dialects is the distinction between eastern and western ones.

In Figure 1 the six most signicant geographical groups of Bulgarian dialects are shown as presented in Stoykov (2002, p. 416).

2.1 Related work P§seni§cnova (1973, 1977) applies statistical methods to the problem of classication of Russian dialects, contrasting these with traditional methods. P§seni§cnova (1973) introduces the notion distance among dialects, which she estimates via the χ2 metric applied to the frequency of phonetic features in dierent varieties. P§seni§cnova (1977) considers 100 dialects using 159 binary features which are not identied concretely, but which are said to be distinctive for the 100 dialects involved in the experiment. One example of a feature is whether the rst vowel in the word for window 'okno' is pronounced as an [a] or as an [o] (akan'e). She measures the similarity of dialects using feature overlap. Wherever two dialects have the same value for a feature, this contributes 1 to their similarity score, and where the features dier, this contributes −1. Uncertain information or inapplicable features contribute 0. This results in a similarity matrix showing the similarity among all pairs of dialects, which she then clusters into 5 groups. P§seni§cnova compares her results with the traditional division, concluding that it is in general the same, noting as well small dierences. She interprets clustering as indicating relative, but not absolute boundaries among dialects, defending a continuum view of the dialect landscape. In contrast to P§seni§cnova we employ a measure of pronunciation dierence below, allowing us to treating pronunciation dierences as numerical data. In addition to clustering, we also deploy multidimensional scaling to analyse the result of the aggregate distance (similarity) measurements. Finally, we focus here entirely on Bulgarian. Dobrev (1981) aims at locating the place of Bulgarian among the Slavic languages. He uses the 60 most frequent morphemes in Slavic texts, taken from a corpus covering 35 Slavic varieties (including dialects, styles, literary

3 4

Figure 1: The map of Bulgarian dialect divisions as presented in Stoykov (2002, p. 416). The vertical lines represent Moesian dialects, the horizontal lines represent Balkan dialects, the broken slanting lines southwestern dialects, the crosses northwestern dialects. The thick broken line represents the `yat' borderline that divides the dialects into two major groups: western and eastern. The nearly horizontal slanting lines on the left side show transitional zones, and the steeply slanting lines at the bottom of the map represent the Rupskian (including Rodopian) dialects. languages) and 700,000 word tokens of texts. Grouping is determined by the calculation of cumulative intersective quotas, which is eectively a clustering technique. Dobrev's results show that the eastern and western Slavic groups are fairly homogeneous, while the South Slavic group is not (at least not in comparison to the others). Our work resumes this line of work of Slavicists in the 1970's and 1980's. We will however, apply the computationally inspired Levenshtein distance to measure dialect dierences, automatically processing entire phonetic transcriptions rather than manually counting individual features. In addition to clustering, employed by both of the authors above, we shall submit our distance measurements to consistency analysis and multidimensional scaling.

3 The Data

We consider in turn the sources of our data and the selection we made, its preparation, and its conversion to digital form.

3.1 Sources The data was digitized from the four volumes of Bulgarian dialect atlases which cover the entire country area: Volume I - southeastern Bulgaria (Stoykov and Bernshteyn 1964), Volume II - northeastern Bulgaria (Stokyov 1966), Volume III - southwestern Bulgaria (Stokyov et al. 1975) and Volume IV - northwestern Bulgaria (Stoykov et al. 1981). The atlas materials appear to have been gathered with an eye toward the identication of the historical roots of Bul- garian. This would explain why the data was gathered only from villages with exclusively Bulgarian populations regardless of geography (this is unlike similar atlases for other languages we have worked with). This also means that the sites are not distributed uniformly within Bulgaria. For example, because most of the purer Bulgarian sites are in the mountains, there are more mountainous sites than one would expect. The atlases consist of two parts: maps and commentary on the maps. The maps present general information, while the commentary focuses on deviations and more elaborate characterizations of pronunciation, etc. The data was collected in some ways that promise an especially faithful characterization of language variation: the researchers that collected the data did not rely on only one in- formant, but instead used several, divided into main informants and additional ones. The researchers' approach to informants was not to ask direct questions, but rather to conduct extensive interviews from which material was selected. The data reects dierent linguistic phenomena, and the phenomena are often instantiated by more than one word. Thus material does not coincide across all atlases, and where it does coincide, it may not be presented as instantiating the same linguistic phenomenon. This required us to study the material closely. The following information is presented in three dierent sorts of maps. The rst sort of information maps a single phonetic or grammatical phenomenon which is not connected to specic lexical material. The second sort characterizes a certain phonetic or grammatical phenomenon within a restricted set of words, and

5 the third and nal sort presents a certain phonetic or grammatical phenomenon exemplied by only one word. The basic principle in the presentation of the information is the default to `accept more general information unless otherwise stated'. Consequently, we sometimes interpreted inexplicitness as conrmation that nothing special needed to be said. In the present paper we extract words from these atlases which we then compare in pronunciation. Our method (described below) relies on transcriptions of entire words, which we took from the atlas as best we could (see discussion above). Where we needed to extrapolate, we always did this conservatively, e.g. using little phonetic detail.

3.2 Data Preparation 3.2.1 The maps For the present experiment we sought the transcribed pronunciations of a common set of words for all Bulgarian areas. In view of this, we proceeded in the following manner: First, we chose 490 dialect sites within Bulgaria with as much geographic distribution as possible. In addition we included the pronunciation Standard Bulgarian.1 The sites were selected with respect to two main criteria: maximally complete coverage of the area covered by the atlas, and a representative number of varieties and sub-varieties. There are altogether 1682 sites in the atlas, so that our selection of 490 constitutes roughly one third of all the sites. See the map in Figure 2.2 We would have preferred using sites selected randomly from a regular grid throughout Bulgaria, but there were no collection sites in large stretches of the country. This explains the patchy impression of the map. Second, we digitized a set of 54 words, which turned out not to be instantiated in every sites, but which includes a subset of 36 words that were instantiated in all the atlases. This dierentiation of two sets arose because, as noted above, the lexical material diers a great deal across the four atlases. But we decided also to view the 36-word and 54-word sets as a test of the consistency of the measurements. So while we concentrate here on reporting the results of the experiment with a common set of 36 words, we recognize that some might consider this to be too small a sample (we would disagree). We shall therefore likewise report on the experiments with 54 words, and in particular on the degree to which the two sets of measurements correlate. Based on our experience using pronunciation distance measures in other areas, we expect a high correlation between the two sets of measurements. The digitization step involved transliterating from a Bulgarian system of phonetic transcription into IPA, which was processed in its computerized form, X-SAMPA. We include several tables to show how we interpreted the Bulgarian phonetic transcription system in terms of equivalents in the International Phon- etic Alphabet (International Phonetic Association 2003). Table 1 presents the

1The standard pronunciations are in accordance with Popov et.al. (1998). 2Note that in all maps presented here, the boundary lines indicate administrative divisions, not dialect areas.

6 Figure 2: The distribution of the Bulgarian sites selected. interpretation of vowel symbols, and Table 2 presents the interpretation of the consonant symbols. The words selected represent various parts of speech (nouns, adverbs, ad- jectives, verbs, numerals, participles) and dierent word forms (singular and plural, 3rd person verb forms, past tense forms, etc.). We rst list the words used in phonemic transcription (see Table 3), abstracting away from the phonetic variation found among the dialect sites. At the risk of repetition, we note that the list below is therefore not meant to suggest the range of phonetic variation found in Bulgaria, nor is it meant to provide phonetic detail about any single variety. For example, nal devoicing, which is very common in Bulgarian dialects, is not represented, nor are vowel reductions. In both cases we reason that this level of detail is subphonemic. Table 3 shows the subset of 36 words which were common to all 490 sites. The words in Table 3 represent many of the most important phonetic features of Bulgarian. They reect the following phenomena: 1. the reections of `yat' in dierent phonetic contexts (stressed and unstressed, word-nally, after fricatives, etc.): ["bjala, "bEli, "grESka, mlE"kar, > "v@trE, vEn"tSilo] 2. the reections of the etymological `ja': ["jazdi] or ["jezdi], [po"ljana] or [po"lena], [gu"ljaj] or [gu"lej] 3. palatal-nonpalatal-semipalatal distinction word-nally: [sol] or [solj], [p@t] or [p@tj], [kon] or [konj], [dEn] or [denj] > > > 4. the realizations of `schwa' under stress: ["b@tSva] or ["botSva] or ["batSva] etc. The same for other words: ["z@lva, s@n, "t@nko, oti"S@l, do"S@l]

7 Bulgarian IPA Bulgarian IPA Bulgarian IPA a a a 2 a 5 ˙ ¤a A ¤a A: ©a a: ¦a A` o O o 0 ˙ o U ¢o 6 ż @ ż¢ 9 ż @ ż 7 ˙ fi u u y u ¤y u fi – ÿ ü e˙ ¢ea E æfl ¢e æ e 8 ©e E: ¦e e e œ Æ Efl ˚. ˙ e £E e 8: è i ¤è Y è i ©è i: ˙ fi û 1 §è §

Table 1: The vowel symbols used in the Bulgarian dialect atlases (Stoykov and Bershteyn 1964, etc.) together with the IPA equivalents they were translated into.

5. the realizations of the nasal vowel: [z@b] or [zob] or [zab] etc. Similarly for other words: ["k@Sta, "s@bota]

6. the metatheses `@l-l@' and `@r-r@': [gr@b] or [g@rb], [Z@lt] or [Zl@t]. > > 7. the realizations of various vowels in dierent contexts: ["ovtSe] or ["ovtSo], > > [kljutS] or [klitS]

8. the reduction of the open vowels in unstressed position: [mlE"kar] or > > [mli"kar], [vEn"tSilo] or [vin"tSilo] etc.

The full list of 54 words includes the 36 words above plus the 18 words shown in Table 4, which were not common to all sites. More precisely, they were available only for the sites in the rst and the second atlases. The phonetic information that this set adds is the behavior of the fricative `x' in some of the words. For example, the presence or absence of `x' is indicated in ["xladno] and its alternative ["ladno]:

3.2.2 Digitization When we started the task, the data was available only in printed form. Since we wished to analyze the data computationally, digitization became a very important subtask. The data was converted to the X-SAMPA (Wells 2003) representation of IPA (the alphabet of the International Phonetic Association, see Handbook IPA 2003), because X-SAMPA is used within the Levenshtein-based toolkit we used (www.let.rug.nl/kleiweg/L04/). X-SAMPA encodes IPA, but

8 Bulgarian IPA Bulgarian IPA Bulgarian IPA á b á' bj á> bj â v â' vj â> vj ã g ã' gj ã> gj ä d ä' dj ä> dj æ Z æ' Zj ç z ç' zj ç> zj ê k ê' kj ê> kj ë l ë' lj ë> lj ì m ì' mj ì> mj í n í' nj í> nj ï p ï' pj ï> pj ð r ð' rj ð> rj ñ s ñ' sj ñ> sj ò t ò' tj ò> tj ô f ô' fj ô> fj ö ts> ö' ts>j ö> ts>j > > ø S ÷ tS ÷' tSj > > u dZ u' dZj " " s dz> s' dz> j ϕ F w B h è x x v G X X ð r ë l l L ˚ " ˚ " ˚ " ë Ï l L n N é j

Table 2: The consonants symbols used in the Bulgarian dialect atlases (Stoykov and Bershteyn 1964, etc.) together with the IPA equivalents they were translated into. The rst 19 lines (above the double horizontal lines) are the nonpalatal/palatal series of consonants. Note that the atlas distinguished three levels nonpalatalized, mildly palatalized and very palatizedwhich we also retain in processing. IPA does not distinguish more than a single degree of palatalization. The unpaired consonants are at the bottom of the table, below the double horizontal rule.

9 > áú÷âà /"b@tSva/, `barrel' çúëâà /"z@lva/, `sister-in-law' äîøúë /do"S@l/, `has come-he' æúëò /Z@lt/, `yellow' çúá /z@b/, `tooth' ñúáîòà /"s@bota/, `Saturday' êúùà /"k@Sta/, `house' áÿëà /"bjala/, `white'- fem áåëè /"bEli/, `white'- pl ÿçäè /"jazdi/, `ride'-3per íåäåëÿ /nE"dElja/, `Sunday' ìëåêàð /mlE"kar/, `milkman' > ãðåøêà /"grESka/, `mistake' âåí÷èëî /vEn"tSilo/, `married life' > > êëþ÷ /kljutS/, `key' ÷àøà /"tSaSa/, `glass; cup' ïúò /p@t/, `road' æàáè /"Zabi/, `frogs' íîùâè /"noStvi/, `hutch' ïîëÿíà /po"ljana/, `glade' > îâ÷å /"ovtSe/, `sheep's' òúíêî /"t@nko/, `narrow-neut' > ãóëÿé /gu"ljaj/, `feast' îâ÷àð /ov"tSar/, `shepherd' êîí /kon/, `horse' ñúí /s@n/, `dream' îòèøúë /oti"S@l/, `has gone-he' âúòðå /"v@trE/, `inside' òåíäæåðà /"tEndZEra/, `pot' äæîá /dZob/, `pocket' > íÿìà /"njama/, `there is no' ÷åðåøà /tSE"rESa/, `cherry' ãðúá /gr@b/, `back' æèâÿ /Zi"vja/, `lived' ñîë /sol/, `salt' äåí /dEn/, `day'

Table 3: The thirty-six Bulgarian words which formed the base of the study in phonemic transcription. All of 490 sites used in this study included phonetic transcriptions of these thirty-six words.

øåïà /"SEpa/, `handful' äâå /dve/, `two' ÿñëà /"jasla/, `manger' øàïêà /"Sapka/, `hat' > æèâ /Ziv/, `alive' ïî÷èâàì /po"tSivam/, `rest-I' øèðîê /Si"rok/, `wide' çåò /zEt/, `son-in-law' åçèê /E"zik/, `tongue' äîáúð /do"b@r/, `good' îòèâàì /o"tivam/, `go-I' ìúðçåëèâ /m@rze"liv/, `lazy' ÿáúëêà /"jab@lka/, `apple' õàéäå /"xajdE/, `let's' õëàäíî /"xladno/, `cold' ñíàõà /sna"xa/, `daughter-in-law' áåðÿõà /bE"rjaxa/, `were picking up-they' äàäîõ /"dadox/, `gave-I'

Table 4: The eighteen additional words which were used whenever sites included phonetic transcriptions for their pronunciations. Some sites lacked transcriptions for them.

10 is ASCII-based, and therefore easily processed by a virtually all software on all platforms. We also provided a more permanent form in an XML format using the facilities of the CLaRK System (Simov, Simov, Ganev, Ivanova & Grigorov 2004), and using Unicode and attempting to standardize geographic references. We hope to report on this separately.

4 Measuring Pronunciation Distance

4.1 Method Levenshtein distance is a technique to compare a pair of strings (words) and to assay their distance from each other. Two dialects are compared by comparing the pronunciation of the same words in the two dialects and then averaging the distances of the pairs of words. An eective way to understand Levenshtein distance is to consider how one pronunciation may be transformed into the other by means of inserting, deleting or substituting individual sounds (symbols). Weights are assigned to these three operations. In the simplest form of the algorithm, all operations have the same cost. We illustrate this with an example of two varieties of a word pronunciation > in northwestern dialects. ×åðåøà (`cherry') is pronounced3 as ["tsrESnja] in some dialects in the most northwestern part (near Kula, Belogradchik) and the > transitional area (around Tryn), and as [tSE"rESa] in other dialects such as those near Soa. Changing one pronunciation into the other may proceed as follows (ignoring suprasegmentals and diacritics for this moment):

> > > tsrESna subst. ts by tS 1 > tSrESna insert E 1 > tSErESna delete n 1 > tSErESa 3

> > In fact many sequence operations map [tsrESna] to [tSErESa]. The power of the Levenshtein algorithm is that it always nds the cost of the cheapest mapping. Levenshtein distance is then the distance assigned by the Levenshtein algorithm, the cost of the least expensive means of mapping one string to another. Comparing pronunciations in this way, the distance between longer pronunciations will generally be greater than the distance between shorter pronunciations. The longer the word, the greater the chance for dierences with respect to the corresponding word in another variety. In order not to overemphasize the importance of longer words, we normalize the raw Levenshtein distance using word length. Thus, the sum of the operations is divided by the length of the longest alignment which gives the minimum cost. The longest alignment has the greatest number of matches. In our example we have the following alignment:

3The encodings in this paper are in IPA (Handbook of the International Phonetic Associ- ation 2003)

11 > ts ∅ r ES n a > tSE r ES ∅ a 1 1 1

The total cost of 3 (1+1+1) is now divided by the length of 7. This gives a word distance of 0.43 or 43%. The simplest versions of this method are based on a notion of phonetic distance in which phonetic overlap is binary: non-identical phones contribute to phonetic distance, identical ones do not. In this simplest version the pair [I,6] counts as dierent to the same degree as [I,e]. In more sensitive versions gradual segment distances are used as weights. The segment distances are based on comparison of feature values or acoustic measurements. In fact, Heeringa (2004) shows that the phone-based methods outperform most of the methods which are using gradual segment distances as operation weights (see p. 186 and p. 194). We use this simple version of the Levenshtein algorithm in this paper.

4.2 Results Given a sample of word pronunciations from two sites, we average the distances of all pairs of corresponding words to obtain an estimate of the aggregate pronunciation dierence between any two sites. We repeat this for all (490 ∗ (490 − 1))/2(= 119, 805) pairs of sites. An important issue is whether the data sample is large enough for us to extract a reliable signal. As a measure of reliability we used Cronbach's α method (for details see Heeringa (2004, pp. 170173)), for which a widely accepted threshold is 0.70. Our results show a value of 0.84 for the set of 36 words. We therefore view the data sample large enough to provide a reliable view of pronunciation dierences. This is incident- ally the justifcation for the remark above ( 3.2.1) that we regard the 36-word sample as large enough. Conceptually, our application of the pronunciation dierence measurements results in a 490 × 490 table of pronunciation dierences we have measured using the Levenshtein algorithm. Naturally, we may restrict our attention to half of the table, since the distances in it are symmetrical. Once we know the (pronunciation) distance from Plovdiv to Soa, we know the distance in return from Soa to Plovdiv. In Figure 3 the connections between the dialects are shown based on the Levenshtein distances for 36 words. Darker lines mark pairs of sites which are more similar to each other. Even at this level, without any clustering or further analysis, the important division between western and eastern dialects emerges clearly in the map. it also suggests that the western dialect groups are closer and more coherent than the eastern ones. Further, it shows some close varieties in the East: the southern Rodopian group, the central Balkan group and the central part of the Moesian dialects. In addition, this map shows moderately strong connections between the Balkan and Moesian dialects, which appear here as merged. There is little discussion in the literature as to whether the western dialects are phonetically more cohesive than the eastern ones as a whole. On the one

12 Figure 3: The average Levenshtein distances between 490 Bulgarian dialects are shown for 36 words. Darker lines indicate close varieties, while lighter lines indicate more remote ones. hand, in the West there have not been as many large migrations as in the East, which would promote coherence. On the other hand, the East is supposed to be more coherent with respect to the `yat' realizations (but less coherent with respect to developments of the Old Bulgarian nasal). The aggregate pronunciation dierences highlighted by the average Levenshtein distance demonstrate the West to be a dialectally cohesive area.

5 Analyzing Dialect Distances

From the line map, discussed above, we have already obtained some useful information about the dialects, but it is not sucient for comparison to the results of traditional maps. For that reason, we continue our analysis by exploring the results further. On the one hand we will use cluster analysis in order to divide the dialects into similarity classes. The resulting classication allows us to compare our ndings with the results of dialectological scholarship, which has focused on the identication of dialect areas. We additionally examine the result of the pronunciation dierences using a composite clustering techniques, intended to mitigate chance eects in clustering. We also analyze the pronunciation distance table using multidimensional scaling, which complements clustering nicely and provides a view of the dialect landscape as a continuum. These techniques are presented in more detail below.

13 5.1 Cluster analysis 5.1.1 Method Clustering is commonly applied in many disciplines as an exploratory technique in data analysis, one intended to expose natural classes of similar instances. The result of clustering is a dendrogram, a tree incorporating all the input elements, and in which more similar elements are grouped lower in the tree, i.e., closer to the input items, the leaves (the nest divisions). See Fig. 4. Clustering is most easily understood procedurally. At each step of the procedure we select the shortest distance in the matrix of Levenshtein distances obtained above, and we fuse the two data points which gave rise to it. Since we wish to iterate the procedure, we have to assign a distance from the newly formed cluster to all remaining points. Although there are many so-called matrix- updating algorithms, (Jain and Dubes, 1988) we may be content with simply averaging the lengths from the two points being fused to other elements being clustered.4 We can measure the quality of a clustering result by comparing the distances in the dendrogram (e.g. the number of nodes which have to be traversed to move from one site to another in the dendrogram) to the distances in the pronunciation distance table via a cophenetic correlation coefficient. This is simply the Pearson correlation coecient of the two distances. The clustering technique we used obtained a cophenetic correlation coecient (0.71). The dendogram obtained with this method explains (0.71)2 ×100 = 50.4% of the variance of the original Levenshtein distances. A disadvantage of clustering is that it is statistically not stable: small dif- ferences in input values may lead to substantial dierences in dendrograms. We use it nevertheless for its value in allowing us to compare our results to those of traditional scholarship.

5.1.2 Results Here we present the results from the clustering in two views: in a dendrogram (Fig. 4) and in a classication map (Fig. 5). The numbers which indicate clusters in the dendrogram correspond to the numbered areas in the area map to facilitate the comparison. As the distances among a large number of varieties (490) were calculated, only one general dendrogram is presented here. The dendrogram in Fig. 4 shows ve clusters, which represent the same divisions as in the dialect area map (see Fig. 5 for comparison). Classication Map We likewise present a classication map, showing the varieties assigned to dialect groups. It represents the projection of the dendrogram in Fig. 4 onto Bulgarian geography. In particular we project the ve-way division from the dendrogram in Fig. 4 in the map in Figure 4.

4For those wanting technical detail, we note that we used the `weighted pair group method using arithmetic averages' (UPGMA). See Jain and Dubes (1988).Nerbonne and Siedle (2005, p. 9) argue that this technique is less sensitive to irregularities in geographic sampling.

14 1 Ruse

5 Pleven

Shumen Varna Lovech 4

3 Teteven

Sofia

Burgas 15 Plovdiv 2 Blagoevgrad Malko Tyrnovo Razlog 1 2 Smolyan 3 4 5 0.0 0.001 0.002 0.003 0.004

Figure 4: The dendrogram on the left shows the ve most signicant clusters of Bulgarian dialects. The leaves on the far left correspond to the individual sites in the sample, which are gradually fused into subclusters as one follows the diagram to the root on the right (WPGMA clustering was used, see text). The scale distance shows average Levenshtein distances as a fraction. The tree structure explains 50.4 % of the average Levenshtein distances. Five areas are distinguished in the map on the right, roughly the following: 1. some Rodopi dialects, 2. eastern dialects, 3. other Rodopi dialects, 4. the Teteven region and the northeastern Balkan region and 5. western and south central dialects. The map in Fig. 5 shows that the Rodopian group is not dialectally uniform. The `yat' border is also not represented to the north of Teteven, and as a consequence, the northwestern dialects appear to be rather indistinct from the eastern ones. Although we need to exercise caution in interpreting the results of clustering, since real distinctions may not be reected well, still we note one possible explanation for this: in the sample the similarity of the northwestern with the eastern dialects is reected extensively in the development of the stressed schwa. Northwestern and eastern dialects have a schwa in contrast to southwestern dialects, where we nd an [a]. The map emphasizes a northeastern Balkan dialect area near Varna and Dobrich. This detail deserves attention, because this group behaves distinctly in the continuum map, too (see below). It also distinguishes one of the central Balkan sub-dialects of the eastern part, namely the Teteven dialect. Recall that this dialect was noted as one of the divisions with respect to the `big nosovka' development the open /e/ vs. /æ/. Fi- nally, the western region around Soa is distinguished, which surprisingly shows similarities with eastern areas. Composite Cluster Map The composite cluster map (Fig. 6) is obtained by repeatedly clustering while adding random small amounts of noise to the input distance table. By repeatedly clustering using noise, we overcome the instability inherent in clustering (see above). See Kleiweg, Nerbonne and Bosveld (2004) for details. The repetitions of the clustering procedure result in maps which are superimposed on one another in Fig. 6. The darker the line, the more frequently the boundary appeared in one of the repetitions of clustering. A composite cluster map for Bulgarian was created with the same clustering method used to obtain the dendogram in Fig. 4. This map shows the most signicant divisions of the groups, where the most signicant border is again the one that divides the dialects into western and eastern. Its shape is certainly similar to the one on the traditional map in Fig. 1. However, the contrast of the northern part of the `yat' border near Pleven is interestingly suppressed. At the same time, the map shows that the most heterogenous group is Rodopi. Additionally, this map distinguishes the northeastern Balkan dialect area near Varna and Dobrich, the Teteven dialect and the Plovdiv dialect.

5.2 Multidimensional scaling 5.2.1 Method On the basis of geographic coordinates one can derive distances between any two locations. The reverse is also possible: on the basis of distances, an coordinate system can be proposed within which the data points may be located. The last is realized by a technique known as `multidimensional scaling' (MDS), which, unlike clustering is statistically stable. In an MDS plot, strongly related dialects are close to each other, while very dierent dialects are located far away from each other (Kruskal and Wish, 1984). As input each dialect is viewed as located by a set of distances, namely the distance to itself and the distances to other dialects. The other sites correspond

16 Ruse

Pleven

Shumen Varna Lovech

Teteven

Sofia

Burgas

Plovdiv Blagoevgrad Malko Tyrnovo Razlog

Smolyan

Figure 5: Composite cluster map for 36 words. The darker lines indicate more likely linguistic borders. thus to dimensions, so that if we have 490 dialects, we get 490 dimensions. Using MDS these 490 dimensions can be reduced to one, two, three or more. Each site is then identied via coordinates in 2-, 3-, or 4-dimensional space. Notice that we can then calculate the distance between any two points in the dimensions provided by the MDS solution: this is simply the Euclidean distance in the space of the n dimensions. We then obtain an idea of the quality of the MDS data reduction by examining how well the original distances correlate with the distances dened implicitly in the given MDS solution. MDS also allows the denition of a color map where the resulting MDS dimensions are assigned colors, and the intensity of a color at a given point reects the MDS coordinates assigned to that site. If our 490 dimensions are scaled to 3 dimensions, and we let the 3 dimensions be the intensities of respectively red, green and blue, then each dialect site is assigned a mixture of these three colors. The result is a map which reects the gradual changes in dialect characteristics (see Figure 7). The map is not constrained to reect a dialect continuum, which is why the polygons surrounding a given site may change color abruptly. To the degree that the color shifts are not abrupt, we are dealing with a continuum.

5.2.2 Results The rst and second dimensions of the MDS solution account correlate with the original distance matrix to almost the same degree (r = 0.625, r = 0.620, respectively), while the third dimension correlates much more weakly (r = 0.205). They thus account independently for 39.0%, 38.5%, and 4.2% of the variation in the data. Since there is some collinearity, the total explained variance is less than the sum of that of the three dimensions, but still high: 88.3%. We used

17 Ruse

Pleven

Shumen Varna Lovech

Teteven

Sofia

Burgas

Plovdiv Blagoevgrad Malko Tyrnovo Razlog

Smolyan

Figure 6: This continuum map highlights two subgroups in western Bulgaria (northeastern and south western) as well as the most northeastern subgroup. this three-dimension reduction, mapping each of the dimensions to a color, as explained above. The map in Fig. 7, which does not depend at all on clustering, represents the dialects as forming a continuum, a view which most dialectologists nd con- genial. The map suggests that eastern Bulgaria is a rather uniform area, where only the far eastern part around Varna and Shumen is distinguished. Also, the western part is clearly divided into two parts: northwestern and southwestern. The dialectally least uniform area remains Rodopi.

5.3 Comparison In general, the important border between western and eastern dialects is present in all the maps. Within the western group two varieties are regularly detected: northwest and southwest groups, in conformity with the traditional division. In eastern Bulgaria the Moesian and Balkan dialects are not distinguished sharply, which may reect the fact that a number of migrations have taken place in this area. Of course, we should bear in mind the fact that our procedure is based exclusively on pronunciation, and not all on morphology or the lexicon. The Balkan group around Varna and Shumen is partly distinguished, which may be due to the irregularity in the grid of data collection sites. We regard our results as partial when compared to the view embodied in the map on Fig. 1. The most heterogenous groups seem to be the Rodopian group, and the eastern Balkan group. We examine these separately below. If we compare the maps we have derived to the maps in the generalized volume of Bulgarian dialects (Kochev 1988), then we can see the following: our maps reect e.g. the `yat' border between western and eastern dialects (see map

18 2/15) and also the borders in the Atlas map 2/2, which distinguishes the Balkan area around Varna and Shumen and the varieties within the Rodopian group and the neighboring southwestern groups. The features which were taken into account in these two traditional dialect maps were the oppositions: palatal vs. non-palatal consonants and closed vs. open front vowels. These were reected in our data as well.

6 Individual Pronunciation Distributions

This section shows and briey discusses the distributions of three words, whose variants correlate well with the main dimensions of variation we found in the multi-dimensional scaling analysis of the aggregate data. Figure 7 illustrates some of the variation on which the aggregate analysis is based. The top map shows examples of the distribution of áåëè /"bEli/, `white(pl.)', which indicates part of the `yat' alternation. This word correlates most strongly with the second dimension of the MDS analysis (r = 0.62). In Standard Bulgarian the pronunciation is [bEli]. The western dialects are quite homogenous with respect to this feature which shows various along the southern and eastern borders. Thus, [æ] appears for /E/ in the Northeast with concur- rent palatalization of the preceding consonant [bjæli]. But note that the [æ] is independently realized, i.e. without the palatalization, in the South. Both these variants however are greatly outnumbered by the otherwise completely consistent of [E] throughout the rest of Bulgaria. The second, middle map in Figure 7 illustrates the greater variety of pronunciations found for æàáè /"Zabi/, `frogs'. This word, which correlates most strongly with the rst dimension of the MDS analysis (r = 0.53), exemplies the alternation [a-E] after alveolar fricatives and aricates. The western dialects do not alternate as a rule, i.e. they keep the Standard pronunciation form with [a]. On the other hand, the eastern dialects show alternation, which varies between [E] and [æ]. This alternation resembles the `yat' alternation of /"bEli/, but note that the distribution is not the same: the South and Southeast show [æ], again without preceding palatalization, just as in the case of /"bEli/, but the Northeast and the rest of the eastern half of the country pronounces [E] in this word. The third and nal map in Figure 7 illustrates some of the variation found in the realizations of the nasal sound `big novoska'. Note that in contrast to the two previous maps, here the variety applies in higher degree to the western dialects. The word êúùà /"k@Sta/, `house', which is nearly the strongest correlate of the third MDS dimension (r = 0.19), corresponds with three pronunciations of the rst vowel, which however is not the exhaustive list of the possible alternations. The majority pronunciation is [@], which is standard and is found everywhere except in the central West and Southwest, where [a] is frequently found, and on the northeastern coast of the Black Sea and in south Rodopi, where [æ] is present. êúùà /"k@Sta/, `house', shown in the bottom map of Fig. 7, is also interesting because it is the single best diagnostic of dialect area in all of the 54

5We re-use the map numbering used in the atlas.

19 words examined. We determined this as follows: when we obtained 54 distance matrices based on each of the 54 words in the sample, it turns out the single word which correlates most highly with the matrix of aggregate distances the matrix generated from êúùà /"k@Sta/, `house' (r = 0.62). All the maps shown here accord with traditional dialect divisions, which are, of course, based on the same and similar phenomena. There may appear to be minor deviations from the traditional maps, however, which have several causes. For example, our sample did not include all the varieties (see the æàáè `frogs' map). Further, the maps illustrate only part of the sample, leaving further variation unmentioned (the plural `yat' form only in the áåëè `white' map).

7 The Rodopian and Eastern Balkan groups

Before seeking groups each of 490 sites was labeled with the name of the main dialect to which it is supposed to belong. The labels are meant to facilitate the systematic examination of the groupsnothing more. The main dialects are abbreviated in the labels, and we used the labels in the dendrograms presented in this section. The dialects are as follows: Northeastern dialects: Moesian (M); Balkan (Ba); South Balkan (Pba) Southeastern dialects: Thracian (Thr); West Rupskian (Zrup); Rodopian (Ro) Northwestern dialects: West Moesian (ZM); Northwestern borders (SzP) Southwestern dialects: Central Southwestern (CYuz); Southwestern borders (YuzP); Standard (St) We focus now in more detail on two interesting and distinct groups, namely the eastern Balkan group and the Rodopian groups. We will in particular examine them via the dendrograms resulting from clustering. The eastern Balkan group around Varna and Shumen is a compact group that is distant from the surrounding dialects.6 This group seems to behave dierently with respect to one very decisive phonetic feature, namely open /e/ vs. /æ/, which is the realization of `yat' before syllables with a closed vowel, of schwa in a closed syllable, of small er, and of the big nosovka. This group is shown in Fig. 8. Note that the closest neighboring dialects are the so called `erkech' dialects (Kazichino and Golitsa), and the Teteven dialects, which are subdialects of the central Balkan dialect. Recall that this connection is present in the area map (Fig. 5). Traditionally, the southeastern group includes: eastern Rupskian dialects (which are presented as Rodopian and Thracian here), central Rupskian dialects (presented as Rodopian here), and western Rupskian dialects (presented as western Rupskian here as well). The maps show that eastern Rupskian and western Rupskian (with some small exceptions) are more closely unied than

6Note that it is marked as Balkan in the traditional map, but we view it as a group of Moesian dialects, because the Moesian and Balkan dialects merge in this area allowing us to reserve the term `Balkan' for the central Balkan area around the Central Old mountains and the Sredna Gora mountain. We nd this clearer, at least in drawing comparisons.

20 ⎆bϯli ⎆bϯli

⎆bϯli ⎆bϯli ⎆bܕæli ⎆bϯli ⎆bϯli ⎆bϯli ⎆bϯli ⎆bϯli ⎆bϯli ⎆bϯli ⎆bϯli ⎆bæli

⎆bæli

⎆Ћabi ⎆Ћϯbi

⎆Ћabi ⎆Ћϯbi ⎆Ћϯbi ⎆Ћabi ⎆Ћϯbi ⎆Ћabi ⎆Ћabi ⎆Ћϯbi ⎆Ћϯbi ⎆Ћϯbi ⎆Ћabi ⎆Ћæbi

⎆Ћæbi

⎆kђЀta ⎆kђЀta

⎆kђЀta ⎆kђЀta ⎆kђЀta / ⎆kæЀta ⎆kaЀta ⎆kђЀta ⎆kaЀta ⎆kaЀta ⎆kђЀta ⎆kђЀta ⎆kђЀta ⎆kђЀta ⎆kђЀta

⎆kæЀta

Figure 7: The variants of three fairly diagnostic words, /"bEli/ `white', /"k@Sta/ `house', and /"Zabi/ `frogs'. The borders shown are administrative, not linguistic, and have been added to aid in orientation. See text for discussion.

21 Ablanitsa (Zrup) Godyashevo (Zrup) Tuhovishta (Zrup) Ribnovo (Zrup) Borovina (Ro) Levochevo (Ro) Pisanitsa (Ro) Korovo (Ro) Rakitovo (Ro) Kostandovo (Ro) Momchilovtsi (Ro) Skrebatno (Zrup) Lehovo (CYuz) Beden (Ro) Breze (Ro) Babyak (Zrup) Bansko (Zrup) Dobyrsko (Zrup) Godlevo (Zrup) Razlog (Zrup) Bunovo (Ba) Chelopech (Ba) Pirdop (Ba) Smolsko (Ba) Tserovo (Pba) Pravets (CYuzS) Batuliya (CYuzS) Belchin (CYuz) Relyovo (CYuz) Kondofrey (CYuz) Kovachevtsi (CYuz) Yarlovo (CYuz) Podgumer (CYuzS) Chepintsi2 (CYuzS) Lyutibrod (CYuzS) Stolnik (CYuz) Oselna (CYuzS) Yablanitsa (CYuzS) Zanoge (ZM) Zanozhene (CYuzS) Seslavtsi (CYuzS) Leskovdol (CYuzS) Sestrimo (Pba) Bistrets (CYuzS) Chelopek (CYuzS) Dyrmantsi (CYuzS) Boboshevo (CYuz) Frolosh (CYuz) Eremiya (CYuz) Drumohar (CYuz) Smolichano (CYuz) Padala (CYuz) Bosnek (CYuz) Chupetlovo (CYuz) Shipochan (CYuz) Kremenik (CYuz) Dospey (CYuz) Madjare (CYuz) Govedartsi (CYuz) Raduil (CYuz) Gabra (CYuz) Kostentsi (CYuzS) Gutsal (CYuz) Zhivkovo (CYuz) Vrachesh (CYuzS) Belitsa3 (Zrup) Yakoruda (Zrup) Osenovo (CYuz) Belasitsa (CYuz) Gabrene (CYuz) Bobeshino (YuzP) Golema Fucha (CYuz) Palatovo (CYuz) Lokvata (CYuz) Mursalevo (CYuz) Novoselyane (CYuz) Ryzhdavitsa (YuzP) Gramazhdeno (CYuz) Lozno (CYuz) Brestovo (CYuz) Leshko (CYuz) Padezh (CYuz) Brest (YuzP) Kryshalevo (YuzP) Lomnitsa (YuzP) Zlogosh (YuzP) Saparevo (CYuz) Brezhani (CYuz) Gorno Osenovo (CYuz) Debrene (CYuz) Stozha (CYuz) Kresna (CYuz) Mechkul (CYuz) Vlahi (CYuz) Moravska (CYuz) Velyushtets (CYuz) Oshtava (CYuz) Churichene (CYuz) Gega (CYuz) Krynalovo (CYuz) Kyrlanovo (CYuz) Lozenitsa (CYuz) Rozhen (CYuz) Lobosh (YuzP) Vyrba (YuzP) Bohova (SzP) Zabel (SzP) Miloslavtsi (SzP) Slivovnik (SzP) Protopopintsi (SzP) Boynitsa (SzP) Bukorovtsi (SzP) Struindol (SzP) Chuprene (SzP) Kopilovtsi (SzP) Chepyrlintsi (SzP) Stanintsi (SzP) Komshtitsa (SzP) Rebro (SzP) Gorni Romantsi (SzP) Velkovtsi (YuzP) Izvor2 (YuzP) Selishten dol (YuzP) Divlya (YuzP) Shipkovtsi (YuzP) Treklyano (YuzP) Vidrar (YuzP) Metohiya (YuzP) Stnyovtsi (YuzP) Gorni Korten (YuzP) Bachkovo (Ro) Kribul (Zrup) Dorkovo (Ro) Mugla (Ro) Selcha (Ro) Strizhba (Ro) Tsarino (Ro) Zhyltusha (Ro) Botevo (M) Levski (M) Kalimantsi (M) Momino (M) Ravna gora (M) Rudnik (M) Sadovo (M) Izgrev (M) Nikolaevka (M) Zasmyano (M) Voditsa (M) Shtipsko (M) Kazichino (Ba) Golitsa (Ba) Cherni Vit (Ba) Ribaritsa (Ba) Golyama Zhelyazna (Ba) Ahryane (Ro) Bogutevo (Ro) Figure 8: The MoesianOsikovo (Ro) dialects group very uniformly together. Recall that M Petvar (Ro) stands for MoesianSitovo and (Pba) Ba for Balkan. Borikovo (Ro) Mogilitsa (Ro) central RupskianChepintsi (Rodopian). (Ro) In this group, the diering features are more prominent thanKozhari the unifying (Ro) ones. The dendrograms show the distinct sub- Vodnipad (Ro) dialect groups: dialectsTrigrad (Ro) of Smolyan, of Zlatograd, of Velingrad, etc. The rst RodopianKremene (Ro) cluster (from Ahryane to Tihomir in Fig. 9) is rather com- Vitina (Ro) pact in the dendrogram.Kremene2 (Ro) The main dialect is Smolyan. We can also distinguish some of the otherSmilyan neighboring (Ro) sub-dialects, such as the `Shiroka lyka' dialect Raykovo (Ro) (Shiroka lyka, Stoykite,Smolyan (Ro) etc.). Its closest non-Rodopian clusters are Balkan and Moesian. Raykovo2 (Ro) The secondSmolyan2 cluster (Ro) (see Fig. 10) presents the Chepinski Rodopian sub- Ustovo (Ro) dialect. It is closelyGela related (Ro) to west Rupskian dialects. The thirdShiroka Rodopian Lyka (Ro) cluster (see Fig. 11) is closer to the southwestern bor- Stoykite (Ro) dering dialects andPetkovo Moesia. (Ro) The fourth RodopianTyryn (Ro) cluster (see Fig. 12) shows the Hvoyna Rodopian sub- dialect, whichChokmanovo is interestingly (Ro) separated from the other Rodopian sub-dialects Vievo (Ro) present in the rstBeslen dendrogram: (Zrup) Smolyan, Shiroka lyka, Chepintsi. It is closer to west RupskianLydja dialects (Ro) and Balkan dialects. Devisilitsa (Ro) The fthMalyk Rodopian Devisil (Ro) cluster (see Figure 13) shows the Zlatograd Rodopian sub-dialectGolyam and itsDevisil relatedness (Ro) to some of the dialects around Smolyan. It shows similarities to theTihomir Thracian (Ro) dialects as well. Agatovo (Ba) To sum up,Byala these reka subclusters (M) emphasize the integrity of smaller dialect areas even in those areasKoevtsi which (Ba) show little overall cohesiveness. The traditional groups Djulyunitsa2 (M) do not always clusterPisarevo (M) uniformly and exhaustively together. On the contrary, they occasionallyNikyup merge, (Ba) even while displaying many connections among sub- Paskalevets (M) dialects. Kramolin (Ba) Debeltsovo (Ba) Draganovtsi2 (Ba) 8 The StandardShilkovtsi (Ba) Mindya (Ba) Yakovsti (Ba) The Bulgarian standardKipilovo (Ba) language is regarded as close to dialects in the Northeast (central BalkanStara dialect), reka (Pba) (Stoykov 2002, p. 108). At the same time, some sim- Obnova (M) Pyrvomaitsi (M) Hlebarovo (M) Tetovo (M) 22 Belyakovets (Ba) Samovodene (Ba) Dichin (Ba) Dragizhevo (Ba) Prisovo (Ba) Pirgovo (M) Gagalya (M) Ryahovo (M) Alfatar (M) Asparuhovo (M) Garvan (M) Pozharevo (M) Vetren (M) Kalipetrovo (M) Katselovo (M) Dyadovo (Pba) Lyubenovo (Pba) Kozarevo (Pba) Tvyrditsa (Pba) Asenovtsi (Ba) Letnitsa (Ba) Butovo (M) Bylgarsko Slivovo (M) Tsarevets (M) Kozlovets (M) Morava (M) Dolna Studena (M) Hadjidimitrovo (M) Krivina (M) Novograd (M) Styrmen (M) Beltsov (M) Djulyunitsa (M) Dolna Lipnitsa (M) Lesicheri (M) Patresh (M) Bjala Cherkova (M) Vishovgrad (M) Musina (M) Karaisen (M) Sredna kula (M) Cherven (M) Gabrovtsi (Ba) Kladni dyal (Ba) Golemanite (Ba) Kulina voda (M) Osym (M) Debrene (M) Opanets (M) Avren (Ro) Planinets (Ro) Belitsa2 (Thr) Mezek (Thr) Trakiets (Thr) Bryagovets (Ro) Svirkovo (Thr) Bylgarin (Thr) Babek (Pba) Belozem (Pba) Voysil (Pba) Rupkite (Pba) Bebrovo (Ba) Palitsi (Ba) Kostel (Ba) Bujnovtsi (Ba) Lazartsi (Ba) Todyuvtsi (Ba) Gradets (Ba) Zheravna (Ba) Medven (Ba) Chiflika (Ba) Drenta (Ba) Lipov ryt (Ba) Karaivantsi (Ba) Dryanovets (M) Svalenik (M) Sadina (M) Krivnya (M) Senovo (M) Orizare (Pba) Brezovo (Pba) Malina (Thr) Bylgarska polyana (Thr) Poroy (Pba) Srem (Thr) Byal kladenets (Pba) Malomir (Pba) Gavrailovo (Pba) Korten (Pba) Rakovo (Pba) Klikach (Pba) Ognen (Pba) Dinya (Pba) Pshenichevo (Pba) Dybrava (Pba) Kaloyanovets (Pba) Kazanka (Pba) Zlatari (Pba) Rosen (Thr) Rusokastro (Thr) Swetlina (Thr) Enina (Ba) Chernokonevo (Thr) Voden (Thr) Yabylkovo (Thr) Debyr (Thr) Karadjalovo (Thr) Aleksandrovo (ZM) Vylchi tryn (ZM) Gorno Pavlikene (Ba) Bylgarene (MZ) Malka Zhelyazna (Ba) Ugyrchin (Ba) Galata (Ba) Golyam Chardak (Pba) Vidrare (Ba) Golets (ZM) Oreshak (Ba) Gumoshtnik (Ba) Grivitsa (ZM) Vybel2 (ZM) Lisets (ZM) Slavyanovo (ZM) Balanite (Ba) Zhyltesh (Ba) Bozhentsite (Ba) Sabotkovtsi (Ba) Todorcheta (Ba) Chernovryh (Ba) Pyrsha (Ba) Gybene (Ba) Toplesh (Ba) Ganchovets (Ba) Belitsa (Ba) Mryzetsi (Ba) Gostilitsa (Ba) Draganovtsi (Ba) Popovtsi (Ba) Kozirog (Ba) Fyrevtsi (Ba) Batoshevo (Ba) Kupen (Ba) Dushevsko (Ba) Kryvenik (Ba) Dedina (Ba) Krushare (Pba) Kaloyanovo2 (Pba) Syrnegor (Ba) Voynyagovo (Ba) Belyovo (CYuz) Goleshovo (CYuz) Dospat (Ro) Syrnitsa (Ro) Kasyka (Ro) Gyrmen (Zrup) Nastan (Ro) Hvoyna (Ro) Orehovo (Ro) Gaytaninovo (Zrup) Nova Lovcha (Zrup) Paril (Zrup) Teshovo (Zrup) Lyki (Zrup) Rozovo (Ro) Bezden (CYuzS) Opitsvet (CYuzS) Bogdanov dol (YuzP) Myrchaevo (CYuz) Gradets2 (CYuzS) Breze2 (CYuzS) Shipot (ZM) Sumer (ZM) Gintsi (SzP) Gurmazovo (CYuzS) Ivanyane (YuzP) Slatina (YuzP) Verdikal (YuzP) Kokalyane (CYuz) Lozen (CYuz) Dondukovo (ZM) Razgrad mahala (ZM) Slavotin (ZM) Momin brod (ZM) Vladichentsi (ZM) Chiren (ZM) Galiche (ZM) Dinkovitsa (ZM) Plakuder (ZM) Inovo (ZM) Yanovets (SzP) Gorno Tserovene (ZM) Vladimirovo (ZM) Rakevo (ZM) Selanovtsi (ZM) Kosta Perchevo (ZM) Osen (ZM) Manastirishte (ZM) Zheglitsa (ZM) Furen (ZM) Gabare (ZM) Sokolare (ZM) Stavertsi (ZM) Pisarovo (ZM) Kalen (ZM) Varvara (Pba) Opanets2 (ZM) Vyrbitsa (CYuzS) Poibrene (Ba) Ovchepoltsi (Pba) Stroevo (Pba) Miromir (Pba) Debrashtitsa (Pba) Gradina (Thr) Seltsi (Ro) Konush (Ro) Dolnoslav (Ro) Gornoslav (Ro) Oreshets (Ro) Gorno Lukovo (Ro) Meden buk (Ro) Tselina (Pba) Momkovo (Thr) Yerusalimovo (Thr) Oryahovo (Thr) Konstantinovo (Thr) Izvor (Thr) Nedelino (Ro) Bryshlyan (Thr) Byala voda (Thr) Moryane (Thr) Zvezdets (Thr) Fakiya (Thr) Zabernovo (Thr) Bylgari (Thr) Kondolovo (Thr) Rezovo (Thr) Gramatikovo (Thr) Stoilovo (Thr) Gospodintsi (Zrup) Cheshnegirovo (Ro) Cherkovna (M) Ravna (M) Nenovo (M) Krivnya2 (M) Kosten (Pba) Podvis (Pba) Krumovo (M) Vybel (M) Prilep (Pba) Rosenovo (Thr) Kirilovo (Thr) Dragoevo (M) Kochovo (M) Kyulevcha (M) Bukovo (Ro) Erma reka (Ro) Davidkovo (Ro) Vyrbina (Ro) Kraynovo (Thr) Vylchan dol (Ro) Zagrazhden (Ro) Zlatograd (Ro) Belene (M) Dolen (Zrup) Pletena (Zrup) Oryahovets (Ro) Oresh (M) Kutela (Ro) Belentsi (ZM) Dermantsi (ZM) Gigen (ZM) Sinyo byrdo (ZM) Gulyantsi (ZM) Radovene (ZM) Tsarevets2 (ZM) Bryshlyanitsa (ZM) Nikolaevo (ZM) Duvanli (Pba) Kaloyanovo (Pba) General Nikolaevo (Pba) Sekirovo (Pba)

0.0 0.001 0.002 0.003 0.004 Ablanitsa (Zrup) Godyashevo (Zrup) Tuhovishta (Zrup) Ribnovo (Zrup) Borovina (Ro) Levochevo (Ro) Pisanitsa (Ro) Korovo (Ro) Rakitovo (Ro) Kostandovo (Ro) Momchilovtsi (Ro) Skrebatno (Zrup) Lehovo (CYuz) Beden (Ro) Breze (Ro) Babyak (Zrup) Bansko (Zrup) Dobyrsko (Zrup) Godlevo (Zrup) Razlog (Zrup) Bunovo (Ba) Chelopech (Ba) Pirdop (Ba) Smolsko (Ba) Tserovo (Pba) Pravets (CYuzS) Batuliya (CYuzS) Belchin (CYuz) Relyovo (CYuz) Kondofrey (CYuz) Kovachevtsi (CYuz) Yarlovo (CYuz) Podgumer (CYuzS) Chepintsi2 (CYuzS) Lyutibrod (CYuzS) Stolnik (CYuz) Oselna (CYuzS) Yablanitsa (CYuzS) Zanoge (ZM) Zanozhene (CYuzS) Seslavtsi (CYuzS) Leskovdol (CYuzS) Sestrimo (Pba) Bistrets (CYuzS) Chelopek (CYuzS) Dyrmantsi (CYuzS) Boboshevo (CYuz) Frolosh (CYuz) Eremiya (CYuz) Drumohar (CYuz) Smolichano (CYuz) Padala (CYuz) Bosnek (CYuz) Chupetlovo (CYuz) Shipochan (CYuz) Kremenik (CYuz) Dospey (CYuz) Madjare (CYuz) Govedartsi (CYuz) Raduil (CYuz) Gabra (CYuz) Kostentsi (CYuzS) Gutsal (CYuz) Zhivkovo (CYuz) Vrachesh (CYuzS) Belitsa3 (Zrup) Yakoruda (Zrup) Osenovo (CYuz) Belasitsa (CYuz) Gabrene (CYuz) Bobeshino (YuzP) Golema Fucha (CYuz) Palatovo (CYuz) Lokvata (CYuz) Mursalevo (CYuz) Novoselyane (CYuz) Ryzhdavitsa (YuzP) Gramazhdeno (CYuz) LoznoAblanitsa (CYuz) (Zrup) BrestovoGodyashevo (CYuz) (Zrup) LeshkoTuhovishta (CYuz) (Zrup) PadezhRibnovo (CYuz) (Zrup) BrestBorovina (YuzP) (Ro) KryshalevoLevochevo (YuzP) (Ro) LomnitsaPisanitsa (YuzP) (Ro) ZlogoshKorovo (YuzP) (Ro) SaparevoRakitovo (CYuz) (Ro) BrezhaniKostandovo (CYuz) (Ro) Gorno OsenovoMomchilovtsi (CYuz) (Ro) DebreneSkrebatno (CYuz) (Zrup) StozhaLehovo (CYuz) (CYuz) KresnaBeden (CYuz) (Ro) MechkulBreze (CYuz) (Ro) VlahiBabyak (CYuz) (Zrup) MoravskaBansko (CYuz) (Zrup) VelyushtetsDobyrsko (CYuz) (Zrup) OshtavaGodlevo (CYuz) (Zrup) ChuricheneRazlog (CYuz) (Zrup) GegaBunovo (CYuz) (Ba) KrynalovoChelopech (CYuz) (Ba) KyrlanovoPirdop (CYuz) (Ba) LozenitsaSmolsko (CYuz) (Ba) RozhenTserovo (CYuz) (Pba) LoboshPravets (YuzP) (CYuzS) BatuliyaVyrba (YuzP) (CYuzS) BohovaBelchin (SzP) (CYuz) RelyovoZabel (SzP) (CYuz) MiloslavtsiKondofrey (SzP) (CYuz) KovachevtsiSlivovnik (SzP) (CYuz) ProtopopintsiYarlovo (SzP) (CYuz) PodgumerBoynitsa (SzP)(CYuzS) BukorovtsiChepintsi2 (SzP)(CYuzS) LyutibrodStruindol (SzP)(CYuzS) ChupreneStolnik (SzP) (CYuz) KopilovtsiOselna (SzP)(CYuzS) ChepyrlintsiYablanitsa (SzP)(CYuzS) StanintsiZanoge (SzP) (ZM) KomshtitsaZanozhene (SzP)(CYuzS) SeslavtsiRebro (SzP)(CYuzS) Gorni LeskovdolRomantsi (SzP)(CYuzS) VelkovtsiSestrimo (YuzP) (Pba) Izvor2Bistrets (YuzP) (CYuzS) SelishtenChelopek dol (YuzP) (CYuzS) DyrmantsiDivlya (YuzP) (CYuzS) ShipkovtsiBoboshevo (YuzP) (CYuz) TreklyanoFrolosh (YuzP) (CYuz) VidrarEremiya (YuzP) (CYuz) MetohiyaDrumohar (YuzP) (CYuz) StnyovtsiSmolichano (YuzP) (CYuz) Gorni KortenPadala (YuzP) (CYuz) BachkovoBosnek (Ro)(CYuz) ChupetlovoKribul (Zrup) (CYuz) ShipochanDorkovo (Ro)(CYuz) KremenikMugla (Ro)(CYuz) DospeySelcha (Ro)(CYuz) StrizhbaMadjare (Ro)(CYuz) GovedartsiTsarino (Ro)(CYuz) ZhyltushaRaduil (Ro)(CYuz) BotevoGabra (CYuz)(M) KostentsiLevski (CYuzS) (M) KalimantsiGutsal (CYuz)(M) ZhivkovoMomino (CYuz)(M) RavnaVrachesh gora (CYuzS) (M) Belitsa3Rudnik (M)(Zrup) YakorudaSadovo (M)(Zrup) OsenovoIzgrev (CYuz)(M) NikolaevkaBelasitsa (CYuz)(M) ZasmyanoGabrene (CYuz)(M) BobeshinoVoditsa (YuzP)(M) GolemaShtipsko Fucha (CYuz)(M) KazichinoPalatovo (Ba)(CYuz) LokvataGolitsa (Ba)(CYuz) MursalevoCherni Vit (Ba)(CYuz) NovoselyaneRibaritsa (Ba)(CYuz) GolyamaRyzhdavitsa Zhelyazna (Ba)(YuzP) GramazhdenoAhryane (Ro)(CYuz) BogutevoLozno (Ro)(CYuz) BrestovoOsikovo (Ro)(CYuz) LeshkoPetvar (Ro)(CYuz) SitovoPadezh (Pba) (CYuz) BorikovoBrest (Ro) (YuzP) KryshalevoMogilitsa (Ro) (YuzP) ChepintsiLomnitsa (Ro) (YuzP) KozhariZlogosh (Ro) (YuzP) VodnipadSaparevo (Ro)(CYuz) BrezhaniTrigrad (Ro)(CYuz) GornoKremene Osenovo (Ro)(CYuz) DebreneVitina (Ro)(CYuz) Kremene2Stozha (Ro)(CYuz) SmilyanKresna (Ro)(CYuz) RaykovoMechkul (Ro)(CYuz) SmolyanVlahi (Ro)(CYuz) Raykovo2Moravska (Ro)(CYuz) VelyushtetsSmolyan2 (Ro)(CYuz) OshtavaUstovo (Ro)(CYuz) ChuricheneGela (Ro)(CYuz) Shiroka GegaLyka (Ro)(CYuz) KrynalovoStoykite (Ro)(CYuz) KyrlanovoPetkovo (Ro)(CYuz) LozenitsaTyryn (Ro)(CYuz) ChokmanovoRozhen (Ro)(CYuz) LoboshVievo (Ro) (YuzP) BeslenVyrba (Zrup) (YuzP) BohovaLydja (Ro) (SzP) DevisilitsaZabel (Ro) (SzP) MalykMiloslavtsi Devisil (Ro) (SzP) Golyam SlivovnikDevisil (Ro) (SzP) ProtopopintsiTihomir (Ro) (SzP) AgatovoBoynitsa (Ba) (SzP) ByalaBukorovtsi reka (M) (SzP) KoevtsiStruindol (Ba) (SzP) Djulyunitsa2Chuprene (M) (SzP) PisarevoKopilovtsi (M) (SzP) ChepyrlintsiNikyup (Ba) (SzP) PaskalevetsStanintsi (M) (SzP) KomshtitsaKramolin (Ba) (SzP) DebeltsovoRebro (Ba) (SzP) GorniDraganovtsi2 Romantsi (Ba) (SzP) ShilkovtsiVelkovtsi (Ba)(YuzP) MindyaIzvor2 (Ba)(YuzP) SelishtenYakovsti dol (Ba)(YuzP) KipilovoDivlya (Ba)(YuzP) StaraShipkovtsi reka (Pba) (YuzP) TreklyanoObnova (YuzP)(M) PyrvomaitsiVidrar (YuzP)(M) HlebarovoMetohiya (YuzP)(M) StnyovtsiTetovo (YuzP)(M) GorniBelyakovets Korten (Ba)(YuzP) SamovodeneBachkovo (Ba) (Ro) DichinKribul (Ba) (Zrup) DragizhevoDorkovo (Ba) (Ro) PrisovoMugla (Ba) (Ro) PirgovoSelcha (M) (Ro) GagalyaStrizhba (M) (Ro) RyahovoTsarino (M) (Ro) ZhyltushaAlfatar (M) (Ro) AsparuhovoBotevo (M) (M) GarvanLevski (M) (M) PozharevoKalimantsi (M) (M) VetrenMomino (M) (M) KalipetrovoRavna gora (M) (M) KatselovoRudnik (M) (M) DyadovoSadovo (Pba) (M) LyubenovoIzgrev (Pba) (M) KozarevoNikolaevka (Pba) (M) TvyrditsaZasmyano (Pba) (M) AsenovtsiVoditsa (Ba) (M) LetnitsaShtipsko (Ba) (M) KazichinoButovo (M) (Ba) Bylgarsko SlivovoGolitsa (M) (Ba) TsarevetsCherni Vit (M) (Ba) KozlovetsRibaritsa (M) (Ba) Golyama ZhelyaznaMorava (M) (Ba) Dolna StudenaAhryane (M) (Ro) HadjidimitrovoBogutevo (M) (Ro) KrivinaOsikovo (M) (Ro) NovogradPetvar (M) (Ro) StyrmenSitovo (M) (Pba) BeltsovBorikovo (M) (Ro) DjulyunitsaMogilitsa (M) (Ro) Dolna LipnitsaChepintsi (M) (Ro) LesicheriKozhari (M) (Ro) PatreshVodnipad (M) (Ro) Bjala CherkovaTrigrad (M) (Ro) VishovgradKremene (M) (Ro) MusinaVitina (M) (Ro) KaraisenKremene2 (M) (Ro) SrednaSmilyan kula (M) (Ro) ChervenRaykovo (M) (Ro) GabrovtsiSmolyan (Ba) (Ro) KladniRaykovo2 dyal (Ba) (Ro) GolemaniteSmolyan2 (Ba) (Ro) Kulina vodaUstovo (M) (Ro) OsymGela (M) (Ro) ShirokaDebrene Lyka (M) (Ro) OpanetsStoykite (M) (Ro) AvrenPetkovo (Ro) (Ro) PlaninetsTyryn (Ro) (Ro) ChokmanovoBelitsa2 (Thr) (Ro) MezekVievo (Thr) (Ro) TrakietsBeslen (Thr) (Zrup) BryagovetsLydja (Ro) (Ro) SvirkovoDevisilitsa (Thr) (Ro) MalykBylgarin Devisil (Thr) (Ro) GolyamBabek Devisil (Pba) (Ro) BelozemTihomir (Pba) (Ro) VoysilAgatovo (Pba) (Ba) RupkiteByala (Pba)reka (M) BebrovoKoevtsi (Ba) (Ba) Djulyunitsa2Palitsi (Ba) (M) KostelPisarevo (Ba) (M) BujnovtsiNikyup (Ba) (Ba) PaskalevetsLazartsi (Ba) (M) TodyuvtsiKramolin (Ba) (Ba) DebeltsovoGradets (Ba) (Ba) Draganovtsi2Zheravna (Ba) (Ba) MedvenShilkovtsi (Ba) (Ba) ChiflikaMindya (Ba) (Ba) DrentaYakovsti (Ba) (Ba) LipovKipilovo ryt (Ba) (Ba) KaraivantsiStara reka (Ba) (Pba) DryanovetsObnova (M) (M) SvalenikPyrvomaitsi (M) (M) HlebarovoSadina (M) (M) Ablanitsa (Zrup) KrivnyaTetovo (M) (M) Godyashevo (Zrup)BelyakovetsSenovo (M) (Ba) Tuhovishta (Zrup)SamovodeneOrizare (Pba) (Ba) Ribnovo (Zrup)BrezovoDichin (Pba) (Ba) Borovina (Ro)DragizhevoMalina (Thr) (Ba) LevochevoBylgarska (Ro) polyanaPrisovo (Thr) (Ba) Pisanitsa (Ro) PoroyPirgovo (Pba) (M) Korovo (Ro) SremGagalya (Thr) (M) RakitovoByal (Ro) kladenetsRyahovo (Pba) (M) Kostandovo (Ro)MalomirAlfatar (Pba) (M) Momchilovtsi (Ro)GavrailovoAsparuhovo (Pba) (M) Skrebatno (Zrup)KortenGarvan (Pba) (M) Lehovo (CYuz)RakovoPozharevo (Pba) (M) Beden (Ro)KlikachVetren (Pba) (M) Breze (Ro)OgnenKalipetrovo (Pba) (M) Babyak (Zrup) DinyaKatselovo (Pba) (M) Bansko Pshenichevo(Zrup) Dyadovo (Pba) (Pba) Dobyrsko (Zrup)DybravaLyubenovo (Pba) (Pba) GodlevoKaloyanovets (Zrup)Kozarevo (Pba) (Pba) Razlog (Zrup)KazankaTvyrditsa (Pba) (Pba) Bunovo (Ba)ZlatariAsenovtsi (Pba) (Ba) Chelopech (Ba) RosenLetnitsa (Thr) (Ba) Pirdop Rusokastro(Ba) Butovo (Thr) (M) SmolskoBylgarsko (Ba)Swetlina Slivovo (Thr) (M) Tserovo (Pba) TsarevetsEnina (Ba) (M) Pravets (CYuzS)ChernokonevoKozlovets (Thr) (M) Batuliya (CYuzS) VodenMorava (Thr) (M) Belchin (CYuz)YabylkovoDolna Studena (Thr) (M) Relyovo (CYuz)HadjidimitrovoDebyr (Thr) (M) Kondofrey (CYuz)KaradjalovoKrivina (Thr) (M) Kovachevtsi (CYuz)AleksandrovoNovograd (ZM) (M) Yarlovo (CYuz)Vylchi Styrmentryn (ZM) (M) Podgumer Gorno(CYuzS) PavlikeneBeltsov (Ba) (M) Chepintsi2 (CYuzS)BylgareneDjulyunitsa (MZ) (M) LyutibrodMalka (CYuzS) ZhelyaznaDolna Lipnitsa (Ba) (M) Stolnik (CYuz)UgyrchinLesicheri (Ba) (M) Oselna (CYuzS) GalataPatresh (Ba) (M) YablanitsaGolyam (CYuzS)Bjala Chardak Cherkova (Pba) (M) Zanoge (ZM) VidrareVishovgrad (Ba) (M) Zanozhene (CYuzS) GoletsMusina (ZM) (M) Seslavtsi (CYuzS)OreshakKaraisen (Ba) (M) Leskovdol (CYuzS)GumoshtnikSredna kula (Ba) (M) Sestrimo (Pba)GrivitsaCherven (ZM) (M) Bistrets (CYuzS)Vybel2Gabrovtsi (ZM) (Ba) Chelopek (CYuzS)KladniLisets dyal (ZM) (Ba) Dyrmantsi (CYuzS)SlavyanovoGolemanite (ZM) (Ba) Boboshevo (CYuz)BalaniteKulina voda (Ba) (M) Frolosh (CYuz)ZhylteshOsym (Ba) (M) Eremiya (CYuz)BozhentsiteDebrene (Ba) (M) Drumohar (CYuz)SabotkovtsiOpanets (Ba) (M) Smolichano (CYuz)TodorchetaAvren (Ba) (Ro) Padala (CYuz)ChernovryhPlaninets (Ba) (Ro) Bosnek (CYuz) PyrshaBelitsa2 (Ba) (Thr) Chupetlovo (CYuz)GybeneMezek (Ba) (Thr) Shipochan (CYuz)TopleshTrakiets (Ba) (Thr) Kremenik (CYuz)GanchovetsBryagovets (Ba) (Ro) Dospey (CYuz) BelitsaSvirkovo (Ba) (Thr) Madjare (CYuz)MryzetsiBylgarin (Ba) (Thr) Govedartsi (CYuz)GostilitsaBabek (Ba) (Pba) Raduil (CYuz)DraganovtsiBelozem (Ba) (Pba) Gabra (CYuz)PopovtsiVoysil (Ba) (Pba) Kostentsi (CYuzS)KozirogRupkite (Ba) (Pba) Gutsal (CYuz)FyrevtsiBebrovo (Ba) (Ba) Zhivkovo (CYuz)BatoshevoPalitsi (Ba) (Ba) Vrachesh (CYuzS) KupenKostel (Ba) (Ba) Belitsa3 (Zrup)DushevskoBujnovtsi (Ba) (Ba) Yakoruda (Zrup)KryvenikLazartsi (Ba) (Ba) Osenovo (CYuz) DedinaTodyuvtsi (Ba) (Ba) Belasitsa (CYuz)KrushareGradets (Pba) (Ba) Gabrene (CYuz)Kaloyanovo2Zheravna (Pba) (Ba) Bobeshino (YuzP)SyrnegorMedven (Ba) (Ba) Golema Fucha (CYuz)VoynyagovoChiflika (Ba) (Ba) Palatovo (CYuz)BelyovoDrenta (CYuz) (Ba) Lokvata (CYuz)GoleshovoLipov (CYuz) ryt (Ba) Mursalevo (CYuz)KaraivantsiDospat (Ro) (Ba) Novoselyane (CYuz)SyrnitsaDryanovets (Ro) (M) Ryzhdavitsa (YuzP)KasykaSvalenik (Ro) (M) Gramazhdeno (CYuz)GyrmenSadina (Zrup) (M) Lozno (CYuz) NastanKrivnya (Ro) (M) Brestovo (CYuz)HvoynaSenovo (Ro) (M) Leshko (CYuz)OrehovoOrizare (Ro) (Pba) Padezh (CYuz)GaytaninovoBrezovo (Zrup) (Pba) Brest Nova(YuzP) LovchaMalina (Zrup) (Thr) KryshalevoBylgarska (YuzP) Paril polyana (Zrup) (Thr) Lomnitsa (YuzP)TeshovoPoroy (Zrup) (Pba) Zlogosh (YuzP) LykiSrem (Zrup) (Thr) Saparevo (CYuz)Byal Rozovokladenets (Ro) (Pba) Brezhani (CYuz)BezdenMalomir (CYuzS) (Pba) Gorno Osenovo (CYuz)OpitsvetGavrailovo (CYuzS) (Pba) DebreneBogdanov (CYuz) dolKorten (YuzP) (Pba) Stozha (CYuz)MyrchaevoRakovo (CYuz) (Pba) Kresna (CYuz)Gradets2Klikach (CYuzS) (Pba) Mechkul (CYuz)Breze2Ognen (CYuzS) (Pba) Vlahi (CYuz) ShipotDinya (ZM) (Pba) Moravska (CYuz)PshenichevoSumer (ZM) (Pba) Velyushtets (CYuz) GintsiDybrava (SzP) (Pba) OshtavaGurmazovo (CYuz)Kaloyanovets (CYuzS) (Pba) Churichene (CYuz)IvanyaneKazanka (YuzP) (Pba) Gega (CYuz)SlatinaZlatari (YuzP) (Pba) Krynalovo (CYuz)VerdikalRosen (YuzP) (Thr) Kyrlanovo (CYuz)KokalyaneRusokastro (CYuz) (Thr) Lozenitsa (CYuz)LozenSwetlina (CYuz) (Thr) Rozhen (CYuz)DondukovoEnina (ZM) (Ba) LoboshRazgrad (YuzP)Chernokonevo mahala (ZM) (Thr) Vyrba (YuzP)SlavotinVoden (ZM) (Thr) Bohova (SzP)MominYabylkovo brod (ZM) (Thr) Zabel (SzP)VladichentsiDebyr (ZM) (Thr) Miloslavtsi (SzP)KaradjalovoChiren (ZM) (Thr) Slivovnik (SzP)AleksandrovoGaliche (ZM) (ZM) Protopopintsi (SzP)DinkovitsaVylchi tryn (ZM) (ZM) Boynitsa Gorno(SzP)Plakuder Pavlikene (ZM) (Ba) Bukorovtsi (SzP)BylgareneInovo (ZM) (MZ) Struindol Malka(SzP)Yanovets Zhelyazna (SzP) (Ba) ChupreneGorno (SzP) TseroveneUgyrchin (ZM) (Ba) Kopilovtsi (SzP)VladimirovoGalata (ZM) (Ba) ChepyrlintsiGolyam (SzP)Rakevo Chardak (ZM) (Pba) Stanintsi (SzP)SelanovtsiVidrare (ZM) (Ba) KomshtitsaKosta (SzP) PerchevoGolets (ZM) (ZM) Rebro (SzP) OsenOreshak (ZM) (Ba) Gorni Romantsi Manastirishte(SzP)Gumoshtnik (ZM) (Ba) Velkovtsi (YuzP)ZheglitsaGrivitsa (ZM) (ZM) Izvor2 (YuzP) FurenVybel2 (ZM) (ZM) Selishten dol (YuzP)GabareLisets (ZM) (ZM) Divlya (YuzP)SokolareSlavyanovo (ZM) (ZM) Shipkovtsi (YuzP)StavertsiBalanite (ZM) (Ba) Treklyano (YuzP)PisarovoZhyltesh (ZM) (Ba) Vidrar (YuzP)BozhentsiteKalen (ZM) (Ba) Metohiya (YuzP)VarvaraSabotkovtsi (Pba) (Ba) Stnyovtsi (YuzP)Opanets2Todorcheta (ZM) (Ba) Gorni Korten (YuzP)VyrbitsaChernovryh (CYuzS) (Ba) Bachkovo (Ro)PoibrenePyrsha (Ba) (Ba) Kribul (Zrup)OvchepoltsiGybene (Pba) (Ba) Dorkovo (Ro)StroevoToplesh (Pba) (Ba) Mugla (Ro)GanchovetsMiromir (Pba) (Ba) SelchaDebrashtitsa (Ro) Belitsa (Pba) (Ba) Strizhba (Ro)GradinaMryzetsi (Thr) (Ba) Tsarino (Ro) GostilitsaSeltsi (Ro) (Ba) Zhyltusha (Ro)DraganovtsiKonush (Ro) (Ba) Botevo (M)DolnoslavPopovtsi (Ro) (Ba) Levski (M)GornoslavKozirog (Ro) (Ba) Figure 11:Kalimantsi Rodopian (M)Oreshets dialectsFyrevtsi (Ro) (Ba) close to bordering southwestern dialects. MominoGorno (M) BatoshevoLukovo (Ro) (Ba) Ravna gora Meden(M) bukKupen (Ro) (Ba) Rudnik (M)TselinaDushevsko (Pba) (Ba) Sadovo (M)MomkovoKryvenik (Thr) (Ba) IzgrevYerusalimovo (M) Dedina (Thr) (Ba) Nikolaevka (M)OryahovoKrushare (Thr) (Pba) ZasmyanoKonstantinovo Kaloyanovo2(M) (Thr) (Pba) Voditsa (M) SyrnegorIzvor (Thr) (Ba) Shtipsko (M)VoynyagovoNedelino (Ro) (Ba) Kazichino (Ba)BryshlyanBelyovo (Thr) (CYuz) Golitsa (Ba)ByalaGoleshovo voda (Thr) (CYuz) Cherni Vit (Ba)MoryaneDospat (Thr) (Ro) Ribaritsa (Ba)ZvezdetsSyrnitsa (Thr) (Ro) Golyama Zhelyazna (Ba) FakiyaKasyka (Thr) (Ro) Ahryane (Ro)ZabernovoGyrmen (Thr) (Zrup) Bogutevo (Ro)BylgariNastan (Thr) (Ro) Osikovo (Ro)KondolovoHvoyna (Thr) (Ro) Petvar (Ro)RezovoOrehovo (Thr) (Ro) Sitovo Gramatikovo(Pba)Gaytaninovo (Thr) (Zrup) Borikovo (Ro)NovaStoilovo Lovcha (Thr) (Zrup) MogilitsaGospodintsi (Ro) Paril (Zrup) (Zrup) ChepintsiCheshnegirovo (Ro) Teshovo (Ro) (Zrup) Kozhari (Ro)CherkovnaLyki (M)(Zrup) Vodnipad (Ro) RavnaRozovo (M) (Ro) Trigrad (Ro)BezdenNenovo (CYuzS) (M) Kremene (Ro)OpitsvetKrivnya2 (CYuzS) (M) FigureVitinaBogdanov (Ro) 12:Kosten The dol (Pba) (YuzP) Hvoyna Rodopian sub-dialect. Kremene2 (Ro)MyrchaevoPodvis (Pba) (CYuz) Smilyan (Ro)Gradets2Krumovo (CYuzS) (M) Raykovo (Ro)Breze2Vybel (CYuzS) (M) Smolyan (Ro) PrilepShipot (Pba) (ZM) Raykovo2 (Ro)RosenovoSumer (Thr) (ZM) Smolyan2 (Ro)KirilovoGintsi (Thr) (SzP) UstovoGurmazovo (Ro)Dragoevo (CYuzS) (M) Gela (Ro)IvanyaneKochovo (YuzP)(M) Shiroka Lyka (Ro)KyulevchaSlatina (YuzP)(M) Stoykite (Ro)BukovoVerdikal (Ro) (YuzP) Petkovo (Ro)ErmaKokalyane reka (Ro)(CYuz) Tyryn (Ro)DavidkovoLozen (Ro)(CYuz) Chokmanovo (Ro)DondukovoVyrbina (Ro) (ZM) VievoRazgrad (Ro)Kraynovo mahala (Thr) (ZM) Beslen (Zrup)VylchanSlavotin dol (Ro) (ZM) Lydja Zagrazhden(Ro)Momin brod (Ro) (ZM) Devisilitsa (Ro)ZlatogradVladichentsi (Ro) (ZM) Malyk Devisil (Ro) BeleneChiren (M) (ZM) Golyam Devisil (Ro)DolenGaliche (Zrup) (ZM) FigureTihomir (Ro) 13:PletenaDinkovitsa The (Zrup) Zlatograd (ZM) Rodopian sub-dialect. Agatovo (Ba)OryahovetsPlakuder (Ro) (ZM) Byala reka (M) OreshInovo (M) (ZM) Koevtsi (Ba)YanovetsKutela (Ro) (SzP) Djulyunitsa2Gorno (M) BelentsiTserovene (ZM) (ZM) Pisarevo (M)DermantsiVladimirovo (ZM) (ZM) Nikyup (Ba) GigenRakevo (ZM) (ZM) PaskalevetsSinyo (M)Selanovtsi byrdo (ZM) (ZM) 24 KramolinKosta (Ba)Gulyantsi Perchevo (ZM) (ZM) Debeltsovo (Ba)RadoveneOsen (ZM) (ZM) Draganovtsi2 (Ba)Tsarevets2Manastirishte (ZM) (ZM) ShilkovtsiBryshlyanitsa (Ba) Zheglitsa (ZM) (ZM) Mindya (Ba)NikolaevoFuren (ZM) (ZM) Yakovsti (Ba)DuvanliGabare (Pba) (ZM) KipilovoKaloyanovo (Ba) Sokolare (Pba) (ZM) StaraGeneral reka (Pba) NikolaevoStavertsi (Pba) (ZM) Obnova (M)SekirovoPisarovo (Pba) (ZM) Pyrvomaitsi (M) Kalen (ZM) Varvara (Pba) Hlebarovo (M) 0.0 0.001 0.002 0.003 0.004 Tetovo (M) Opanets2 (ZM) Belyakovets (Ba)Vyrbitsa (CYuzS) Samovodene (Ba) Poibrene (Ba) Dichin (Ba)Ovchepoltsi (Pba) Dragizhevo (Ba) Stroevo (Pba) Prisovo (Ba) Miromir (Pba) Pirgovo (M)Debrashtitsa (Pba) Gagalya (M) Gradina (Thr) Ryahovo (M) Seltsi (Ro) Alfatar (M) Konush (Ro) Asparuhovo (M) Dolnoslav (Ro) Garvan (M) Gornoslav (Ro) Pozharevo (M) Oreshets (Ro) Vetren Gorno(M) Lukovo (Ro) Kalipetrovo (M)Meden buk (Ro) Katselovo (M) Tselina (Pba) Dyadovo (Pba) Momkovo (Thr) Lyubenovo (Pba)Yerusalimovo (Thr) Kozarevo (Pba) Oryahovo (Thr) Tvyrditsa (Pba)Konstantinovo (Thr) Asenovtsi (Ba) Izvor (Thr) Letnitsa (Ba) Nedelino (Ro) Butovo (M) Bryshlyan (Thr) Bylgarsko Slivovo (M)Byala voda (Thr) Tsarevets (M) Moryane (Thr) Kozlovets (M) Zvezdets (Thr) Morava (M) Fakiya (Thr) Dolna Studena (M)Zabernovo (Thr) Hadjidimitrovo (M) Bylgari (Thr) Krivina (M)Kondolovo (Thr) Novograd (M) Rezovo (Thr) Styrmen (M)Gramatikovo (Thr) Beltsov (M) Stoilovo (Thr) Djulyunitsa (M)Gospodintsi (Zrup) Dolna LipnitsaCheshnegirovo (M) (Ro) Lesicheri (M) Cherkovna (M) Patresh (M) Ravna (M) Bjala Cherkova (M) Nenovo (M) Vishovgrad (M) Krivnya2 (M) Musina (M) Kosten (Pba) Karaisen (M) Podvis (Pba) Sredna kula (M) Krumovo (M) Cherven (M) Vybel (M) Gabrovtsi (Ba) Prilep (Pba) Kladni dyal (Ba)Rosenovo (Thr) Golemanite (Ba) Kirilovo (Thr) Kulina voda (M) Dragoevo (M) Osym (M) Kochovo (M) Debrene (M) Kyulevcha (M) Opanets (M) Bukovo (Ro) Avren (Ro) Erma reka (Ro) Planinets (Ro)Davidkovo (Ro) Belitsa2 (Thr) Vyrbina (Ro) Mezek (Thr) Kraynovo (Thr) Trakiets (Thr)Vylchan dol (Ro) Bryagovets (Ro)Zagrazhden (Ro) Svirkovo (Thr) Zlatograd (Ro) Bylgarin (Thr) Belene (M) Babek (Pba) Dolen (Zrup) Belozem (Pba) Pletena (Zrup) Voysil (Pba)Oryahovets (Ro) Rupkite (Pba) Oresh (M) Bebrovo (Ba) Kutela (Ro) Palitsi (Ba) Belentsi (ZM) Kostel (Ba)Dermantsi (ZM) Bujnovtsi (Ba) Gigen (ZM) Lazartsi (Ba)Sinyo byrdo (ZM) Todyuvtsi (Ba) Gulyantsi (ZM) Gradets (Ba)Radovene (ZM) Zheravna (Ba)Tsarevets2 (ZM) Medven (Ba)Bryshlyanitsa (ZM) Chiflika (Ba) Nikolaevo (ZM) Drenta (Ba) Duvanli (Pba) Lipov ryt (Ba)Kaloyanovo (Pba) KaraivantsiGeneral (Ba) Nikolaevo (Pba) Dryanovets (M) Sekirovo (Pba) Svalenik (M) Sadina (M) 0.0 0.001 0.002 0.003 0.004 Krivnya (M) Senovo (M) Orizare (Pba) Brezovo (Pba) Malina (Thr) Bylgarska polyana (Thr) Poroy (Pba) Srem (Thr) Byal kladenets (Pba) Malomir (Pba) Gavrailovo (Pba) Korten (Pba) Rakovo (Pba) Klikach (Pba) Ognen (Pba) Dinya (Pba) Pshenichevo (Pba) Dybrava (Pba) Kaloyanovets (Pba) Kazanka (Pba) Zlatari (Pba) Rosen (Thr) Rusokastro (Thr) Swetlina (Thr) Enina (Ba) Chernokonevo (Thr) Voden (Thr) Yabylkovo (Thr) Debyr (Thr) Karadjalovo (Thr) Aleksandrovo (ZM) Vylchi tryn (ZM) Gorno Pavlikene (Ba) Bylgarene (MZ) Malka Zhelyazna (Ba) Ugyrchin (Ba) Galata (Ba) Golyam Chardak (Pba) Vidrare (Ba) Golets (ZM) Oreshak (Ba) Gumoshtnik (Ba) Grivitsa (ZM) Vybel2 (ZM) Lisets (ZM) Slavyanovo (ZM) Balanite (Ba) Zhyltesh (Ba) Bozhentsite (Ba) Sabotkovtsi (Ba) Todorcheta (Ba) Chernovryh (Ba) Pyrsha (Ba) Gybene (Ba) Toplesh (Ba) Ganchovets (Ba) Belitsa (Ba) Mryzetsi (Ba) Gostilitsa (Ba) Draganovtsi (Ba) Popovtsi (Ba) Kozirog (Ba) Fyrevtsi (Ba) Batoshevo (Ba) Kupen (Ba) Dushevsko (Ba) Kryvenik (Ba) Dedina (Ba) Krushare (Pba) Kaloyanovo2 (Pba) Syrnegor (Ba) Voynyagovo (Ba) Belyovo (CYuz) Goleshovo (CYuz) Dospat (Ro) Syrnitsa (Ro) Kasyka (Ro) Gyrmen (Zrup) Nastan (Ro) Hvoyna (Ro) Orehovo (Ro) Gaytaninovo (Zrup) Nova Lovcha (Zrup) Paril (Zrup) Teshovo (Zrup) Lyki (Zrup) Rozovo (Ro) Bezden (CYuzS) Opitsvet (CYuzS) Bogdanov dol (YuzP) Myrchaevo (CYuz) Gradets2 (CYuzS) Breze2 (CYuzS) Shipot (ZM) Sumer (ZM) Gintsi (SzP) Gurmazovo (CYuzS) Ivanyane (YuzP) Slatina (YuzP) Verdikal (YuzP) Kokalyane (CYuz) Lozen (CYuz) Dondukovo (ZM) Razgrad mahala (ZM) Slavotin (ZM) Momin brod (ZM) Vladichentsi (ZM) Chiren (ZM) Galiche (ZM) Dinkovitsa (ZM) Plakuder (ZM) Inovo (ZM) Yanovets (SzP) Gorno Tserovene (ZM) Vladimirovo (ZM) Rakevo (ZM) Selanovtsi (ZM) Kosta Perchevo (ZM) Osen (ZM) Manastirishte (ZM) Zheglitsa (ZM) Furen (ZM) Gabare (ZM) Sokolare (ZM) Stavertsi (ZM) Pisarovo (ZM) Kalen (ZM) Varvara (Pba) Opanets2 (ZM) Vyrbitsa (CYuzS) Poibrene (Ba) Ovchepoltsi (Pba) Stroevo (Pba) Miromir (Pba) Debrashtitsa (Pba) Gradina (Thr) Seltsi (Ro) Konush (Ro) Dolnoslav (Ro) Gornoslav (Ro) Oreshets (Ro) Gorno Lukovo (Ro) Meden buk (Ro) Tselina (Pba) Momkovo (Thr) Yerusalimovo (Thr) Oryahovo (Thr) Konstantinovo (Thr) Izvor (Thr) Nedelino (Ro) Bryshlyan (Thr) Byala voda (Thr) Moryane (Thr) Zvezdets (Thr) Fakiya (Thr) Zabernovo (Thr) Bylgari (Thr) Kondolovo (Thr) Rezovo (Thr) Gramatikovo (Thr) Stoilovo (Thr) Gospodintsi (Zrup) Cheshnegirovo (Ro) Cherkovna (M) Ravna (M) Nenovo (M) Krivnya2 (M) Kosten (Pba) Podvis (Pba) Krumovo (M) Vybel (M) Prilep (Pba) Rosenovo (Thr) Kirilovo (Thr) Dragoevo (M) Kochovo (M) Kyulevcha (M) Bukovo (Ro) Erma reka (Ro) Davidkovo (Ro) Vyrbina (Ro) Kraynovo (Thr) Vylchan dol (Ro) Zagrazhden (Ro) Zlatograd (Ro) Belene (M) Dolen (Zrup) Pletena (Zrup) Oryahovets (Ro) Oresh (M) Kutela (Ro) Belentsi (ZM) Dermantsi (ZM) Gigen (ZM) Sinyo byrdo (ZM) Gulyantsi (ZM) Radovene (ZM) Tsarevets2 (ZM) Bryshlyanitsa (ZM) Nikolaevo (ZM) Duvanli (Pba) Kaloyanovo (Pba) General Nikolaevo (Pba) Sekirovo (Pba)

0.0 0.001 0.002 0.003 0.004 Ruse

Pleven

Shumen Varna Lovech

Teteven

Sofia

Burgas

Plovdiv Blagoevgrad Malko Tyrnovo Razlog

Smolyan

Figure 14: The lightest areas are most similar to Standard Bulgarian, while the darkest ones are least similar. ilarity to the southeastern sites would not be surprising either, rst, because it would conrm the unity of eastern dialects with respect to the important feature of `yat' realizations, and second, because there are phonetic features that are typical for the Northeast, but not accepted in the standard language, including the reduction of open vowels and the existence of nal palatal consonants. We measured the distance from the standard to each of the 490 varieties in our sample using the same techniques described above. We illustrate the results of these measurements in Fig. 14, displaying the distances of the varieties to the standard. Several comments are in order. The general picture suggests that the most similar dialects are the central, southeastern and northwestern ones. The western and Rodopian dialects are least similar, which is not surprising bearing in mind the source of the standard language. Note that the Rodopian area (Smo- lyan) as well the west Rupskian area (Razlog) are not homogeneously similar to the standard. At the same time, it is interesting to notice that the most similar and the least similar areas are often geographically close. For example, the Teteven region shows dissimilarity near the `yat' border and similarity in the East. At the same time the Shumen region and the entire Northeast seem to be unlike the standard. Perhaps these eects are caused by the transitional dialect zones (Pleven, Teteven, Razlog). It is said that the Shumen dialect is very close to the one in Thessaloniki, Greece.

25 9 Discussion of the Larger Sample

For the sake of completeness, and in order to put our claim of consistency to the test, we briey discuss the similarities and dierences between results based on the 54-word set and those based on the 36-word set (discussed above). No new maps need to be presented, but we will wish to refer to the ones above. First, we note that the correlation between distances assigned on the basis of the 36-word set and those assigned on the basis of the 54-word set is very high indeed (r = 0.97). The high correlation indicates that the 36 common words are representative of the data in the atlases. For example, the line map for 54-word set (similar to Fig. 3), the MDS results (similar to Fig. 7), and the composite cluster maps (similar to Fig. 6) all highly resemble the results for the 36-word set. Nonetheless, the more unstable clustering gives somewhat diering results, which may be taken into consideration as well. We shall examine these results briey, even at the risk of overemphasizing less stable aspects of the analysis. The dendrogram for 54 words and 5 groups highlights the western, eastern, Rodopian and northeastern Balkan groups. The 36-word dendrogram seems to recognize more distinct areas in Bulgaria (Fig. 4). As a result, the classication map for 54 words (not shown) stresses the east-west distinction in a way that is closer to the traditional division. In the 54-word set some more dialect features are included (the ones not relevant throughout the four atlases) and these justify the east-west division more substantially. One such feature is, for example, the presence/omission of the fricative [x] in dierent positions in the word. Some of the west Moesian dialects are distinguishedin the area of Pleven and Lovech. This is again not very surprising as these dialects are at the border of the `yat' division, and they exhibit similarities to western as well as eastern dialects. The other distinct area in 54-word set is the dialect of Pavlikjans, in the Plovdiv region. It contains a lot of archaic Rodopian features. Furthermore, the clusters for Rodopi dier in number. For the 54-word set it shows only 3 compact groups of dialects: Smolyan, Zlatograd and the varieties closer to west Rupskian.

10 Conclusions and prospects

In this paper we applied dialectological techniques which had been successfully applied to Germanic and Romance language to the Slavic language, Bulgarian. The application involved including a novel set of sounds and a dialect area which is famous for its contact phenomena. Nonetheless, the techniques could be applied straightforwardly, indicating that they are more general than had been shown until now. We analyzed the relations among Bulgarian dialects based on the pronunciation distances calculated between all pairs of 490 Bulgarian varieties, and including additionally the standard language. The data was selected and digitized from the four-volume set of Atlases of Bulgarian Dialects. The relations among Bulgarian dialects turned out to be quite complex. The comparison with

26 the earlier maps and their dialect divisions showed that the most signicant border remains the `yat' realization border. However, our research also conrms Balan's claims (Teodorov-Balan 1904) about the need to introduce a third major dialect area, namely: Rodopi. With respect to the four-way distinction in the atlases, our results suggest that the western dialects are more cohesive than the eastern ones, and that within the western dialects, the North and South are additionally distinguished. The East shows weaker cohesion among dialects, and its northern and southern areas are more uniform. One explanation could be that more migrations have taken place in the East than in the West. At the same time, several special subclusters are also distinguished: the Rodopian region, the northeastern Balkan region and some interesting areas near the `yat' border. We also explored the relation of the dialects to the standard language. We found that the standard is most similar to south Balkan dialects from the North- east, which was expected. We also observed interestingly that most similar and least similar varieties need not be geographically near one another. In our opinion, this work paints a faithful picture of Bulgarian dialects even though it is based on a limited number of word pronunciations. In addition, our more inclusive data set conrmed the results of the smaller one, even while hinting at potentially dierent dialect subgroups. We see the future work in several directions. We should like to examine dierent dialect data, and in particular data collected from sites that were not selected for being purely Bulgarian. It would be important to identify the regular aspects of the distinctions at the base of the analysis here, i.e. the linguistic basis of the aggregate analysis, and, in fact, we have initiated that work in collaboration with a Ph.D. student. It would be interesting to include lexical variation (Nerbonne and Kleiweg, 2003) in a parallel analysis, and to examine the degree to which lexical dierences correlate with dierences in pronunciations. Finally, Bulgarian and the Balkan are most famous linguistically for the extensive language contact which has developed there, and it would be fas- cinating to modify and apply the quantitative techniques used here in order to explore and analyze language contact.

11 Acknowledgments

The authors would like to thank Kiril Simov for his kind help in the digitization of the data, Luchia Antonova for her valuable comments on the X-SAMPA conversion, for the selection of the Bulgarian sites and her recommendations on Bulgarian dialects in the process of our work, Christine Siedle for her kind help with the geographical coordinates and the maps, Peter Kleiweg for the software facilities and his quick reactions on software questions. We also owe thanks to Renee van Bezooijen, Charlotte Gooskens, Marco Spruit, Bill Kretzschmar and the editors of Zeitschrift f¤urslavische Philologie for very useful discussion and suggestions on a preliminary version of this paper.

27 References

Bolognesi, Roberto & Wilbert Heeringa (2002), `De invloed van dominante talen op het lexicon en de fonologie van sardische dialecten', Gramma/TTT: Tijdschrift voor Taalwetenschap 9(1), 4584. URL: http://www. let.rug.nl/˜heeringa/dialectology/ papers/ Dobrev, Ivan (1981), `Ìÿñòîòî íà áúëãàðñêèÿ åçèê ñðåä äðóãèòå ñëàâÿíñêè åçèöè', Ñúïîñòàâèòåëíî åçèêîçíàíèå IV, 314.

Gooskens, Charlotte & Wilbert Heeringa (2004), `Perceptive evaluation of levenshtein dialect distance measurements using norwegian dialect data', Lan- guage Variation and Change 16(2), 189207.

Heeringa, Wilbert (2004), Measuring Dialect Pronunciation Dierences using Levenshtein Distance, PhD thesis, University of Groningen, Groningen. URL: http://www.let.rug.nl/˜heeringa/dialectology/thesis/ Hoppenbrouwers, Cor & Geer Hoppenbrouwers (2001), De indeling van de Nederlandse streektalen: Dialekten van 156 steden en dorpen geklasseerd volgens de FFM, Koninkljike Van Gorcum B.V., Assen.

International Phonetic Association (2003), Handbook of the IPA. A Guide to the Use of the International Phonetic Alphabet, Cambridge University Press, Cambridge.

Jain, Anil K. & Richard C. Dubes (1988), Algorithms for Clustering Data, Pren- tice Hall, Englewood Clis, New Jersey.

Kleiweg, Peter, John Nerbonne & Leonie Bosveld (2004), Geographic projection of cluster composites, in A.Blackwell, K.Marriott & A.Shimojima, eds, `Diagramatic Representation and Inference. Third International Con- ference, Diagrams. Cambridge, UK, March 2004 Lecture Notes on Articial Intelligence 2980', Springer, Berlin, pp. 392394.

Ko§cev,Ivan (1988), Áúëãàðñêè äèàëåêòåí àòëàñ, Vol. 2 of Îáîáùàâàù òîì, Publisher of the Bulgarian Academy of Science, Soa.

Kruskal, Joseph & Myron Wish (1978), Multidimensional scaling, Number 07- 011 in Sage University Paper Series on Quantative Applications in the Social Sciences, Sage Publications, Newbury Park.

Nerbonne, John & Christine Siedle (2005), `Dialektklassikation auf der Grundlage aggregierter Ausspracheunterschiede', Zeitschrift f¤urDialekto- logie und Linguistik 72(2), 129147.

Nerbonne, John & Peter Kleiweg (2003), `Lexical distance in LAMSAS', Com- putational Methods in Dialectometry. Special issue of Computers and the Humanities 37(3), 339357.

28 Nerbonne, John & Wilbert Heeringa (2001), `Computational comparison and classication of dialects', Dialectologia et Geolinguistica 9, 6983. Spec. iss. ed. by Jaap van Marle and Jan Berens w. selections from 2nd Int'l Congress of Dialectologists and Geolinguists, Amsterdam, 1997.

Nerbonne, John, Wilbert Heeringa, Erik van den Hout, Peter van der Kooi, Simone Otten & Willem van de Vis (1996), Phonetic distance between dutch dialects, in G.Durieux, W.Daelemans & S.Gillis, eds, `CLIN VI, Papers from the sixth CLIN meeting', University of Antwerp, Center for Dutch Language and Speech, Antwerpen, pp. 185202. URL: http://www.let.rug.nl/˜heeringa/dialectology/ papers/ Popov, Dimit¢ar, Kiril Simov & Svetlomira Vidinska (1998), Ðå÷íèê çà ïðàâîãîâîð, ïðàâîïèñ, ïóíêòóàöèÿ, Atlantis, Soa.

P§seni§cnova, N. N. (1973), `Ñòàòèñòè÷åñêèé àíàëèç ìíîãî÷ëåííîãî äèàëåêòíîãî ñîîòâåòñòâèÿ (Ãëóõèå äîëãèå øèïÿùèå â ðóññêèõ ãîâîðàõ)', Âîïðîñû ÿçûêîçíàíèÿ 6, 8091.

P§seni§cnova, N. N. (1977), Ïðèìåíåíèå ìåòîäà òàêñîíîìè÷åñêîãî àíàëèçà ê êëàñèôèêàöèè ãîâîðîâ, in `Äèàëåêòè÷åñêèå èññëåäîâàíèÿ ïî ðóññêîìó ÿçûêó', Moscow, pp. 314.

Simov, Kiril, Alexand¢arSimov, Xristo Ganev, Krasimira Ivanova & Ilko Grig- orov (2004), The CLaRK system: XML-based corpora development system for rapid prototyping, in `Proceedings of LREC 2004', Lisbon, Portugal, pp. 235238.

Stojkov, Stojko (1966), Áúëãàðñêè äèàëåêòåí àòëàñ: Ñåâåðîèçòî÷íà Áúëãàðèÿ, Vol. 2, Publisher of the Bulgarian Academy of Science, Soa.

Stojkov, Stojko (2002), Áúëãàðñêà äèàëåêòîëîãèÿ, Publisher of the Bulgarian Academy of Science, Soa. ÷åòâúðòî èçäàíèå.

Stojkov, Stojko, Ivan Ko§cev& Maksim Mladenov (1981), Áúëãàðñêè äèàëåêòåí àòëàñ: Ñåâåðîèçòî÷íà Áúëãàðèÿ, Vol. 4, Publisher of the Bulgarian Academy of Science, Soa.

Stojkov, Stojko, Kiril Mir§cev, Ivan Ko§cev & Maksim Mladenov (1975), Áúëãàðñêè äèàëåêòåí àòëàñ: Þãîçàïàäíà Áúëãàðèÿ, Vol. 3, Publisher of the Bulgarian Academy of Science, Soa.

Stojkov, Stojko & S. B. Bernstein (1964), Áúëãàðñêè äèàëåêòåí àòëàñ: Þãîèçòî÷íà Áúëãàðèÿ, Vol. 1, Publisher of the Bulgarian Academy of Science, Soa.

Teodorov-Balan, Aleksand¢ar(1904), Ðîäîïñêîòî íàðå÷èå, in `Ñáîðíèê ñòàòèé ïî ñëàâÿíîâåäåíèþ, ïîñâåùåííûõ ïðîô. Ì. Ñ. Äðèíîâó', Õàðêîâ, pp. 111127.

29 Wells, John (2003), `Computer-coding the IPA: a proposed extension of SAMPA'. URL: http://www.phon.ucl.ac.uk/home/sampa/x-sampa.htm