fix(waterdata): small coherence cleanups (annotations, column order, defensiveness)

Five small, low-risk fixes surfaced by the package review:

1. get_latest_continuous / get_latest_daily: `value` was annotated `int`, but
   every other getter (and the docstrings) use `str | Iterable[str]`; the `int`
   hint also rejected the multi-value list filtering the others advertise.
2. get_time_series_metadata: `thresholds` was annotated `int`, vs
   `float | list[float]` on get_combined_metadata for the same queryable.
3. _arrange_cols: the "move the synthetic per-record id column to the end" set
   omitted peak_id, channel_measurements_id, combined_meta_id, and
   field_series_id, so those four getters left their id at the front instead of
   the end like daily_id. Added them for consistent column layout.
4. _next_req_url: returned a falsy `href` ("") instead of None, contradicting
   its Optional[str] contract. Return None.
5. _get_resp_data (geopandas branch): mirror the non-geopandas branch's
   `f.get("id")` so a feature missing a top-level id yields None rather than a
   KeyError.

Co-Authored-By: Claude Opus 4.8 (1M context) <noreply@anthropic.com>

Author: thodson-usgs

dataretrieval/waterdata/api.py

@@ -768,7 +768,7 @@ def get_time_series_metadata(
     unit_of_measure: str | Iterable[str] | None = None,
     computation_period_identifier: str | Iterable[str] | None = None,
     computation_identifier: str | Iterable[str] | None = None,
-    thresholds: int | None = None,
+    thresholds: float | list[float] | None = None,
     sublocation_identifier: str | Iterable[str] | None = None,
     primary: str | Iterable[str] | None = None,
     parent_time_series_id: str | Iterable[str] | None = None,
@@ -1213,7 +1213,7 @@ def get_latest_continuous(
     approval_status: str | Iterable[str] | None = None,
     unit_of_measure: str | Iterable[str] | None = None,
     qualifier: str | Iterable[str] | None = None,
-    value: int | None = None,
+    value: str | Iterable[str] | None = None,
     last_modified: str | Iterable[str] | None = None,
     skip_geometry: bool | None = None,
     time: str | Iterable[str] | None = None,
@@ -1407,7 +1407,7 @@ def get_latest_daily(
     approval_status: str | Iterable[str] | None = None,
     unit_of_measure: str | Iterable[str] | None = None,
     qualifier: str | Iterable[str] | None = None,
-    value: int | None = None,
+    value: str | Iterable[str] | None = None,
     last_modified: str | Iterable[str] | None = None,
     skip_geometry: bool | None = None,
     time: str | Iterable[str] | None = None,

dataretrieval/waterdata/utils.py

@@ -803,7 +803,7 @@ def _next_req_url(
             continue
         href = link.get("href")
         if not href:
-            return href
+            return None
         # Refuse to follow a next-page link to a different host —
         # the request's headers/auth were minted for the original
         # host and shouldn't leak to whatever a poisoned response
@@ -905,7 +905,9 @@ def _get_resp_data(
 
     # Organize json into geodataframe and make sure id column comes along.
     df = gpd.GeoDataFrame.from_features(features)
-    df["id"] = pd.json_normalize(features)["id"].values
+    # Mirror the non-geopandas branch's defensive ``f.get("id")`` so a feature
+    # missing a top-level ``id`` yields None rather than a KeyError.
+    df["id"] = [f.get("id") for f in features]
     df = df[["id"] + [col for col in df.columns if col != "id"]]
 
     # If no geometry present, then return pandas dataframe. A geodataframe
@@ -1299,6 +1301,10 @@ def _arrange_cols(
             "daily_id",
             "continuous_id",
             "field_measurement_id",
+            "field_series_id",
+            "peak_id",
+            "channel_measurements_id",
+            "combined_meta_id",
         }
     )