Utrzymanie uczestników i profil badań terenowych
Podczas formalnych prac terenowych w grotach Longmen uzyskano łącznie 480 kwestionariuszy. Po weryfikacji kryteriów kwalifikowalności i kontroli jakości danych, do analizy zakwalifikowano 430 kwestionariuszy, co dało wskaźnik poprawnych odpowiedzi na poziomie 89,6%. Ostateczna liczebność poprawnej próby przekroczyła wcześniej określone minimalne wymaganie 40 poprawnych kwestionariuszy. Wykluczone przypadki obejmowały brak spełnienia kryteriów kwalifikowalności, nieprawidłowe wyniki testów uwagi, zbyt szybkie wypełnienie, odpowiedzi liniowe lub niemal niezmienne wzorce odpowiedzi, powtórzone odpowiedzi oraz inne problemy związane z pracami terenowymi. W przypadku, gdy jeden kwestionariusz spełniał więcej niż jeden warunek wykluczenia, jako główną kategorię wykluczenia przyjęto pierwszą zastosowaną regułę protokołu. Schemat selekcji od 480 surowych kwestionariuszy do 430 poprawnych przypadków przedstawiono na Rysunku 2.
Zachowane przypadki rozdzielono pomiędzy trzy planowane punkty badania ankietowego. Obszar wyjścia dla zwiedzających dostarczył 139 prawidłowych przypadków, publiczna strefa odpoczynku poza główną trasą zwiedzania dostarczyła 146 przypadków, a obszar obsługi odwiedzających w pobliżu punktu przesiadkowego dostarczył 145 przypadków. Osoby odwiedzające miejsce po raz pierwszy stanowiły 5,6% prawidłowych przypadków, a osoby powracające 4,4%. Uczestnicy zorganizowanych wycieczek grupowych stanowili 23,0%, co pozostało poniżej uprzednio określonej granicy monitorowania. Rozkład prac terenowych oraz profil odwiedzających przedstawiono w Tabeli 6.
Rozkład konstruktów i diagnostyka pomiarowa
Wyniki dla konstruktów wykazały wystarczający rozrzut do przeprowadzenia analizy. Średnie wyniki mieściły się w zakresie od 3,258 do 3,415, a odchylenia standardowe od 0,837 do 0,892. Żaden konstrukt nie wykazał silnego efektu sufitu w oparciu o wcześniej określone weryfikacje rozrzutu. Zmienna postrzeganego tłoku została zachowana w oryginalnym kierunku, zatem wyższy wynik wskazywał na silniejsze poczucie tłoku, a nie na bardziej pozytywne doświadczenia odwiedzających.
Współczynnik alfa Cronbacha wynosił od 0,69 dla postrzeganego tłoku do 0,812 dla deklarowanych odpowiedzialnych zachowań turystycznych. Wszystkie konstrukty spełniły albo dopuszczalny próg niezawodności, albo regułę niezawodności marginalnej. Postrzegany tłok uznano za marginalnie niezawodny, ponieważ jego wartość alfa była nieco poniżej 0,70. Został on zachowany, gdyż skorygowane korelacje element-całość spełniały próg protokołu, a sam konstrukt był niezbędny do reprezentacji warunków presji odwiedzających w obszarze dziedzictwa o dużym natężeniu ruchu.
Dziewięcioczynnikowa analiza czynnikowa konfirmacyjna potwierdziła planowaną strukturę pomiarową. Wskaźniki dopasowania modelu były akceptowalne: χ2 (398) = 712.4, CFI = 0.932, TLI = 0.918, RMSEA = 0.043 oraz SRMR = 0.052. Standaryzowane ładunki czynnikowe mieściły się w zakresie od 0.56 do 0.84. Wartości złożonej rzetelności mieściły się w zakresie od 0.702 do 0.842. Średnia wariancja wyekstrahowana mieściła się w zakresie od 0.453 do 0.575, przy czym najniższa wartość wystąpiła dla postrzeganego zatłoczenia.
Weryfikacja trafności dyskryminacyjnej nie wykazała poważnego nakładania się konstruktów. Maksymalna wartość HTMT wyniosła 0,642, co znajduje się poniżej wcześniej określonego progu weryfikacyjnego. Postawa wobec ochrony przyrody, chęć ochrony przyrody oraz samoopisane odpowiedzialne zachowania podróżnicze pozostały rozróżnialne w modelu pomiarowym. Pierwszy czynnik nierotowany wyjaśnił 18,7% całkowitej wariancji, co mieści się poniżej 40-procentowego progu przesiewowego dla błędów metodologicznych (common-method bias). Wynik ten potraktowano jedynie jako wynik przesiewowy, a nie jako dowód na brak błędów metodologicznych. Najsilniejsza korelacja między konstruktami wystąpiła pomiędzy postawą wobec ochrony przyrody a chęcią ochrony przyrody, r = 0,432, a żadna para konstruktów nie przekroczyła progu redundancji wynoszącego 0,70. Główne parametry diagnostyczne pomiaru podsumowano w Tabeli 7. Pełne wyniki niezawodności na poziomie pozycji, ładunki CFA , CR, AVE, HTMT oraz korelacje między konstruktami przedstawiono w Tabeli uzupełniającej 1, Tabeli uzupełniającej 2, oraz Tabeli uzupełniającej 3.
Wyniki modeli i kontrole czułości
Protokół pozwolił na uzyskanie interpretowalnych wyników regresji oraz modeli ścieżkowych bez wystąpienia poważnej współliniowości. Trzy modele regresji wykazały umiarkowaną moc wyjaśniającą. Model postawy wobec ochrony przyrody uzyskał R2 = 0,267 oraz skorygowane R2 = 0,24. Model gotowości do ochrony przyrody uzyskał R2 = 0,33 oraz skorygowane R2 = 0,292. Model deklarowanego odpowiedzialnego zachowania podczas podróży uzyskał R2 = 0,275 oraz skorygowane R2 = 0,28. Wartości te mieściły się w wcześniej określonym prawdopodobnym zakresie dla przekrojowych danych dotyczących zachowań odwiedzających.
Model ścieżkowy odpowiadał zaplanowanej sekwencji: predyktory związane z lokalizacją → postawa proekologiczna → gotowość do ochrony przyrody → deklarowane odpowiedzialne zachowania podczas podróżowania. Dopasowanie modelu było akceptowalne: CFI = 0.927, TLI = 0.904, RMSEA = 0.049 oraz SRMR = 0.041. Główne szacunki ścieżek były interpretowalne. Postawa proekologiczna była dodatnio skorelowana z gotowością do ochrony przyrody, β = 0.276. Postawa proekologiczna oraz gotowość do ochrony przyrody były dodatnio skorelowane z deklarowanymi odpowiedzialnymi zachowaniami podczas podróżowania, odpowiednio β = 0.232 i β = 0.271. Efekt pośredni od postawy proekologicznej do deklarowanych odpowiedzialnych zachowań podczas podróżowania poprzez gotowość do ochrony przyrody wyniósł 0.075, przy 95% przedziale ufności bootstrap 0.036–0.12. Pełne współczynniki regresji, szacunki ścieżek, efekty pośrednie, efekty całkowite oraz pliki z wynikami modelu znajdują się w Pliku uzupełniającym 3, Tabeli uzupełniającej 4 oraz Tabeli uzupełniającej 5.
Diagnostyka regresji nie wskazała na występowanie poważnej współliniowości. Wartości VIF mieściły się w zakresie od 1,012 do 1,42, czyli poniżej określonego progu 3,0. Diagnostyka przypadków wpływowych nie uzasadniła usunięcia żadnych przypadków z modeli głównych. Przypadki oznaczone za pomocą dystansu Cooka zostały zachowane w analizie głównej i ocenione w ramach analizy wrażliwości.
Analizy wrażliwości potwierdziły stabilność wyników protokołu. Żaden z prawidłowo wypełnionych kwestionariuszy nie przekroczył progu czasu wypełniania wynoszącego 1 80 s. Kierunki głównych współczynników pozostały niezmienne po usunięciu uczestników zorganizowanych wycieczek grupowych, zastosowaniu punktacji konstruktów opartej na medianie oraz wykluczeniu przypadków oznaczonych przy użyciu odległości Cooka. W ramach różnych specyfikacji analizy wrażliwości, standaryzowany współczynnik dla postawy wobec ochrony → chęci ochrony mieścił się w przedziale od 0,261 do 0,290, dla postawy wobec ochrony → deklarowanych odpowiedzialnych zachowań podczas podróżowania w przedziale od 0,198 do 0,260, a dla chęci ochrony → deklarowanych odpowiedzialnych zachowań podczas podróżowania w przedziale od 0,243 do 0,298. Pełne wyniki analizy wrażliwości przedstawiono w Tabeli uzupełniającej 6.
Podsumowanie wyników protokołu
Przepływ pracy pozwolił na uzyskanie poprawnej próbki analitycznej przekraczającej wcześniej określone minimalne wymagania dotyczące wielkości próby, z wypełnionymi kwestionariuszami rozłożonymi na wszystkie planowane punkty badawcze. Zachowana próbka obejmowała zróżnicowanie pod kątem pochodzenia odwiedzających oraz sposobu podróżowania, a proporcja zorganizowanych wycieczek grupowych pozostała w granicach monitoringu. Kontrole pomiarowe potwierdziły zasadność zastosowania wyników na poziomie konstruktu: większość konstruktów spełniła próg niezawodności, postrzegane zatłoczenie uznano za marginalnie niezawodne, dziewięcioczynnikowa CFA wykazała akceptowalne dopasowanie, a testy trafności różnicowej potwierdziły rozróżnienie między postawą proekologiczną, gotowością do ochrony przyrody a deklarowanym odpowiedzialnym zachowaniem podczas podróży.
Wyniki modelu posłużyły do wykazania, że protokół pozwala na wygenerowanie dowodów zachowań zwiedzających, które podlegają analizie i audytowi. Wyniki regresji i modelu ścieżkowego były interpretowalne, test efektów pośrednich był dostępny dla planowanego ciągu postawa–chęć–zachowanie, a testy wrażliwości nie zmieniły kierunku podstawowych powiązań. Wyniki te wspierają proceduralne zastosowanie tego przepływu pracy w ankietach wśród zwiedzających obiekty dziedzictwa w celu zminimalizowania zakłóceń, a nie służą jako bezpośrednie dowody obserwacyjne zachowań zwiedzających. Zanonimizowany zestaw danych, kwestionariusze, pomocnicze dokumenty badania, skrypt analizy, pełne wyniki modelu oraz dane źródłowe dla wszystkich rysunków i tabel zostały zdeponowane w Zenodo (DOI: https://doi.org/10.5281/zenodo.216406).

Rycina 1Przebieg procedury badania ankietowego odwiedzających Groty Longmen. Schemat przedstawia sekwencję działań zastosowaną w niniejszym protokole: od opracowania kwestionariusza, testów pilotażowych, przechwytywania odwiedzających w przestrzeni publicznej, przeprowadzania ankiet i blokowania surowych danych, aż po kontrolę jakości, punktowanie konstruktów, walidację pomiarów, estymację modelu, analizy wrażliwości i archiwizację plików badania. Proszę kliknąć tutaj, aby wyświetlić powiększoną wersję tej figury.

Rycina 2Schemat selekcji ankietowej dla końcowej próby analitycznej. Schemat przedstawia przejście od 480 surowych kwestionariuszy do 430 prawidłowych kwestionariuszy po przeprowadzeniu weryfikacji kryteriów kwalifikowalności, przeglądu pytań kontrolnych (attention-check), weryfikacji czasu wypełniania, wykryciu odpowiedzi jednolitych (straight-lining), kontroli duplikatów odpowiedzi oraz kontroli jakości związanych z pracami terenowymi. W przypadku wystąpienia więcej niż jednego warunku ostrzegawczego, wykluczenia zostały policzone według głównego powodu wykluczenia. Proszę kliknąć tutaj, aby wyświetlić powiększoną wersję tej ryciny.
Tabela 1: Mapa konstruktów i struktura kwestionariusza dla badania odwiedzających Groty Longmen.W tabeli podsumowano definicje konstruktów, skróty, kody pozycji, liczbę pozycji, kierunek punktacji, pozycję analityczną oraz źródło lub podstawę adaptacji wykorzystaną do pomiaru gotowości odwiedzających do konserwacji oraz deklarowanego odpowiedzialnego zachowania podczas podróży. Pełna treść pytań znajduje się w Pliku uzupełniającym 1. Wszystkie pytania merytoryczne były punktowane w skali od 1 = zdecydowanie się nie zgadzam do 5 = zdecydowanie się zgadzam. PC zachowano w oryginalnym kierunku; wyższy wynik PC wskazuje na silniejsze postrzegane zatłoczenie. Kliknij tutaj, aby pobrać tę tabelę.
Tabela 2: Harmonogram prac terenowych i systematyczny plan próbkowania metodą przechwytywania.Tabela ta przedstawia formalny okres zbierania danych, punkty badań publicznych, dzienne okna czasowe, interwał próbkowania, docelową liczbę próbek surowych, minimalne wymaganie dotyczące poprawnych próbek, zakresy monitorowania składu odwiedzających, zasady zawieszenia oraz zmienne rejestrowane podczas prac terenowych. Prosimy kliknąć tutaj, aby pobrać tę tabelę.
Tabela 3: Reguły kwalifikowalności, kontroli uwagi i screeningu jakości. Tabela ta przedstawia kryteria kwalifikowalności, reguły kontroli uwagi, progi czasu wykonania, reguły dotyczące powtarzalności odpowiedzi (straight-lining), sprawdzanie duplikatów odpowiedzi, kody wykluczenia oraz decyzje dotyczące zachowania danych, zastosowane przed punktowaniem konstruktów. Prosimy kliknąć tutaj, aby pobrać tę tabelę.
Tabela 4: Plan kodowania zmiennych i obliczania wyników konstruktu. Tabela ta definiuje reguły kodowania identyfikatora respondenta, zmiennych demograficznych, zmiennych związanych z podróżami, rekordów badania, zmiennych kontroli uwagi, elementów skali Likerta, wyników konstruktu oraz statusu wykluczenia. Wszystkie wyniki konstruktów pozostały w oryginalnej skali 1–5. Nie zastosowano imputacji brakujących wartości na poziomie elementów, ponieważ wszystkie merytoryczne pytania w kwestionariuszu były polami obowiązkowymi. Kliknij tutaj, aby pobrać tę tabelę.
Tabela 5: Wstępnie określone progi diagnostyczne stosowane przed interpretacją modelu.W tabeli przedstawiono progi dla niezawodności, retencji pozycji, analizy czynnikowej potwierdzającej, trafności dyskryminacyjnej, przesiewowej analizy wspólnej metody, kontroli rozkładu, przeglądu korelacji, współliniowości, diagnostyki przypadków wpływowych, wiarygodności modelu, testowania istotności oraz interpretacji efektów pośrednich. CFA = analiza czynnikowa potwierdzająca; CFI = wskaźnik dopasowania porównawczego;
TLI = wskaźnik Tuckera-Lewisa; RMSEA = pierwiastek błędu średniokwadratowego przybliżenia; SRMR = ustandaryzowany pierwiastek średniokwadratowego błędu resztowego; CR = niezawodność złożona; AVE = średnia wariancja wyodrębniona; HTMT = stosunek heterotrait-monotrait; VIF = współczynnik inflacji wariancji. Kliknij tutaj, aby pobrać tę tabelę.
Tabela 6: Rozkład prac terenowych i charakterystyka demograficzna reprezentatywnej próby. Tabela ta podsumowuje rozkład prawidłowo wypełnionych kwestionariuszy według punktu badawczego, okna czasowego, okresu badania, płci, grupy wiekowej, wykształcenia, miesięcznego dochodu osobistego, statusu wizyty oraz środka transportu. Procenty obliczono dla końcowej próby analitycznej, n = 430. Suma procentów może nie wynosić 10,0 ze względu na zaokrąglenia. Kliknij tutaj, aby pobrać tę tabelę.
Tabela 7: Podsumowanie rozkładu konstruktów, diagnostyki pomiarowej, wyników modelu i analizy wrażliwości.Tabela ta zawiera podsumowanie statystyk opisowych na poziomie konstruktów, wartości alfa Cronbacha, indeksy dopasowania CFA, rzetelność złożoną, średnią wariancję wyekstrahowaną, HTMT, wyniki badania błędów wspólnej metody, dopasowanie modelu regresji, dopasowanie modelu ścieżek, testowanie efektów pośrednich, diagnostykę współliniowości oraz stabilność analizy wrażliwości. Pełne wyniki rzetelności na poziomie pozycji, ładunki CFA, macierz HTMT, korelacje między konstruktami, współczynniki regresji, estymacje ścieżek, efekty pośrednie, efekty całkowite oraz wyniki analizy wrażliwości znajdują się w Pliku Uzupełniającym 3. HK = wiedza o dziedzictwie; PHV = postrzegana wartość dziedzictwa; PA = przywiązanie do miejsca; PMQ = jakość proekologicznego przekazu; PC = postrzegane zatłoczenie; FSW = wsparcie obiektów i usług; CA = postawa konserwatorska; CW = gotowość do ochrony; RTB = samoopisowe odpowiedzialne zachowania podróżnicze. PC zachowano w pierwotnym kierunku; wyższa wartość PC wskazuje na silniejsze postrzegane zatłoczenie. Kliknij tutaj, aby pobrać tę tabelę.
Plik uzupełniający 1: Kwestionariusz i materiały dotyczące zgody.Plik ten zawiera końcowy kwestionariusz, stronę świadomej zgody, pytania sprawdzające kryteria kwalifikacyjności, pytania kontrolne uwagi, zapis wyświetlania kwestionariusza na platformie oraz zapis przeglądu pilotażowego.Prosimy kliknąć tutaj, aby pobrać ten plik.
Plik uzupełniający 2: Rejestry prac terenowych i zarządzania danymi. Plik ten zawiera formularz dziennika prac terenowych, rejestr przydziału ról, dziennik zarządzania danymi, formularz rejestru wykluczeń oraz strukturę książki kodowej.Proszę kliknąć tutaj, aby pobrać ten plik.
Plik uzupełniający 3: Materiały dotyczące analizy i powtarzalności.Plik ten zawiera zanonimizowany zbiór danych analitycznych, skrypt analizy, wyjście funkcji sessionInfo(), wyniki rzetelności na poziomie pozycji, wyniki CFA, macierz HTMT, macierz korelacji między konstruktami, wyniki regresji, wyniki modelu ścieżkowego, wyniki efektów pośrednich, wyniki kontroli wrażliwości, kompletny arkusz wyników modelu oraz listę kontrolną danych źródłowych dla rycin i tabel.Prosimy kliknąć tutaj, aby pobrać ten plik.
Tabela uzupełniająca 1: Statystyki opisowe na poziomie pozycji i wyniki rzetelności. Wszystkie skorygowane korelacje pozycja-całość wynosiły ≥0,30. PC = postrzegane zatłoczenie (perceived crowding).Prosimy kliknąć tutaj, aby pobrać ten plik.
Tabela uzupełniająca 2: Standaryzowane ładunki CFA, CR i AVE. CFA = analiza czynnikowa potwierdzająca; CR = niezawodność złożona; AVE = średnia wariancja wyjaśniona.Prosimy kliknąć tutaj, aby pobrać ten plik.
Tabela uzupełniająca 3: Macierz korelacji między konstruktami. HK = wiedza o dziedzictwie (heritage knowledge); PHV = postrzegana wartość dziedzictwa (perceived heritage value); PA = przywiązanie do miejsca (place attachment); PMQ = jakość proekologicznego przekazu (pro-environmental message quality); PC = postrzegane zatłoczenie (perceived crowding); FSW = wsparcie w zakresie infrastruktury i usług (facility and service support); CA = postawa konserwatorska (conservation attitude); CW = gotowość do ochrony (conservation willingness); RTB = deklarowane odpowiedzialne zachowania podróżnych (self-reported responsible travel behavior).Prosimy kliknąć tutaj, aby pobrać ten plik.
Tabela uzupełniająca 4: Pełne wyniki regresji. Uwzględniono zmienne kontrolne: płeć, grupę wiekową, wykształcenie, miesięczny dochód osobisty, status pierwszej wizyty oraz sposób dotarcia. VIF = czynnik inflacji wariancji.Proszę kliknąć tutaj, aby pobrać ten plik.
Tabela uzupełniająca 5: Wyniki modelu ścieżkowego i efektów pośrednich.CA = postawa proekologiczna; CW = chęć ochrony przyrody; RTB = samoopisowe odpowiedzialne zachowania podczas podróżowania.Prosimy kliknąć tutaj, aby pobrać ten plik.
Tabela uzupełniająca 6: Analizy wrażliwości dla głównych wyników modelu. CA = postawa wobec ochrony przyrody; CW = chęć ochrony przyrody; RTB = deklarowane odpowiedzialne zachowania podczas podróży.Prosimy kliknąć tutaj, aby pobrać ten plik.
Tabela uzupełniająca 7: Lista kontrolna danych źródłowych dla rysunków i tabel. Prosimy kliknąć tutaj, aby pobrać ten plik.