NGC 3227 Blinks: An Extreme Changing-Look Event in a Seyfert Galaxy

Introduction

Changing-look active galactic nuclei (CLAGNs) are AGNs that exhibit dramatic changes in their spectral characteristics (C. Ricci & B. Trakhtenbrot [40]), such as the appearance or disappearance of broad emission lines or changing X-ray spectral slopes. Changing-look events provide direct probes of changes in the inner accretion flow, the formation and destruction of circum-nuclear structures, and the gaseous environment surrounding AGNs.

These transitions are generally attributed to two distinct physical mechanisms (C. Ricci & B. Trakhtenbrot [40]). In the first case, changes in the intrinsic accretion flow onto the supermassive black hole (SMBH) lead to the strengthening or weakening of the ionizing continuum, resulting in the appearance or disappearance of broad emission lines ([42, 18]). In the second, known as changing obscuration, variations in the line-of-sight column density modulate the observed emission, producing pronounced spectral changes ([41, 48]). Such obscuration-driven events display X-ray spectral hardening, often accompanied by suppressed soft X-ray emission and, in some cases, blue-shifted broad absorption-line troughs in the ultraviolet ([26, 23]).

While most optical CLAGN events are attributed to intrinsic accretion-rate changes ([40, 22 M. J. Temple et al.]

Email: [email protected] 2022), direct evidence that obscuration alone can drive large-amplitude optical continuum and broad emission-line variability has been reported in only a few cases [52, 35, 53].

NGC 3227 is a nearby, well-studied Seyfert galaxy ($z = 0.003859$) at a distance of 23.7 Mpc [2], with a black hole mass of $M_{\mathrm{BH}} \sim1 \times10^{7}\ M_{\odot}$ from reverberation mapping [4, 2] and dynamical modeling [2, 1]. It is in an early-stage interaction with its companion, NGC 3226. This interaction has resulted in substantial internal reddening in the host galaxy of NGC 3227, as revealed by H i radio mapping [34] and by the reddened optical and ultraviolet spectral shape [6, 32]. We account for this reddening in our spectral decomposition (§ 3.1).

NGC 3227 has a long history of X-ray variability and obscuration events, with episodes of X-ray spectral hardening reported as early as 2000 using RXTE observations [28] (A. G. Markowitz et al. 2014). Subsequent X-ray monitoring has revealed additional hardening events [3, 49, 51] (A. G. Markowitz et al. 2014), suggesting that such obscuration events recur. However, none of these episodes was accompanied by high-cadence optical spectroscopic monitoring.

In November 2019, another X-ray hardening event in NGC 3227 was detected by Swift [32], triggering intensive multi-wavelength follow-up observations with HST, XMM-Newton, and NuSTAR. During this campaign, the nucleus exhibited a pronounced decline in flux from the UV to the X-ray bands, with the suppression strongest at soft X-ray energies, consistent with the appearance of new obscuring gas along the line of sight [32]. Between the two target-of-opportunity epochs, the obscurer column density increased from $N_{\mathrm{H}} \sim10^{22}\ \mathrm{cm}^{-2}$ in mid-November to $\sim10^{23}\ \mathrm{cm}^{-2}$ by early December, tracking the deepening X-ray suppression [31].

In this Letter, we present contemporaneous high-cadence optical spectroscopic monitoring of NGC 3227 that covers this event. We describe the observations and data reduction in § 2, the spectral analysis in § 3, and our results in § 4, comprising the changing-look event itself, a damped random walk null test, the reverberation analysis, and photoionization modeling of the obscurer. We discuss the implied obscuration geometry in § 5 and summarize in § 6.

Observations and Data Reduction

WIRO spectroscopic monitoring

We monitored NGC 3227 between 2019 November 14 and 2020 May 29 with the Wyoming Infrared Observatory (WIRO) 2.3 m telescope using the Long Slit Spectrograph in remote-operation mode [13], following the procedures established for the Monitoring AGNs with H$\beta$ Asymmetry (MAHA) campaign (P. Du et al. 2018). The 900 line mm$^{-1}$ grating provided a dispersion of 1.49 Å pixel$^{-1}$ over 4000–7000 Å. Observations were carried out with a 5″ wide slit fixed in the north–south orientation. Each epoch typically consisted of three 5 minute exposures. At least one spectrophotometric standard star (usually Feige 34, occasionally G191B2B) was observed per night, along with CuAr lamps for wavelength calibration and dome flats from a uniformly illuminated white card.

Data were reduced following the standard MAHA pipeline (P. Du et al. 2018), which employs IRAF v2.16 [46, 47] for the basic reduction steps. One-dimensional spectra were extracted with apertures of \pm6\farcs8 and background windows spanning 7\farcs6–15\farcs1 on both sides. Flux calibration and extinction correction were applied using spectrophotometric standards, but additional corrections were required to account for variable atmospheric transparency and telescope tracking errors.

We adopted the [O III] λ5007 narrow-line emission as an internal flux calibrator, since it originates from an extended narrow-line region that is stable on timescales of years [37]. Variable seeing and tracking errors between epochs introduce small shifts in the [O III] line centroid, which would propagate into the flux calibration. To mitigate this, all spectra were convolved to a common spectral resolution before calibration. This reduces the effective spectral resolution but does not affect the broad-line measurements that are the focus of this work. Relative spectrophotometric calibration was then performed using the algorithm of E. van Groningen & I. Wanders (1992), with minor modifications (P. Du et al. 2018). We measured the [O III] fluxes over 5005–5050 Å from spectra obtained under photometric conditions. Final flux-calibrated spectra were obtained by noise-weighted averaging of exposures from the same night.

Archival ZTF photometry

To supplement the continuum light curve and increase the sampling cadence, we used archival $g$-band photometry of NGC 3227 from the Zwicky Transient Facility (ZTF; [21]). The data were extracted from a region centered on the nucleus, with only measurements having catflag = 0 retained for quality control. All data were drawn from object IDs 571109400017782, 572112300010962, and 1617102200002807, covering the heliocentric Julian date (HJD) range 2458800 to 2459010.

The inter-calibration between the ZTF $g$-band light curve and the WIRO spectroscopic photometry at 5100 \AA{} was performed using PyCALI (Y.-R. Li et al. 2014). Before calibration, the ZTF observations were rebinned by averaging all measurements within a single day, and three binned epochs exceeding $3\sigma$ from the reconstructed light curve were removed to improve data quality.

Spectral Analysis

Host decomposition and reddening correction

The reduced spectra of NGC 3227 were corrected for Milky Way extinction using the G23 extinction curve of K. D. Gordon et al. [16] with $R_V = 3.1$ and $E(B-V) = 0.02$ mag, then shifted to the rest frame. The AGN and host galaxy contributions were separated using spectral decomposition with the galaxy templates of A. L. Kinney et al. [27]. The final host spectrum was constructed as a linear combination of Sa, Sb, and bulge templates with scaling factors of $1/10.9$, $1/8.3$, and $1/100$, respectively. These factors were fixed by two joint requirements: (i) the residual AGN continuum remains non-negative in the lowest-flux state, and (ii) the host stellar absorption features—in particular the Na i D doublet near 5900 \AA{}—are removed in the residual AGN spectra. The host contribution is dominated by the Sa and Sb components, with a negligible bulge contribution. This host spectrum, assumed constant across epochs, was subtracted from all spectra.

NGC 3227 is known to suffer significant internal reddening [6, 34]. Following M. Mehdipour et al. [32], we applied the SMC Bar extinction curve of K. D. Gordon et al. [17] with $R_V = 4.0$ and $E(B-V) = 0.45$ mag, which they found to best characterize the dust extinction in the nuclear region. The AGN-only spectra were dereddened using this extinction curve to recover the intrinsic AGN continuum and emission-line fluxes.

Mean and rms spectra

We constructed mean and root-mean-square (rms) spectra from the reddening-corrected AGN spectra to characterize the average emission profile and its variability. Integration windows were then defined for the $\mathrm{H}\beta$ emission line and the [O iii] $\lambda\lambda4959$, 5007 doublet, with boundaries set from the extent of the line emission in the mean and rms spectra. The $\mathrm{H}\beta$ profile was measured over 4800–4929 \AA{}, the [O iii] doublet over 4929–5040 \AA{}, and the 5100 \AA{} continuum over 5075–5125 \AA{}. The minimal [O iii] in the rms spectrum confirms the stability of the flux calibration. The mean and rms spectra, together with the wavelength windows defined above, are shown in Figure 1.

Recovery of the broad H component

To remove the local AGN continuum around the $\mathrm{H}\beta$ line, we defined continuum windows on either side of the $\mathrm{H}\beta$+[O iii] complex, specifically 4750–4790 \AA{} and 5075–5125 \AA{}. These windows were chosen to coincide with line-free regions in the mean and rms spectra, as identified by visual inspection. The local continuum was then modeled by linear interpolation between the two windows, and subtracted to isolate the emission-line flux. The continuum-subtracted spectra showing the broad $\mathrm{H}\beta$ variability and the stable [O iii] doublet across all monitoring epochs are shown in Figure 2.

!Continuum-subtracted Hβ and [O III] emission-line profiles bbox=0.13,0.49,0.46,0.85

The broad $\mathrm{H}\beta$ component was deblended by subtracting the narrow-line contribution from the continuum-subtracted $\mathrm{H}\beta$ emission region. The narrow $\mathrm{H}\beta$ flux and profile were assumed to be stable during the monitoring campaign, as indicated by the constancy of the [O iii] emission line. To determine the narrow $\mathrm{H}\beta$ profile, we used a low-state spectrum (2019 December 19) to minimize contamination from the broad $\mathrm{H}\beta$ component and residual continuum. Emission-line profiles were modeled using a Markov chain Monte Carlo (MCMC) approach to fully explore the parameter space of the broad and narrow components. The fitting was restricted to a narrow wavelength interval of 4840–4880 \AA{}. This restriction is required because the broad $\mathrm{H}\beta$ emission exhibits complex substructure in the lowest state that cannot be adequately represented by a single Gaussian; fitting the full profile would otherwise bias the decomposition by allowing the broad component to dominate the $\chi^2$ minimization. The posterior distributions for the narrow $\mathrm{H}\beta$ component yield amplitude $A = 2.75^{+0.21}_{-0.24} \times10^{-14}$ erg cm$^{-2}$ s$^{-1}$ \AA$^{-1}$, mean $\mu= 4860.62^{+0.24}_{-0.23}$ \AA{}, and width $\sigma= 8.74^{+0.69}_{-0.62}$ \AA{}. The broad $\mathrm{H}\beta$ profiles were then recovered by subtracting the best-fitting narrow component from the observed spectra. The resulting narrow-line flux ratio of [O iii] $\lambda5007/\mathrm{H}\beta\approx9$ is consistent with the typical value of $\sim10$ observed among AGNs [25, 43], supporting the robustness of the decomposition.

Light curves

Given the complex structure of the $\mathrm{H}\beta$ profile, the broad $\mathrm{H}\beta$ flux in the WIRO spectra was measured by direct integration over the recovered broad-line spectrum within the H$\beta$ window. The 5100 Å continuum flux was estimated as the mean flux density within the 5075–5125 Å continuum window. The total uncertainties consist of two components: Poisson noise from the observations and systematic uncertainties. Poisson noise was propagated throughout the data analysis. Systematic uncertainties were estimated using two complementary approaches, combined in quadrature, following standard MAHA practices ([4] (P. Du et al. 2018)). The first quantifies within-night variations: for epochs with multiple exposures, the scatter among individual flux measurements captures the effect of changing atmospheric conditions and tracking errors. The second captures short-term variability across epochs: the light curve is smoothed with a five-point median filter ([8, 9]), and the standard deviation of the residuals provides an estimate of rapid intrinsic fluctuations. The resulting broad H$\beta$ and 5100 Å continuum light curves are shown in Figure 3.

Results

An extreme optical changing-look event

Our optical spectroscopic monitoring with WIRO, which overlaps with the 2019 November multi-wavelength campaign, shows a contemporaneous decline in optical flux that evolves on weekly timescales (Figure 3). The AGN optical continuum at 5100 Å and the broad H$\beta$ light curve show substantial declines of approximately 70% and 80%, respectively. The broad H$\beta$ profile nearly disappears in the low state (2019 December 19) and recovers by the high state (2020 January 31), while the narrow H$\beta$ emission remains unchanged throughout. The Swift/XRT hardening began before the first WIRO epoch and continued through the early optical drop. During the period of joint coverage, both the optical continuum and the Swift/XRT count rates declined, though the X-ray sampling is sparse and ends prior to the optical minimum. The overlapping X-ray hardening supports a physical association with the optical decline, but we cannot conclusively demonstrate that both behaviors were produced by the exact same obscuring structure.

Given this rapid decline in continuum flux, stochastic accretion variability alone is an unlikely explanation: only ∼0.2% of damped random walk simulations calibrated to NGC 3227’s normal-state variability reproduce a decline this steep (§ 4.3). Combined with the X-ray hardening at the onset of the optical decline, we argue that this event is driven by obscuration. Reverberation mapping through the dim state (§ 4.2) recovers a robust 6–8 day lag between the continuum and broad Hβ, showing that the broad-line region continues to reverberate and placing the obscuring material interior to the broad-line region.

Reverberation lags

We divided the light curves into two temporal segments. The obscuration phase (2019 November 14 to 2020 January 31) captures the sharp decline in Hβ flux, following an X-ray hardening event, and its subsequent recovery. The normal phase (2020 February 2 to 2020 May 29) reflects typical AGN stochastic variability.

We employed the interpolated cross-correlation function (ICCF; C. M. Gaskell & L. S. Sparke [14]), implemented with the PyCCF package (M. Sun et al. [44]), to estimate the time lag between the ionizing continuum and the Hβ emission line, thereby inferring the broad-line region (BLR) size via the light-travel distance. To complement this classical approach, we also applied Multiple and Inhomogeneous Component Analysis (MICA; Y.-R. Li et al. (2016)), a Bayesian reverberation-mapping framework. MICA simultaneously models the stochastic continuum and the line response using a flexible transfer-function parameterization.

To improve coverage during the obscuration phase, we augmented the 5100 Å continuum with inter-calibrated ZTF data, as sparse sampling degrades both ICCF and MICA results more than the modest systematics introduced by inter-calibration. In contrast, ZTF data were excluded from the normal phase to preserve the homogeneity and precision of our primary dataset, though including the inter-calibrated ZTF data in the normal phase yields an essentially identical lag. ICCF and MICA analyses were performed independently for each phase, enabling a direct comparison of the delay behavior.

During the obscuration phase (Figure 4), the ICCF centroid distribution gives a lag of $6.13^{+1.44}_{-0.99}$ days, while MICA yields a broader posterior at $7.72^{+2.09}_{-2.10}$ days. During the normal phase (Figure 5), ICCF analysis yields a lag of $4.99^{+1.53}_{-1.43}$ days and MICA gives $5.38^{+0.87}_{-0.88}$ days; both methods converge on a consistent lag of $\sim5$ days, indicating a relatively stable BLR response under unobscured conditions and consistent with previous reverberation mapping results for NGC 3227 ([2, 4]).

Robustness of the obscuration-phase lag

To assess whether the deep-valley shape and reduced cadence of the obscuration phase could bias the recovered lag, we performed mock recovery tests using ICCF. We treated the MICA-reconstructed mean continuum light curve as the deterministic “true” continuum, convolved it with a Gaussian transfer function of standard deviation $\sigma\approx1.67$ days to generate an H$\beta$ light curve, interpolated both at the actual observation epochs, and assigned per-epoch uncertainties by rescaling the observed flux errors so that each mock point preserves the fractional uncertainty of the corresponding real measurement. We tested two input lags: $\tau_0 = 0$ (null) and $\tau_0 = 6$ days (recovery). ICCF recovered both inputs within their quoted uncertainties ($\tau= 0.11^{+0.94}_{-0.53}$ d for $\tau_0 = 0$ and $\tau= 5.50^{+0.98}_{-0.97}$ d for $\tau_0 = 6$ d), demonstrating that the obscuration-phase lag is not a spurious artifact of the sampling. Drawing pairs from the observed and mock zero-lag cross-correlation centroid distributions, we find $p(\tau_{\mathrm{null}} \ge\tau_{\mathrm{obs}}) \approx6.8 \times10^{-5}$, equivalent to a $3.8\sigma$ rejection of the null hypothesis $\tau= 0$; the observed lag distribution is therefore inconsistent with arising from a zero-lag light curve under the actual sampling and noise.

Damped random walk null test

To assess whether the observed $\sim70\%$ optical continuum decline could arise from stochastic accretion fluctuations rather than an external obscuration event, we performed a null-hypothesis test using damped random walk (DRW) simulations.

We used five seasons of MAHA optical monitoring of NGC 3227 obtained outside the obscuration episode: the 2017 season ([4]) and four subsequent seasons spanning 2019–2023 (Zhang et al. 2027, MAHA VI, in preparation). A monotonically declining segment of the 2023 season was excluded to ensure that the DRW parameters reflect stochastic variability. For each season, we independently inferred the DRW parameters, $\sigma$ (variability amplitude) and $\tau$ (characteristic timescale), and combined the posterior samples to determine representative median values. We obtained $\tau= 9.66^{+8.89}_{-3.69}$ days and $\sigma= 3.63^{+1.87}_{-1.67} \times10^{-15}$ erg cm$^{-2}$ s$^{-1}$ Å$^{-1}$.

Using these median parameters, we generated 100,000 mock light curves drawn from a DRW Gaussian process, adopting the same temporal sampling and photometric uncertainties as during the monitoring campaign. For each realization, we searched for events that exhibited a $\ge70\%$ decline in continuum flux on timescales shorter than 40 days. Such events occur in only 0.2% of the simulated light curves, and only a similarly small fraction of light curves would recover to a high state on the same timescale. This low occurrence rate indicates that the observed continuum drop and subsequent rise is very unlikely to arise from stationary stochastic accretion variability alone.

Photoionization modeling of the obscurer

The host-galaxy subtracted rest-frame 5100 Å continuum flux drops from $2.24 \times10^{-14}$ to $6.68 \times10^{-15}$ erg cm$^{-2}$ s$^{-1}$ Å$^{-1}$ between the high and low states, a decline of $\sim70\%$. To interpret this attenuation, we performed photoionization modeling with Cloudy C25 ([19]) to identify the obscurer column density required to reproduce the observed low-state 5100 Å flux. The ionizing source was modeled with the intrinsic AGN spectral energy distribution of NGC 5548 ([33]), normalized to the high-state luminosity $\lambda L_{5100\,\text{\AA}} \approx7.76 \times10^{42}$ erg s$^{-1}$. The obscurer was modeled as a plane-parallel slab of solar-abundance gas.

The exact location and density of the obscurer are not directly constrained by our data; therefore, we choose sampling ranges that are physically plausible for an obscurer interior to the 5100 Å continuum-emitting region. Following the disk size–luminosity relation of C. J. Burke et al. [5], the characteristic 2500 Å emitting radius of NGC 3227 is $\sim0.09$ light-days, and the empirical continuum reverberation scaling $\tau\propto\lambda^{4/3}$ ([11, 12]) extrapolates this to $\sim$0.23 light-days at 5100 Å. We therefore sampled obscurer inner radii of $r = 0.1$–$0.2$ light-days, just interior to the 5100 Å emitting region, and hydrogen number densities $\log(n_{\mathrm{H}}/\mathrm{cm}^{-3}) = 9$–$10$, typical of gas at BLR scales ([36]).

We performed grid photoionization modeling over this parameter space, with column densities spanning $\log(N_{\mathrm{H}}/\mathrm{cm}^{-2}) = 23$–$24.5$. Across this grid the ionization parameter spans $\log U \approx3$–$4.5$, leaving the obscurer fully ionized in all cases. The transmitted 5100 Å flux therefore depends only on $N_{\mathrm{H}}$, with negligible sensitivity to $n_{\mathrm{H}}$ or $r$. The Cloudy results coincide with the analytic Thomson prediction $T = \exp(-1.17 N_{\mathrm{H}} \sigma_{\mathrm{T}})$ for fully ionized hydrogen + helium gas at solar abundance, confirming that bound-free and line opacities are negligible at optical wavelengths and that electron scattering accounts for the entire attenuation. Reproducing the observed flux decrease requires a Compton-thick column density of $N_{\mathrm{H}} \simeq1.5 \times10^{24}$ cm$^{-2}$. Across the explored Cloudy grid, energy is conserved to within $\sim$3%; at the best-fit column the obscurer’s diffuse outward continuum at 5100 Å remains below $\sim4 \times10^{-4}$ of the transmitted flux, confirming that the fully ionized cloud does not self-illuminate the optical band probed in this work.

mitted flux confirming that the fully ionized cloud does not self-illuminate the optical band probed in this work.

Discussion

The detection of a coherent reverberation lag between the continuum and broad H$\beta$ during the obscuration phase localizes the obscuring material: the H$\beta$ decline is the BLR’s light-travel response to an attenuated ionizing continuum, not a direct extinction of the line-emitting gas along our sightline, so the BLR continues to reprocess the continuum it receives. This excludes obscuration exterior to the BLR. Such an obscurer dims H$\beta$ only by attenuating the line emission, so suppressing the broad H$\beta$ flux by $\sim$80% requires it to cover most of the spatially extended BLR — a structure comparable in scale to the BLR itself. If the obscurer covers the continuum source and BLR simultaneously (for example, an in-situ rise in opacity along the sightline), the two decline together, producing correlated, zero-lag drops.

If instead the obscurer transits the system, it occults the continuum source and the surrounding BLR at different times and smears the H$\beta$ decline over the crossing — in neither case reproducing the coherent $\sim$6–8 day delay we observe. A well-defined positive lag, in which the entire H$\beta$ light curve follows the continuum after a fixed delay, is instead the signature of light-travel reprocessing and requires the obscurer to lie interior to the BLR, where it attenuates the continuum before it reaches the line-emitting gas. This inner location is consistent with the X-ray analysis of the 2019 event by J. Mao et al. [31], and with this system’s earlier X-ray occultations, whose absorbing clouds were also placed within or at the inner edge of the BLR (T. J. Turner et al. [49]; T. Beuchert et al. [3]).

Photoionization modeling further shows that reproducing the observed attenuation of the 5100 $\mathrm{\AA}$ continuum requires a fully ionized Compton-thick column density of $N_{\mathrm{H}} \sim1.5 \times10^{24}\ \mathrm{cm}^{-2}$ along the line of sight. Such a large column density is qualitatively consistent with the pronounced soft X-ray suppression and X-ray spectral hardening observed at the onset of this event, although the X-ray-derived column density is one to two orders of magnitude lower. This difference likely reflects both temporal and geometric factors: the X-ray column was likely still evolving when last measured, increasing from $\sim10^{22}$ to $\sim10^{23}\ \mathrm{cm}^{-2}$ over the weeks before the optical minimum (J. Mao et al. [31]), and no X-ray observation samples the epoch of deepest optical suppression (2019 December 19, both later and deeper than the final X-ray epoch), precluding a direct comparison at the same evolutionary stage. Furthermore, the X-ray column density traces a narrow sightline to the compact corona, whereas the optical continuum sees a sightline to the continuum source at a larger radius, where the obscurer’s column density may differ.

To illustrate the observational results, we present a proposed obscuration geometry in Figure 6, which is similar to that previously offered for Mrk 817 (E. Kara et al. [24]), where a comparable inner obscurer left the broad UV lines decoupled from the observed continuum (a “holiday”; Y. Homayouni et al. [20]). In NGC 3227 the covering instead suppresses the continuum reaching the H$\beta$-emitting BLR, producing a full H$\beta$ changing-look event. Unlike the Mrk 817 holiday, however, the persistent $\sim6$–8 day lag shows that H$\beta$ still reverberates in response to the attenuated continuum. The physical nature of the obscurer — whether a disk wind or another circumnuclear structure — cannot be determined from the present data and remains an open question. We note, however, that a time-dependent, line-driven accretion-disk wind is one plausible candidate: simulations show that such winds can become Compton-thick along equatorial sightlines (D. Proga et al. [39]; D. Proga & T. R. Kallman [38]).

The event in NGC 3227 is currently unique in exhibiting X-ray obscuration accompanied by an optical changing-look transition, in which the broad emission lines nearly vanished alongside the continuum. In contrast, most X-ray changing-obscuration events do not show a clear correlation between the X-ray and UV/optical changing-look behavior. Certain X-ray changing-obscuration events (e.g., Mrk 817 and NGC 5548) instead exhibit (a) blueshifted absorption-line troughs rather than the near-complete continuum and emission-line suppression observed here, and (b) a BLR “holiday,” in which the broad lines decouple from the ionizing-continuum variability [15, 7, 24, 20]. Such holidays may depend on the obscurer’s column density ($N_{\mathrm{H}}$): if the material is transparent to specific bands of the spectral energy distribution, only partial suppression occurs. Because Compton-thick, near-total-coverage obscuration events are likely far rarer than ordinary obscuration, NGC 3227 remains the only unambiguous case identified to date.

Summary

We have presented high-cadence optical spectroscopic monitoring of an extreme changing-look event in NGC 3227 and shown that: (1) the 5100 Å continuum and broad H$\beta$ flux fell by $\sim70\%$ and $\sim80\%$ and recovered within two months, a decline that damped random walk simulations disfavor as stochastic accretion variability ($p = 0.002$) and that coincides with Swift X-ray hardening; (2) a persistent 6–8 day continuum-to-H$\beta$ lag through the dim state shows the broad-line region still reverberates, placing the obscurer interior to the BLR; and (3) reproducing the optical suppression requires a fully ionized, Compton-thick column, $N_{\mathrm{H}} \simeq1.5 \times10^{24}\ \mathrm{cm}^{-2}$. Together these establish that compact, high-column-density obscuration interior to the BLR can drive an extreme optical changing-look event on week timescales. Such events are hard to characterize from either band alone — X-ray and UV data fix the absorber’s column and ionization state, high-cadence optical spectroscopy localizes it through reverberation mapping — so simultaneous coverage is needed to constrain the obscurer’s geometry. Wide-field time-domain surveys with the Rubin Observatory should uncover many more and motivate the coordinated follow-up such events require.

Acknowledgments

We thank the WIRO engineering staff James Weger and Conrad Vogel for their dedicated technical support in maintaining the telescope and instrumentation that made these observations possible. Based on observations obtained with the Samuel Oschin Telescope 48 inch and the 60 inch Telescope at the Palomar Observatory as part of the Zwicky Transient Facility project. ZTF is supported by the National Science Foundation under Grants No. AST-1440341 and AST-2034437 and a collaboration including current partners Caltech, IPAC, the Weizmann Institute for Science, the Oskar Klein Center at Stockholm University, the University of Maryland, Deutsches Elektronen-Synchrotron and Humboldt University, the TANGO Consortium of Taiwan, the University of Wisconsin at Milwaukee, Trinity College Dublin, Lawrence Livermore National Laboratories, IN2P3, University of Warwick, Ruhr University Bochum, Northwestern University and former partners the University of Washington, Los Alamos National Laboratories, and Lawrence Berkeley National Laboratories. Operations are conducted by COO, IPAC, and UW.

This work was supported by the U.S. National Science Foundation through grant AST-2511198. P.D. acknowledges financial support from the National Key R&D Program of China (2023YFA1607903 and 2021YFA1600404), from NSFC grants NSFC-12521005, -12022301 and -11991051, the National Science and Technology Major Project (2024ZD03000303), and by the China Manned Space Project with no. CMS-CSST-2025-A07.

Facilities: WIRO (Wyoming Infrared Observatory), ZTF (Zwicky Transient Facility), Swift

Software: Cloudy [19], IRAF [46, 47], MICA (Y.-R. Li et al. 2016), PyCALI (Y.-R. Li et al. 2014), PyCCF [44]

Data Availability

The underlying data, light curves, spectra, and processing/analysis scripts used to produce the figures and results presented in this manuscript are available via DOI: 10.5281/zenodo.23191065.