Abstract
The maximal covering location problem was first introduced by Church and ReVelle in 1973 at the North American Regional Science Council meetings and subsequently published in Papers in Regional Science (formerly Papers of the Regional Science Association) in 1974. It has proven to be a seminal contribution to location analysis and modeling, in terms of both technical merit and practical significance. With more than 1,500 citations in the academic literature, it has truly stood the test of time and may actually be as relevant or more so today than when it was first presented/published. Not only is it the subject of broad application and extension, but it has also been integrated in a number of geographic information system–based commercial software packages, including ArcGIS and TransCAD, for general use. This article provides an overview of the maximal covering location problem, highlighting the use, application, solution, evolution, and extension of this important location analytic approach.
Keywords
Introduction
Location analysis and modeling has a long history of theoretical and methodological development that has produced a wide array of methods to support analysis, planning, and decision making in public- and private-sector contexts. Referred to as location theory historically and location analytics in more contemporary contexts, location analysis and modeling has been used both to support prescriptive planning efforts seeking to site service facilities and to provide descriptive insights regarding service performance of existing operations. In both descriptive and prescriptive capacities, location theory has had a major role in regional science, with a long history founded upon the work of von Thunen, Weber, Christallar, Losch, Hotelling, among others (Ghosh and Rushton 1987; Isard 2003; Murray 2010). Overviews of location modeling can be found in Owen and Daskin (1998), Hamacher and Drezner (2002), Church and Murray (2009), Murray (2010), and Eiselt and Marianov (2011).
A primary category of location analytic approaches involves coverage, where a facility provides service based on spatial proximity. Facilities might correspond to fire stations, health care clinics, cellular antennae, or a chain of restaurants. Each facility then has a spatial footprint that reflects the area that it primarily “serves.” This may be based on travel time, as is the case for police, ambulance, and fire response, or could be based on line of sight or audibility factors. Such spatial footprints reflect service and may be regular or irregular in shape and contiguous or fragmented in terms of areal extent. One of the most highly cited location modeling approaches is to optimize coverage, the maximal covering location problem (MCLP), detailed in Church and ReVelle (1974). A concise problem statement is: Maximize coverage of demand within a desired service standard S by locating a fixed number of facilities. This approach has proven to be a seminal contribution in regional science and beyond, in terms of both technical merit and practical significance (see Isserman 2003; Ratick, Osleeb, and Si 2016; Snyder and Haight 2016; Wei 2016). As a result, it has been included in the Modern Classics in Regional Science series edited by Haynes et al. (1996). With more than 1,550 citations to date in the academic literature, it has truly stood the test of time and may actually be as relevant or more so today than when it was first published. Not only is it the subject of broad application and extension, but it has also been integrated in a number of geographic information system (GIS)-based commercial software packages, including ArcGIS and TransCAD, given its general utility.
The assumption in the original specification of the model is that a network of nodes and arcs represent geographic space, where nodes are concentrations of anticipated service demand and/or potential facility sites and arcs reflect travel between nodes. With this in mind, notation for use in mathematical specification is now introduced:
i = index of demand nodes (entire set I),
j = index of potential facility sites (entire set J),
dij
= shortest distance or travel time from demand node i to potential facility j,
ai
= population to be served at demand node i,
p = number of facilities to be located,
Ni
= set of potential facilities that suitably cover demand node i,
S = service distance or time standard.
Through evaluation and assessment, it is possible to identify a priori which potential facility sites suitably serve/cover any particular demand node i, the result of which is reflected in the set Ni . Such assessment involves determining whether a potential facility site j is within the service distance/time standard S of node i, or formally if dij is less than S. Worth noting is that coverage may be nonspatial as well, derived based upon other criteria (see Church, Stoms, and Davis 1996; Ando et al. 1998). In addition to the above notation, decision variables are needed to track or account for facility siting decisions as well as to assess whether demand nodes are suitably served. These variables are as follows:
With notation and decision variables, the mathematical formulation of the MCLP is:
The objective, (1), is to maximize total demand in the region suitably served or covered. The Yi
decision variable accounts for whether demand node i is served by a sited facility, so the goal is a maximization of the demand weighted sum of those nodes served by sited facilities,
Unknown for a particular application of this model is which potential facility sites to select,
The abstraction of geographic space, reducing the complexities of the real world to something more simplified and manageable, has been important for the MCLP, and location theory and regional science more generally. Because demand is conceived to be at a discrete and finite number of points (nodes), with travel-along arcs of the network, this makes assessment and interpretation of coverage rather straightforward. A facility node j is either able to serve a demand node i, dij ≤ S, or it is not. Yet, there is evidence that point simplification and/or manipulation introduces unintended measurement and interpretation errors (Daskin et al. 1989; Current and Schilling 1990; Tong and Murray 2009). Church (1999) and Murray (2010) highlight, in fact, that demand and facility sites could be any sort of object (e.g., point, line, polygon, etc.). Recent research demonstrates specific issues associated with the application of the MCLP. Worth noting in particular is the work by Tong and Murray (2009), Alexandris and Giannikos (2010), Cromley, Lin, and Merwin (2012), and Yin and Mu (2015) suggesting that non-point-based contexts cannot be sufficiently addressed under classic assumptions through the direct use of the MCLP.
The significance and implications of advancements in addressing coverage optimization for regional science are important and the motivation for this article. Enhanced computing enables faster system response, larger study regions examined, and more variables to be considered. GIS proliferation and widespread adoption mean that there now exist better, more detailed data. These data are based on less simplification and more accurate abstractions of reality. As a result, for many planning and policy contexts, there is a better understanding of spatiotemporal, behavioral, and other kinds of relationships and political processes. From a modeling perspective, this has meant that analytical tools must evolve, and the MCLP is a good example of this. Simple point-based representation of geographic space has proven to be insufficient, leading to errors and uncertainty in analysis (see Tong and Murray 2009; Cromley, Lin, and Merwin 2012). Nevertheless, there are ways to address the fundamental goals of regional coverage optimization, essentially involving the use of better geographic information as well as taking advantage of derived spatial knowledge/insights. This is no doubt important given the practical significance of the MCLP as it is capable of encapsulating analysis, planning, and decision-making contexts, but also that the MCLP is reflective of other commonly relied upon regional science methods.
This article provides an overview of an important location modeling and analysis approach, the MCLP. The next section reviews the wide range of application contexts that have been addressed using this model. This is followed by a description of solution approaches that have been relied upon. The evolution of model and solution capabilities attributable to advances in GIScience are then examined. Noteworthy extensions are then detailed. This article ends with discussion and concluding comments.
Application
The MCLP has been particularly remarkable for the breadth of planning, management, and decision-making contexts to which it has been applied. Chung (1986) and Schilling, Jayaraman, and Barkhi (1993) summarize some of the various planning milieus that the MCLP has been utilized to address. Considerably more work has appeared since. A select summary of work relying on the MCLP is given in Table 1 in order to highlight the array of applications, beginning in the 1970s to now.
Application Evolution of Reported MCLP Usage for Planning and Analysis.
Note: MCLP = maximal covering location problem.
Table 1 summary identifies the general application context (Application), the facility being sited (Facility), and key reference(s) reporting this application. Again, the table is not meant to be comprehensive in terms of representing all of the various application areas utilizing the MCLP nor does it include all references in the literature where similar studies may have been undertaken. The intent, rather, is to merely give an overview that highlights the range and diversity of MCLP application. Historically, there has been much reliance on the MCLP to support emergency response contexts, such as the positioning and dispatching of ambulances, siting fire stations, policing, and so on. In fact, Church and ReVelle (1974) intimate that the MCLP would be particularly valuable in such areas when it was originally proposed and formulated. Of course, application has extended well beyond, with the MCLP being used to support air quality and natural resource management, advertising, fabrication, image processing, and even dentistry, among others. Another point to be taken from Table 1 is the continued relevance of the MCLP over the past four decades, with new applications continuing to emerge.
Figure 1 displays the historical trend of citations to Church and ReVelle (1974) over the years as reported by Google Scholar, with 1,534 total citations reported through 2014 (as of May 11, 2015). Already in 2015, there are 47 citations, but this is not included in Figure 1 as the year is not complete. To give a sense of how the MCLP continues to be used, an examination of the 152 reported studies in 2014 was carried out. A coarse characterization is that emergency services were some 28 percent of the contexts. Theory-oriented work extending the MCLP, relaxing assumptions, establishing properties, and so on, represented 18.4 percent of studies. Service facility siting, such as retail stores, was the theme of 13.6 percent of the cases. In 7.5 percent of the studies, public facilities were of interest, like libraries, parks, and so on. The remaining 32 percent was broken down in the following categories: clinics (6.7 percent), transportation (5.4 percent), fueling/refueling stations (4.7 percent), algorithm/heuristic development (4.7 percent), sensor siting (4 percent), nature reserve/wildlife protection (3.3 percent) and other (districting, telecommunication, sirens, and disease surveillance; 3.3 percent).

Citations to Church and ReVelle (1974). Source: Google Scholar May 11, 2015.
In order to provide more depth about the way in which the MCLP is being utilized to support planning and decision making, details of select problem applications are now provided. The three illustrative studies focus on telecommunications, surveillance, and emergency services, respectively.
An evaluation of telecommunications was examined in Grubesic and Murray (2002), specifically broadband access provided through digital subscriber lines, xDSL. Of particular interest was emerging issues of an Internet digital divide between those with quality high-speed access and those with comparatively poor-quality access. The focus was on which telephone company central offices in an urban region providers would be equipped with xDSL, given limited capital to invest in equipment. A major implication of this investment is that only those areas within 18,000 feet of an equipped central office would have quality access to xDSL. The MCLP was applied to assess potential inequities that would emerge under different development conditions and strategies.
Policing was considered in Murray and Tong (2007) with monitoring sensors utilized for surveillance of an area or region. The sensors record information and feed this to a control room where it is processed and evaluated by human personnel. In this case, the police department of a university campus monitors a wall display panel with live feeds from positioned video cameras. The visibility possible from a camera is a function of type, height, and physical obstructions. The MCLP was used to select sites for a limited number of sensors, video cameras, in order to maximize surveillance potential.
The evaluation, assessment, and planning of fire stations was reported in Murray, Tong, and Grubesic (2012). A community experiencing continuing population growth and development sought to make future station location decisions given a limited budget. They needed to expand their current fire response system and wanted to be able to answer a call for service within five minutes to the greatest extent possible. To accomplish this, the MCLP was used to identify siting configurations satisfying future needs.
The Essential Air Service program was examined in Grubesic, Murray, and Matisziw (2013). This program was a response to the Airline Deregulation Act of 1978, effectively providing airports in rural/remote communities a financial bridge to support continued commercial airline service. Fiscal deficits and tough economic times prompted planning for decreased services. Since the Essential Air Service program seeks to serve people and businesses within seventy miles of rural airport, the MCLP was used to assess system efficiency and identify effective strategies for continued support of select airports, given limited funds.
This section provided an overview of the various application contexts that the MCLP has been utilized to address. Clearly, the areas are growing and diversifying, reflecting the continued relevance and significance of this particular analytical approach. Citation trends suggest this will continue over the coming decades.
Solution
With such a broad range of problem contexts being addressed, there continues to be a need for better, faster, and more efficient techniques for solving the MCLP. Exact and heuristic solution techniques for general coverage models are discussed in Schilling, Jayaraman, and Barkhi (1993) and Farahani et al. (2012), with specific mention of some that have been used to solve the MCLP. ReVelle, Scholssberg, and Williams (2008) provide one of the more recent reviews of methods that expressly solve the MCLP. A summary of many of the methods that have been employed are included in Table 2, along with associated references reporting the developed approaches. The noteworthy point to be taken from Table 2 is that the application of the MCLP has required the development of specific approaches tailored to more efficiently solve problem instances. This is not particularly surprising because the MCLP is NP-complete, meaning that a given application instance is likely to be difficult to solve, optimally or approximately.
Solution Approaches for the MCLP.
Note: MCLP = maximal covering location problem.
A standard exact approach for deriving solutions of an application instance of the MCLP has been to rely on commercial general purpose optimization software, such as Gurobi Optimizer, IBM ILOG CPLEX, FICO Xpress, and so on. Similarly, open source software exists as well, like LP_Solve. Using these approaches effectively involves setting up the problem using a scripting or programming language, then call a solver library associated with the software package. Typically, problems solved along these lines would make use of linear programing or linear-integer programming (branch and bound). Assuming convergence, the problem solutions would be exact, or provably optimal.
Another approach for solving the MCLP is to use a heuristic. While not provably optimal, heuristics offer approximate solutions that are hoped to be of high quality. Big data and real-time planning environments have meant that finding good solutions fast is of great practical value. For this reason, there remains a need for improved solution techniques, reflected in the continued development reported in Table 2.
It is worth mentioning that access to and solution of the MCLP is also possible in commercial GIS software packages. For example, structuring and solving an MCLP application instance in ArcGIS and TransCAD, among others, is easily accomplished. In ArcGIS, there are location-allocation functions in network analyst, and in particular “maximize coverage” would enable an MCLP to be set up as the problem type. Similarly, in TransCAD, one would use “facility location analysis” to accomplish this. There are in fact a number of reported application studies that have relied on this basic approach using GIS, such as Gerrard et al. (1997) and Garcia-Palomares, Gutiérrez, and Latorre (2012). It is important to note that GIS packages facilitating access to the MCLP typically solve problems using a heuristic (see Murray 2010).
Evolution
The evolving use and application of the MCLP is particularly evident in Tables 1 and 2 as well as in Figure 1. In many ways, this has happened in concert with the development, growth, and maturation of GIS. Church and Murray (2009) and Longley et al. (2015), among others, discuss that GIS is a combination of hardware, software, and procedures enabling data management and spatial query. Key components of GIS include data capture, modeling, management, manipulation, analysis, and display. What is noteworthy about GIS beyond the basics of generating and dealing with spatial information are the capabilities that it provides for supporting a range of analytical functions, many of which facilitate MCLP specification and application. Highlighted in Murray (2010) are capabilities in GIS to undertake suitability analysis as well as operations associated with containment, overlay, distance, buffering, and spatial interpolation, among others. This enables feasible potential facility sites, J, to be identified. As well, potential demand (a i ), proximity (dij ), and identification of potential facility sites capably of serving an area (Ni ) can also be readily derived.
Much spatial information now exists, expediting coverage location analysis as well. While a discussion of geographic information and how to obtain it can be found in GIS texts, such as Church and Murray (2009) and Longley et al. (2015), geolibraries or geoportals often curate federal, state, and/or local government produced spatial information. One example is DATA.GOV. Beyond this, GPS, satellites, aircraft, and drones represent a rich source of data. Vendors too have a variety of detailed data. Finally, there are much individual user-generated spatial data, often referred to volunteered geographic information. The point is that far more geographic data are accessible than it was when the MCLP was first conceived in the early 1970s.
Enhanced computing capabilities, GIS, and a wealth of spatial information combine to facilitate carrying out analysis in real time, enabling the MCLP to be applied in creative ways. Add to this a better understanding of practical problem contexts with less abstraction/simplification and the stage has been set for broad utilization of the MCLP. It is an important modeling approach that reflects planning and management nuances and can be supported by data and analytical tools.
Another facet of the evolution of the MCLP is a recognition of potential issues associated with its application, specifically errors that may arise due to the representation of geographic space. Early research by Daskin et al. (1989) and Current and Schilling (1990) recognized that results and insights derived using the MCLP may be impacted by spatial aggregation of input data (demand areas), a common data management and simplification approach. Tong and Murray (2009), Alexandris and Giannikos (2010), Tong and Church (2012), Cromley, Lin, and Merwin (2012), and Yin and Mu (2015) demonstrated that the MCLP is sensitive to how space is abstracted, either as points or as polygons. Such findings associated with aggregation and/or abstraction are suggestive of modifiable areal unit problem (MAUP) impacts (see Murray 2005). That is, the analytical method can be manipulated or varies based upon a change in how geographic space is represented. What this means is that error/uncertainty is likely to be introduced by the use of a digital approximation of geographic space when the MCLP is applied.
In many ways, the MCLP is an attempt to discretize a planning context that may be viewed as a continuous space problem. MAUP and abstraction issues, therefore, are not too surprising. To address this, Murray et al. (2008) and Wei and Murray (2015) detail the continuous space MCLP. A formulation is as follows:
where
fi
( ) = function indicating the amount of demand i served by a given set of sited facilities, Δ = set of sited facilities (to be selected by the model), (ϕ
j
, λ
j
) = decision variables indicating the locational coordinates of sited facility j,
Ai
= decision variable indicating how much demand in area i is suitably served.
This formulation essentially reflects what is mathematically structured in the original specification of the MCLP using discrete points as demand and discrete points as potential facility locations. It also formalizes the description of a continuous space MCLP offered in Watson-Gandy (1982), though it was referred to as the “m-partial cover problem.” What differs is that now demand is non-point based and is represented as a continuous subarea of a region. Further, potential facility sites may be anywhere in the region, not limited to predefined points. The continuous space surface in Figure 2 illustrates this, where demand varies and feasible sites are anywhere in the region. With this in mind, the objective, (5), is to cover as much potential demand as possible. Constraint (6) accounts for how much demand in each area i is suitably served, derived using integration. The number of facilities to be sited is limited to p in constraint (7). Worth noting is that | | indicates the number of members in a set. Finally, conditions on decision variables are specified in constraints (8) and (9).

Continuous space region (darker color indicates higher demand density).
The formulation of the continuous space MCLP, (5)–(9), is inherently nonlinear. It is therefore very challenging to solve, heuristically or exactly. As a result, many special cases of this general model have been examined and solved in the literature. Friedman (1976), Watson-Gandy (1982), Mehrez and Stulman (1982), and Church (1984) consider the case where facilities may be sited anywhere in continuous space, but demand is discrete and point based. Further, the service coverage of a facility is a circle of given radius. Murray and Tong (2007) extend this case to account for demand as any spatial object, point, line, or polygon, but also allow for the service coverage of a facility to be of any shape. Murray et al. (2008) and Matisziw and Murray (2009) examine the continuous space MCLP assuming that demand is uniformly distributed. Cromley, Lin, and Merwin (2012) and Wei and Murray (2014) study the case when demand is continuously distributed and potential facility sites are discrete, finite, and known a priori.
What is noteworthy about the special cases considered by Mehrez and Stulman (1982), Church (1984), Murray and Tong (2007), Murray et al. (2008), Matisziw and Murray (2009), Cromley, Lin, and Merwin (2012), and Wei and Murray (2014) is that they reflect application contexts of importance but also are amenable to solution approaches exploiting problem features, assumptions, and so on. Interestingly, this exploitation is largely facilitated by the use of GIS in some way. Three GIS functions in particular stand out: overlay, finite dominating set (FDS), and skeleton.
Mentioned earlier was the work of Cromley, Lin, and Merwin (2012) and Wei and Murray (2014) considering continuously distributed demand and potential facility sites that are given. Through the use of polygon overlay, Cromley, Lin, and Merwin (2012) and Wei and Murray (2014) are able to discretize continuous space while also accounting for demand heterogeneity. Discussion and description of overlay can be found in de Berg et al. (2000), Church and Murray (2009), and Longley et al. (2015). A formal description is given two layers of nonoverlapping objects (polygons), L
1 and L
2, then O(L
1
, L
2) is the overlay operator producing a new layer, L, of faces f ∊ L that are a maximally connected subset

Overlay process.
As detailed previously, Mehrez and Stulman (1982), Church (1984), and Murray and Tong (2007) consider the case where facilities may be sited anywhere in continuous space, but demand is represented as discrete objects. While this technically means that facilities could be sited at an infinite number of locations, Mehrez and Stulman (1982), Church (1984), and Murray and Tong (2007) are able to exploit spatial knowledge about areas that offer enhanced coverage over other areas. This is effectively an FDS, a discrete set of locations that can be proven to contain an optimal solution to the continuous-space MCLP. Murray and Tong outline the general process for identification of the FDS, where first one must identify demand objects that are to be served (given), then derive the coverage boundary for each demand object, and finally identify a point in the intersection of coverage boundaries. The intersection of coverage boundaries is the by-product of overlay, in fact. This set derived from overlay faces represents those facility sites that only need to be considered, reducing infinite continuous space to a finite set of points. A representative point for each face is shown in Figure 4 and would constitute an FDS based upon the abovementioned conditions. With the FDS, it is then possible to apply the MCLP as originally formulated, (1)–(4). Effectively, the continuous-space MCLP is reduced to an MCLP through the identification and use of the FDS.

Representative point for each finite dominating set (FDS) polygon/face.
Recall that Murray et al. (2008) and Matisziw and Murray (2009) examined the case where demand is uniformly distributed. Observing that the most coverage would be achieved by siting in central locations of a region, Murray et al. (2008) and Matisziw and Murray (2009) were able to prove under certain conditions that the optimal single facility location was along the skeleton (or medial axis) of a region. The skeleton is the locus of points that are the centers of all disks maximally inscribed within a region, where inscription corresponds to touching the region boundary at two or more locations. Derivation of a polygon skeleton is a standard GIS function. The significance of this is that infinite space is reduced to line segments. The skeleton for a region is given in Figure 5. Combined with intelligent line search techniques, this enables an optimal solution to be obtained in some cases (single facility) and good solutions found using heuristic techniques in other cases (p facilities).

Polygon skeleton (medial axis).
These examples highlight the evolution of the MCLP from a discrete problem specification to a more generalized continuous space MCLP. Because of properties, conditions, rich spatial data, and derived knowledge from geographic information, a number of techniques have been developed and employed to solve special cases of the continuous-space MCLP, both exact and heuristic. It is clear that much more work will continue along these lines because there are limited capabilities for solving the general case problem, the continuous-space MCLP.
Extensions
As noted previously in this article, there are many coverage model reviews in the literature, such as Daskin, Hogan, and ReVelle (1988), Schilling, Jayaraman, and Barkhi (1993), Erkut et al. (2009), and Farahani et al. (2012), among others. What can be inferred is that there is both generalization of the MCLP and extension. Accounting for continuous space, in terms of either potential demand or potential siting of facilities, could be considered both generalization and extension. Above all, the distinction of discrete and continuous space was recognized as part of the evolution of the MCLP, given data and computing advances. Admittedly, this could be cast in other ways. In this section, a discussion is offered to support that there are most certainly other types of generalization and/or extension that can be recognized associated with the MCLP. This arises for many reasons, and could be due to theory, application need, model specification, solution advancement, and so on. Specifically highlighted here are six broad categories of extension. This categorization ranges from parameter and constraint changes to the basic MCLP to the more technically substantial involving probabilistic functions of certainty/uncertainty. It is worth pointing out, however, that even relatively minor modification can be both interesting and computationally significant, necessitating further research and specialized solution approaches.
Parameter and constraint modification
One category of MCLP extension has involved subtle and not so subtle modifications to the MCLP. A combination of the MCLP with the location set covering problem, referred to in Church and ReVelle (1974) as the MCLP with mandatory closeness constraints, is an example of this. In essence, this involves the addition of a set of constraints to the MCLP, (1)–(4), where each demand must be within a designated service standard of a sited facility. Another example is the addition of service facility capacity constraints, as reviewed in Pirkul and Schilling (1991). Noteworthy in this case is that incorporation of capacity constraints on facilities and/or workload has proven to make problem solution much more difficult. Different measures and metrics representing demand have also been considered, such as the criticality index in Oztekin et al. (2010) and data envelopment analysis–based efficiency scores in Grubesic et al. (2016). Finally, recent exploration of temporal dimensions of service and decision making can be found in Zarandi, Davari, and Sisakht (2013).
Multiple objectives
There has been considerable focus in MCLP extension on models that incorporate different objectives, possibly with associated constraints. Examples include the two different demand objectives in Schilling et al. (1980), the weighted benefit MCLP discussed in Church and Roberts (1983), backup models reviewed in Hogan and ReVelle (1986), the combined primary/backup coverage and capacity model of Araz, Selim, and Ozkarahan (2007), and the incorporation of alternative demand metrics in Oztekin et al. (2010), among others. These modeling extensions build on the MCLP, often by adding one or more coverage attributes to consider. The most significant issue is that the problem now becomes more complicated as many objectives must be considered. How to do this is not simple and remains an active area of research. A host of multiple objective methods for identifying trade-off solutions exist, with the weighting method and constraint method being among the most common. However, the associated problem instances may be difficult to solve, and understanding trade-off solution space is not trivial.
Probabilistic certainty
Much of the extension of MCLP has been devoted to accounting for aspects of probabilistic certainty. This ranges from demand for service being uncertain to possibilities that a service facility may not be available when needed. As an example, an ambulance is out on an emergency call when another call for service arrives reflects such a scenario. To address these types of situations, extensions of the MCLP have attempted to account for various aspects of uncertainty, including the maximum expected and maximum availability approaches reviewed in Daskin, Hogan, and ReVelle (1988). Queuing-based approach has also been developed (see Marianov and Serra 1998). A comparative study of select approaches in this area is offered in Erkut et al. (2009). Other forms of accounting for uncertainty also include the notion of hedging, as structured within an MCLP framework by Ratick, Osleeb, and Si (2016).
Service coverage
Observed behavior and performance suggests that in some application contexts, service is not uniform up to the prespecified distance or time standard S, but rather could degrade in some manner, be amplified, or could otherwise be impacted by local conditions. Church and Roberts (1983) focused on one aspect of this for the MCLP by varying coverage levels, where different time/distance intervals account for different service quality characteristics. A more generalized overview of this as well as associated research can be found in Berman, Drezner, and Krass (2010), highlighting gradual, cooperative, and variable radius approaches for extending the MCLP. Davari, Zarandi, and Turksen (2013) and Averbakh et al. (2014) offer approaches for solving one or more combinations of these cases.
Path structure
An important class of MCLP extension involves simultaneous routing and siting of facilities that provide coverage. Current, ReVelle, and Cohon (1985) sought to minimize route path length and maximize potential demand covered by facilities located along the route. Multiple path coverage was considered in Wu and Murray (2005) while route extension along these lines was structured in Matisziw, Murray, and Kim (2006). A timber harvesting context was evaluated in Bont, Heinimann, and Church (2015). Finally, Lin and Wong (2014) incorporate equity and coverage considerations in bus feeder design. Certainly, a challenge is combining routing with coverage planning as both are difficult spatial optimization problems in their own right.
Interdiction
While probabilistic certainty recognizes that a facility might be busy or in use, and therefore would delay response, interdiction reflects that a facility is no longer able to provide any service due to failure, damage, sabotage, and so on. Church, Scaparra, and Middleton (2004) proposed an extension of the MCLP to identify a subset of facilities that degrade system coverage the most. O’Hanley and Church (2011) extend this work to explicitly account for pre- and post-interdiction levels of coverage and detail solution algorithms. Lei, Tong, and Church (2014) extend the MCLP to account for spatially varying facility interdiction/failure probabilities. The goal in interdiction modeling is to identify weakness in the system in some cases while in others, it is to plan for service coverage patterns with damage or sabotage in mind. In either case, the intent is facility hardening and/or design that is more resilient to failure/damage.
While this review could continue with many additional categories and distinguishing characteristics of MCLP extensions, due to space restrictions the categories will be limited to those mentioned earlier. Recent reviews on coverage modeling offer other interpretations and additional insights.
Conclusions
The MCLP introduced in Church and ReVelle (1974) has proven to be a seminal contribution to location analysis and modeling. Subsequent work has applied the MCLP to a range of problem contexts, developed better solution techniques, relaxed assumptions and established additional theoretical foundations, and extended the MCLP in various ways. Arguably, enhanced computing and improved spatial data have contributed to a better understanding of problems and as a result have led to extensions of the MCLP. What has not changed, however, is the importance of coverage goals underlying many problem contexts when service facilities are being sited and/or evaluated. Providing service to the greatest amount of potential demand with limited resources resonates. The areas of MCLP application are growing and diversifying, reflecting the continued relevance and significance of this particular coverage model. Supporting model extensions and solution development remain strong as well. Citation trends suggest this will continue over the coming decades with many exciting developments associated with the MCLP lying ahead.
Footnotes
Declaration of Conflicting Interests
The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.
Funding
The author(s) received no financial support for the research, authorship, and/or publication of this article.
