Modeling Spatial Determinants of Blue School Certification: A Maxent Approach in Mallorca

1. Introduction

Blue Schools are educational institutions committed to advancing ocean literacy (OL) and promoting sustainable practices related to both marine and freshwater environments. They play a key role in teaching students about the importance of aquatic ecosystems while encouraging active participation in conservation and sustainability initiatives. According to the European Commission and UNESCO, a Blue School “brings the ocean into the classroom” through project-based learning, fostering environmental awareness and a sense of responsibility toward marine life and human–ocean interactions [1,2]. This means that everyday subjects, from science and geography to art and history, can incorporate ocean-related topics, helping students understand our deep connection with ocean and water systems.

The Blue Schools concept builds upon earlier efforts to embed ocean themes in education. A notable foundation was the EU-funded Sea Change project (2015–2018), which positioned OL as a strategic educational objective in Europe. Through public engagement campaigns and school-based pilots, Sea Change demonstrated the potential of marine-themed learning in formal education and laid the groundwork for what would later become national and international Blue School initiatives [3].

A pivotal development in the Blue Schools movement occurred in 2017, when Portugal launched “Escola Azul”, a national Blue School program coordinated by the Directorate-General for Maritime Policy (DGPM) under the Ministry of the Sea. This initiative showcased how OL could be embedded across school curricula at scale. Built on partnerships with universities, municipalities, NGOs, and the private sector, Escola Azul created a network of stakeholders supporting teachers and students in ocean-themed projects [4]. It also proactively addressed systemic challenges that often hinder innovation in schools, such as rigid curricula and limited teacher training on new topics. Through an interdisciplinary and intercultural model, Escola Azul has contributed to improvements in OL, sustainability education, scientific and civic engagement among students [5].

Inspired in part by the Portuguese experience, several other Blue School networks and programs have emerged in recent years globally, such as the All-Atlantic Blue Schools Network[6], and more recently the UNESCO Global Blue Schools Network[7].

1.1. The European Blue Schools Initiative

In Europe, the Blue Schools initiative was officially launched on June 8, 2020, during World Oceans Day, as part of the EU4Ocean Coalition [8]. Supported by the European Commission, this coalition brings together educators, organizations, youth groups, policymakers, and individuals to deepen the understanding of the ocean’s critical role in sustaining life on Earth. Schools participating in this network engage students in hands-on projects that connect them to the ocean, fostering sustainability and marine awareness. Schools that successfully integrate OL into their school activities receive the European Blue School Certification, recognizing their commitment to marine education.

To obtain EU Blue School Certification, schools must meet five key requirements [9]: 1. Schools must develop a well-defined, tangible outcome related to OL, such as an educational campaign, research project, awareness initiative, or exhibition showcasing students’ learning experiences. 2. They are required to implement structured projects that interconnect various marine and maritime topics, ensuring that students engage with the ocean in a comprehensive manner. 3. EU Blue School projects must also be inclusive, involving all students across different levels and subject areas to make OL accessible to the entire school community. 4. Additionally, collaboration with external organizations, including marine research institutions, environmental NGOs, or maritime businesses, is essential to providing students with real-world insights. 5. Finally, schools must effectively communicate project results, sharing their findings with the local community and beyond to foster broader engagement and awareness in OL efforts.

Figure 1. EU Blue Schools criteria. Source: Network of EU Blue Schools.

Figure 1. EU Blue Schools criteria. Source: Network of EU Blue Schools.

1.2. The Case of Mallorca

While the European Blue Schools framework presents a promising approach, its adoption remained limited, with only 40 schools [10] applying for certification by the end of 2021 under a top-down strategy. This prompted the authors to explore whether bottom-up approaches might be more effective for educational institutions, given their daily operational commitments. To test this hypothesis, a grassroots strategy was implemented in a localized setting, with the island Mallorca chosen as the pilot case. With support from the Erasmus+ program, the authors, along with additional two partners [11], initiated a pilot project aiming to certify 100 schools, equivalent to one in three schools on the island. This strategy was structured into key phases: conducting background research on schools through websites and social media, engaging in one-on-one outreach, facilitating school visits, developing five model Blue School projects as multipliers, supporting teachers throughout the certification process, and organizing certification events to recognize achievements.

Figure 2. We are Blue Schools Project implementation approach.

Figure 2. We are Blue Schools Project implementation approach.

The project achieved notable success, receiving formal recognition from the European Commission’s Maritime Forum [12] and being awarded “Best Project” by the WestMED Initiative [13] for its innovative and inclusive approach to fostering OL at the local level. As a result, Mallorca now accounts for approximately one in five certified European Blue Schools. However, despite the success and a consistent and equal outreach strategy applied across the island, a substantial number of schools chose not to pursue certification. This uneven uptake highlights an important research question: “What territorial and institutional factors influence a school’s decision to engage with the Blue Schools program?”

Figure 3. 100 Certified EU Blue Schools (in blue) and non-certified schools (in red).

Figure 3. 100 Certified EU Blue Schools (in blue) and non-certified schools (in red).

1.3. Exploring Factors for Blue School Certification: Towards a Place-Based Understanding

Limited Research on Territorial Factors

Despite the growing visibility of the Blue Schools movement, the territorial factors influencing certification remain insufficiently studied. While the initiative promotes inclusive OL across diverse contexts, from coastal to inland schools [9], existing research has largely focused on school-level dynamics, such as teacher motivation [14,15], student [16] and teachers OL knowledge level, barriers to teach OL in schools [17], or curricular integration [18,19].

A case study from the Azores highlights Blue School participation but primarily focuses on teachers’ awareness and the limited classroom integration of ocean literacy [20]. Another study examining ocean literacy levels among students in three Mediterrananean countries, notes the presence of Blue School projects, yet the emphasis remains on student knowledge gaps rather than spatial determinants of school participation [21].

The ProBleu methodology supports school engagement through citizen science and peer-based mentoring strategies, promoting inclusivity and capacity-building within the network [22]. Yet these internal factors do not operate in isolation. The territorial context in which a school is embedded, its proximity to the sea, access to marine institutions, infrastructural connectivity, or rural isolation, can shape both the feasibility and form of marine education projects.

A Hypothesis Supported by Programmatic Insights

Although the academic literature on territorial drivers of Blue School certification is still emerging, preliminary evidence from national networks suggests that territorial context does, in fact, play a significant role. For instance, in the Portuguese Blue School Network [5], schools operate across a variety of settings including rural areas, fishing villages, urban centers, and islands, each with distinct conditions that shape project implementation. While island schools benefit from direct access to marine environments and urban schools often enjoy better infrastructure and partnerships with science centers, rural schools face logistical challenges despite strong community cohesion. These variations underscore the importance of context-sensitive approaches in advancing OL. In the broader Mediterranean region, frameworks like the Mediterranean Sea Literacy (MSL) initiative [23] also emphasize aligning educational content with local ecological and cultural realities.

Beyond the scope of Blue Schools specifically, previous research has highlighted the influence of human geography variables on OL more broadly. Research has shown that younger students tend to be more receptive to OL initiatives, while female learners consistently demonstrate higher levels of environmental concern and a greater willingness to engage in marine conservation efforts, however male show more levels of knowledge [24]. Among adults, factors such as level of education, income, and residential proximity to the coast have been positively associated with both ocean knowledge and concern for marine issues [25]. Furthermore, evidence from the UK–Ireland context shows that coastal residents are significantly more likely than inland residents to engage in ocean-related safety behaviors such as tide-checking, illustrating how spatial proximity to the sea influences practical OL [26].

1.4. Modeling Spatial Potential for Blue School Certification: A MaxEnt Approach

Understanding how educational innovations like Blue School certification expand across space calls for modeling approaches that can work with the kind of data we often have in the real world, cases where we only know where the program exists. Maximum Entropy (MaxEnt) modeling offers a useful solution in this context. It’s a presence-only method that estimates the most likely distribution of a phenomenon based on the conditions observed at known presence locations, without requiring data on absences [27].

While MaxEnt was originally developed for ecological applications, particularly species distribution modeling [27], it has become widely used across disciplines thanks to some practical advantages. It performs well even with small datasets, can handle correlations among predictor variables, and models nonlinear relationships effectively [28]. These strengths make it appealing for use in education and social research, where clean, large, and balanced datasets are often unavailable.

Over the past decade, MaxEnt has been applied to an increasingly broad range of human and environmental challenges. In environmental sciences, for example, it has been used to map flood-prone areas [29], forecast fire risk [30,31,32], and identify hotspots of human–wildlife conflict [33]. In each case, presence data, whether flood events, fire occurrences, or conflict incidents, are modeled against a set of environmental and social predictors to estimate future risk or likelihood.

More recently, researchers have extended the application of MaxEnt to model subjective or behavior-driven phenomena. In the context of land-use planning, MaxEnt has been used to visualize the spatial distribution of cultural ecosystem services by integrating participatory mapping with biophysical and landscape variables [34]. Similarly, other studies have employed MaxEnt to predict tourist concentrations by coupling the model with infrastructure, land use data, and geotagged social media check-ins[35].

Building on this interdisciplinary potential, we apply MaxEnt to model the spatial distribution of Blue School certification in Mallorca. Using the locations of certified schools as presence data, and drawing on a range of socio-economic, environmental, and infrastructural variables, we estimate the probability of certification across the island. This helps us not only understand the factors that influence certification, but also identify areas with potential for future expansion, contributing both to methodological development and more informed educational policy.

1.5. Research Aims and Hypotheses

Based on the conceptual framework of the Blue Schools and the spatial dynamics of educational participation, this study explores the territorial conditions that enable or constrain Blue School certification. The research focuses on three interrelated dimensions:

(1) the spatial distribution of certified schools in relation to environmental, socio-economic, and infrastructural predictors;

(2) the explanatory power of Maxent modeling in identifying high-suitability areas for certification using presence-only data; and

(3) the potential for using spatial prediction to guide strategic expansion of the Blue Schools programs.

From these objectives, the following hypotheses are proposed:

H1.

The spatial distribution of Blue School certifications in Mallorca is not explained solely by coastal proximity, but by a wider set of territorial variables including transport accessibility, presence of marine institutions, and economic structure.

H2.

Maxent modeling can reliably predict the suitability of school locations for Blue School certification using presence-only data, achieving high classification accuracy through a combination of geospatial predictors.

H3.

Spatial modeling will reveal non-certified schools located in areas with high territorial suitability, particularly in inland or underserved municipalities, challenging assumptions that Blue School certification is primarily a coastal phenomenon.

2. Materials and Methods

Mallorca, the largest island in Spain’s Balearic archipelago, features a diverse geography: the Serra de Tramuntana mountains in the northwest, the fertile Es Pla plain in the centre, and a heavily urbanized southern coastline anchored by Palma. As of the first quarter of 2025, continuous census data indicate Mallorca's population stood at approximately 967,322 inhabitants, marking a 0.7 % growth compared to the previous year [36,37]. The island encompasses 53 municipalities, with 27 coastal municipalities, providing varied environments for marine education and experiences.

Figure 4. Methodology.

Figure 4. Methodology.

2.1. Study Area

Although tourism remains the dominant economic sector [38], its influence is equally evident in the island’s training and education systems [39] which continue to prioritize the needs of the tourism industry [40,41]. However, regional policy documents for 2023–2030 outline a strategic shift toward Blue Economy development, emphasizing marine conservation, sustainable fisheries, renewable ocean energy, and maritime research [42]. This shift is reinforced by Mallorca’s expanding maritime infrastructure, which includes over 40 ports and marinas [43,44], alongside a growing ecosystem of marine-focused NGOs, maritime museums, research organisations, and educational institutions. Together, these developments create favorable conditions for advancing OL initiatives across the island.

Figure 5. Study area: Mallorca.

Figure 5. Study area: Mallorca.

2.2. Data and Variables

In order to identify which geographical factors influence Blue School certification and should therefore be included in the MaxEnt modeling framework, all certified projects from the We Are Blue Schools database were reviewed. The analysis focused on extracting keywords that reflect spatially relevant themes. These were then compared with factors previously highlighted in academic literature, such as distance to the coast or income [19] and the urban–rural divide [5] to ensure both empirical grounding and theoretical alignment.

From this keyword analysis, several territorial themes emerged recurrently, including Beach clean up areas (43 mentions), Posidonia oceanica meadows (22 mentions), marine biodiversity and coastal ecosystems (12), natural parks and marine protected areas (8), freshwater systems and waterways (7), and contrasts between virgin and tourist zones (5). These project-level insights supported the design of a spatial modeling framework grounded in real educational content.

Table 1. Keywords extracted from the 100 projects of the We are Blue Schools project database.

#	Times mentioned	Keyword (s)
1	43	Beach cleanup (1)
2	26	Environmental conservation (1)
3	22	Posidonia oceanica meadows (1)
4	18	Waste management and recycling (1)
5	15	plastic pollution & reduction (1)
6	14	Microplastics (1)
7	12	Coastal ecosystem knowledge & protection; Marine biodiversity (2)
8	9	Sailing and water sports (1)
9	8	Marine ecosystem; natural park (2)
10	7	Marine environment; Water and waterways; dunes; learning camp (4)
11	6	Art and upcycling solutions; Human actions impact in marine environment (4)
12	5	Virgin areas vs. tourist areas; Mediterranean Sea (5)

Table 1. Keywords extracted from the 100 projects of the We are Blue Schools project database.

#	Times mentioned	Keyword (s)
1	43	Beach cleanup (1)
2	26	Environmental conservation (1)
3	22	Posidonia oceanica meadows (1)
4	18	Waste management and recycling (1)
5	15	plastic pollution & reduction (1)
6	14	Microplastics (1)
7	12	Coastal ecosystem knowledge & protection; Marine biodiversity (2)
8	9	Sailing and water sports (1)
9	8	Marine ecosystem; natural park (2)
10	7	Marine environment; Water and waterways; dunes; learning camp (4)
11	6	Art and upcycling solutions; Human actions impact in marine environment (4)
12	5	Virgin areas vs. tourist areas; Mediterranean Sea (5)

A total of 16 spatial variables were compiled and grouped into environmental, socio-economic, and infrastructural categories. All layers were converted into raster format at a spatial resolution of 1 km² and harmonized to a common projection (EPSG:25831) for compatibility with Maxent modeling. The selected variables were as follows:

Table 2. Variables selected and data sources.

Variable	Description	Data Source
Distance to coast	Straight-line distance to nearest coastline	Own calculation (GIS)
Elevation	Elevation above sea level (EU-DEM)	Copernicus EU-DEM
Distance to Forest	Proximity to forested areas	CORINE Land Cover
Natura 2000 areas	Distance to Natura 2000 areas	IDEIB
Hidrological network	Length of water streams per grid cell	IDEIB
Freshwater vulnerability	Index of aquifer vulnerability	IDEIB
Population working in Primary sector	% of employed population in agriculture/fisheries	INE
Population working in Secondary sector	% in industry/construction	INE
Population working in Terciary sector	% in services	INE
Young population (under 16yrs)	% of population under 19	INE
GDP per capita	Average income per person	INE
Time spent to reach the coast by public transport	Time to reach coast using public transport in min	Transport Illes Balears
Time spent to reach the coast by private transport	Time to reach coast using private transport in min	Google Maps
Distance to Blue Economy organizations	Distance to nearest Blue Economy actor	PortsIB, APEAM, CAIB
Distance to cleanup sites	Distance to known cleanup activity locations	Let's Clean Up Europe, MARNOBA
Distance to rural	Distance to rural areas	CORINE Land Cover

Table 2. Variables selected and data sources.

Variable	Description	Data Source
Distance to coast	Straight-line distance to nearest coastline	Own calculation (GIS)
Elevation	Elevation above sea level (EU-DEM)	Copernicus EU-DEM
Distance to Forest	Proximity to forested areas	CORINE Land Cover
Natura 2000 areas	Distance to Natura 2000 areas	IDEIB
Hidrological network	Length of water streams per grid cell	IDEIB
Freshwater vulnerability	Index of aquifer vulnerability	IDEIB
Population working in Primary sector	% of employed population in agriculture/fisheries	INE
Population working in Secondary sector	% in industry/construction	INE
Population working in Terciary sector	% in services	INE
Young population (under 16yrs)	% of population under 19	INE
GDP per capita	Average income per person	INE
Time spent to reach the coast by public transport	Time to reach coast using public transport in min	Transport Illes Balears
Time spent to reach the coast by private transport	Time to reach coast using private transport in min	Google Maps
Distance to Blue Economy organizations	Distance to nearest Blue Economy actor	PortsIB, APEAM, CAIB
Distance to cleanup sites	Distance to known cleanup activity locations	Let's Clean Up Europe, MARNOBA
Distance to rural	Distance to rural areas	CORINE Land Cover

The spatial and socio-economic variables used in the Maxent model were derived from a combination of publicly available datasets and original geospatial analyses. Proximity to Blue Economy organizations was calculated through custom mapping using data provided by Ports de les Illes Balears [44] and the Balearic Association of Marine Enterprises [45]. Travel time to the coast via public and private transport was modeled using data from Transport Illes Balears [46], while straight-line distance to the coast was calculated through original GIS processing.

Cleanup site locations were compiled from the Let's Clean Up Europe initiative [47] and the MARNOBA citizen science platform [48]. Demographic and economic indicators, including employment by sector, youth population (under 19 years), and GDP per capita were obtained from the Spanish National Statistics Institute [49].

Land classification data for identifying rural and urban areas, as well as forest coverage, were sourced from the CORINE Land Cover inventory [50]. Natura 2000 protected areas, the hydrological network, and freshwater vulnerability indicators were acquired from the Balearic Spatial Data Infrastructure [51]. Elevation data were derived from the Copernicus Digital Elevation Model [52].

2.3. The Method: Maxent Modeling

We employed Maxent version 3.4.4 to model the spatial distribution of Blue School certifications using presence-only data. The model was trained on a 75% random sample of the data and tested on the remaining 25%, with five cross-validation replicates to ensure stability and reduce overfitting. The logistic output format was used to produce continuous suitability values ranging from 0 to 1, representing the estimated relative probability of Blue School certification under the constraints imposed by the input variables.

Maxent operates by estimating the probability distribution of maximum entropy (i.e., closest to uniform), constrained by the empirical average of each environmental variable across presence locations. Given a set of features fj(x) the algorithm seeks the probability distribution P(x) that maximizes entropy:

where:

• P(x) is the estimated suitability at location x,
• fj(x) are the environmental or socio-economic features,
• λj are the weights assigned to each feature j, learned during model training
• and Z is a normalizing constant ensuring that the total probability sums to one.

This formulation ensures that the expected value of each feature under the estimated distribution matches its empirical average over the observed presences, while otherwise assuming as little as possible about unknown data.

Model performance was evaluated using the Area Under the Curve (AUC) of the Receiver Operating Characteristic (ROC), alongside omission rate plots to assess calibration. To determine the relative contribution of each variable, we analyzed percent contribution (based on model training), permutation importance (based on reduction in model performance when variables are randomized), and jackknife tests (measuring gain from each variable alone or excluded). This multi-metric approach provides a robust assessment of which spatial factors most influence the likelihood of Blue School certification.

2.4. School-Level Suitability Extraction

Following the generation of the median prediction surface, we extracted the predicted suitability value at each known school location in Mallorca. This allowed us to assign a spatial probability score, ranging from 0 to 1, to both certified and non-certified schools across the island. By juxtaposing the full population of schools with the prediction map, we were able to identify high-suitability schools that have not yet been certified. This spatial overlay serves as a practical tool for prioritizing outreach and capacity-building efforts, particularly in inland or underserved areas that exhibit strong alignment with the territorial drivers of certification identified by the model.

3. Results

3.1. Model Performance

The MaxEnt model exhibited strong predictive performance, achieving an average test AUC (Area Under the Curve) of 0.840 across five cross-validated replicates, with a standard deviation of ±0.044 (Figure 6). This high AUC value indicates excellent discriminatory power in differentiating between locations with and without Blue School certification. In practical terms, it reflects a well-calibrated balance between sensitivity (true positive rate) and specificity (true negative rate), confirming the model’s robustness in identifying areas of high suitability based on the selected geospatial predictors.

Additional evidence of strong model calibration is provided by the omission versus predicted area curve (Figure 5). The observed omission rate (cyan line) closely tracks the predicted omission rate (black line) across cumulative thresholds, remaining largely within the orange confidence envelope (± one standard deviation). This alignment suggests that the predicted probabilities correspond well with the actual spatial distribution of certified schools, reinforcing the model’s internal consistency.

Figure 7. Average Omission and Predicted Area for MaxEnt Model.

Figure 7. Average Omission and Predicted Area for MaxEnt Model.

Figure 8. Receiver Operating Characteristic (ROC) Curve and AUC for MaxEnt Model.

Figure 8. Receiver Operating Characteristic (ROC) Curve and AUC for MaxEnt Model.

3.2. Variable Importance and Spatial Drivers

Among the predictor variables, distance to Blue Economy organisations emerged as the most influential according to percent contribution (30.2%), and it yielded the highest training gain when evaluated independently in the jackknife analysis. This underscores the pivotal role of proximity to maritime enterprises, research institutions, and NGOs in facilitating school engagement with OL. However, its low permutation importance (4.4%) suggests that its predictive influence may be shared with correlated variables, reducing its unique contribution when others are present.

In contrast, the hydrological network surfaced as the most informative and non-redundant variable. Its omission resulted in the greatest drop in model performance, with a permutation importance of 23.5%, indicating that freshwater systems, such as torrents (dry rivers), reservoirs, and wetlands are critical territorial enablers of aquatic environmental education.

Table 3. Importance and contribution of the variables to the model.

Variable	Percent contribution	Permutation importance
Distance to blue economy organisations	30.2	4.4
Hidrological network	18.4	23.5
Population working in Primary sector	12.9	21.2
Time spent to reach the coast by public transport	11.9	25.6
Natura 2000 areas	10	7.9
Distance to Forest	6.9	1.8
Distance to Rural	2.7	0.8
GDP per capita	1.9	3
Distance to beach and river cleanups	1.5	0.2
Elevation	1.1	2.8
Population working in Terciary sector	0.9	2.3
Freshwater vulnerability	0.6	0.7
Population working in Secondary sector	0.5	2.1
Distance to coast	0.4	0.9
Time spent to reach the coast by private transport	0.2	2
Young population (under 16yrs)	0	0.7

Table 3. Importance and contribution of the variables to the model.

Variable	Percent contribution	Permutation importance
Distance to blue economy organisations	30.2	4.4
Hidrological network	18.4	23.5
Population working in Primary sector	12.9	21.2
Time spent to reach the coast by public transport	11.9	25.6
Natura 2000 areas	10	7.9
Distance to Forest	6.9	1.8
Distance to Rural	2.7	0.8
GDP per capita	1.9	3
Distance to beach and river cleanups	1.5	0.2
Elevation	1.1	2.8
Population working in Terciary sector	0.9	2.3
Freshwater vulnerability	0.6	0.7
Population working in Secondary sector	0.5	2.1
Distance to coast	0.4	0.9
Time spent to reach the coast by private transport	0.2	2
Young population (under 16yrs)	0	0.7

Other high-ranking variables include the percentage of the population employed in the primary sector (21.2%) and travel time to the coast by public transport (25.6%), both of which point to contextual factors that facilitate marine engagement. The former reflects the relevance of traditional economies, such as agriculture and fisheries, in cultivating environmental awareness and stewardship. The latter highlights the enabling role of sustainable accessibility, schools with efficient public transport connections to coastal zones (e.g. train and direct public lines) are better positioned to carry out field-based marine learning activities.

Moderate contributions were also observed for proximity to Natura 2000 areas (10%) and distance to forested zones (6.9%), suggesting that natural protected areas and transitional ecosystems support place-based learning and curricular innovation. These landscapes likely serve as focal points for school projects that align with biodiversity and conservation goals.

By contrast, socio-economic indicators such as GDP per capita (3%), employment in the tertiary sector (2.3%), and youth population share (0.7%) played marginal roles in the model. This suggests that structural and geographic accessibility to marine or freshwater environments is more decisive than general economic wealth or demographic composition in driving school participation.

Together, these results indicate that the spatial distribution of Blue School certification is not solely determined by proximity to the coast, but rather by a complex interplay of access to marine-related institutions, engagement with water systems, and community-based drivers rooted in economic activity and environmental values.

Figure 9. MaxEnt response curves of the studied variables.

Figure 9. MaxEnt response curves of the studied variables.

3.3. Suitability Mapping and School-Level Predictions

The median prediction surface generated by the MaxEnt model reveals distinct spatial patterns in the likelihood of Blue School certification across the island of Mallorca. Zones of high suitability are clearly concentrated along the southeastern and eastern coastal corridors, as well as within the Palma metropolitan area. These areas exhibit a strong convergence of favorable territorial conditions, most notably, close proximity to Blue Economy organizations, efficient public transport access to the coast, and dense hydrological networks, which were among the highest-ranking variables in the model. Their spatial overlap highlights the synergistic effect of infrastructure, ecosystem connectivity, and institutional presence in fostering OL engagement.

In particular, the Palma region stands out not only due to urban density but also because of its clustering of marine research institutions, NGOs, port authorities, and educational hubs, which together create an ecosystem of support for schools seeking certification. Several schools within Palma’s urban fabric scored above 0.94 in predicted suitability, despite not yet being certified. This suggests a significant latent potential in urban environments where functional and institutional access compensates for distance to the physical coastline.

Figure 10. New Blue Schools suitability map.

Figure 10. New Blue Schools suitability map.

Similarly, the eastern and southeastern coastal belts align with productive landscapes and Natura 2000 protected areas, where environmental awareness may be higher and partnerships with conservation actors more feasible. The predictive contribution of the primary sector employment variable (notably high in rural coastal areas) may explain the engagement of schools located in small fishing villages and agricultural zones, where communities are more attuned to marine and aquatic resource sustainability.

In contrast, areas of low suitability are predominantly concentrated in the central interior (Pla de Mallorca) and the northwestern mountainous zones (Serra de Tramuntana). These regions are typically more rural, dispersed, and infrastructurally isolated, with fewer marine institutions, weaker transport connectivity, and lower exposure to both marine and freshwater ecosystems. While some of these areas host schools with strong environmental traditions, the absence of structural enablers, such as public transport links to the coast or nearby marine stakeholders, limits their baseline potential for certification under current conditions. This highlights how territorial disadvantage may stem not from distance per se, but from multi-dimensional access barriers.

Moreover, the model’s surface reveals nuanced cases where inland schools score relatively high in suitability, particularly where freshwater systems intersect with educational centers, access to public transport, and Natura 2000 areas. These "outliers" reinforce the importance of freshwater engagement and protected area proximity as alternative pathways for schools not located near the sea to participate in the Blue School program.

As a next step, the prediction surface was overlaid with the complete set of school locations in Mallorca, enabling the extraction of site-specific suitability scores for both certified and non-certified institutions. This spatial overlay revealed that more than 30 non-certified schools exhibit predicted suitability values exceeding 0.85, indicating a high latent potential for engagement with the Blue School programme. Notably, several urban schools in Palma emerged as particularly promising candidates, with predicted scores above 0.94. These include IES Joan Alcover, IES Ramon Llull, IES Politècnic, and CEIP Aina Moll i Marquès.

Table 4. Schools with higher probability to become Blue.

name	Value	name	Value
IES Joan Alcover	0.943	EOI Andreu Crespí Plaza	0.926
IES Ramon Llull	0.943	Escola Jafudà Cresques	0.926
IES Politècnic	0.943	Col·legi de les Escolàpies	0.921
CEIP Aina Moll i Marquès	0.943	CEIP Son Pisà	0.921
Col·legi Lluís Vives	0.943	escola Simó Ballester	0.907
Escola S. de Disseny IB	0.943	Col·legi de la Puresa de Maria	0.907
Escola Joan Capó	0.936	col·legi de Sant Vicenç de Paül	0.907
Col·legi de l'Arcàngel Sant Rafel	0.926	escola de Sa Torre	0.907

Table 4. Schools with higher probability to become Blue.

name	Value	name	Value
IES Joan Alcover	0.943	EOI Andreu Crespí Plaza	0.926
IES Ramon Llull	0.943	Escola Jafudà Cresques	0.926
IES Politècnic	0.943	Col·legi de les Escolàpies	0.921
CEIP Aina Moll i Marquès	0.943	CEIP Son Pisà	0.921
Col·legi Lluís Vives	0.943	escola Simó Ballester	0.907
Escola S. de Disseny IB	0.943	Col·legi de la Puresa de Maria	0.907
Escola Joan Capó	0.936	col·legi de Sant Vicenç de Paül	0.907
Col·legi de l'Arcàngel Sant Rafel	0.926	escola de Sa Torre	0.907

Their elevated suitability values, despite their current uncertified status, highlight clear opportunities for targeted outreach and capacity-building, especially in metropolitan areas where key enabling factors such as institutional proximity, infrastructure, and access to Blue Economy actors are already in place. These findings provide actionable insight into how the Blue School program can expand strategically. Rather than focusing solely on coastal proximity, the suitability map underscores the importance of multi-dimensional territorial factors, ranging from ecological networks to transport infrastructure and economic structure, in enabling educational engagement with the marine environment. In practical terms, the results support a data-driven approach to scaling the Blue Schools program by identifying and prioritizing schools already situated within favorable territorial conditions, yet not yet part of the certification network.

Figure 11. Probability map of new schools becoming "Blue".

Figure 11. Probability map of new schools becoming "Blue".

4. Discussion

4.1. Spatial Context and Blue School Engagement

This study’s findings underscore that territorial context is pivotal in shaping Blue School certification outcomes. The MaxEnt model revealed that factors like spatial accessibility to Blue Economy institutions, the extent of hydrological networks, and availability of sustainable transport infrastructure are among the most influential predictors of school participation. These results challenge any assumption that coastal proximity or economic prosperity alone drive OL engagement. Instead, it is the combination of functional access, ecosystem connectivity, and community relevance that enables schools to develop meaningful marine education projects. This aligns with the Blue Schools philosophy that successful OL initiatives must connect learning to local communities and real-world issues [2,4,5].

In fact, OL frameworks are increasingly being adapted to fit regional and local contexts, for example, the Mediterranean Sea Literacy principles were developed to tailor OL to the specific ecological and cultural setting of the Mediterranean region [23]. Such efforts reinforce that place-based context and relevance are key to engaging students, whether their schools are on the coast or farther inland.

4.2. Inland Schools and Hydrological Connectivity

An intriguing insight is the high importance of inland hydrological networks (rivers, wetlands, etc.) in the model. This reinforces that OL is not confined to strictly coastal settings. Schools situated away from the coast but embedded in rich freshwater ecosystems or longstanding place-based environmental education traditions may be equally well-positioned to join the Blue Schools network. Indeed, the Blue Schools initiative itself emphasizes inclusivity across diverse geographies, aiming to reach “children everywhere” as ocean stewards [9].

By analogy, high potential for Blue School certification can emerge in areas where multiple enabling conditions, such as a strong environmental context, community partnerships, and access to resources, converge, even if those areas are not located directly along the coast. Inland schools with robust water-related environmental programs and active community engagement can thus become vibrant nodes of OL, extending the Blue Schools movement beyond the shoreline. For example, the NOAA-supported Ocean Acidification Education Partnership has demonstrated that landlocked schools can effectively engage in OL through watershed-based approaches and inter-school collaborations with coastal institutions [53]. Similarly, ocean education initiatives in Canada emphasize the pressing need to improve marine education for the millions of citizens living inland, showing that physical proximity to the ocean is not a prerequisite for meaningful engagement [54]. Inland schools with robust water-related environmental programs or community engagement can thus become vibrant nodes of OL, extending the Blue Schools movement beyond the shoreline.

4.3. Identifying Opportunity Zones for Expansion

One of the most impactful outcomes of this spatial analysis is the identification of 33 non-certified schools with very high predicted suitability scores (>0.85) for the Blue School program. Notably, many of these schools are clustered in and around the city of Palma, the densely populated capital of Mallorca. This confirms that urban density and institutional proximity can provide fertile ground for Blue School initiatives.

Urban centers like Palma host numerous educational institutions and nearby marine or science organizations (museums, aquariums, universities, NGOs), creating a concentration of resources and potential partners to support Blue School projects. These results offer a clear roadmap for strategically expanding the program through targeted outreach. Schools that show this kind of territorial readiness, where the geographic, social, and infrastructural conditions are ripe, can be prioritized for certification support.

From a policy standpoint, this approach enables more strategic resource allocation by identifying “opportunity zones” with the highest potential for successful Blue School certification. The predictive map can serve as a practical tool for guiding mentorship, capacity-building, and funding efforts toward schools where OL initiatives are most likely to flourish, thereby enhancing overall impact.

This territorial analysis adds a valuable layer to existing engagement models like ProBleu [22], which promote inclusivity and peer learning among schools. Incorporating geographic suitability into such initiatives could help tailor outreach and support efforts to those schools where the conditions for ocean literacy projects are already favorable.

4.4. Extending MaxEnt to Educational Context

Beyond the substantive findings on Blue Schools identification, this research showcases an innovative application of the MaxEnt modeling technique outside of its traditional ecological realm. MaxEnt was originally developed for species distribution modeling [27], but here we adapted it to presence-only data on educational program participation. By doing so, we open new possibilities for predicting where educational or social innovations are most likely to succeed based on spatial context.

This study contributes to the expanding use of MaxEnt in modeling complex social phenomena by demonstrating its relevance in the education sector. Specifically, it shows how MaxEnt can effectively identify spatial opportunities to strengthen OL in schools. By applying a tool traditionally used in ecological modeling, this work highlights the method's potential for social and educational planning, particularly in contexts where conventional data, such as indicators of non-engagement, are limited or unavailable.

4.5. Transferability and Policy Implications

While our case study focused on the island of Mallorca, the approach we developed is highly transferable to other regions. The predictive modeling framework can be readily replicated by education authorities or OL networks in other parts of Spain, through the Spanish Network of Blue Schools [55], and Europe. By adjusting the input layers to local data, such as each region’s Blue Economy infrastructure, hydrological features, and community indices, stakeholders can identify under-engaged but spatially promising schools in their own territories. In this sense, MaxEnt modelg can serve as both an analytical and strategic tool for scaling Blue Schools equitably across diverse geographies.

Our findings build on this momentum by offering a data-driven method to guide where the next waves of Blue Schools could emerge. Policymakers and program coordinators can use these insights to proactively target outreach, funding, and partnerships in areas with high latent potential for OL education. Ultimately, aligning expansion efforts with the geography of opportunity will help ensure that the growth of the Blue School network is both efficient and inclusive, bringing the benefits of OL to coastal and inland communities alike, and empowering the next generation of ocean-literate citizens across Europe

5. Conclusions

This study demonstrates the value of spatial modeling in guiding the equitable expansion of the Blue Schools network. By applying MaxEnt to presence-only data, we identified key territorial factors that shape school engagement with OL.

Importantly, the model revealed high potential in several non-certified urban and inland schools, challenging the assumption that proximity to the coast is a prerequisite for participation. These insights offer a practical foundation for strategic outreach and policy planning, ensuring that OL can reach diverse educational contexts.

Beyond its specific findings, this research highlights the broader applicability of MaxEnt for educational planning in data-scarce environments. As EU initiatives continue to promote OL, tools like this can support targeted investment and inclusive program growth across regions.

Future research could refine and extend this approach by integrating dynamic or qualitative variables, such as institutional readiness, teacher training capacity, or levels of student environmental interest, or by applying the model longitudinally to evaluate impact over time. Additionally, replicating this analysis in other European regions would help validate the model’s generalizability and support transnational coordination across the Network of European Blue Schools.