skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Coupling agent-based simulation and spatial optimization models to understand spatially complex and co-evolutionary behavior of cocaine trafficking networks and counterdrug interdiction
Despite more than 40 years of counterdrug interdiction efforts in the Western Hemisphere, cocaine trafficking, or ‘narco-trafficking’, networks continue to evolve and increase their global reach. Counterdrug interdiction continues to fall short of performance targets due to the adaptability of narco-trafficking networks and spatially complex constraints on interdiction operations (e.g., resources, jurisdictional). Due to these dynamics, current modeling approaches offer limited strategic insights into time-varying, spatially optimal allocation of counterdrug interdiction assets. This study presents coupled agent-based and spatial optimization models to investigate the co-evolution of counterdrug interdiction deployment and narco-trafficking networks’ adaptive responses. Increased spatially optimized interdiction assets were found to increase seizure volumes. However, the value per seized shipment concurrently decreased and the number of active nodes increased or was unchanged. Narco-trafficking networks adaptively responded to increased interdiction pressure by spatially diversifying routes and dispersing shipment volumes. Thus, increased interdiction pressure had the unintended effect of expanding the spatial footprint of narcotrafficking networks. This coupled modeling approach enabled the study of narco-trafficking network evolution while being subjected to varying interdiction pressure as a spatially complex adaptive system. Capturing such co-evolution dynamics is essential for simulating traffickers’ realistic adaptive responses to a wide range of interdiction scenarios.  more » « less
Award ID(s):
1837698 2039975
PAR ID:
10358669
Author(s) / Creator(s):
; ; ; ; ;
Date Published:
Journal Name:
IISE Transactions
ISSN:
2472-5854
Page Range / eLocation ID:
1 to 23
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Counterdrug interdiction efforts designed to seize or disrupt cocaine shipments between South American source zones and US markets remain a core US “supply side” drug policy and national security strategy. However, despite a long history of US-led interdiction efforts in the Western Hemisphere, cocaine movements to the United States through Central America, or “narco-trafficking,” continue to rise. Here, we developed a spatially explicit agent-based model (ABM), called “NarcoLogic,” of narco-trafficker operational decision making in response to interdiction forces to investigate the root causes of interdiction ineffectiveness across space and time. The central premise tested was that spatial proliferation and resiliency of narco-trafficking are not a consequence of ineffective interdiction, but rather part and natural consequence of interdiction itself. Model development relied on multiple theoretical perspectives, empirical studies, media reports, and the authors’ own years of field research in the region. Parameterization and validation used the best available, authoritative data source for illicit cocaine flows. Despite inherently biased, unreliable, and/or incomplete data of a clandestine phenomenon, the model compellingly reproduced the “cat-and-mouse” dynamic between narco-traffickers and interdiction forces others have qualitatively described. The model produced qualitatively accurate and quantitatively realistic spatial and temporal patterns of cocaine trafficking in response to interdiction events. The NarcoLogic model offers a much-needed, evidence-based tool for the robust assessment of different drug policy scenarios, and their likely impact on trafficker behavior and the many collateral damages associated with the militarized war on drugs. 
    more » « less
  2. We consider a new class of multi-period network interdiction problems, where interdiction and restructuring decisions are decided upon before the network is operated and implemented throughout the time horizon.We discuss how we apply this new problem to disrupting domestic sex trafficking networks, and introduce a variant where a second cooperating attacker has the ability to interdict victims and prevent the recruitment of prospective victims. This problem is modeled as a bilevel mixed integer linear program (BMILP), and is solved using column-and-constraint generation with partial information. We also simplify the BMILP when all interdictions are implemented before the network is operated. Modeling-based augmentations are proposed to significantly improve the solution time in a majority of instances tested. We apply our method to synthetic domestic sex trafficking networks, and discuss policy implications from our model. In particular, we show how preventing the recruitment of prospective victims may be as essential to disrupting sex trafficking as interdicting existing participants. 
    more » « less
  3. Illicit Wildlife Trade (IWT) is a serious global crime that negatively impacts biodiversity, human health, national security, and economic development. Many flora and fauna are trafficked in different product forms. We investigate a network interdiction problem for wildlife trafficking and introduce a new model to tackle key challenges associated with IWT. Our model captures the interdiction problem faced by law enforcement impeding IWT on flight networks, though it can be extended to other types of transportation networks. We incorporate vital issues unique to IWT, including the need for training and difficulty recognizing illicit wildlife products, the impact of charismatic species and geopolitical differences, and the varying amounts of information and objectives traffickers may use when choosing transit routes. Additionally, we incorporate different detection probabilities at nodes and along arcs depending on law enforcement’s interdiction and training actions. We present solutions for several key IWT supply chains using realistic data from conservation research, seizure databases, and international reports. We compare our model to two benchmark models and highlight key features of the interdiction strategy. We discuss the implications of our models for combating IWT in practice and highlight critical areas of concern for stakeholders. 
    more » « less
  4. Despite significant progress, the treatment of estrogen receptor-positive (ER+) breast cancer remains clinically challenging due to reversible drug resistance and immune evasion. Drug resistance often arises as cells undergo a dynamic epithelial-to-mesenchymal transition (EMT), while elevated PD-L1 levels contribute to immune escape. While these phenotypic features can variably co-occur, the impact of co-occurrence on the availability of synergistic treatment strategies remains unknown. To investigate their interplay, we constructed an ER-EMT-PD-L1 gene regulatory network and simulated these networks as coupled ordinary differential equations with biologically informed parameters, to generate steady-state expression profiles. Our study revealed that the relevant overarching network generated antagonistic epithelial and mesenchymal modules, capable of producing monostable, bistable, and tristable dynamics. We further examined the link between phenotypes and immune evasion by quantifying average PD-L1 expression, and found that epithelial-sensitive states consistently exhibited low PD-L1. In contrast, hybrid- and mesenchymal-resistant states were associated with a non-linear, stepwise increase in PD-L1, highlighting a strong coupling between EMT, resistance, and immune evasion. Extending on these network-level insights, we further used a spatially explicit agent-based model seeded with GRN-derived phenotypes to probe tumor behavior under therapeutic pressure. Simulations revealed that sustained tumor expansion occurred only when resistance, motility, and immune evasion traits co-existed, and this requirement remained robust across GRN landscapes with differing stability. Plasticity and multistability increased the accessible phenotypic state-space and accelerated shifts toward high-fitness resistant states. We further identified combination therapies that significantly reduced phenotypic diversification and improved immune infiltration in silico. Taken together, our modeling work links regulatory dynamics with tumor-level adaptation and highlights strategies to reprogram resistant cell states toward sensitivity, which are difficult to infer from bulk or cross-sectional data alone. In addition, it provides a controllable in silico testbed to systematically evaluate candidate treatment combinations and their effects on tumor phenotypic transitions and spatial T cell access, thereby helping to prioritize experimental regimens for follow-up. 
    more » « less
  5. null (Ed.)
    Abstract. Important uncertainties remain in our understanding of the spatial andtemporal variability of atmospheric hydroxyl radical concentration ([OH]).Carbon-14-containing carbon monoxide (14CO) is a useful tracer that canhelp in the characterization of [OH] variability. Prior measurements ofatmospheric 14CO concentration ([14CO] are limited in both theirspatial and temporal extent, partly due to the very large air sample volumes that have been required for measurements (500–1000 L at standardtemperature and pressure, L STP) and the difficulty and expense associatedwith the collection, shipment, and processing of such samples. Here wepresent a new method that reduces the air sample volume requirement to≈90 L STP while allowing for [14CO] measurement uncertainties that are on par with or better than prior work (≈3 % or better, 1σ). The method also for the first time includes accurate characterization of the overall procedural [14CO] blank associated with individual samples, which is a key improvement over prior atmospheric 14CO work. The method was used to make measurements of [14CO] at the NOAA Mauna Loa Observatory, Hawaii, USA, between November 2017 and November 2018. The measurements show the expected [14CO] seasonal cycle (lowest in summer)and are in good agreement with prior [14CO] results from anotherlow-latitude site in the Northern Hemisphere. The lowest overall [14CO]uncertainties (2.1 %, 1σ) are achieved for samples that aredirectly accompanied by procedural blanks and whose mass is increased to≈50 µgC (micrograms of carbon) prior to the 14Cmeasurement via dilution with a high-CO 14C-depleted gas. 
    more » « less