TY - JOUR AU - Liu, Tongjing AB - Introduction In 2011, the world witnessed the Fukushima disaster, marking the most severe nuclear catastrophe since the Chernobyl disaster, leaving an indelible imprint on human history. In the initial six weeks following the disaster, the amount of radioactive materials released from Fukushima was detected to be substantial, with the release of cesium amounting to 42% of the total release in the Chernobyl disaster, and the release of xenon reaching an unprecedented level [1]. In response to this crisis, the Japanese government adopted water injection measures to slow down the melting speed of the fuel rods. However, this measure led to another crisis—large amounts of nuclear-contaminated water needed to be sealed and stored, posing not only a high cost but also a potential permanent burden (International Atomic Energy Agency, IAEA, 2021) [2]. To address this pressing issue, the Japanese government initiated a series of intense discussions and evaluations, exploring various possible solutions. Among numerous solutions, discharging wastewater into the sea was considered a relatively practical approach, albeit filled with controversy and risk. After a series of careful considerations and weighing of pros and cons, on April 13, 2021, the Japanese government officially decided to discharge 1.2 million tons of Fukushima wastewater into the sea [3], a decision that immediately attracted widespread attention worldwide. Despite the Japanese government’s repeated assurances that the discharged polluted water had been treated and the radiation values met safety standards, this action still provoked strong opposition from the international community [4], as reflected in the diverse responses captured in media representations across Japan and neighboring countries, discussed in detail by Gong et al. (2024) [5]. On August 24, 2023, the Japanese government officially commenced the discharge of nuclear wastewater from the Fukushima Daiichi Nuclear Power Station into the ocean, once again triggering extensive international attention and controversy [6]. Although the International Atomic Energy Agency (IAEA) deems Japan’s plan to discharge treated water from Fukushima into the sea to be in compliance with international safety standards [7], many countries have expressed strong opposition due to the uncertainty of the potential dangers and long-term impacts brought by such astonishing discharge. Among them, China firmly opposes Japan’s unilateral decision affecting the entire international community and condemns this action [8]. Japan’s nuclear wastewater discharge plan is essentially resolving its own nuclear safety threat by sacrificing the safety of the human marine environment [9]. In this contentious context, Ng Kwan Hoong and colleagues (2022) offer critical perspectives on the possible avenues through which the Japanese government could engage more effectively with the international community to alleviate concerns and build consensus around their wastewater management strategies [10]. As early as the end of 2018, Tokyo Electric Power Company had already admitted that about 80%, i.e., 890,000 tons, of the then stored 1.1 million tons of treated water still contained 90Sr, 60Co, 106Ru, and many other radioactive nuclides exceeding the limit values, and the purification system failed to adequately remove these radioactive nuclides [11]. With the progress of Japan’s nuclear wastewater discharge plan, the world is on the verge of experiencing the most severe nuclear pollution. The seas near Fukushima Prefecture are not only the economic source for coastal residents but also an integral part of the Pacific and the world’s oceans. The substantial amount of radioactive materials will have immeasurable impacts on marine life, the natural environment, and human health. The act of discharging has already occurred, and the rational handling of Fukushima’s nuclear wastewater is a critical issue that the discharging country and even all stakeholders should focus on. Finding effective approaches is key to resolving governance dilemmas and achieving sustainable development of the oceans. In this context, exploring the long-term relationships between the Japanese government, other countries, and the Japanese Fisheries Association is necessary. Viewing these three parties as stakeholders, the dynamic game among them encompasses the outcomes of handling nuclear contaminated water. We utilize an evolutionary game model to discern the pivotal influencing factors and conduct an in-depth study on the impact of various parameters on the evolutionary process under the condition of Japan refraining from discharging nuclear wastewater, aspiring to contribute to the conservation of the living environment for all of humanity. In Section 2 of this paper, we review literature pertinent to our study. Section 3 provides a detailed explanation of each parameter used in the model, offering a clear theoretical basis for subsequent analysis and discussion. In Section 4, we propose a set of initial assumptions aimed at more accurately constructing the model of benefits and costs for the three parties—Japan, other countries, and the Fisheries Association. Based on these assumptions, we further develop the replicator dynamics equations in the evolutionary game model to reveal the strategic changes and choices of the parties during the game process. Section 5 delves deeply into the behavioral stability of the three parties under different conditions. We not only analyze various types of asymptotic stability in detail but also explore the necessary conditions and possibilities for achieving such stability, providing important theoretical support for the refinement of the theoretical model and empirical analysis. In Section 6, we validate the feasibility and accuracy of the model through numerical simulations, revealing evolutionary trends under different stable conditions. Specifically, we focus on the scenario where Japan decides to cease ocean discharge, analyzing the impact of key parameters on the evolutionary process under this scenario. This not only provides theoretical references for the three parties to seek optimal strategies in practical problems but also helps in more comprehensively understanding how each party responds to different strategies and environmental variables, thereby affecting the final outcome of the evolutionary game. In the wake of the 2011 Fukushima disaster, significant attention has been devoted to the containment and mitigation of ongoing environmental threats posed by nuclear wastewater. While considerable literature has focused on the immediate chemical hazards and environmental impacts of such discharges, there remains a substantial gap in understanding the strategic interactions among key stakeholders involved in these critical decisions. Our research aims to fill this gap by applying evolutionary game theory to explore the strategic dynamics and decision-making processes involved in the discharge of Fukushima’s nuclear wastewater. This approach allows us to identify and analyze the factors influencing these decisions and to propose evolutionarily stable strategies that could lead to more sustainable and acceptable outcomes. By doing so, we contribute a new perspective to the discourse on nuclear disaster management, emphasizing the need for a strategic approach that balances environmental integrity with political and economic realities. In this study, we utilize key parameters in our evolutionary game theory model to understand the strategic interactions concerning Fukushima’s nuclear wastewater discharge. Each symbol represents a specific variable critical to modeling the decisions and interactions among the Japanese government (J), other countries (C), and the Japanese Fisheries Association (F). These include probabilities such as x and y which represent Japan’s likelihood to continue or cease discharge and other countries’ decisions to sanction, respectively. Costs like CDJ (discharge cost) and CSJ (storage cost) quantify the financial implications of each decision. Additionally, variables like IJ and CIF reflect the international image repercussions for Japan and the Fisheries Association. Further, ERF, TRJ, and CSC assess economic impacts from changes in revenue, tax, and seafood production costs. The inclusion of these parameters provides a comprehensive framework to simulate and analyze the outcomes of different strategic scenarios, as detailed in Table 1. Download: PPT PowerPoint slide PNG larger image TIFF original image Table 1. Description of parameters. https://doi.org/10.1371/journal.pone.0317419.t001 Literature review Chemical hazards of Japan’s nuclear wastewater discharge into the sea On March 11, 2011, a severe earthquake occurred in the northeastern sea area of Japan, triggering a tsunami and causing serious damage to the Fukushima Daiichi Nuclear Power Plant (FDNPP). The failure of the cooling system led to a temperature rise in the reactor, generating hydrogen and other gases, causing explosions and releasing radioactive gases and fragments into the atmosphere. Simultaneously, untreated cooling water was directly discharged into the northeastern sea area of Japan, making this incident one of the most severe uncontrolled inputs of artificial radionuclides into the ocean in history [12]. According to the water treatment facility monitoring reports by Tokyo Electric Power Company Holdings (TEPCO, 2012–2020) and the seawater quality monitoring reports by the Nuclear Regulation Authority (NRA, 2013–2020), the radioactive substances mainly include 3H, 14C, 134Cs, 137Cs, 60Co, 125Sb, 90Sr, 129I, 99Tc, 106Ru, and 238Pu [13]. Although the concentrations of these isotopes are lower than tritium, they are more likely to integrate into marine organisms and seabed sediments (Buesseler, 2020) [14]. If nuclear wastewater is discharged into the Pacific Ocean, potential risks will exist for hundreds or even thousands of years in the future. The marine food web tightly connects marine organisms through predation (Albouy et al., 2019) [15]. Cesium isotopes are the main pollutants in the accumulated seawater inside Japanese nuclear facilities and tanks. Additionally, cesium is a highly soluble radioactive nuclide in seawater, posing long-term radiation risks to the environment and is more likely to integrate into marine organisms or sediments under high concentration factors (Buesseler, 2020) [14]. The harm of cesium to humans is mainly reflected in its radioactive radiation and human absorption. Moreover, the concentration of 14C in treated nuclear wastewater is also high (the Korea Times, 2021b) [16], it releases low-energy β particles and may enter the biosphere and accumulate in marine ecosystems (Williams et al, 2010) [17]. Co is also a radioactive isotope that releases gamma rays, which can penetrate the human body and cause cell damage (Khajeh et al, 2017) [18]. 90Sr can simulate calcium in the human body, significantly increasing the risk of osteosarcoma and leukemia (Khani et al., 2012) [19]. The ocean has a large volume and a complex current system, with strong abilities to dilute and disperse radioactive substances, but long half-life radioactive nuclides will still exist in the marine environment for a long time (Men and Deng et al., 2017) [20], causing significant consequences to the natural environment, marine life, and human health. Evolutionary game theory Game theory was originally designed to analyze economic behavior and, with modifications, can also be applied to evolving populations. The concept of “Evolutionarily Stable Strategy” (ESS) proposed by Smith (1982) reveals the relevance of the optimal behavior of animals or plants under the influence of the behavior of other individuals [21]. In classical game theory, the basic components of a game include players, strategies, and payoff matrices, and the results of the game can be clearly displayed through the payoff matrices. Evolutionary game theory, originating from classical game theory, has gradually expanded its range of application. Smith (1973) regarded evolutionary game theory as a theoretical framework for studying the interaction between individual behavior and adaptation. According to this theory, the evolutionary process is a dynamic strategy selection process, where the fitness of an individual depends not only on the chosen strategy but also on interactions with other individuals [22]. Adami et al. (2016) believed that the main goal in the application of evolutionary game theory is to find suitable strategies or optimal decision sequences to resolve existing conflicts and obtain maximum benefits. However, the competition between different strategies does not occur simultaneously at the same point in time but changes over specific periods. New strategies continuously emerge, compete with existing ones, and may replace the originally dominant strategies by gaining advantages in space and time. Therefore, the success of a strategy depends on its competition with other strategies during a specific period [23]. Evolutionary game models are important tools for exploring social behavior. Predecessors have constructed various types of evolutionary game models through extensive research to explain the emergence and continuation of various social behaviors. It has been widely used in economics, sociology, and statistical physics to analyze various factors affecting the formation of group behavior, producing rich and influential research results (Perc and Jordan et al., 2017; Khoo and Fu et al., 2018; Wang and He et al., 2018)[24–27]. Notably, ESS (Evolutionarily Stable Strategy) is a powerful tool in game theory to find optimal strategies that can prevail in the competition with other strategies [28]. Application of evolutionary game theory in this study This global concern has spurred significant research into the containment and mitigation of ongoing environmental threats posed by such nuclear wastewater. Despite extensive discussions on the immediate chemical hazards and environmental impacts of these discharges, there remains an evident gap in the literature regarding the strategic interactions among the key stakeholders managing these critical decisions. Following Japan’s decision to discharge nuclear wastewater into the sea, researchers have begun to employ evolutionary game theory to study this issue. Zheng (2022) constructed an evolutionary game model involving fishermen, consumers, and the government, analyzing the impact of each entity’s behavior probability on the strategies of other entities and the stability of the entire system. However, this study did not consider the impact of international cooperation and conflicts on strategy evolution [29]. Xin et al. (2022) explored an evolutionary game considering emotional impacts to study the evolutionarily stable strategies of the Japanese government and Japanese fishermen. This study primarily focused on domestic factors and did not fully consider the dynamics and interactions at the international level [30]. Xu et al. (2022) established a tripartite evolutionary game model including the International Atomic Energy Agency (IAEA), the Discharging Country (DC), and the Fisheries Cooperative Association of Japan (FCA), aiming to propose management insights for international cooperation. However, this study did not fully explore the changes and evolution of strategies among the parties in the game, an aspect this paper intends to analyze more comprehensively [31]. Su et al. (2023) constructed a Tripartite Evolutionary Game (TEG) model to analyze the evolutionary stability of strategy choices among fishermen, cooperatives, and government departments, exploring influencing factors and their relationships. The results indicate that the profits of cooperatives and the cooperation costs of fishermen are key factors affecting the strategy choices of the three parties in the game. However, this study did not consider the impact of international factors on strategy choices [32]. Zhang et al. (2023) adopted a complex network evolutionary game based on the WS small-world model, focusing on the fisherman population, to explore the dynamic evolutionary laws of cooperative behavior diffusion in the development of marine carbon sink fisheries. This study did not involve the roles and strategies of other countries and international organizations, an aspect that is supplemented and expanded upon in this paper [33]. Previous research on Japan’s nuclear wastewater discharge was predominantly predictive, given that Japan only commenced the actual discharge on August 24, 2023. Consequently, studies prior to this date were often based on limited information and considerations. In contrast, our research, informed by post-discharge realities, incorporates a broader spectrum of practical factors and employs evolutionary game theory to delve into the strategic decisions and interactions among the stakeholders. Our approach amalgamates perspectives from economics, environmental science, and international relations, with a particular emphasis on real-time international impacts and responses—elements frequently overlooked in earlier studies. Notably, our work pioneers the exploration of the balance of interests among the Japanese government, other nations, and the Japan Fisheries Association. Furthermore, we are the first to apply evolutionary game theory within this specific context. Our study not only broadens the knowledge frontier of evolutionary game models in environmental policy and international relations but also, through a tripartite framework, adeptly addresses the literature gap concerning the tripartite balance of interests in Japan’s marine discharge policy. In essence, our analysis, rooted in an international dimension, offers a comprehensive and realistic perspective on the dynamic games surrounding Japan’s nuclear wastewater issue. Chemical hazards of Japan’s nuclear wastewater discharge into the sea On March 11, 2011, a severe earthquake occurred in the northeastern sea area of Japan, triggering a tsunami and causing serious damage to the Fukushima Daiichi Nuclear Power Plant (FDNPP). The failure of the cooling system led to a temperature rise in the reactor, generating hydrogen and other gases, causing explosions and releasing radioactive gases and fragments into the atmosphere. Simultaneously, untreated cooling water was directly discharged into the northeastern sea area of Japan, making this incident one of the most severe uncontrolled inputs of artificial radionuclides into the ocean in history [12]. According to the water treatment facility monitoring reports by Tokyo Electric Power Company Holdings (TEPCO, 2012–2020) and the seawater quality monitoring reports by the Nuclear Regulation Authority (NRA, 2013–2020), the radioactive substances mainly include 3H, 14C, 134Cs, 137Cs, 60Co, 125Sb, 90Sr, 129I, 99Tc, 106Ru, and 238Pu [13]. Although the concentrations of these isotopes are lower than tritium, they are more likely to integrate into marine organisms and seabed sediments (Buesseler, 2020) [14]. If nuclear wastewater is discharged into the Pacific Ocean, potential risks will exist for hundreds or even thousands of years in the future. The marine food web tightly connects marine organisms through predation (Albouy et al., 2019) [15]. Cesium isotopes are the main pollutants in the accumulated seawater inside Japanese nuclear facilities and tanks. Additionally, cesium is a highly soluble radioactive nuclide in seawater, posing long-term radiation risks to the environment and is more likely to integrate into marine organisms or sediments under high concentration factors (Buesseler, 2020) [14]. The harm of cesium to humans is mainly reflected in its radioactive radiation and human absorption. Moreover, the concentration of 14C in treated nuclear wastewater is also high (the Korea Times, 2021b) [16], it releases low-energy β particles and may enter the biosphere and accumulate in marine ecosystems (Williams et al, 2010) [17]. Co is also a radioactive isotope that releases gamma rays, which can penetrate the human body and cause cell damage (Khajeh et al, 2017) [18]. 90Sr can simulate calcium in the human body, significantly increasing the risk of osteosarcoma and leukemia (Khani et al., 2012) [19]. The ocean has a large volume and a complex current system, with strong abilities to dilute and disperse radioactive substances, but long half-life radioactive nuclides will still exist in the marine environment for a long time (Men and Deng et al., 2017) [20], causing significant consequences to the natural environment, marine life, and human health. Evolutionary game theory Game theory was originally designed to analyze economic behavior and, with modifications, can also be applied to evolving populations. The concept of “Evolutionarily Stable Strategy” (ESS) proposed by Smith (1982) reveals the relevance of the optimal behavior of animals or plants under the influence of the behavior of other individuals [21]. In classical game theory, the basic components of a game include players, strategies, and payoff matrices, and the results of the game can be clearly displayed through the payoff matrices. Evolutionary game theory, originating from classical game theory, has gradually expanded its range of application. Smith (1973) regarded evolutionary game theory as a theoretical framework for studying the interaction between individual behavior and adaptation. According to this theory, the evolutionary process is a dynamic strategy selection process, where the fitness of an individual depends not only on the chosen strategy but also on interactions with other individuals [22]. Adami et al. (2016) believed that the main goal in the application of evolutionary game theory is to find suitable strategies or optimal decision sequences to resolve existing conflicts and obtain maximum benefits. However, the competition between different strategies does not occur simultaneously at the same point in time but changes over specific periods. New strategies continuously emerge, compete with existing ones, and may replace the originally dominant strategies by gaining advantages in space and time. Therefore, the success of a strategy depends on its competition with other strategies during a specific period [23]. Evolutionary game models are important tools for exploring social behavior. Predecessors have constructed various types of evolutionary game models through extensive research to explain the emergence and continuation of various social behaviors. It has been widely used in economics, sociology, and statistical physics to analyze various factors affecting the formation of group behavior, producing rich and influential research results (Perc and Jordan et al., 2017; Khoo and Fu et al., 2018; Wang and He et al., 2018)[24–27]. Notably, ESS (Evolutionarily Stable Strategy) is a powerful tool in game theory to find optimal strategies that can prevail in the competition with other strategies [28]. Application of evolutionary game theory in this study This global concern has spurred significant research into the containment and mitigation of ongoing environmental threats posed by such nuclear wastewater. Despite extensive discussions on the immediate chemical hazards and environmental impacts of these discharges, there remains an evident gap in the literature regarding the strategic interactions among the key stakeholders managing these critical decisions. Following Japan’s decision to discharge nuclear wastewater into the sea, researchers have begun to employ evolutionary game theory to study this issue. Zheng (2022) constructed an evolutionary game model involving fishermen, consumers, and the government, analyzing the impact of each entity’s behavior probability on the strategies of other entities and the stability of the entire system. However, this study did not consider the impact of international cooperation and conflicts on strategy evolution [29]. Xin et al. (2022) explored an evolutionary game considering emotional impacts to study the evolutionarily stable strategies of the Japanese government and Japanese fishermen. This study primarily focused on domestic factors and did not fully consider the dynamics and interactions at the international level [30]. Xu et al. (2022) established a tripartite evolutionary game model including the International Atomic Energy Agency (IAEA), the Discharging Country (DC), and the Fisheries Cooperative Association of Japan (FCA), aiming to propose management insights for international cooperation. However, this study did not fully explore the changes and evolution of strategies among the parties in the game, an aspect this paper intends to analyze more comprehensively [31]. Su et al. (2023) constructed a Tripartite Evolutionary Game (TEG) model to analyze the evolutionary stability of strategy choices among fishermen, cooperatives, and government departments, exploring influencing factors and their relationships. The results indicate that the profits of cooperatives and the cooperation costs of fishermen are key factors affecting the strategy choices of the three parties in the game. However, this study did not consider the impact of international factors on strategy choices [32]. Zhang et al. (2023) adopted a complex network evolutionary game based on the WS small-world model, focusing on the fisherman population, to explore the dynamic evolutionary laws of cooperative behavior diffusion in the development of marine carbon sink fisheries. This study did not involve the roles and strategies of other countries and international organizations, an aspect that is supplemented and expanded upon in this paper [33]. Previous research on Japan’s nuclear wastewater discharge was predominantly predictive, given that Japan only commenced the actual discharge on August 24, 2023. Consequently, studies prior to this date were often based on limited information and considerations. In contrast, our research, informed by post-discharge realities, incorporates a broader spectrum of practical factors and employs evolutionary game theory to delve into the strategic decisions and interactions among the stakeholders. Our approach amalgamates perspectives from economics, environmental science, and international relations, with a particular emphasis on real-time international impacts and responses—elements frequently overlooked in earlier studies. Notably, our work pioneers the exploration of the balance of interests among the Japanese government, other nations, and the Japan Fisheries Association. Furthermore, we are the first to apply evolutionary game theory within this specific context. Our study not only broadens the knowledge frontier of evolutionary game models in environmental policy and international relations but also, through a tripartite framework, adeptly addresses the literature gap concerning the tripartite balance of interests in Japan’s marine discharge policy. In essence, our analysis, rooted in an international dimension, offers a comprehensive and realistic perspective on the dynamic games surrounding Japan’s nuclear wastewater issue. Symbol definitions All symbols and their corresponding meanings are detailed and illustrated in Table 1. Methodology Basic assumptions In this tripartite game, Japan (J), the National Fisheries Association of Japan (F), and other countries (C) are set as agents. Due to limited information and rationality, the decisions of each party are constrained. Employing the method of evolutionary game theory, the aforementioned three parties are considered as stakeholders. The dynamic game among these stakeholders involves their respective pursuits of interest and choices of action. In this game, Japan (J), as the country discharging nuclear-contaminated water, might weigh domestic economic interests against international reputation, deciding whether to continue the discharge of nuclear-contaminated water. The probability of discharging is denoted as x, and conversely, as 1 − x. Other countries (C), as additional stakeholders, may exhibit attitudes of either sanctioning or not sanctioning. Some countries might express dissatisfaction with Japan’s discharge strategy and adopt sanction measures with a probability of y, to protect international environmental and marine resource interests. Conversely, some might choose not to sanction, with a probability of 1 − y. The National Fisheries Association of Japan (F) represents the interests of Japanese fisheries. Their decisions may manifest in two forms: one opposing Japan’s discharge strategy with a probability of z, perceiving it as detrimental to fishery resources; the other, with a probability of 1 − z, might not oppose and thus would endorse Japan’s discharge actions. This study solely analyzes the game strategies of each party from an economic perspective and does not consider the impact of political and other factors on the strategic choices of each participant. Each participant’s occurrence in all possible scenarios under different strategies is illustrated in Fig 1. Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 1. https://doi.org/10.1371/journal.pone.0317419.g001 Evolutionary game model The payoff matrix for the Japanese government, other countries, and the Japanese Fisheries Association is shown in Table 2. Download: PPT PowerPoint slide PNG larger image TIFF original image Table 2. Payoff matrix of Japan, other countries, and japanese fisheries association. https://doi.org/10.1371/journal.pone.0317419.t002 Assume that the utility of the Japanese government choosing to discharge nuclear wastewater is denoted as U11, which encapsulates the benefits and costs associated with this action, including economic gains and potential international backlash. Conversely, the utility of the Japanese government deciding against discharging is denoted as U12. This reflects the potential benefits from maintaining environmental standards and avoiding international sanctions, minus any internal costs of managing stored wastewater. The average expected utility for Japan, taking into account the probability of each strategy being employed, is then represented by . The formula for is given by: (1) where x is the probability of choosing to discharge. (2) (3) The replicator dynamic equation for the Japanese government, denoted as S(x), is: (4) For other countries, the expected utility of sanctioning the Japanese government is denoted as U21. This utility reflects the benefits derived from imposing sanctions aimed at enforcing international environmental norms, potentially including improved global environmental outcomes and international prestige, minus the costs of enforcing these sanctions. The utility of not sanctioning is denoted as U22, which could include maintaining trade relations with Japan or avoiding the diplomatic costs of sanctioning. The average expected utility for other countries, reflecting their collective strategy towards Japan, is represented by . It is calculated as follows: (5) where y is the probability of sanctioning Japan. (6) (7) The replicator dynamic equation for other countries, denoted as G(y), is: (8) Assume that the utility of the Japan Fisheries Association opposing the discharge of nuclear wastewater is represented by U31. This utility might include the benefits of preserving marine life and fisheries resources, which are critical to their economic interests, offset by any potential conflicts with government policies or economic repercussions from reduced national industrial activities. Conversely, the utility of the Fisheries Association supporting the discharge is denoted as U32. This might reflect the potential benefits of aligning with government policies or economic gains from continued industrial operations, minus the environmental costs and potential long-term damage to fisheries. The average expected utility for the Japan Fisheries Association, which considers the probabilities of opposing or supporting the discharge, is calculated as: (9) where z is the probability of the Association opposing the discharge. (10) (11) (12) The replicator dynamic equation for the Japan Fisheries Association, denoted as P(z), is: (13) To further discuss the stable points in evolutionary game theory, we combine the Eqs (4), (8) and (13) to obtain the system of replicator dynamic equations as follows: (14) Basic assumptions In this tripartite game, Japan (J), the National Fisheries Association of Japan (F), and other countries (C) are set as agents. Due to limited information and rationality, the decisions of each party are constrained. Employing the method of evolutionary game theory, the aforementioned three parties are considered as stakeholders. The dynamic game among these stakeholders involves their respective pursuits of interest and choices of action. In this game, Japan (J), as the country discharging nuclear-contaminated water, might weigh domestic economic interests against international reputation, deciding whether to continue the discharge of nuclear-contaminated water. The probability of discharging is denoted as x, and conversely, as 1 − x. Other countries (C), as additional stakeholders, may exhibit attitudes of either sanctioning or not sanctioning. Some countries might express dissatisfaction with Japan’s discharge strategy and adopt sanction measures with a probability of y, to protect international environmental and marine resource interests. Conversely, some might choose not to sanction, with a probability of 1 − y. The National Fisheries Association of Japan (F) represents the interests of Japanese fisheries. Their decisions may manifest in two forms: one opposing Japan’s discharge strategy with a probability of z, perceiving it as detrimental to fishery resources; the other, with a probability of 1 − z, might not oppose and thus would endorse Japan’s discharge actions. This study solely analyzes the game strategies of each party from an economic perspective and does not consider the impact of political and other factors on the strategic choices of each participant. Each participant’s occurrence in all possible scenarios under different strategies is illustrated in Fig 1. Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 1. https://doi.org/10.1371/journal.pone.0317419.g001 Evolutionary game model The payoff matrix for the Japanese government, other countries, and the Japanese Fisheries Association is shown in Table 2. Download: PPT PowerPoint slide PNG larger image TIFF original image Table 2. Payoff matrix of Japan, other countries, and japanese fisheries association. https://doi.org/10.1371/journal.pone.0317419.t002 Assume that the utility of the Japanese government choosing to discharge nuclear wastewater is denoted as U11, which encapsulates the benefits and costs associated with this action, including economic gains and potential international backlash. Conversely, the utility of the Japanese government deciding against discharging is denoted as U12. This reflects the potential benefits from maintaining environmental standards and avoiding international sanctions, minus any internal costs of managing stored wastewater. The average expected utility for Japan, taking into account the probability of each strategy being employed, is then represented by . The formula for is given by: (1) where x is the probability of choosing to discharge. (2) (3) The replicator dynamic equation for the Japanese government, denoted as S(x), is: (4) For other countries, the expected utility of sanctioning the Japanese government is denoted as U21. This utility reflects the benefits derived from imposing sanctions aimed at enforcing international environmental norms, potentially including improved global environmental outcomes and international prestige, minus the costs of enforcing these sanctions. The utility of not sanctioning is denoted as U22, which could include maintaining trade relations with Japan or avoiding the diplomatic costs of sanctioning. The average expected utility for other countries, reflecting their collective strategy towards Japan, is represented by . It is calculated as follows: (5) where y is the probability of sanctioning Japan. (6) (7) The replicator dynamic equation for other countries, denoted as G(y), is: (8) Assume that the utility of the Japan Fisheries Association opposing the discharge of nuclear wastewater is represented by U31. This utility might include the benefits of preserving marine life and fisheries resources, which are critical to their economic interests, offset by any potential conflicts with government policies or economic repercussions from reduced national industrial activities. Conversely, the utility of the Fisheries Association supporting the discharge is denoted as U32. This might reflect the potential benefits of aligning with government policies or economic gains from continued industrial operations, minus the environmental costs and potential long-term damage to fisheries. The average expected utility for the Japan Fisheries Association, which considers the probabilities of opposing or supporting the discharge, is calculated as: (9) where z is the probability of the Association opposing the discharge. (10) (11) (12) The replicator dynamic equation for the Japan Fisheries Association, denoted as P(z), is: (13) To further discuss the stable points in evolutionary game theory, we combine the Eqs (4), (8) and (13) to obtain the system of replicator dynamic equations as follows: (14) Evolutionary stability analysis Asymptotic stability analysis of the three parties Let the dynamic replicator equations (14) be defined as S(x) = G(y) = P(z) = 0. By solving these, we can identify eight pure strategy equilibria and one mixed strategy equilibrium, which are γ1(0, 0, 0)T, γ2(1, 0, 0)T, γ3(0, 1, 0)T, γ4(0, 0, 1)T, γ5(1, 1, 0)T, γ6(1, 0, 1)T, γ7(0, 1, 1)T, γ8(1, 1, 1)T, and γ9(x*, y*, z*)T, where γ9(x*, y*, z*)T is determined by the following equations: (15) In asymmetric games, Evolutionarily Stable Strategies (ESS) are pure strategy equilibria, hence we only consider the asymptotic stability of pure strategies. The stability of equilibrium points in evolutionary games is typically determined by the stability of the Jacobian matrix, as illustrated below. (16) where (17) (18) (19) (20) (21) (22) (23) (24) (25) According to Lyapunov stability theory, the stability of equilibrium points in evolutionary games can be determined by the eigenvalues of the Jacobian matrix [34]. An equilibrium point is asymptotically stable, and the strategy is an Evolutionarily Stable Strategy (ESS), only if all the eigenvalues of the Jacobian matrix are negative. By analyzing the computed results for each point (in Appendix), the stability analysis of these points can be found as shown in Table 3. Download: PPT PowerPoint slide PNG larger image TIFF original image Table 3. Stability analysis of equilibrium point. https://doi.org/10.1371/journal.pone.0317419.t003 Table 3 shows that the equilibrium points γ1(0, 0, 0)T, γ2(1, 0, 0)T, γ3(0, 1, 0)T, γ5(1, 1, 0)T and γ7(0, 1, 1)T all have positive eigenvalues, hence they are unstable points. The equilibrium points γ4(0, 0, 1)T, γ6(1, 0, 1)T and γ8(1, 1, 1)T can all have negative eigenvalues. Therefore, they are stable points. When CSJ < (CLF + CMJ + CDJ), the eigenvalues of the Jacobian matrix corresponding to the evolutionary stable point γ4(0, 0, 1)T are all negative, indicating that γ4(0, 0, 1)T is an evolutionarily stable point. Under this condition, when the cost for Japan to store nuclear wastewater is less than the sum of the costs of discharging it into the sea, marine monitoring fees, and litigation compensation to the Fisheries Association, Japan is more likely to choose the “non-discharge” strategy. By adopting this strategy, Japan can reduce the direct impact on the marine ecosystem, respond to public concerns, protect the image and reputation of the government, potentially gain a competitive advantage in the field of environmental protection technology, and promote international cooperation. Other countries tend to choose the “non-sanction” strategy, as they might worry that sanctioning Japan when it chooses “non-discharge” could have negative impacts on bilateral trade, investment, and economic cooperation, thereby damaging their own economic interests, leading them to choose a more conservative strategy. For the Japanese Fisheries Association, it tends to choose the “oppose” strategy because, even under this condition where the cost of storing nuclear wastewater is less than other related costs, storing nuclear wastewater cannot completely eliminate the potential risks to the marine environment. The association needs to consider that such risks might cause long-term damage to the marine ecosystem and fishery resources, bringing uncertainty to fishing activities and the livelihoods of fishermen. When CSJ > (CLF + CMJ + CDJ) and BSP + CLC < CSC, the eigenvalues of the Jacobian matrix corresponding to the evolutionary stable point γ6(1, 0, 1)T are all negative, indicating that γ6(1, 0, 1)T is an evolutionarily stable point. At this time, the cost for Japan to store nuclear wastewater exceeds the sum of the costs of discharging it into the sea, marine monitoring fees, and litigation compensation to the Japanese Fisheries Association; the potential benefits and litigation compensation obtained by other countries from introducing substitutes for Japanese seafood are less than the additional costs for other countries to develop their own seafood products. In this scenario, Japan is more likely to choose the “discharge” strategy, allowing it to reduce the costs associated with storing nuclear wastewater. Other countries are more inclined to choose the “non-sanction Japan” strategy, as they can fill the domestic supply gap of Japanese seafood by introducing substitutes and increase their own seafood exports, hence they might prefer introducing seafood substitutes over pursuing litigation compensation. The Japanese Fisheries Association may choose to “oppose” due to a sharp decline in seafood income caused by marine pollution. When CSJ > (CDJ + CHJ + CLC + CLF + CMJ + IJ + TRJ) and CSC < CLC + BSP, all the eigenvalues of the Jacobian matrix corresponding to the evolutionary stable point γ8(1, 1, 1)T are negative. This denotes that γ8(1, 1, 1)T is an evolutionarily stable point. Under such circumstances, the cost incurred by Japan to store nuclear wastewater surpasses the aggregate of the costs related to sea discharge, acquiring aid from other nations, compensations from litigations to other nations and the domestic Fisheries Association, marine surveillance fees, impacts on Japan’s global image, and the decrement in export tax revenue due to sea discharge. Concurrently, the supplementary cost for other nations to cultivate their seafood is inferior to the aggregate of the litigation compensations and the prospective gains to their seafood industry from adopting Japanese seafood substitutes. In this context, Japan is inclined to opt for the “discharge” strategy, mitigating the expenses linked to nuclear wastewater storage. The Japanese Fisheries Association is predisposed to “oppose,” as its cardinal objective is to safeguard and advance the domestic fisheries sector’s interests. The discharge of nuclear wastewater into the sea could inflict irreversible damages to the marine ecosystem and pose potential hazards to fisheries resources and the fishing sector. This discharge may contaminate fisheries resources and compromise the quality of the catches, adversely impacting the revenue and living conditions of the fishermen. When the additional cost for other nations to cultivate their seafood is less than the combined litigation compensations and the prospective advantages to their seafood industry from adopting Japanese seafood substitutes, other nations are likely to select the “sanction” strategy. This choice aims to preserve the interests of their seafood sector, secure economic gains from litigation compensations, and sustain a commendable international image and reputation. Asymptotic stability analysis of the three parties Let the dynamic replicator equations (14) be defined as S(x) = G(y) = P(z) = 0. By solving these, we can identify eight pure strategy equilibria and one mixed strategy equilibrium, which are γ1(0, 0, 0)T, γ2(1, 0, 0)T, γ3(0, 1, 0)T, γ4(0, 0, 1)T, γ5(1, 1, 0)T, γ6(1, 0, 1)T, γ7(0, 1, 1)T, γ8(1, 1, 1)T, and γ9(x*, y*, z*)T, where γ9(x*, y*, z*)T is determined by the following equations: (15) In asymmetric games, Evolutionarily Stable Strategies (ESS) are pure strategy equilibria, hence we only consider the asymptotic stability of pure strategies. The stability of equilibrium points in evolutionary games is typically determined by the stability of the Jacobian matrix, as illustrated below. (16) where (17) (18) (19) (20) (21) (22) (23) (24) (25) According to Lyapunov stability theory, the stability of equilibrium points in evolutionary games can be determined by the eigenvalues of the Jacobian matrix [34]. An equilibrium point is asymptotically stable, and the strategy is an Evolutionarily Stable Strategy (ESS), only if all the eigenvalues of the Jacobian matrix are negative. By analyzing the computed results for each point (in Appendix), the stability analysis of these points can be found as shown in Table 3. Download: PPT PowerPoint slide PNG larger image TIFF original image Table 3. Stability analysis of equilibrium point. https://doi.org/10.1371/journal.pone.0317419.t003 Table 3 shows that the equilibrium points γ1(0, 0, 0)T, γ2(1, 0, 0)T, γ3(0, 1, 0)T, γ5(1, 1, 0)T and γ7(0, 1, 1)T all have positive eigenvalues, hence they are unstable points. The equilibrium points γ4(0, 0, 1)T, γ6(1, 0, 1)T and γ8(1, 1, 1)T can all have negative eigenvalues. Therefore, they are stable points. When CSJ < (CLF + CMJ + CDJ), the eigenvalues of the Jacobian matrix corresponding to the evolutionary stable point γ4(0, 0, 1)T are all negative, indicating that γ4(0, 0, 1)T is an evolutionarily stable point. Under this condition, when the cost for Japan to store nuclear wastewater is less than the sum of the costs of discharging it into the sea, marine monitoring fees, and litigation compensation to the Fisheries Association, Japan is more likely to choose the “non-discharge” strategy. By adopting this strategy, Japan can reduce the direct impact on the marine ecosystem, respond to public concerns, protect the image and reputation of the government, potentially gain a competitive advantage in the field of environmental protection technology, and promote international cooperation. Other countries tend to choose the “non-sanction” strategy, as they might worry that sanctioning Japan when it chooses “non-discharge” could have negative impacts on bilateral trade, investment, and economic cooperation, thereby damaging their own economic interests, leading them to choose a more conservative strategy. For the Japanese Fisheries Association, it tends to choose the “oppose” strategy because, even under this condition where the cost of storing nuclear wastewater is less than other related costs, storing nuclear wastewater cannot completely eliminate the potential risks to the marine environment. The association needs to consider that such risks might cause long-term damage to the marine ecosystem and fishery resources, bringing uncertainty to fishing activities and the livelihoods of fishermen. When CSJ > (CLF + CMJ + CDJ) and BSP + CLC < CSC, the eigenvalues of the Jacobian matrix corresponding to the evolutionary stable point γ6(1, 0, 1)T are all negative, indicating that γ6(1, 0, 1)T is an evolutionarily stable point. At this time, the cost for Japan to store nuclear wastewater exceeds the sum of the costs of discharging it into the sea, marine monitoring fees, and litigation compensation to the Japanese Fisheries Association; the potential benefits and litigation compensation obtained by other countries from introducing substitutes for Japanese seafood are less than the additional costs for other countries to develop their own seafood products. In this scenario, Japan is more likely to choose the “discharge” strategy, allowing it to reduce the costs associated with storing nuclear wastewater. Other countries are more inclined to choose the “non-sanction Japan” strategy, as they can fill the domestic supply gap of Japanese seafood by introducing substitutes and increase their own seafood exports, hence they might prefer introducing seafood substitutes over pursuing litigation compensation. The Japanese Fisheries Association may choose to “oppose” due to a sharp decline in seafood income caused by marine pollution. When CSJ > (CDJ + CHJ + CLC + CLF + CMJ + IJ + TRJ) and CSC < CLC + BSP, all the eigenvalues of the Jacobian matrix corresponding to the evolutionary stable point γ8(1, 1, 1)T are negative. This denotes that γ8(1, 1, 1)T is an evolutionarily stable point. Under such circumstances, the cost incurred by Japan to store nuclear wastewater surpasses the aggregate of the costs related to sea discharge, acquiring aid from other nations, compensations from litigations to other nations and the domestic Fisheries Association, marine surveillance fees, impacts on Japan’s global image, and the decrement in export tax revenue due to sea discharge. Concurrently, the supplementary cost for other nations to cultivate their seafood is inferior to the aggregate of the litigation compensations and the prospective gains to their seafood industry from adopting Japanese seafood substitutes. In this context, Japan is inclined to opt for the “discharge” strategy, mitigating the expenses linked to nuclear wastewater storage. The Japanese Fisheries Association is predisposed to “oppose,” as its cardinal objective is to safeguard and advance the domestic fisheries sector’s interests. The discharge of nuclear wastewater into the sea could inflict irreversible damages to the marine ecosystem and pose potential hazards to fisheries resources and the fishing sector. This discharge may contaminate fisheries resources and compromise the quality of the catches, adversely impacting the revenue and living conditions of the fishermen. When the additional cost for other nations to cultivate their seafood is less than the combined litigation compensations and the prospective advantages to their seafood industry from adopting Japanese seafood substitutes, other nations are likely to select the “sanction” strategy. This choice aims to preserve the interests of their seafood sector, secure economic gains from litigation compensations, and sustain a commendable international image and reputation. Numerical simulation Numerical simulation of different stable points In the aforementioned tripartite evolutionary game model, the strategy equilibrium of each game participant is influenced by the strategy choices of other participants. To more intuitively explore the strategy selection process of the three stable equilibrium points mentioned above, we set different parameter values according to different scenarios for numerical simulation to analyze the evolutionary trajectories of stable strategies of the Japanese government, other countries, and the Japanese Fisheries Association. As can be seen from Table 3, the stable points are γ4(0, 0, 1)T, γ6(1, 0, 1)T, and γ8(1, 1, 1)T. To delve deeper into the dynamic characteristics and evolutionary processes of these stable points, we simulate the evolutionary trajectories of these three stable points. We have set a group of parameter values that meet the stable conditions for each stable point. These parameter values will be used to simulate the evolutionary trajectories of these stable points to understand their dynamic behavior more accurately. The specific parameter values are shown in Table 4, with the initial points all set to [0.5, 0.5, 0.5]. Download: PPT PowerPoint slide PNG larger image TIFF original image Table 4. Basic parameter settings for three stable equilibrium points. https://doi.org/10.1371/journal.pone.0317419.t004 Evolutionary trajectory under Condition 1. When the system satisfies Condition 1, i.e., CSJ < (CLF + CMJ + CDJ), we refer to the parameters of Condition 1 in Table 4 for analysis. This condition implies that when the cost for Japan to store nuclear wastewater is less than the total of other related costs (these costs include the cost of discharging into the sea, marine monitoring fees, and litigation compensation to the Japanese domestic Fisheries Association), the Japanese government is more likely to choose the “non-discharge” strategy. This is because, under such circumstances, choosing to store nuclear wastewater instead of discharging it into the sea is more economical for Japan from a financial perspective. The results of numerical simulation further confirm this viewpoint. The simulation results show that under this condition, the stable point of the system is (0, 0, 1). Fig 2a displays the evolution of the probability of each decision-maker over time, clearly showing how each entity gradually tends towards this stable point over time. Fig 2b illustrates the overall evolutionary trajectory under Condition 1, providing us with a macroscopic perspective to observe the dynamic behavior of the system. Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 2. https://doi.org/10.1371/journal.pone.0317419.g002 Evolutionary trajectory under Condition 2. When the system satisfies Condition 2, i.e., CSJ > (CLF + CMJ + CDJ) and BSP + CLC < CSC, we refer to the parameters of Condition 3 in Table 4 for analysis. This condition implies that the cost for Japan to store nuclear wastewater has exceeded the costs of discharging into the sea, marine monitoring fees, and litigation compensation to the Japanese domestic Fisheries Association; the potential benefits and litigation compensation obtained by other countries through the introduction of Japanese seafood substitutes are less than the additional costs of developing their own seafood. Under these circumstances, Japan is more likely to choose the “discharge” strategy. This decision involves considerations from multiple aspects, including economic costs, environmental impacts, social responsibility, and public opinion. The results of numerical simulation further confirm this viewpoint. The simulation results show that under this condition, the stable point of the system is (1, 0, 1). Fig 3a displays the evolution of the probability of each decision-maker over time, clearly showing how each entity gradually tends towards this stable point over time. Fig 3b illustrates the overall evolutionary trajectory under Condition 2, providing us with a macroscopic perspective to observe the dynamic behavior of the system. Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 3. https://doi.org/10.1371/journal.pone.0317419.g003 Evolutionary trajectory under Condition 3. When the system satisfies Condition 3, i.e., CSJ > (CDJ + CHJ + CLC + CLF + CMJ + IJ + TRJ) and CSC < CLC + BSP, we refer to the parameters of Condition 3 in Table 4 for analysis. The cost for Japan to store nuclear wastewater has already exceeded the costs of discharging into the sea, marine monitoring fees, loss of international image, reduction in export taxes, and litigation compensation to the Japanese domestic Fisheries Association and other countries; the additional cost for other countries to develop their own seafood is less than the potential benefits and litigation compensation obtained by introducing Japanese seafood substitutes. Under these circumstances, Japan is more likely to choose the “discharge” strategy. The results of numerical simulation further confirm this viewpoint. The simulation results show that under this condition, the stable point of the system is (1, 0, 1). Fig 4a displays the evolution of the probability of each decision-maker over time, clearly showing how each entity gradually tends towards this stable point over time. Fig 4b illustrates the overall evolutionary trajectory under Condition 3, providing us with a macroscopic perspective to observe the dynamic behavior of the system. Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 4. https://doi.org/10.1371/journal.pone.0317419.g004 Impact of key parameters on evolutionary trajectories Our goal is for Japan to ultimately choose the strategy to stop discharging into the sea, i.e., Condition 1. To delve deeper into the key parameters under this condition, such as the cost of Japan discharging into the sea (CDJ), the cost of Japan storing nuclear wastewater (CSJ), litigation compensation from other countries (CLC), litigation compensation from the Fisheries Association (CLF), aid received by Japan from other countries (CHJ), and the initial strategy choice probabilities of Japan, other countries, and the Japanese Fisheries Association (x0, y0, z0) on the Evolutionarily Stable Strategy (ESS) of the three parties, this paper has conducted a series of numerical simulations under Condition 1. Under these settings, we explore the impact of these parameters on the evolution results and trajectories of the discharge strategy of the Japanese government, the sanction strategy of other countries, and the opposition strategy of the Japanese Fisheries Association, with specific parameter values as shown in Table 5. Our simulation results reveal how these key parameters influence the strategy choices and evolutionary dynamics of each party in the game, providing in-depth insights for our understanding of the formulation and implementation of international environmental policies. Download: PPT PowerPoint slide PNG larger image TIFF original image Table 5. Initial parameter settings under Condition 1. https://doi.org/10.1371/journal.pone.0317419.t005 Cost of Japan’s discharge into the sea CDJ. We explore the cost of Japan’s discharge into the sea, denoted as CDJ, on the evolutionary outcomes and trajectories of the three game entities. We assume CDJ = 1.0, 2.0, 3.0, 4.0, 5.0, and 6.0, while maintaining other parameters as shown in Table 5. The specific impacts of varying CDJ are illustrated in Fig 5. As illustrated in Fig 5a, with the increase in CDJ, the evolutionary stable point of the Japanese government, γ4(0, 0, 1)T, remains constant, but the evolutionary speed rapidly accelerates, meaning the Japanese government reaches the equilibrium point more swiftly. This might imply that higher discharge costs prompt the Japanese government to quickly solidify its strategy to mitigate additional costs induced by uncertainties. As depicted in Fig 5b, with the increase in CDJ, the evolutionary stable point of other countries, γ4(0, 0, 1)T, also remains constant, and the evolutionary speed does not significantly change. This might be attributed to the decision-making of other countries being predominantly influenced by their domestic factors and international political elements, having little association with Japan’s discharge costs. As shown in Fig 5c, with the increase in CDJ, the evolutionary stable point of the Japanese Fisheries Association, γ4(0, 0, 1)T, likewise remains constant, but the evolutionary speed decelerates, meaning the Japanese Fisheries Association reaches the equilibrium point more slowly. This might be because the Japanese Fisheries Association requires more time to assess and adapt to higher discharge costs and the potential impacts of these costs on fisheries and related industries. Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 5. https://doi.org/10.1371/journal.pone.0317419.g005 In conclusion, higher costs of processing nuclear wastewater not only augment the probability of Japan adopting the discharge strategy but also expedite its evolutionary speed. The evolutionary speed and strategy selection of other countries and the Japanese Fisheries Association, on the other hand, are more influenced by their respective internal and external factors. Cost of storing nuclear wastewater in Japan CSJ. We explore the cost of storing nuclear wastewater in Japan, denoted as CSJ, on the evolutionary outcomes and trajectories of the three game entities. We assume CSJ = 25, 27, 29, 31, 33, and 35, while maintaining other parameters as shown in Table 5. The specific impacts of varying CSJ are illustrated in Fig 6. As CSJ increases, the evolutionary stable point of the Japanese government, γ4(0, 0, 1)T, remains constant, but the evolutionary speed decreases, implying that the Japanese government requires more time to reach the equilibrium point. This may suggest that higher costs of storing nuclear wastewater necessitate a longer deliberation period for the Japanese government to formulate its strategy, aiming to mitigate the additional costs brought about by uncertainties. As depicted in Fig 6b, with the increase in CSJ, the evolutionary stable point of other countries, γ4(0, 0, 1)T, also remains unchanged, and the evolutionary speed is essentially constant. This could be attributed to the fact that the decisions of other countries are predominantly influenced by their domestic considerations and international political factors, having little correlation with Japan’s cost of discharging into the sea. As shown in Fig 6c, with the increase in CSJ, the evolutionary stable point of the Japanese Fisheries Association, γ4(0, 0, 1)T, also remains constant, but the evolutionary speed consistently accelerates, meaning the Japanese Fisheries Association reaches the equilibrium point more rapidly. This indicates that the Japanese Fisheries Association places a higher emphasis on environmental issues and desires the Japanese government to adopt stricter environmental protection measures promptly. Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 6. https://doi.org/10.1371/journal.pone.0317419.g006 In conclusion, the varying costs of storing nuclear wastewater in Japan not only influence the strategic deliberation and evolutionary pace of the Japanese government but also reflect the environmental priorities and expectations of the Japanese Fisheries Association, while other countries remain relatively unaffected, maintaining their strategic stances based on their internal and international considerations. Litigation compensation to other countries CLC. We analyze the impact of litigation compensation to other countries, denoted as CLC, on the evolutionary outcomes and trajectories of the three game entities. We assume CLC = 1.0, 5.0, 15.0, 20.0, 25.0, and 30.0, while keeping other parameters constant as shown in Table 5. The specific impacts of varying CLC are illustrated in Fig 7. Fig 7a demonstrates that as CLC increases, the evolutionary stable point of the Japanese government, γ4(0, 0, 1)T, remains constant, but the evolutionary speed slows down, indicating that the Japanese government takes a longer time to reach the equilibrium point. This could signify that higher litigation costs from other countries impose greater economic pressure and uncertainty on Japan, necessitating a more cautious and deliberate approach in strategy formulation. As depicted in Fig 7b, with the increase in CLC, the evolutionary stable point of other countries, γ4(0, 0, 1)T, also remains constant, and the evolutionary speed gradually decelerates. This could be interpreted as other countries facing an increased number of legal disputes and litigations, requiring more resources and time to address these legal matters, thereby inducing greater economic pressure and uncertainty on their own economies. Fig 7c illustrates that with the increase in CLC, the evolutionary stable point of the Japanese Fisheries Association, γ4(0, 0, 1)T, also remains constant, but the evolutionary speed progressively slows down, meaning the Japanese Fisheries Association reaches the equilibrium point more slowly. This could be attributed to the Japanese Fisheries Association needing more time to assess and adapt to the potential impacts of higher litigation compensations paid by the Japanese government to other countries on the fisheries and related industries in Japan. Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 7. https://doi.org/10.1371/journal.pone.0317419.g007 In conclusion, the varying levels of litigation compensation to other countries not only influence the strategic deliberation and evolutionary pace of the Japanese government and other countries but also reflect the adaptive responses and strategic considerations of the Japanese Fisheries Association in the face of increased economic pressures and uncertainties. Litigation compensation to fisheries association CLF. We analyze the impact of litigation compensation to the Fisheries Association, denoted as CLF, on the evolutionary outcomes and trajectories of the three game entities. We assume CLF = 30, 32, 34, 36, 38, and 40, while keeping other parameters constant as shown in Table 5. The specific impacts of varying CLF are illustrated in Fig 8. Fig 8a demonstrates that as CLF increases, the evolutionary stable point of the Japanese government, γ4(0, 0, 1)T, remains constant, but the evolutionary speed rapidly accelerates, meaning the Japanese government reaches the equilibrium point more swiftly. This could signify that higher litigation costs to the Fisheries Association prompt the Japanese government to quickly solidify its strategy to mitigate additional costs induced by uncertainties, emphasizing the government’s strategic adaptability in response to internal compensatory obligations. As depicted in Fig 8b, with the increase in CLF, the evolutionary stable point of other countries, γ4(0, 0, 1)T, also remains constant, and the evolutionary speed does not significantly change. This could be interpreted as other countries’ decision-making being predominantly influenced by their domestic factors and international political elements, having little association with the compensations paid by the Japanese government to its domestic Fisheries Association, reflecting the relative independence of international actors from Japan’s internal compensatory dynamics. Fig 8c illustrates that with the increase in CLF, the evolutionary stable point of the Japanese Fisheries Association, γ4(0, 0, 1)T, also remains constant, but the evolutionary speed progressively slows down, meaning the Japanese Fisheries Association reaches the equilibrium point more slowly. This could be attributed to the Japanese Fisheries Association needing more time to assess and adapt to the potential impacts of receiving higher litigation compensations. Despite the potential influx of compensatory funds, the association must weigh the long-term considerations for the fishermen and the future development of fisheries livelihoods, highlighting the intricate balance between immediate financial relief and long-term sustainable development. Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 8. https://doi.org/10.1371/journal.pone.0317419.g008 In conclusion, the varying levels of litigation compensation to the Fisheries Association not only influence the strategic deliberation and evolutionary pace of the Japanese government and the Japanese Fisheries Association but also underscore the nuanced interplay between economic compensations and sustainable strategic considerations in the context of environmental disputes. Aid received by Japan from other countries CHJ. We explore the implications of the aid received by Japan from other countries, denoted as CHJ, on the evolutionary outcomes and trajectories of the three game entities. We assume CHJ = 5, 9, 13, and 17, while maintaining other parameters as shown in Table 5. The specific impacts of varying CHJ are illustrated in Fig 9. Fig 9a demonstrates that as CHJ increases, the evolutionary stable point of the Japanese government, γ4(0, 0, 1)T, remains constant, but the evolutionary speed progressively decelerates, meaning the Japanese government reaches the equilibrium point more slowly. This suggests that the aid from other countries may induce the Japanese government to deliberate more meticulously on its strategies, with the slower evolutionary speed aimed at mitigating additional costs induced by uncertainties. Enhanced international aid from other countries acts as a deterrent to Japan’s discharge strategy and prolongs the time to reach the evolutionary stable point, serving as an incentivizing measure for Japan to refrain from discharging. As depicted in Fig 9b, with the increase in CHJ, the evolutionary stable point of other countries, γ4(0, 0, 1)T, also remains constant, but the evolutionary speed accelerates. This could be interpreted as the larger the international aid provided by other countries, the slower the progression of the crisis, and when improvements are observed, other countries also accelerate their evolution, reflecting the adaptive nature of international actors in response to the evolving crisis dynamics. Fig 9c illustrates that with the increase in CHJ, the evolutionary stable point of the Japanese Fisheries Association, γ4(0, 0, 1)T, also remains constant, but the evolutionary speed accelerates, meaning the Japanese Fisheries Association reaches the equilibrium point more swiftly. This could be attributed to the Japanese Fisheries Association perceiving the aid from other countries as an incentive that could potentially foster the development of Japanese fisheries and secure the income of fishermen, highlighting the perceived positive impact of international aid on domestic industries. Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 9. https://doi.org/10.1371/journal.pone.0317419.g009 In conclusion, the varying levels of international aid received by Japan not only influence the strategic deliberation and evolutionary pace of the involved entities but also underscore the intricate interplay between international cooperation and strategic environmental decision-making in the context of international environmental disputes. Initial strategy selection probabilities of Japan, other countries, and the japanese fisheries Association: x0, y0, z0. Assuming (x0, y0, z0) = (0.5, 0.5, 0.5), (0.8, 0.1, 0.1), (0.2, 0.7, 0.1), (0.7, 0.2, 0.1), while maintaining other parameters as shown in Table 5, we can observe the evolutionary outcomes and trajectories of the three game entities under varying initial strategy selection probabilities, as depicted in Fig 10. From the figure, it is evident that when initial strategy selection probabilities vary, the evolutionary stable point γ4(0, 0, 1)T remains constant, indicating that it does not affect the evolutionary outcomes but does influence the evolutionary trajectories. For Japan, a stronger initial strategy selection probability, x0, will delay Japan’s discharge strategy and decelerate its evolutionary speed, as illustrated in Fig 10a. For other countries, the evolutionary stable point γ4(0, 0, 1)T also remains constant when initial strategy selection probabilities vary, signifying that it does not impact the evolutionary outcomes but does have certain effects on the evolutionary trajectories. A stronger initial strategy selection probability, y0, will delay the opposition strategy of the Fisheries Association and slow down its evolutionary speed, as shown in Fig 10b. Lastly, for the Fisheries Association, the evolutionary stable point γ4(0, 0, 1)T still remains constant when initial strategy selection probabilities change, suggesting that it does not alter the evolutionary outcomes but does affect the evolutionary trajectories. A stronger initial strategy selection probability, z0, will delay the opposition strategy of the Fisheries Association and decelerate its evolutionary speed, as depicted in Fig 10c. Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 10. https://doi.org/10.1371/journal.pone.0317419.g010 Numerical simulation of different stable points In the aforementioned tripartite evolutionary game model, the strategy equilibrium of each game participant is influenced by the strategy choices of other participants. To more intuitively explore the strategy selection process of the three stable equilibrium points mentioned above, we set different parameter values according to different scenarios for numerical simulation to analyze the evolutionary trajectories of stable strategies of the Japanese government, other countries, and the Japanese Fisheries Association. As can be seen from Table 3, the stable points are γ4(0, 0, 1)T, γ6(1, 0, 1)T, and γ8(1, 1, 1)T. To delve deeper into the dynamic characteristics and evolutionary processes of these stable points, we simulate the evolutionary trajectories of these three stable points. We have set a group of parameter values that meet the stable conditions for each stable point. These parameter values will be used to simulate the evolutionary trajectories of these stable points to understand their dynamic behavior more accurately. The specific parameter values are shown in Table 4, with the initial points all set to [0.5, 0.5, 0.5]. Download: PPT PowerPoint slide PNG larger image TIFF original image Table 4. Basic parameter settings for three stable equilibrium points. https://doi.org/10.1371/journal.pone.0317419.t004 Evolutionary trajectory under Condition 1. When the system satisfies Condition 1, i.e., CSJ < (CLF + CMJ + CDJ), we refer to the parameters of Condition 1 in Table 4 for analysis. This condition implies that when the cost for Japan to store nuclear wastewater is less than the total of other related costs (these costs include the cost of discharging into the sea, marine monitoring fees, and litigation compensation to the Japanese domestic Fisheries Association), the Japanese government is more likely to choose the “non-discharge” strategy. This is because, under such circumstances, choosing to store nuclear wastewater instead of discharging it into the sea is more economical for Japan from a financial perspective. The results of numerical simulation further confirm this viewpoint. The simulation results show that under this condition, the stable point of the system is (0, 0, 1). Fig 2a displays the evolution of the probability of each decision-maker over time, clearly showing how each entity gradually tends towards this stable point over time. Fig 2b illustrates the overall evolutionary trajectory under Condition 1, providing us with a macroscopic perspective to observe the dynamic behavior of the system. Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 2. https://doi.org/10.1371/journal.pone.0317419.g002 Evolutionary trajectory under Condition 2. When the system satisfies Condition 2, i.e., CSJ > (CLF + CMJ + CDJ) and BSP + CLC < CSC, we refer to the parameters of Condition 3 in Table 4 for analysis. This condition implies that the cost for Japan to store nuclear wastewater has exceeded the costs of discharging into the sea, marine monitoring fees, and litigation compensation to the Japanese domestic Fisheries Association; the potential benefits and litigation compensation obtained by other countries through the introduction of Japanese seafood substitutes are less than the additional costs of developing their own seafood. Under these circumstances, Japan is more likely to choose the “discharge” strategy. This decision involves considerations from multiple aspects, including economic costs, environmental impacts, social responsibility, and public opinion. The results of numerical simulation further confirm this viewpoint. The simulation results show that under this condition, the stable point of the system is (1, 0, 1). Fig 3a displays the evolution of the probability of each decision-maker over time, clearly showing how each entity gradually tends towards this stable point over time. Fig 3b illustrates the overall evolutionary trajectory under Condition 2, providing us with a macroscopic perspective to observe the dynamic behavior of the system. Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 3. https://doi.org/10.1371/journal.pone.0317419.g003 Evolutionary trajectory under Condition 3. When the system satisfies Condition 3, i.e., CSJ > (CDJ + CHJ + CLC + CLF + CMJ + IJ + TRJ) and CSC < CLC + BSP, we refer to the parameters of Condition 3 in Table 4 for analysis. The cost for Japan to store nuclear wastewater has already exceeded the costs of discharging into the sea, marine monitoring fees, loss of international image, reduction in export taxes, and litigation compensation to the Japanese domestic Fisheries Association and other countries; the additional cost for other countries to develop their own seafood is less than the potential benefits and litigation compensation obtained by introducing Japanese seafood substitutes. Under these circumstances, Japan is more likely to choose the “discharge” strategy. The results of numerical simulation further confirm this viewpoint. The simulation results show that under this condition, the stable point of the system is (1, 0, 1). Fig 4a displays the evolution of the probability of each decision-maker over time, clearly showing how each entity gradually tends towards this stable point over time. Fig 4b illustrates the overall evolutionary trajectory under Condition 3, providing us with a macroscopic perspective to observe the dynamic behavior of the system. Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 4. https://doi.org/10.1371/journal.pone.0317419.g004 Evolutionary trajectory under Condition 1. When the system satisfies Condition 1, i.e., CSJ < (CLF + CMJ + CDJ), we refer to the parameters of Condition 1 in Table 4 for analysis. This condition implies that when the cost for Japan to store nuclear wastewater is less than the total of other related costs (these costs include the cost of discharging into the sea, marine monitoring fees, and litigation compensation to the Japanese domestic Fisheries Association), the Japanese government is more likely to choose the “non-discharge” strategy. This is because, under such circumstances, choosing to store nuclear wastewater instead of discharging it into the sea is more economical for Japan from a financial perspective. The results of numerical simulation further confirm this viewpoint. The simulation results show that under this condition, the stable point of the system is (0, 0, 1). Fig 2a displays the evolution of the probability of each decision-maker over time, clearly showing how each entity gradually tends towards this stable point over time. Fig 2b illustrates the overall evolutionary trajectory under Condition 1, providing us with a macroscopic perspective to observe the dynamic behavior of the system. Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 2. https://doi.org/10.1371/journal.pone.0317419.g002 Evolutionary trajectory under Condition 2. When the system satisfies Condition 2, i.e., CSJ > (CLF + CMJ + CDJ) and BSP + CLC < CSC, we refer to the parameters of Condition 3 in Table 4 for analysis. This condition implies that the cost for Japan to store nuclear wastewater has exceeded the costs of discharging into the sea, marine monitoring fees, and litigation compensation to the Japanese domestic Fisheries Association; the potential benefits and litigation compensation obtained by other countries through the introduction of Japanese seafood substitutes are less than the additional costs of developing their own seafood. Under these circumstances, Japan is more likely to choose the “discharge” strategy. This decision involves considerations from multiple aspects, including economic costs, environmental impacts, social responsibility, and public opinion. The results of numerical simulation further confirm this viewpoint. The simulation results show that under this condition, the stable point of the system is (1, 0, 1). Fig 3a displays the evolution of the probability of each decision-maker over time, clearly showing how each entity gradually tends towards this stable point over time. Fig 3b illustrates the overall evolutionary trajectory under Condition 2, providing us with a macroscopic perspective to observe the dynamic behavior of the system. Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 3. https://doi.org/10.1371/journal.pone.0317419.g003 Evolutionary trajectory under Condition 3. When the system satisfies Condition 3, i.e., CSJ > (CDJ + CHJ + CLC + CLF + CMJ + IJ + TRJ) and CSC < CLC + BSP, we refer to the parameters of Condition 3 in Table 4 for analysis. The cost for Japan to store nuclear wastewater has already exceeded the costs of discharging into the sea, marine monitoring fees, loss of international image, reduction in export taxes, and litigation compensation to the Japanese domestic Fisheries Association and other countries; the additional cost for other countries to develop their own seafood is less than the potential benefits and litigation compensation obtained by introducing Japanese seafood substitutes. Under these circumstances, Japan is more likely to choose the “discharge” strategy. The results of numerical simulation further confirm this viewpoint. The simulation results show that under this condition, the stable point of the system is (1, 0, 1). Fig 4a displays the evolution of the probability of each decision-maker over time, clearly showing how each entity gradually tends towards this stable point over time. Fig 4b illustrates the overall evolutionary trajectory under Condition 3, providing us with a macroscopic perspective to observe the dynamic behavior of the system. Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 4. https://doi.org/10.1371/journal.pone.0317419.g004 Impact of key parameters on evolutionary trajectories Our goal is for Japan to ultimately choose the strategy to stop discharging into the sea, i.e., Condition 1. To delve deeper into the key parameters under this condition, such as the cost of Japan discharging into the sea (CDJ), the cost of Japan storing nuclear wastewater (CSJ), litigation compensation from other countries (CLC), litigation compensation from the Fisheries Association (CLF), aid received by Japan from other countries (CHJ), and the initial strategy choice probabilities of Japan, other countries, and the Japanese Fisheries Association (x0, y0, z0) on the Evolutionarily Stable Strategy (ESS) of the three parties, this paper has conducted a series of numerical simulations under Condition 1. Under these settings, we explore the impact of these parameters on the evolution results and trajectories of the discharge strategy of the Japanese government, the sanction strategy of other countries, and the opposition strategy of the Japanese Fisheries Association, with specific parameter values as shown in Table 5. Our simulation results reveal how these key parameters influence the strategy choices and evolutionary dynamics of each party in the game, providing in-depth insights for our understanding of the formulation and implementation of international environmental policies. Download: PPT PowerPoint slide PNG larger image TIFF original image Table 5. Initial parameter settings under Condition 1. https://doi.org/10.1371/journal.pone.0317419.t005 Cost of Japan’s discharge into the sea CDJ. We explore the cost of Japan’s discharge into the sea, denoted as CDJ, on the evolutionary outcomes and trajectories of the three game entities. We assume CDJ = 1.0, 2.0, 3.0, 4.0, 5.0, and 6.0, while maintaining other parameters as shown in Table 5. The specific impacts of varying CDJ are illustrated in Fig 5. As illustrated in Fig 5a, with the increase in CDJ, the evolutionary stable point of the Japanese government, γ4(0, 0, 1)T, remains constant, but the evolutionary speed rapidly accelerates, meaning the Japanese government reaches the equilibrium point more swiftly. This might imply that higher discharge costs prompt the Japanese government to quickly solidify its strategy to mitigate additional costs induced by uncertainties. As depicted in Fig 5b, with the increase in CDJ, the evolutionary stable point of other countries, γ4(0, 0, 1)T, also remains constant, and the evolutionary speed does not significantly change. This might be attributed to the decision-making of other countries being predominantly influenced by their domestic factors and international political elements, having little association with Japan’s discharge costs. As shown in Fig 5c, with the increase in CDJ, the evolutionary stable point of the Japanese Fisheries Association, γ4(0, 0, 1)T, likewise remains constant, but the evolutionary speed decelerates, meaning the Japanese Fisheries Association reaches the equilibrium point more slowly. This might be because the Japanese Fisheries Association requires more time to assess and adapt to higher discharge costs and the potential impacts of these costs on fisheries and related industries. Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 5. https://doi.org/10.1371/journal.pone.0317419.g005 In conclusion, higher costs of processing nuclear wastewater not only augment the probability of Japan adopting the discharge strategy but also expedite its evolutionary speed. The evolutionary speed and strategy selection of other countries and the Japanese Fisheries Association, on the other hand, are more influenced by their respective internal and external factors. Cost of storing nuclear wastewater in Japan CSJ. We explore the cost of storing nuclear wastewater in Japan, denoted as CSJ, on the evolutionary outcomes and trajectories of the three game entities. We assume CSJ = 25, 27, 29, 31, 33, and 35, while maintaining other parameters as shown in Table 5. The specific impacts of varying CSJ are illustrated in Fig 6. As CSJ increases, the evolutionary stable point of the Japanese government, γ4(0, 0, 1)T, remains constant, but the evolutionary speed decreases, implying that the Japanese government requires more time to reach the equilibrium point. This may suggest that higher costs of storing nuclear wastewater necessitate a longer deliberation period for the Japanese government to formulate its strategy, aiming to mitigate the additional costs brought about by uncertainties. As depicted in Fig 6b, with the increase in CSJ, the evolutionary stable point of other countries, γ4(0, 0, 1)T, also remains unchanged, and the evolutionary speed is essentially constant. This could be attributed to the fact that the decisions of other countries are predominantly influenced by their domestic considerations and international political factors, having little correlation with Japan’s cost of discharging into the sea. As shown in Fig 6c, with the increase in CSJ, the evolutionary stable point of the Japanese Fisheries Association, γ4(0, 0, 1)T, also remains constant, but the evolutionary speed consistently accelerates, meaning the Japanese Fisheries Association reaches the equilibrium point more rapidly. This indicates that the Japanese Fisheries Association places a higher emphasis on environmental issues and desires the Japanese government to adopt stricter environmental protection measures promptly. Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 6. https://doi.org/10.1371/journal.pone.0317419.g006 In conclusion, the varying costs of storing nuclear wastewater in Japan not only influence the strategic deliberation and evolutionary pace of the Japanese government but also reflect the environmental priorities and expectations of the Japanese Fisheries Association, while other countries remain relatively unaffected, maintaining their strategic stances based on their internal and international considerations. Litigation compensation to other countries CLC. We analyze the impact of litigation compensation to other countries, denoted as CLC, on the evolutionary outcomes and trajectories of the three game entities. We assume CLC = 1.0, 5.0, 15.0, 20.0, 25.0, and 30.0, while keeping other parameters constant as shown in Table 5. The specific impacts of varying CLC are illustrated in Fig 7. Fig 7a demonstrates that as CLC increases, the evolutionary stable point of the Japanese government, γ4(0, 0, 1)T, remains constant, but the evolutionary speed slows down, indicating that the Japanese government takes a longer time to reach the equilibrium point. This could signify that higher litigation costs from other countries impose greater economic pressure and uncertainty on Japan, necessitating a more cautious and deliberate approach in strategy formulation. As depicted in Fig 7b, with the increase in CLC, the evolutionary stable point of other countries, γ4(0, 0, 1)T, also remains constant, and the evolutionary speed gradually decelerates. This could be interpreted as other countries facing an increased number of legal disputes and litigations, requiring more resources and time to address these legal matters, thereby inducing greater economic pressure and uncertainty on their own economies. Fig 7c illustrates that with the increase in CLC, the evolutionary stable point of the Japanese Fisheries Association, γ4(0, 0, 1)T, also remains constant, but the evolutionary speed progressively slows down, meaning the Japanese Fisheries Association reaches the equilibrium point more slowly. This could be attributed to the Japanese Fisheries Association needing more time to assess and adapt to the potential impacts of higher litigation compensations paid by the Japanese government to other countries on the fisheries and related industries in Japan. Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 7. https://doi.org/10.1371/journal.pone.0317419.g007 In conclusion, the varying levels of litigation compensation to other countries not only influence the strategic deliberation and evolutionary pace of the Japanese government and other countries but also reflect the adaptive responses and strategic considerations of the Japanese Fisheries Association in the face of increased economic pressures and uncertainties. Litigation compensation to fisheries association CLF. We analyze the impact of litigation compensation to the Fisheries Association, denoted as CLF, on the evolutionary outcomes and trajectories of the three game entities. We assume CLF = 30, 32, 34, 36, 38, and 40, while keeping other parameters constant as shown in Table 5. The specific impacts of varying CLF are illustrated in Fig 8. Fig 8a demonstrates that as CLF increases, the evolutionary stable point of the Japanese government, γ4(0, 0, 1)T, remains constant, but the evolutionary speed rapidly accelerates, meaning the Japanese government reaches the equilibrium point more swiftly. This could signify that higher litigation costs to the Fisheries Association prompt the Japanese government to quickly solidify its strategy to mitigate additional costs induced by uncertainties, emphasizing the government’s strategic adaptability in response to internal compensatory obligations. As depicted in Fig 8b, with the increase in CLF, the evolutionary stable point of other countries, γ4(0, 0, 1)T, also remains constant, and the evolutionary speed does not significantly change. This could be interpreted as other countries’ decision-making being predominantly influenced by their domestic factors and international political elements, having little association with the compensations paid by the Japanese government to its domestic Fisheries Association, reflecting the relative independence of international actors from Japan’s internal compensatory dynamics. Fig 8c illustrates that with the increase in CLF, the evolutionary stable point of the Japanese Fisheries Association, γ4(0, 0, 1)T, also remains constant, but the evolutionary speed progressively slows down, meaning the Japanese Fisheries Association reaches the equilibrium point more slowly. This could be attributed to the Japanese Fisheries Association needing more time to assess and adapt to the potential impacts of receiving higher litigation compensations. Despite the potential influx of compensatory funds, the association must weigh the long-term considerations for the fishermen and the future development of fisheries livelihoods, highlighting the intricate balance between immediate financial relief and long-term sustainable development. Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 8. https://doi.org/10.1371/journal.pone.0317419.g008 In conclusion, the varying levels of litigation compensation to the Fisheries Association not only influence the strategic deliberation and evolutionary pace of the Japanese government and the Japanese Fisheries Association but also underscore the nuanced interplay between economic compensations and sustainable strategic considerations in the context of environmental disputes. Aid received by Japan from other countries CHJ. We explore the implications of the aid received by Japan from other countries, denoted as CHJ, on the evolutionary outcomes and trajectories of the three game entities. We assume CHJ = 5, 9, 13, and 17, while maintaining other parameters as shown in Table 5. The specific impacts of varying CHJ are illustrated in Fig 9. Fig 9a demonstrates that as CHJ increases, the evolutionary stable point of the Japanese government, γ4(0, 0, 1)T, remains constant, but the evolutionary speed progressively decelerates, meaning the Japanese government reaches the equilibrium point more slowly. This suggests that the aid from other countries may induce the Japanese government to deliberate more meticulously on its strategies, with the slower evolutionary speed aimed at mitigating additional costs induced by uncertainties. Enhanced international aid from other countries acts as a deterrent to Japan’s discharge strategy and prolongs the time to reach the evolutionary stable point, serving as an incentivizing measure for Japan to refrain from discharging. As depicted in Fig 9b, with the increase in CHJ, the evolutionary stable point of other countries, γ4(0, 0, 1)T, also remains constant, but the evolutionary speed accelerates. This could be interpreted as the larger the international aid provided by other countries, the slower the progression of the crisis, and when improvements are observed, other countries also accelerate their evolution, reflecting the adaptive nature of international actors in response to the evolving crisis dynamics. Fig 9c illustrates that with the increase in CHJ, the evolutionary stable point of the Japanese Fisheries Association, γ4(0, 0, 1)T, also remains constant, but the evolutionary speed accelerates, meaning the Japanese Fisheries Association reaches the equilibrium point more swiftly. This could be attributed to the Japanese Fisheries Association perceiving the aid from other countries as an incentive that could potentially foster the development of Japanese fisheries and secure the income of fishermen, highlighting the perceived positive impact of international aid on domestic industries. Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 9. https://doi.org/10.1371/journal.pone.0317419.g009 In conclusion, the varying levels of international aid received by Japan not only influence the strategic deliberation and evolutionary pace of the involved entities but also underscore the intricate interplay between international cooperation and strategic environmental decision-making in the context of international environmental disputes. Initial strategy selection probabilities of Japan, other countries, and the japanese fisheries Association: x0, y0, z0. Assuming (x0, y0, z0) = (0.5, 0.5, 0.5), (0.8, 0.1, 0.1), (0.2, 0.7, 0.1), (0.7, 0.2, 0.1), while maintaining other parameters as shown in Table 5, we can observe the evolutionary outcomes and trajectories of the three game entities under varying initial strategy selection probabilities, as depicted in Fig 10. From the figure, it is evident that when initial strategy selection probabilities vary, the evolutionary stable point γ4(0, 0, 1)T remains constant, indicating that it does not affect the evolutionary outcomes but does influence the evolutionary trajectories. For Japan, a stronger initial strategy selection probability, x0, will delay Japan’s discharge strategy and decelerate its evolutionary speed, as illustrated in Fig 10a. For other countries, the evolutionary stable point γ4(0, 0, 1)T also remains constant when initial strategy selection probabilities vary, signifying that it does not impact the evolutionary outcomes but does have certain effects on the evolutionary trajectories. A stronger initial strategy selection probability, y0, will delay the opposition strategy of the Fisheries Association and slow down its evolutionary speed, as shown in Fig 10b. Lastly, for the Fisheries Association, the evolutionary stable point γ4(0, 0, 1)T still remains constant when initial strategy selection probabilities change, suggesting that it does not alter the evolutionary outcomes but does affect the evolutionary trajectories. A stronger initial strategy selection probability, z0, will delay the opposition strategy of the Fisheries Association and decelerate its evolutionary speed, as depicted in Fig 10c. Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 10. https://doi.org/10.1371/journal.pone.0317419.g010 Cost of Japan’s discharge into the sea CDJ. We explore the cost of Japan’s discharge into the sea, denoted as CDJ, on the evolutionary outcomes and trajectories of the three game entities. We assume CDJ = 1.0, 2.0, 3.0, 4.0, 5.0, and 6.0, while maintaining other parameters as shown in Table 5. The specific impacts of varying CDJ are illustrated in Fig 5. As illustrated in Fig 5a, with the increase in CDJ, the evolutionary stable point of the Japanese government, γ4(0, 0, 1)T, remains constant, but the evolutionary speed rapidly accelerates, meaning the Japanese government reaches the equilibrium point more swiftly. This might imply that higher discharge costs prompt the Japanese government to quickly solidify its strategy to mitigate additional costs induced by uncertainties. As depicted in Fig 5b, with the increase in CDJ, the evolutionary stable point of other countries, γ4(0, 0, 1)T, also remains constant, and the evolutionary speed does not significantly change. This might be attributed to the decision-making of other countries being predominantly influenced by their domestic factors and international political elements, having little association with Japan’s discharge costs. As shown in Fig 5c, with the increase in CDJ, the evolutionary stable point of the Japanese Fisheries Association, γ4(0, 0, 1)T, likewise remains constant, but the evolutionary speed decelerates, meaning the Japanese Fisheries Association reaches the equilibrium point more slowly. This might be because the Japanese Fisheries Association requires more time to assess and adapt to higher discharge costs and the potential impacts of these costs on fisheries and related industries. Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 5. https://doi.org/10.1371/journal.pone.0317419.g005 In conclusion, higher costs of processing nuclear wastewater not only augment the probability of Japan adopting the discharge strategy but also expedite its evolutionary speed. The evolutionary speed and strategy selection of other countries and the Japanese Fisheries Association, on the other hand, are more influenced by their respective internal and external factors. Cost of storing nuclear wastewater in Japan CSJ. We explore the cost of storing nuclear wastewater in Japan, denoted as CSJ, on the evolutionary outcomes and trajectories of the three game entities. We assume CSJ = 25, 27, 29, 31, 33, and 35, while maintaining other parameters as shown in Table 5. The specific impacts of varying CSJ are illustrated in Fig 6. As CSJ increases, the evolutionary stable point of the Japanese government, γ4(0, 0, 1)T, remains constant, but the evolutionary speed decreases, implying that the Japanese government requires more time to reach the equilibrium point. This may suggest that higher costs of storing nuclear wastewater necessitate a longer deliberation period for the Japanese government to formulate its strategy, aiming to mitigate the additional costs brought about by uncertainties. As depicted in Fig 6b, with the increase in CSJ, the evolutionary stable point of other countries, γ4(0, 0, 1)T, also remains unchanged, and the evolutionary speed is essentially constant. This could be attributed to the fact that the decisions of other countries are predominantly influenced by their domestic considerations and international political factors, having little correlation with Japan’s cost of discharging into the sea. As shown in Fig 6c, with the increase in CSJ, the evolutionary stable point of the Japanese Fisheries Association, γ4(0, 0, 1)T, also remains constant, but the evolutionary speed consistently accelerates, meaning the Japanese Fisheries Association reaches the equilibrium point more rapidly. This indicates that the Japanese Fisheries Association places a higher emphasis on environmental issues and desires the Japanese government to adopt stricter environmental protection measures promptly. Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 6. https://doi.org/10.1371/journal.pone.0317419.g006 In conclusion, the varying costs of storing nuclear wastewater in Japan not only influence the strategic deliberation and evolutionary pace of the Japanese government but also reflect the environmental priorities and expectations of the Japanese Fisheries Association, while other countries remain relatively unaffected, maintaining their strategic stances based on their internal and international considerations. Litigation compensation to other countries CLC. We analyze the impact of litigation compensation to other countries, denoted as CLC, on the evolutionary outcomes and trajectories of the three game entities. We assume CLC = 1.0, 5.0, 15.0, 20.0, 25.0, and 30.0, while keeping other parameters constant as shown in Table 5. The specific impacts of varying CLC are illustrated in Fig 7. Fig 7a demonstrates that as CLC increases, the evolutionary stable point of the Japanese government, γ4(0, 0, 1)T, remains constant, but the evolutionary speed slows down, indicating that the Japanese government takes a longer time to reach the equilibrium point. This could signify that higher litigation costs from other countries impose greater economic pressure and uncertainty on Japan, necessitating a more cautious and deliberate approach in strategy formulation. As depicted in Fig 7b, with the increase in CLC, the evolutionary stable point of other countries, γ4(0, 0, 1)T, also remains constant, and the evolutionary speed gradually decelerates. This could be interpreted as other countries facing an increased number of legal disputes and litigations, requiring more resources and time to address these legal matters, thereby inducing greater economic pressure and uncertainty on their own economies. Fig 7c illustrates that with the increase in CLC, the evolutionary stable point of the Japanese Fisheries Association, γ4(0, 0, 1)T, also remains constant, but the evolutionary speed progressively slows down, meaning the Japanese Fisheries Association reaches the equilibrium point more slowly. This could be attributed to the Japanese Fisheries Association needing more time to assess and adapt to the potential impacts of higher litigation compensations paid by the Japanese government to other countries on the fisheries and related industries in Japan. Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 7. https://doi.org/10.1371/journal.pone.0317419.g007 In conclusion, the varying levels of litigation compensation to other countries not only influence the strategic deliberation and evolutionary pace of the Japanese government and other countries but also reflect the adaptive responses and strategic considerations of the Japanese Fisheries Association in the face of increased economic pressures and uncertainties. Litigation compensation to fisheries association CLF. We analyze the impact of litigation compensation to the Fisheries Association, denoted as CLF, on the evolutionary outcomes and trajectories of the three game entities. We assume CLF = 30, 32, 34, 36, 38, and 40, while keeping other parameters constant as shown in Table 5. The specific impacts of varying CLF are illustrated in Fig 8. Fig 8a demonstrates that as CLF increases, the evolutionary stable point of the Japanese government, γ4(0, 0, 1)T, remains constant, but the evolutionary speed rapidly accelerates, meaning the Japanese government reaches the equilibrium point more swiftly. This could signify that higher litigation costs to the Fisheries Association prompt the Japanese government to quickly solidify its strategy to mitigate additional costs induced by uncertainties, emphasizing the government’s strategic adaptability in response to internal compensatory obligations. As depicted in Fig 8b, with the increase in CLF, the evolutionary stable point of other countries, γ4(0, 0, 1)T, also remains constant, and the evolutionary speed does not significantly change. This could be interpreted as other countries’ decision-making being predominantly influenced by their domestic factors and international political elements, having little association with the compensations paid by the Japanese government to its domestic Fisheries Association, reflecting the relative independence of international actors from Japan’s internal compensatory dynamics. Fig 8c illustrates that with the increase in CLF, the evolutionary stable point of the Japanese Fisheries Association, γ4(0, 0, 1)T, also remains constant, but the evolutionary speed progressively slows down, meaning the Japanese Fisheries Association reaches the equilibrium point more slowly. This could be attributed to the Japanese Fisheries Association needing more time to assess and adapt to the potential impacts of receiving higher litigation compensations. Despite the potential influx of compensatory funds, the association must weigh the long-term considerations for the fishermen and the future development of fisheries livelihoods, highlighting the intricate balance between immediate financial relief and long-term sustainable development. Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 8. https://doi.org/10.1371/journal.pone.0317419.g008 In conclusion, the varying levels of litigation compensation to the Fisheries Association not only influence the strategic deliberation and evolutionary pace of the Japanese government and the Japanese Fisheries Association but also underscore the nuanced interplay between economic compensations and sustainable strategic considerations in the context of environmental disputes. Aid received by Japan from other countries CHJ. We explore the implications of the aid received by Japan from other countries, denoted as CHJ, on the evolutionary outcomes and trajectories of the three game entities. We assume CHJ = 5, 9, 13, and 17, while maintaining other parameters as shown in Table 5. The specific impacts of varying CHJ are illustrated in Fig 9. Fig 9a demonstrates that as CHJ increases, the evolutionary stable point of the Japanese government, γ4(0, 0, 1)T, remains constant, but the evolutionary speed progressively decelerates, meaning the Japanese government reaches the equilibrium point more slowly. This suggests that the aid from other countries may induce the Japanese government to deliberate more meticulously on its strategies, with the slower evolutionary speed aimed at mitigating additional costs induced by uncertainties. Enhanced international aid from other countries acts as a deterrent to Japan’s discharge strategy and prolongs the time to reach the evolutionary stable point, serving as an incentivizing measure for Japan to refrain from discharging. As depicted in Fig 9b, with the increase in CHJ, the evolutionary stable point of other countries, γ4(0, 0, 1)T, also remains constant, but the evolutionary speed accelerates. This could be interpreted as the larger the international aid provided by other countries, the slower the progression of the crisis, and when improvements are observed, other countries also accelerate their evolution, reflecting the adaptive nature of international actors in response to the evolving crisis dynamics. Fig 9c illustrates that with the increase in CHJ, the evolutionary stable point of the Japanese Fisheries Association, γ4(0, 0, 1)T, also remains constant, but the evolutionary speed accelerates, meaning the Japanese Fisheries Association reaches the equilibrium point more swiftly. This could be attributed to the Japanese Fisheries Association perceiving the aid from other countries as an incentive that could potentially foster the development of Japanese fisheries and secure the income of fishermen, highlighting the perceived positive impact of international aid on domestic industries. Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 9. https://doi.org/10.1371/journal.pone.0317419.g009 In conclusion, the varying levels of international aid received by Japan not only influence the strategic deliberation and evolutionary pace of the involved entities but also underscore the intricate interplay between international cooperation and strategic environmental decision-making in the context of international environmental disputes. Initial strategy selection probabilities of Japan, other countries, and the japanese fisheries Association: x0, y0, z0. Assuming (x0, y0, z0) = (0.5, 0.5, 0.5), (0.8, 0.1, 0.1), (0.2, 0.7, 0.1), (0.7, 0.2, 0.1), while maintaining other parameters as shown in Table 5, we can observe the evolutionary outcomes and trajectories of the three game entities under varying initial strategy selection probabilities, as depicted in Fig 10. From the figure, it is evident that when initial strategy selection probabilities vary, the evolutionary stable point γ4(0, 0, 1)T remains constant, indicating that it does not affect the evolutionary outcomes but does influence the evolutionary trajectories. For Japan, a stronger initial strategy selection probability, x0, will delay Japan’s discharge strategy and decelerate its evolutionary speed, as illustrated in Fig 10a. For other countries, the evolutionary stable point γ4(0, 0, 1)T also remains constant when initial strategy selection probabilities vary, signifying that it does not impact the evolutionary outcomes but does have certain effects on the evolutionary trajectories. A stronger initial strategy selection probability, y0, will delay the opposition strategy of the Fisheries Association and slow down its evolutionary speed, as shown in Fig 10b. Lastly, for the Fisheries Association, the evolutionary stable point γ4(0, 0, 1)T still remains constant when initial strategy selection probabilities change, suggesting that it does not alter the evolutionary outcomes but does affect the evolutionary trajectories. A stronger initial strategy selection probability, z0, will delay the opposition strategy of the Fisheries Association and decelerate its evolutionary speed, as depicted in Fig 10c. Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 10. https://doi.org/10.1371/journal.pone.0317419.g010 Discussion and conclusion This study employed evolutionary game theory to explore the strategic interactions and evolutionary dynamics among Japan, other countries, and the Japanese Fisheries Association in the context of nuclear wastewater discharge. The need to understand these potential decision-making processes under varying environmental and political pressures is critical, and our research aimed to address this by delineating three evolutionarily stable strategies (ESSs) and their corresponding stability conditions. These findings reveal how different scenarios could potentially stabilize based on the strategic decisions of the involved parties. Additionally, the study highlighted the profound impact of several parameters, including Japan’s international image and economic repercussions such as litigation and export tax reductions, alongside environmental costs on the strategic behaviors of the stakeholders. The existing literature often focuses on isolated aspects of environmental policy or employs static models that do not account for dynamic strategic interactions. By integrating evolutionary game theory, our study offers a dynamic perspective that encapsulates the complex interplay of economic, environmental, and social factors influencing policy decisions related to nuclear wastewater management. This provides a novel framework for predicting and analyzing the outcomes of such contentious environmental issues, enriching the theoretical foundations of environmental policy studies and offering a new lens through which to view the strategic decisions of countries and organizations dealing with global environmental crises. The findings suggest that policymakers need to consider a wider array of factors, including potential international backlash and internal economic pressures, when deciding on environmental policies. We recommend that Japan reassess its approach to nuclear wastewater management, considering not only the economic but also the long-term environmental and social impacts. Furthermore, other countries and the Fisheries Association should strengthen their collaborative efforts to apply pressure and advocate for sustainable environmental practices, ensuring that their strategies are informed by both scientific research and public opinion. In conclusion, our research illuminates the complexities involved in managing environmental crises through the lens of evolutionary game theory, providing critical insights into the potential evolutionary trajectories of policy decisions regarding nuclear wastewater discharge in Japan. The identified ESSs and the significant role of key parameters offer policymakers and stakeholders a grounded basis for reevaluating their strategies to ensure they align with both domestic needs and international expectations. This framework not only aids in understanding the current situation but also assists in forecasting and preparing for future challenges in environmental management. Our recommendations aim to foster more informed, responsible, and cooperative environmental governance, essential for addressing the multifaceted challenges posed by nuclear wastewater and other similar global environmental issues. Appendix Jacobian matrices for different points (26) (27) (28) (29) (30) (31) (32) (33) Eigenvalues for different points (34) (35) (36) (37) (38) (39) (40) (41) Jacobian matrices for different points (26) (27) (28) (29) (30) (31) (32) (33) Eigenvalues for different points (34) (35) (36) (37) (38) (39) (40) (41) Supporting information S1 File. https://doi.org/10.1371/journal.pone.0317419.s001 (DOCX) TI - An evolutionary game theory analysis on the environmental impact of discharging Fukushima’s nuclear wastewater: International stakeholders and strategic dynamics JO - PLoS ONE DO - 10.1371/journal.pone.0317419 DA - 2025-01-27 UR - https://www.deepdyve.com/lp/public-library-of-science-plos-journal/an-evolutionary-game-theory-analysis-on-the-environmental-impact-of-oqxC0lIYEp SP - e0317419 VL - 20 IS - 1 DP - DeepDyve ER -